SPOJ PHRASES Relevant Phrases of Annihilation

这道题注意要是不重叠的,一开始这里WA了一次

其他的思路应该挺简单的

#include<stdio.h>
#include<algorithm>
#include<string.h>
using namespace std;
#define N 100100
int r[N];
char s[11][10005];
int wa[N],wb[N],wv[N],ws[N];
int sa[N],rank[N],height[N];
int m[N];
int cmp(int *r,int a,int b,int l){
	return r[a]==r[b]&&r[a+l]==r[b+l];
}
void da(int *r,int *sa,int n,int m){
	int i,j,p,*x=wa,*y=wb;
	for(i=0;i<m;i++) ws[i]=0;
	for(i=0;i<n;i++) ws[x[i]=r[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;m=p,j<<=1){
		for(i=n-j,p=0;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];
		swap(x,y);
		for(i=1,p=1,x[sa[0]]=0;i<n;i++) x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
	}
}
void calheight(int *r,int *sa,int n){
	int i,j,k=0;
	for(i=1;i<=n;i++) rank[sa[i]]=i;
	for(i=0;i<n;height[rank[i++]]=k)
		for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
	//for(i=1;i<=n;i++) printf("%d ",height[i]);
}
int vis[15],mi[15],ma[15];
int judge(int k,int l,int n){//注意要是不重叠的
	int i,j,size=0;
	memset(vis,0,sizeof(vis));
	for(i=1;i<=l;i++){
		if (height[i]>=k){
			if(vis[m[sa[i]]]==0) {
				mi[m[sa[i]]]=sa[i];
				ma[m[sa[i]]]=sa[i]+k-1;
				vis[m[sa[i]]]++;
			}
			else if(vis[m[sa[i]]]==1){
				if(sa[i]>ma[m[sa[i]]]||sa[i]+k-1<mi[m[sa[i]]]) {
					vis[m[sa[i]]]++;
					size++;
				}
				else{
					mi[m[sa[i]]]=max(sa[i],mi[m[sa[i]]]);
					ma[m[sa[i]]]=min(sa[i]+k-1,ma[m[sa[i]]]);
				}
			}
		}
		else{
			size=0;
			memset(vis,0,sizeof(vis));
			if(m[sa[i]]>=0){
				vis[m[sa[i]]]++;
				mi[m[sa[i]]]=sa[i];
				ma[m[sa[i]]]=sa[i]+k-1;
			}
		}
		if(size>=n) return 1;
	}
	return 0;
}
int main(){
	#ifndef ONLINE_JUDGE
	freopen("in.txt","r",stdin);
	#endif
	int T,n;
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);
		memset(m,-1,sizeof(m));
		int i,j,l=0,max=129;
		for(i=0;i<n;i++) scanf("%s",s[i]);
		for(i=0;i<n;i++){
			int len=strlen(s[i]);
			for(j=0;j<len;j++){
				m[l]=i;
				r[l++]=s[i][j];
			}
			r[l++]=max++;
		}
		r[l]=0;
		da(r,sa,l+1,max);
		calheight(r,sa,l);
		int low=0,high=10005,mid;
		while(low+1<high){
			mid=(low+high)>>1;
			if(judge(mid,l,n)) low=mid;
			else high=mid;
		}
		printf("%d\n",low);
	}
	return 0;
}
时间: 2024-08-03 09:28:15

SPOJ PHRASES Relevant Phrases of Annihilation的相关文章

SPOJ 220 Relevant Phrases of Annihilation (后缀数组)

题目大意: 求在m个串中同时出现两次以上且不覆盖的子串的长度. 思路分析: 二分答案,然后check是否满足,判断不覆盖的方法就是用up down 来处理边界. #include <cstdio> #include <iostream> #include <algorithm> #include <cstring> #include <map> #include <string> #define maxn 110005 using n

SPOJ 220. Relevant Phrases of Annihilation(后缀数组多次不重叠子串)

