POJ2774 Long Long Message 【后缀数组lcp】

长长的消息

时间限制： 4000MS 内存限制： 131072K提交总数： 32393 接受： 13079案件时间限制： 1000MS

描述

小猫在拜特兰的首府物理专业。最近有一个不幸的消息传来，他的母亲生病了。担心在火车票上花费太多（拜特兰是一个如此大的国家，他不得不在火车上花费16个小时到家乡），他决定只和母亲发短信。 

小猫生活在一个非常富有的家庭，所以他经常到移动服务中心去查看他在短信上花了多少钱。昨天，服务中心的电脑坏了，打印了两条很长的信息。很快就发现了这只聪明的小猫： 

1.所有消息中的字符都是小写的拉丁字母，没有标点符号和空格。 
2.所有短信都被相互追加 - （i + 1）短信直接在第i个短信之后 - 这就是为什么这两个短信很长。 
3.他自己的短信已经被附加在一起，但是由于破损的计算机，可能有很多冗余字符向左和向右出现。 
例如：如果他的短信是“motheriloveyou”，那么由该机器打印的长信息可能是“hahamotheriloveyou”，“motheriloveyoureally”，“motheriloveyouornot”，“bbbmotheriloveyouaaa”等之一 
。4.对于这些破碎的问题，猫已经两次打印他的原始文本（所以出现两个很长的消息）。即使两个打印的消息中的原始文本保持不变，两边的冗余字符也可能不同。 

给你这两条很长的信息，你必须输出由小猫写的最长可能的原始文本的长度。 

背景： 
Byterland移动服务中的短信以每字节为单位收费。这就是为什么这只小猫正在担心最长的原始文本能持续多久。 

为什么要你写程序？有四个方面的问题： 
1，这些小猫在物理课上很忙碌; 
小猫想保留他对母亲所说的话。 
POJ是一个很棒的在线评委。 
来自“现代汉英综合大词典”小猫想从POJ那里挣点钱，并试图说服他的母亲去看医生： 

输入

在两个输入行上分别输入两个小写字母的字符串。每个字符的数量不会超过100000。

产量

单个整数的单行 - 小猫写的原始文本的最大长度是多少？

示例输入

yeshowmuchiloveyoumydearmotherreallyicannotbelieveityeaphowmuchiloveyoumydearmother

示例输出

27

经典的双串匹配最长公共长度问题

我们在第一个串末尾末尾加一个'#'这样的东西，在把第二个串连上去，这样子#一定不会参与匹配而且我们还把两个串分隔开了

接下来就是求lcp了，我们只需找到跨越#的最大的height就可以了

#include<iostream>#include<cstdio>#include<cstring>#include<algorithm>#define LL long long int#define REP(i,n) for (int i = 1; i <= (n); i++)#define Redge(u) for (int k = head[u]; k != -1; k = edge[k].next)using namespace std;const int maxn = 200005,maxm = 100005,INF = 1000000000;inline int RD(){int out = 0,flag = 1; char c = getchar();while (c < 48 || c > 57) {if (c == '-') flag = -1; c = getchar();}while (c >= 48 && c <= 57) {out = (out << 1) + (out << 3) + c - '0'; c = getchar();}return out * flag;}char s[maxn],B[maxn];int la,lb,n,m = 256,sa[maxn],rank[maxn],height[maxn],t1[maxn],t2[maxn],c[maxn];void SA(){int *x = t1,*y = t2;for (int i = 0; i <= m; i++) c[i] = 0;for (int i = 1; i <= n; i++) c[x[i] = s[i]]++;for (int i = 1; i <= m; i++) c[i] += c[i - 1];for (int i = n; i >= 1; i--) sa[c[x[i]]--] = i;for (int k = 1; k <= n; k <<= 1){int p = 0;for (int i = n - k + 1; i <= n; i++) y[++p] = i;for (int i = 1; i <= n; i++) if (sa[i] - k > 0) y[++p] = sa[i] - k;for (int i = 0; i <= m; i++) c[i] = 0;for (int i = 1; i <= n; i++) c[x[y[i]]]++;for (int i = 1; i <= m; i++) c[i] += c[i - 1];for (int i = n; i >= 1; i--) sa[c[x[y[i]]]--] = y[i];swap(x,y);p = 1; x[sa[1]] = 1;for (int i = 2; i <= n; i++)x[sa[i]] = (y[sa[i]] == y[sa[i - 1]] && y[sa[i] + k] == y[sa[i - 1] + k]) ? p : ++p;if (p >= n) break;m = p;}for (int i = 1; i <= n; i++) rank[sa[i]] = i;int k = 0;for (int i = 1; i <= n; i++){if (k) k--;int j = sa[rank[i] - 1];while (s[i + k] == s[j + k]) k++;height[rank[i]] = k;}}void solve(){int ans = 0;for (int i = 2; i <= n; i++){if (sa[i] <= la && sa[i - 1] > la || sa[i] > la && sa[i - 1] <= la){ans = max(ans,height[i]);}}cout<<ans<<endl;}int main(){scanf("%s",s + 1); la = strlen(s + 1);scanf("%s",B + 1); lb = strlen(B + 1);strcat(s + 1,"#"); strcat(s + 1,B + 1);n = strlen(s + 1);SA();solve();return 0;}