V.UVA11107 Life Forms
这题同上题类似，只不过把“在全部串中出现”变成了“在超过一半（即\(\left\lfloor\dfrac{n}{2}\right\rfloor+1\)）个串中出现”。
这题中我的方法是上题中提到的“two-pointers+单调队列”算法。第一遍跑求出所有满足“出现次数”为\(\left\lfloor\dfrac{n}{2}\right\rfloor+1\)的极小区间，因为小区间的\(\min ht_i\)肯定大于大区间的。
在单调队列求出公共子串长度的最大值后，需要通过关于\(<ans\)的位置分段（即二分中check的方法）来输出答案——因为按照上面的单调队列算法某些子串有可能会被重复输出。
复杂度\(O(n)\)，假如你用DC3算法的话可惜我不会。
代码：
#include<bits/stdc++.h>
using namespace std;
int all,n,m,x[200100],y[200100],buc[200100],sa[200100],ht[200100],rk[200100],id[200100],s[200100],occ[200100];
char str[200100];
void SA(){
for(int i=0;i<=m;i++)buc[i]=0;
for(int i=0;i<n;i++)buc[x[i]=s[i]]++;
for(int i=1;i<=m;i++)buc[i]+=buc[i-1];
for(int i=n-1;i>=0;i--)sa[--buc[x[i]]]=i;
for(int k=1;k<n;k<<=1){
int num=0;
for(int i=n-k;i<n;i++)y[num++]=i;
for(int i=0;i<n;i++)if(sa[i]>=k)y[num++]=sa[i]-k;
for(int i=0;i<=m;i++)buc[i]=0;
for(int i=0;i<n;i++)buc[x[y[i]]]++;
for(int i=1;i<=m;i++)buc[i]+=buc[i-1];
for(int i=n-1;i>=0;i--)sa[--buc[x[y[i]]]]=y[i];
swap(x,y);
x[sa[0]]=num=0;
for(int i=1;i<n;i++)x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k])?num:++num;
m=num;
}
for(int i=0;i<n;i++)rk[sa[i]]=i;
for(int i=0,k=0;i<n;i++){
if(!rk[i])continue;
if(k)k--;
int j=sa[rk[i]-1];
while(i+k<n&&j+k<n&&s[i+k]==s[j+k])k++;
ht[rk[i]]=k;
}
}
deque<int>q;
int main(){
scanf("%d",&all);
while(all){
n=0;
for(int i=1;i<=all;i++){
scanf("%s",str);
m=strlen(str);
for(int j=0;j<m;j++)id[n]=i,s[n]=(int)str[j]+5,n++;
id[n]=-1,s[n]=i,n++;
}
m='z'+5;
SA();
// for(int i=0;i<n;i++)printf("%2d ",id[sa[i]]);puts("");
// for(int i=0;i<n;i++)printf("%2d ",ht[i]);puts("");
int res=0;
for(int i=0,j=0,k=0;i<n;i++){
while(j<n&&k<all/2+1){
if(id[sa[j]]!=-1)k+=!occ[id[sa[j]]]++;
while(!q.empty()&&ht[q.back()]>=ht[j])q.pop_back();
q.push_back(j++);
}
while(!q.empty()&&q.front()<=i)q.pop_front();
if(k>=all/2+1)res=max(res,ht[q.front()]);
if(id[sa[i]]!=-1)k-=!--occ[id[sa[i]]];
}
if(!res)puts("?");
else for(int i=0,j=0;i<n;i=j){
if(id[sa[i]]==-1){j++;continue;}
int k=0;
k+=!occ[id[sa[i]]],occ[id[sa[i]]]=true;
for(j=i+1;j<n&&ht[j]>=res;j++)k+=!occ[id[sa[j]]],occ[id[sa[j]]]=true;
for(int l=i;l<j;l++)occ[id[sa[l]]]=false;
if(k>=all/2+1){for(int l=0;l<res;l++)printf("%c",s[sa[i]+l]-5);puts("");}
}
for(int i=0;i<n;i++)s[i]=sa[i]=ht[i]=rk[i]=x[i]=y[i]=id[i]=0;
scanf("%d",&all);
if(all)puts("");
}
return 0;
}