题目大意:给定N个串,求每个串至少出现两次的最长子串. 解题思路:每个字符串至少出现两次且不可重叠的最长子串:二分枚举长度后在同一分组中对每一个字符串保留一个最小的位置和一个最大的位置,最后查看是否每个串在同一组中都有至少两个后缀,并且后缀的坐标差大于枚举的长度. POJ Problem Set (classical) 220. Relevant Phrases of Annihilation Problem code: PHRASES You are the King of Byteland.

SPOJ 220 Relevant Phrases of Annihilation(后缀数组+二分答案)

[题目链接] http://www.spoj.pl/problems/PHRASES/ [题目大意] 求在每个字符串中出现至少两次的最长的子串 [题解] 注意到这么几个关键点:最长,至少两次,每个字符串. 首先对于最长这个条件,我们可以想到二分答案, 然后利用后缀数组所求得的三个数组判断是否满足条件. 其次是出现两次,每次出现这个条件的时候, 我们就应该要想到这是最大值最小值可以处理的, 将出现在同一个字符串中的每个相同字符串的起始位置保存下来, 如果最小值和最大值的差距超过二分长度L,则表明在

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

You are the King of Byteland. Your agents have just intercepted a batch of encrypted enemy messages concerning the date of the planned attack on your island. You immedietaly send for the Bytelandian Cryptographer, but he is currently busy eating popc

SPOJ220 Relevant Phrases of Annihilation

http://www.spoj.com/problems/PHRASES/ 题意:给n个串,求n个串里面都有2个不重叠的最长的字串长度. 思路:二分答案,然后就可以嘿嘿嘿 PS:辣鸡题目毁我青春,一开始二分的时候ans没有赋初值为0,结果没答案的时候就会输出奇怪的数字T_T,其实主要还是怪我不小心.. 1 #include<cstdio> 2 #include<iostream> 3 #include<cmath> 4 #include<cstring> 5

【SPOJ220】Relevant Phrases of Annihilation (SA)

成功完成3连T!   嗯没错,三道TLE简直爽到不行,于是滚去看是不是模版出问题了..拿了3份其他P党的模版扔上去,嗯继续TLE...蒟蒻表示无能为力了... 思路像论文里面说的,依旧二分长度然后分组...然后记录下每个字符的最大和最小值去判断是否满足全部成立...完事...写起来其实蛮简单的... const maxn=100419; var h,sum,rank,x,y,sa,c,lx,rx,col:array[0..maxn] of longint; n,k,maxlen,t,q:long

SPOJ220 Relevant Phrases of Annihilation(后缀数组)

引用罗穗骞论文中的话: 先将n 个字符串连起来,中间用不相同的且没有出现在字符串中的字符隔开,求后缀数组.然后二分答案,再将后缀分组.判断的时候,要看是否有一组后缀在每个原来的字符串中至少出现两次,并且在每个原来的字符串中,后缀的起始位置的最大值与最小值之差是否不小于当前答案(判断能否做到不重叠,如果题目中没有不重叠的要求,那么不用做此判断).这个做法的时间复杂度为O(nlogn). 二分枚举长度,对每个长度遍历height[]数组,将height[]数组分块,每个块内任意两串的lcp均大于等于

【SPOJ220】Relevant Phrases of Annihilation(后缀数组,二分)

题意: n<=10,len<=1e4 思路: 1 #include<cstdio> 2 #include<cstring> 3 #include<string> 4 #include<cmath> 5 #include<iostream> 6 #include<algorithm> 7 #include<map> 8 #include<set> 9 #include<queue> 10

后缀数组 &amp; 题目

后缀数组被称为字符串处理神器,要解决字符串问题,一定要掌握它.(我这里的下标全部都是从1开始) 首先后缀数组要处理出两个数组,一个是sa[],sa[i]表示排名第i为的后缀的起始位置是什么,rank[i]表示第i个字符为起始点的后缀,它的排名是什么.可以知道sa[rank[i]] = i; rank[sa[i]] = i; 由于每个后缀各不相同,至起码长度不同,所以每个后缀是不可能相等的. 解除一个值,就能在O(n)时间内得到另外一个. 定义:suffix(i)表示从[i, lenstr]这个后