博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
POJ 3729 Facer’s string(后缀数组)
阅读量:6682 次
发布时间:2019-06-25

本文共 1532 字,大约阅读时间需要 5 分钟。

题目链接:http://poj.org/problem?id=3729

题意:给出两个数字串A和B以及K。求A有多少个长度为K的子串同时是B的子串且这个子串加上在A中后面一个数字时不是B的子串?

思路:统计大于等于K的减去大于等于K+1的。

int r[N],sa[N],wa[N],wb[N],wd[N],rank[N],h[N];int cmp(int *r,int a,int b,int len){    return r[a]==r[b]&&r[a+len]==r[b+len];}void da(int *r,int *sa,int n,int m){    int i,j,p,*x=wa,*y=wb,*t;    FOR0(i,m) wd[i]=0;    FOR0(i,n) wd[x[i]=r[i]]++;    FOR1(i,m-1) wd[i]+=wd[i-1];    FORL0(i,n-1) sa[--wd[x[i]]]=i;    for(j=1,p=1;p
<<=1,m=p) { p=0; FOR(i,n-j,n-1) y[p++]=i; FOR0(i,n) if(sa[i]>=j) y[p++]=sa[i]-j; FOR0(i,m) wd[i]=0; FOR0(i,n) wd[x[i]]++; FOR1(i,m-1) wd[i]+=wd[i-1]; FORL0(i,n-1) sa[--wd[x[y[i]]]]=y[i]; t=x;x=y;y=t;p=1;x[sa[0]]=0; FOR1(i,n-1) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++; }}void calHeight(int *r,int *sa,int n){ int i,j,k=0; FOR1(i,n) rank[sa[i]]=i; FOR0(i,n) { if(k) k--; j=sa[rank[i]-1]; while(i+k
M) break; R=L; while(R<=M&&h[R]>=x) R++; k=t=0; FOR(j,L-1,R-1) { if(sa[j]
n) t=1; } ans+=(i64)k*t; i=R; } return ans;}int main(){ while(scanf("%d%d%d",&n,&m,&K)!=-1) { int i,j,L=0; FOR0(i,n) RD(r[L]),r[L]++,L++; r[L++]=10002; FOR0(i,m) RD(r[L]),r[L]++,L++; r[L]=0; da(r,sa,L+1,10003); calHeight(r,sa,L); M=L; i64 ans=cal(K)-cal(K+1); PR(ans); } return 0;}

  

转载地址:http://wqxao.baihongyu.com/

你可能感兴趣的文章
windows7系统笔记本设置成虚拟WiFi热点(即“无线路由器”)
查看>>
IE 的浏览器模式和文本模式(一)
查看>>
StaticFilesServer静态文件服务器
查看>>
对SpringAop的思考之基于jdk的动态代理
查看>>
openstack学习笔记五 多节点部署之 rabbitmq信息中枢与元数据
查看>>
count(*),count(1)和count(主键)的区别
查看>>
揭秘设计模式:适配器模式(Adapter)
查看>>
centos救援模式修改root密码
查看>>
sysprep重置windows,封装系统
查看>>
我的友情链接
查看>>
冒泡排序
查看>>
sed学习笔记-2
查看>>
Linux系统的启动和修复模式
查看>>
Citrix ICA协议简要介绍
查看>>
软件发布版本区别介绍
查看>>
python操作selenium的基本操作
查看>>
kvm虚拟机迁移
查看>>
Docker 修改docker容器内部时间
查看>>
解决windows下redis狂占C盘内存
查看>>
yii2高级模板添加新增模块
查看>>