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那里挣点钱，并试图说服他的母亲去看医生：
小猫生活在一个非常富有的家庭，所以他经常到移动服务中心去查看他在短信上花了多少钱。昨天，服务中心的电脑坏了，打印了两条很长的信息。很快就发现了这只聪明的小猫：
1.所有消息中的字符都是小写的拉丁字母，没有标点符号和空格。
2.所有短信都被相互追加 - （i + 1）短信直接在第i个短信之后 - 这就是为什么这两个短信很长。
3.他自己的短信已经被附加在一起，但是由于破损的计算机，可能有很多冗余字符向左和向右出现。
例如：如果他的短信是“motheriloveyou”，那么由该机器打印的长信息可能是“hahamotheriloveyou”，“motheriloveyoureally”，“motheriloveyouornot”，“bbbmotheriloveyouaaa”等之一
。4.对于这些破碎的问题，猫已经两次打印他的原始文本（所以出现两个很长的消息）。即使两个打印的消息中的原始文本保持不变，两边的冗余字符也可能不同。
给你这两条很长的信息，你必须输出由小猫写的最长可能的原始文本的长度。
背景：
Byterland移动服务中的短信以每字节为单位收费。这就是为什么这只小猫正在担心最长的原始文本能持续多久。
为什么要你写程序？有四个方面的问题：
1，这些小猫在物理课上很忙碌;
小猫想保留他对母亲所说的话。
POJ是一个很棒的在线评委。
来自“现代汉英综合大词典”小猫想从POJ那里挣点钱，并试图说服他的母亲去看医生：
输入
在两个输入行上分别输入两个小写字母的字符串。每个字符的数量不会超过100000。
产量
单个整数的单行 - 小猫写的原始文本的最大长度是多少？
示例输入
yeshowmuchiloveyoumydearmotherreallyicannotbelieveit yeaphowmuchiloveyoumydearmother
示例输出
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;
}
