hdu 6194
来源:互联网 发布:万网域名怎么转让 编辑:程序博客网 时间:2024/06/06 03:09
题目链接:点击打开链接
题解思路:用后缀数组处理得到的hight数组后,对它再进行线段建树,维护的是区间内最小的一个值,那么我们每次取枚举一个长度为n的区间时,取它的最小值这表示从1-ans_min的字符串长度都是满足的,那么题目还要求正好n次重复,那么我们就再取两次一个的左边界加一的区间值,一个是右边界加一的区间值,取两值的max那么此区间可取范围就应该是(左右边界+1的max,ans_min],时间复杂度就应该为lenlog(len)
代码:
#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int mx = 1e5+10;char str[mx];int c[mx],wa[mx],wb[mx],sa[mx],hight[mx],ran[mx];void sort_sa(int m,int n){ int *x=wa,*y=wb,p=0; for(int i=0;i<m;i++) c[i]=0; for(int i=0;i<n;i++) c[x[i]=str[i]]++; for(int i=1;i<m;i++) c[i]+=c[i-1]; for(int i=n-1;i>=0;i--) sa[--c[x[i]]]=i; for(int k=1;p<n;k<<=1,m=p){ p=0; for(int i=n-k;i<n;i++) y[p++]=i;//按个位排序 for(int i=0;i<n;i++) if(sa[i]>=k) y[p++]=sa[i]-k; for(int i=0;i<m;i++) c[i]=0; for(int i=0;i<n;i++) c[x[y[i]]]++; for(int i=1;i<m;i++) c[i]+=c[i-1]; for(int i=n-1;i>=0;i--) sa[--c[x[y[i]]]]=y[i]; swap(x,y);x[sa[0]] = 0,p = 1; for(int i=1;i<n;i++) if(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]) x[sa[i]]=p-1; else x[sa[i]]=p++; }}void Gethight(){ int k=0,len=strlen(str); for(int i=1;i<=len;i++) ran[sa[i]]=i; for(int i=0;i<len;i++){ if(k) k--; int j=sa[ran[i]-1]; while(str[j+k]==str[i+k]) k++; hight[ran[i]]=k; }}int sum[mx<<2];void build(int l,int r,int rt){ if(l==r){ sum[rt] = hight[l]; return ; } int mid = (l+r)>>1; build(l,mid,rt<<1); build(mid+1,r,rt<<1|1); sum[rt] = min(sum[rt<<1],sum[rt<<1|1]);}int query(int L,int R,int l,int r,int rt){ if(L<=l&&R>=r) return sum[rt]; int mid = (l+r)>>1; int ans = 1e9; if(mid>=L) ans = query(L,R,l,mid,rt<<1); if(mid<R) ans = min(ans,query(L,R,mid+1,r,rt<<1|1)); return ans;}int main(){ int t,n; scanf("%d",&t); while(t--){ scanf("%d",&n); scanf("%s",str); int len=strlen(str),ans = 0; sort_sa(200,len+1); Gethight(); hight[len+1] = 1e9; build(1,len,1); if(n>len) { puts("0"); continue; } if(n==1){ ans += len - sa[1] - hight[2]; ans += len - sa[len] - hight[len]; for(int i=2;i<len;i++) ans += len - sa[i] - max(hight[i],hight[i+1]); printf("%d\n",ans); continue; } int R = query(2,n,1,len,1),L = min(R,hight[n+1]); ans += R-L; for(int i=3;i<=len-n+1;i++){ R = query(i,i+n-2,1,len,1),L = max(min(R,hight[i+n-1]),min(R,hight[i-1])); ans += R-L; } R = query(len-n+2,len,1,len,1),L = min(R,hight[len-n+1]); ans += R-L; printf("%d\n",ans); } return 0;}
阅读全文
0 0
- hdu 6194
- hdu 6194
- HDU 6194
- hdu 6194 后缀数组
- hdu
- hdu
- HDU
- hdu ()
- hdu
- hdu
- HDU
- HDU
- hdu
- hdu
- HDU
- Hdu
- hdu
- hdu-
- http协议无状态中的 "状态" 到底指的是什么?!
- CodeForces 614D 二分+贪心
- Java内存区域与内存溢出异常
- HDOJ HDU 1085 Holding Bin-Laden Captive!
- 高复用率的RTSPClient组件EasyRTSPClient调用说明
- hdu 6194
- 001_JavaWeb中文乱码问题
- http invoker
- 高复用率的RTSPClient组件EasyRTSPClient设计流程
- python 迭代器和生成器
- 九度[1095]-2的幂次方
- hibernate.current_session_context_class
- bzoj3676 回文串 回文树
- Kafka设计解析(一)- Kafka背景及架构介绍