题意:求一个串中不重叠出现两次以上的子串的数量
解法:后缀数组的分组解法实在是经典 但是今天我明显非主流 然后使用了 后缀数组以后没有进行分组 而是rmq与处理后判断所有后缀的height最多可以有多少成为该串的做法 这当然还需要字符串的hash 但是WA了 然后就没有然后了
以下是老老实实的分组做法:
#include<map>
#include<math.h>
#include<cstdio>
#include<string.h>
#include<string>
#include<algorithm>
using namespace std;
#define ll unsigned long long
#define mod 131
#define maxn 11111
char txt[maxn];
int sa[maxn],no_1[maxn],no_2[maxn],ex[maxn],he[maxn],rk[maxn],n,limi;
int mi[25][maxn],lg[maxn],len;
void _sa(char *st){
int i,k,p,*x=no_1,*y=no_2;
for(i=0; i<limi; i++) ex[i]=0;
for(i=0; i<n; i++) ex[x[i]=st[i]]++;
for(i=1; i<limi; i++) ex[i]+=ex[i-1];
for(i=n-1; i>=0; i--)
sa[--ex[x[i]]]=i;
for(k=1,p=1; p<n; k<<=1,limi=p){
for(p=0,i=n-k; i<n; i++) y[p++]=i;
for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
for(i=0; i<limi; i++) ex[i]=0;
for(i=0; i<n; i++) ex[x[y[i]]]++;
for(i=1; i<limi; i++) ex[i]+=ex[i-1];
for(i=n-1; i>=0; i--) sa[--ex[x[y[i]]]]=y[i];
for(swap(x,y),p=1,x[sa[0]]=0,i=1; i<n; i++)
x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
}
}
void _he(char *st){
int i,j,k=0;
for(i=0;i<n;i++) rk[sa[i]]=i;
for(i=0;i<n-1;i++)
{
if(k) k--;
j=sa[rk[i]-1];
while(st[i+k]==st[j+k]) k++;
he[rk[i]]=k;
}
}int q,x,y;
//char now[1111];
void rmq_init(){
for(int i=0;i<n;i++){mi[0][i]=he[i];}
for(int i=1;i<=lg[n];i++)
for(int j=0;j+(1<<i)-1<n;j++){
mi[i][j]=min(mi[i-1][j],mi[i-1][j+(1<<(i-1))]);
}
}
int rmq_min(int l,int r){
if(l>r){
l--;r++;swap(l,r);
}
int tmp=lg[r-l+1];
return min(mi[tmp][l],mi[tmp][r-(1<<tmp)+1]);
}
int abs(int x){return x>0?x:-x;}
map<ll,int>mp;
int main(){
lg[0]=-1;for(int i=1;i<maxn;i++)lg[i]=lg[i>>1]+1;
while(~scanf("%s",txt)){
mp.clear();
limi=255,n=(int)strlen(txt);n++;
_sa(txt);_he(txt);
rmq_init();
n--;int ans=0;
// for(int i=1;i<=n;++i){
//
// printf("%d\n",rmq_min(i,i));
//
// }
// printf("%d\n",rmq_min(2,3));
int mi=0,ma=0,fl=0;
for(int z=1;z<=n;++z){
fl=0;
for(int i=1;i<=n;++i){
// for(int a=sa[i];a<=n;++a)printf("%c",txt[a]);
// printf("\n");
// printf("%d %d\n",sa[i],he[i]);
if(he[i]>=z){
if(!fl){
fl=1;ma=max(sa[i],sa[i-1]);mi=min(sa[i],sa[i-1]);
}else{
ma=max(sa[i],ma);mi=min(sa[i],mi);
}
}else{
if(!fl){
// continue;
}else{
if(ma-mi>=z)ans++;
fl=0;
}
}
// printf("len:%d i:%d %d %d %d %d\n\n",z,i,ma,mi,ans,fl);
}
if(ma-mi>=z&&fl)ans++;
}
printf("%d\n",ans);
}
return 0;
}