[POJ2774]Long Long Message（后缀数组）

题目：

我是超链接

题意：给出两个字符串，求最长公共子串。

题解：

swap的操作最好用algorithm和G++.........

字符串的任何一个子串都是这个字符串某个后缀的前缀

我们把这两个字符串用一个不会出现的字符连接，做一个sa数组，这样当两者不是同一个字符串的后缀时，最大的height就是答案

代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#define N 100005
using namespace std;
char s[N*2],s2[N];
int c[N*2],sa[N*2],y[N*2],x[N*2],height[N*2],rank[N*2],m=200,n;
void make_sa()//sa排第几的是谁，rank你排第几 
{
	int i,k;
	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 (k=1;k<=n;k<<=1)
	{
		int p=0;
		for (i=n-k;i<n;i++) y[p++]=i;
		for (i=0;i<n;i++) if (sa[i]-k>=0) y[p++]=sa[i]-k;
		
		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]] && ((sa[i-1]+k>=n?-1:y[sa[i-1]+k])==(sa[i]+k>=n?-1:y[sa[i]+k]))?p-1:p++;
         //后面千万不要忘了加（） 
		if (p>n) break; 
		m=p; 
	}
}
void make_height()
{
	int i;
	for (i=0;i<n;i++) rank[sa[i]]=i;
	int k=0; height[0]=0;
	for (i=0;i<n;i++)
	{
		if (!rank[i]) continue;
		if (k) k--;
		int j=sa[rank[i]-1];
		while (i+k<n && j+k<n && s[i+k]==s[j+k]) k++;
		height[rank[i]]=k;
	}
}
int main()
{
	int i;
	scanf("%s",s); scanf("%s",s2); 
	int l1=strlen(s); 
	s[l1]='&'; s[l1+1]='\0';
	int l2=strlen(s2);
	strcat(s,s2);
	n=l1+l2+1;
	make_sa();
	make_height();
	int ans=0;
	for (i=1;i<n;i++)
	  if (sa[i-1]>l1 && sa[i]<l1 || sa[i-1]<l1 && sa[i]>l1)
	    ans=max(ans,height[i]);
	printf("%d",ans);
}