后缀数组的题目,把后缀连接起来,这个还是先二分答案,然后选取一段连续的height值,判断这些height代表的后缀有没有覆盖一半以上的字符串。
得出答案的长度之后还要在枚举连续的heigh,判断有没有答案,有的话标记其中一个。
最后再按照sa输出答案。这样就可以保证字典序。
#include <iostream> #include <cstdio> #include <cstring> using namespace std; const int maxn=1e5+10000; int r[maxn]; char a[1111]; int *rank,height[maxn],sa[maxn]; int wx[maxn],wy[maxn],cnt[maxn]; int n,m; inline bool cmp(int *r,int a,int b,int l) { return r[a]==r[b]&&r[a+l]==r[b+l]; } void da(int *r,int n,int m) { r[n+1]=0; int i,l,p,*x=wx,*y=wy,*t; memset(cnt,0,sizeof(int)*(m+1)); for(int i=1;i<=n;i++) cnt[x[i]=r[i]]++; for(int i=1;i<=m;i++) cnt[i]+=cnt[i-1]; for(int i=n;i>=1;i--) sa[cnt[x[i]]--]=i; for(l=1,p=1;p<n;l<<=1,m=p) { for(p=1,i=n-l+1;i<=n;i++) y[p++]=i; for(i=1;i<=n;i++) if(sa[i]>l) y[p++]=sa[i]-l; memset(cnt,0,sizeof(int)*(m+1)); for(i=1;i<=n;i++) cnt[x[i]]++; for(i=1;i<=m;i++) cnt[i]+=cnt[i-1]; for(i=n;i>=1;i--) sa[cnt[x[y[i]]]--]=y[i]; for(t=x,x=y,y=t,p=1,x[sa[1]]=1,i=2;i<=n;i++) x[sa[i]]=cmp(y,sa[i-1],sa[i],l)?p:++p; } rank=x; int j,k=0; for(i=1;i<=n;height[rank[i++]]=k) for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++); return; } bool fla[1111]; int c[maxn]; bool ans[maxn]; bool chk(int ret) { for(int i=2;i<=n;i++) if(height[i]>=ret) { int j=i+1; while(j<=n&&height[j]>=ret) j++; j--; memset(fla,0,sizeof(fla)); for(int k=i-1;k<=j;k++) fla[c[sa[k]]]=1; int cnt=0; for(int k=1;k<=m;k++) cnt+=fla[k]; if(cnt>m/2) return true; i=j; } return false; } int main() { // freopen("in.txt","r",stdin); int ret; int tt=0; while(scanf("%d",&m),m) { if(tt++) printf(" "); n=0; for(int i=1;i<=m;i++) { scanf("%s",a+1); ret=strlen(a+1); for(int j=1;j<=ret;j++) r[n+j]=a[j]; for(int j=n+1;j<=n+ret;j++) c[j]=i; n+=ret; r[++n]=570+i; } da(r,n,680); int st=0,ed=1111,mid; while(st<ed) { mid=st+ed+1>>1; if(chk(mid)) st=mid; else ed=mid-1; } if(m==1) { for(int i=1;i<n;i++) printf("%c",r[i]); printf(" "); } else if(st==0) printf("? "); else { memset(ans,0,sizeof(ans)); for(int i=1;i<=n;i++) if(height[i]>=st) { int j=i+1; while(height[j]>=st) j++; j--; memset(fla,0,sizeof(fla)); for(int k=i-1;k<=j;k++) fla[c[sa[k]]]=1; int cnt=0; for(int k=1;k<=m;k++) cnt+=fla[k]; if(cnt>m/2) { ans[sa[i-1]]=1; } i=j; } for(int i=1;i<=n;i++) if(ans[sa[i]]) { for(int j=1;j<=st;j++) printf("%c",r[j+sa[i]-1]); printf(" "); } } } return 0; }