SPOJ - PHRASES Relevant Phrases of Annihilation(后缀数组+二分)

题目链接：点击查看

题目大意：给出n个字符串，求出n个字符串中的最长公共子串，且该字串必须在每个字符串中出现至少两次且不能重叠

题目分析：结合之前整理过的：

1. 一个字符串中两个不重叠的重复子串的最大长度
2. n个字符串中的出现次数大于等于k次的最长公共子串

这两点的解决方法，一结合就是本题的答案了，虽然说得轻松，但是比较考验代码能力，静下心来慢慢写就好啦

具体解决方法直接看我的模板博客最下方吧：https://blog.csdn.net/qq_45458915/article/details/104026037

代码：

#include<iostream>
#include<cstdio> 
#include<string>
#include<ctime>
#include<cstring>
#include<algorithm>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<cmath>
#include<sstream>
using namespace std;
 
typedef long long LL;
 
const int inf=0x3f3f3f3f;

const int N=1e5+100;
 
char str[N];
 
int sa[N]; //SA数组，表示将S的n个后缀从小到大排序后把排好序的
//的后缀的开头位置顺次放入SA中
int t1[N],t2[N],c[N];
 
int rk[N],height[N],belong[N],len,n,sum;
 
int s[N];

bool vis[15];

int mmax[15],mmin[15],cnt[15];
 
void build_sa(int s[],int n,int m)//n为添加0后的总长
{
    int i,j,p,*x=t1,*y=t2;
    for(i=0;i<m;i++) 
		c[i]=0;
    for(i=0;i<n;i++) 
		c[x[i]=s[i]]++;
    for(i=1;i<m;i++) 
		c[i]+=c[i-1];
    for(i=n-1;i>=0;i--) 
		sa[--c[x[i]]]=i;
    for(j=1;j<=n;j<<=1) 
	{
        p=0;
        for(i=n-j;i<n;i++) 
			y[p++]=i;
        for(i=0;i<n;i++) 
            if(sa[i]>=j) 
                y[p++]=sa[i]-j;
        for(i=0;i<m;i++) 
			c[i]=0;
        for(i=0;i<n;i++) 
			c[x[y[i]]]++;
        for(i=1;i<m;i++) 
			c[i]+=c[i-1];
        for(i=n-1;i>=0;i--) 
			sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=1,x[sa[0]]=0;
        for(i=1;i<n;i++) 
            x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+j]==y[sa[i]+j]?p-1:p++;
        if(p>=n) 
			break;
        m=p;
    }
}
 
void get_height(int s[],int n)//n为添加0后的总长
{
    int i,j,k=0;
    for(i=0;i<=n;i++)
        rk[sa[i]]=i;
    for(i=0;i<n;i++) 
	{
        if(k) 
			k--;
        j=sa[rk[i]-1];
        while(s[i+k]==s[j+k]) 
			k++;
        height[rk[i]]=k;
    }
}
 
void solve(int base=128)
{
	build_sa(s,len+1,base);
	get_height(s,len);
}

void init(int pos)
{
	memset(vis,false,sizeof(vis));
	sum=0;
	for(int i=1;i<=n;i++)
		mmax[i]=-inf,mmin[i]=inf,cnt[i]=0;
	int id=belong[sa[pos]];
	mmax[id]=mmin[id]=sa[pos];
	cnt[id]++;
}

bool check(int mid)
{
	init(1);
	for(int i=2;i<=len;i++)
	{
		if(!belong[sa[i]])
			break;
		if(height[i]<mid)
		{
			init(i);
		}
		else
		{
			int id=belong[sa[i]];
			mmax[id]=max(mmax[id],sa[i]);
			mmin[id]=min(mmin[id],sa[i]);
			cnt[id]++;
			if(cnt[id]>=2&&!vis[id]&&mmax[id]-mmin[id]>=mid)
			{
				sum++;
				vis[id]=true;
			}
		}
		if(sum==n)
			return true;
	}
	return false;
}

int main()
{
//	freopen("input.txt","r",stdin);
//	ios::sync_with_stdio(false);
	int w;
	cin>>w;
	while(w--)
	{
		scanf("%d",&n);
		int pos=30,tot=0;
		for(int i=1;i<=n;i++)
		{
			scanf("%s",str);
			for(int j=0;str[j];j++)
			{
				s[tot]=str[j]-'a'+1;
				belong[tot]=i;
				tot++;
			}
			s[tot++]=pos++;
		}
		s[tot]=0;
		len=tot;
		solve(pos);
		int l=0,r=inf,ans;
		while(l<=r)
		{
			int mid=l+r>>1;
			if(check(mid))
			{
				ans=mid;
				l=mid+1;
			}
			else
				r=mid-1;
		}
		printf("%d\n",ans);
	}
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	
	return 0;
}