POJ 2774 Long Long Message ——后缀数组

【题目分析】

用height数组RMQ的性质去求最长的公共子串。

要求sa[i]和sa[i-1]必须在两个串中,然后取height的MAX。

利用中间的字符来连接两个字符串的思想很巧妙,记得最后还需要空一个位置避免冲突。

【代码】

#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>

#include <map>
#include <set>
#include <queue>
#include <string>
#include <iostream>
#include <algorithm>

using namespace std;

#define maxn 200005
#define inf 0x3f3f3f3f
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)

void Finout()
{
    #ifndef ONLINE_JUDGE
    freopen("in.txt","r",stdin);
//    freopen("out.txt","w",stdout);
    #endif
}

int Getint()
{
    int x=0,f=1; char ch=getchar();
    while (ch<‘0‘||ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();}
    while (ch>=‘0‘&&ch<=‘9‘) {x=x*10+ch-‘0‘; ch=getchar();}
    return x*f;
}

int tot=0,l1,l2;
char s1[maxn],s2[maxn];

struct Suf_Arr{
	int s[maxn<<1],l;
	int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
	int sa[maxn],rank[maxn],height[maxn];
	int cmp(int *r,int a,int b,int l)
	{return r[a]==r[b]&&r[a+l]==r[b+l];}
	void getsa(int n,int m)
	{
		int i,j,p,*x=wa,*y=wb,*t;
		for (i=0;i<m;++i) ws[i]=0;
		for (i=0;i<n;++i) ws[x[i]=s[i]]++;
		for (i=1;i<m;++i) ws[i]+=ws[i-1];
		for (i=n-1;i>=0;--i) sa[--ws[x[i]]]=i;
		for (j=1,p=1;p<n;j*=2,m=p)
		{
			for (p=0,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<n;++i) wv[i]=x[y[i]];
			for (i=0;i<m;++i) ws[i]=0;
			for (i=0;i<n;++i) ws[wv[i]]++;
			for (i=1;i<m;++i) ws[i]+=ws[i-1];
			for (i=n-1;i>=0;--i) sa[--ws[wv[i]]]=y[i];
			for (t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;++i)
				x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
		}
	}
	void gethi(int n)
	{
		int i,j,k=0;
		for (i=1;i<=n;++i) rank[sa[i]]=i;
		for (int i=0;i<n;height[rank[i++]]=k)
			for (k?k--:0,j=sa[rank[i]-1];s[i+k]==s[j+k];k++);
	}
	void build()
	{
		getsa(l+2,30);
//		for (int i=0;i<=l;++i) cout<<sa[i]<<" "; cout<<endl;
		gethi(l+1);
	}
	void solve()
	{
		int maxx=0;
		for(int i=2;i<l;++i)
			if (height[i]>maxx)
			{
				if (0<=sa[i-1]&&sa[i-1]<l1&&l1<sa[i]) maxx=height[i];
				if (0<=sa[i]&&sa[i]<l1&&l1<sa[i-1]) maxx=height[i];
			}
		printf("%d\n",maxx);
	}
}arr;

int main()
{
    Finout();
    scanf("%s",s1); scanf("%s",s2);
    l1=strlen(s1); l2=strlen(s2); arr.l=l1+l2+1;
    for (int i=0;i<l1;++i) arr.s[i]=s1[i]-‘a‘+1;
    for (int i=0;i<l2;++i) arr.s[i+l1+1]=s2[i]-‘a‘+1;
    arr.s[l1+l2+1]=0; arr.s[l1]=28;
//    printf("len is %d\n",arr.l);
//    for (int i=0;i<=arr.l;++i) cout<<(char)(arr.s[i]+‘a‘-1); cout<<endl;
    arr.build();
    arr.solve();
}

  

时间: 2024-10-20 07:46:04

POJ 2774 Long Long Message ——后缀数组的相关文章

poj 2774 Long Long Message(后缀数组入门题)

1 /****************************************************************** 2 题目: Long Long Message(poj 2774) 3 链接: http://poj.org/problem?id=2774 4 题意: 给两个字符串,找最长的公共子串 5 算法: 后缀数组 6 算法思想: 后缀数组就是套模板求先应得数组,这题用到了两个数组,分 7 别是sa[],height[];sa[i]表示所有后缀按字典数排序后以s[i

poj 2774 Long Long Message 后缀数组基础题

Time Limit: 4000MS   Memory Limit: 131072K Total Submissions: 24756   Accepted: 10130 Case Time Limit: 1000MS Description The little cat is majoring in physics in the capital of Byterland. A piece of sad news comes to him these days: his mother is ge

poj 2774 Long Long Message 后缀数组LCP理解

题目链接 题意:给两个长度不超过1e5的字符串,问两个字符串的连续公共子串最大长度为多少? 思路:两个字符串连接之后直接后缀数组+LCP,在height中找出max同时满足一左一右即可: #include<iostream> #include<cstdio> #include<cstring> #include<string.h> #include<algorithm> #include<map> #include<queue&

poj 2774 Long Long Message 后缀数组

点击打开链接题目链接 Long Long Message Time Limit: 4000MS   Memory Limit: 131072K Total Submissions: 23327   Accepted: 9566 Case Time Limit: 1000MS Description The little cat is majoring in physics in the capital of Byterland. A piece of sad news comes to him

PKU 2774 Long Long Message (后缀数组练习模板题)

题意:给你两个字符串,求最长公共字串的长度. by:罗穗骞模板 #include <iostream> #include <stdio.h> #include <string.h> #include <algorithm> using namespace std; #define M 303 #define inf 0x3fffffff #define maxn 500000 #define ws ww #define rank RANK #define F

POJ 2774 Long Long Message (最长公共子串)

Long Long Message Time Limit: 4000MS   Memory Limit: 131072K Total Submissions: 27062   Accepted: 11010 Case Time Limit: 1000MS Description The little cat is majoring in physics in the capital of Byterland. A piece of sad news comes to him these days

POJ 1743 Musical Theme (后缀数组)

题目大意: 刚才上88个键弹出来的音符. 如果出现重复的,或者是高一个音阶的重复的都算. 思路分析: 具体可以参考训练指南222. height数组表示按照排序后的sa最近的两个后缀的最长前缀. 将height 分块.然后二分答案,二分答案之后去判断是否满足. 要考虑到不重合,还有大于5. 所以二分的时候要从5开始,然后判断的时候要加一个 up - down >len #include <cstdio> #include <iostream> #include <alg

HUID 5558 Alice&#39;s Classified Message 后缀数组+单调栈+二分

http://acm.hdu.edu.cn/showproblem.php?pid=5558 对于每个后缀suffix(i),想要在前面i - 1个suffix中找到一个pos,使得LCP最大.这样做O(n^2) 考虑到对于每一个suffix(i),最长的LCP肯定在和他排名相近的地方取得. 按排名大小顺序枚举位置,按位置维护一个递增的单调栈,对于每一个进栈的元素,要算一算栈内元素和他的LCP最大是多少. 如果不需要输出最小的下标,最大的直接是LCP(suffix(st[top]),  suff

POJ 3729 Facer’s string (后缀数组)

题目大意: 串1中有多少个后缀和 串2中的某个后缀 的lcp 为 k 思路分析: 先找出 长度至少为k的对数有多少. 再找出 至少为k+1的有多少 然后相减. #include <cstdio> #include <iostream> #include <algorithm> #include <cstring> #include <map> #include <string> #define maxn 110005 using na