HDU 4691 后缀数组+RMQ
思路:
求一发后缀数组,求个LCP 就好了
注意数字有可能不只一位 (样例2)
//By SiriusRen #include <bits/stdc++.h> using namespace std; const int N=100050; char s[N]; int cntA[N],cntB[N],A[N],B[N],sa[N],rk[N],tsa[N],ht[N],f[N][20]; int n,t,base[N],l[N],r[N]; void SA(){ for(int i=1;i<=256;i++)cntA[i]=0; for(int i=1;i<=n;i++)cntA[s[i]]++; for(int i=1;i<=256;i++)cntA[i]+=cntA[i-1]; for(int i=n;i;i--)sa[cntA[s[i]]--]=i; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(s[sa[i]]!=s[sa[i-1]]); for(int l=1;rk[sa[n]]<n;l<<=1){ memset(cntA,0,sizeof(cntA)); memset(cntB,0,sizeof(cntB)); for(int i=1;i<=n;i++)cntA[A[i]=rk[i]]++,cntB[B[i]=(i+l<=n?rk[i+l]:0)]++; for(int i=1;i<=n;i++)cntA[i]+=cntA[i-1],cntB[i]+=cntB[i-1]; for(int i=n;i;i--)tsa[cntB[B[i]]--]=i; for(int i=n;i;i--)sa[cntA[A[tsa[i]]]--]=tsa[i]; rk[sa[1]]=1; for(int i=2;i<=n;i++)rk[sa[i]]=rk[sa[i-1]]+(A[sa[i]]!=A[sa[i-1]]||B[sa[i]]!=B[sa[i-1]]); } for(int i=1,j=0;i<=n;i++){ j=j?j-1:0; while(s[i+j]==s[sa[rk[i]-1]+j])j++; ht[rk[i]]=j; } } void init_rmq(){ base[0]=-1; for(int i=1;i<=n;i++)f[i][0]=ht[i],base[i]=base[i>>1]+1; for(int j=1;j<=18;j++) for(int i=1;i+(1<<(j-1))<=n;i++) f[i][j]=min(f[i][j-1],f[i+(1<<(j-1))][j-1]); } int LCP(int x,int y){ if(x==y)return n-x+1; x=rk[x],y=rk[y]; if(x>y)swap(x,y);x++; int t=base[y-x+1]; return min(f[x][t],f[y-(1<<t)+1][t]); } int main(){ while(~scanf("%s",s+1)){ long long ans=0;n=strlen(s+1); SA(),init_rmq(); scanf("%d",&t); for(int i=1;i<=t;i++)scanf("%d%d",&l[i],&r[i]),ans+=r[i]-l[i]+1,l[i]++; printf("%lld ",ans); for(int i=2;i<=t;i++){ int tmp=min(LCP(l[i-1],l[i]),min(r[i]-l[i],r[i-1]-l[i-1])+1); ans+=(tmp?(int)log10(tmp):0)-tmp; } printf("%lld\n",ans+t*2); } }