后缀数组(UVa11107)

第一次写后缀数组啦啦啦。。。

对于后缀数组的认知问题,我们还是看看大神罗穗骞的论文吧:《后缀数组--处理字符串的有力工具》


好,那下面我来提一提模板题UVa11107的解法:

题目链接:UVA-11107

题目大意:n个字符串,求长度最大字符串,要求在超过一半的字符串中出现。

题解:

1、先在中间用分隔符把他们拼成一个长字符串,求出他的后缀数组和height数组。

2、二分答案:每次判断是否有一个长度为k的串在超过一半的串中出现

  判断方法:先扫描一遍height数组,分段,每段的height的值都不小于k。

         这样,只要某一段中包含了超过n/2个原串的后缀,k就是满足条件的。

对于后缀数组的题目,使用二分答案和以height分段的方法非常常见,一定要掌握。


下面附上使用倍增法的代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<vector>
using namespace std;

#define maxn 140080
#define inf 0x3f3f3f3f

int key[maxn],belong[maxn],str[maxn];
bool vis[120];
int sa[maxn],t[maxn],t2[maxn],c[maxn];
int height[maxn],Rank[maxn];  //注意尽量不要用rank,我一开始就是这样CE的

void build(int * s,int n,int m) //倍增法构建后缀数组
{
	int i,*x = t,*y = t2;
	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(int 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)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 = 0;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]] && y[sa[i-1] + k] == y[sa[i]+k] ? p-1:p++;
		if(p >= n) break;
		m = p;
	}
}

void getHeight(int *s,int n) //求出height
{
	int i,j,k = 0;
	for(i = 0;i < n;i++) Rank[sa[i]] = i;
	for(i = 0;i < n;i++)
	{
		if(k) k--;
		int j = sa[Rank[i]-1];
		while(s[i+k] == s[j+k]) k++;
		height[Rank[i]] = k;
	}
}

bool judge(int n,int len,int num)
{
	int cnt = 0;
	memset(vis,0,sizeof(vis));
	vis[0] = 1;
	if(!vis[belong[sa[0]]])
	{
		cnt++;
		vis[belong[sa[0]]] = 1;
	}
	for(int i = 1;i < n;i++)
	{
		if(height[i] < len)
		{
			cnt = 0;
			memset(vis,0,sizeof(vis));
			vis[0] = 1;
			if(!vis[belong[sa[i]]])
			{
				vis[belong[sa[i]]] = 1;
				cnt++;
			}
		}
		else
			if(!vis[belong[sa[i]]])
			{
				cnt++;
				vis[belong[sa[i]]] = 1;
			}
		if(cnt >= num) return 1;
	}
	return 0;
}

void print(int n,int len,int num)
{
	int cnt = 0;
	memset(vis,0,sizeof(vis));
	vis[0] = 1;
	if(!vis[belong[sa[0]]]) cnt ++;
	vis[belong[sa[0]]] = 1;
	
	for(int i = 1;i < n;i++)
	{
		if(height[i] < len)
		{
			if(cnt >= num)
			{
				for(int j = sa[i-1];j < sa[i-1]+len;j++)
					printf("%c",str[j] - 18);
				printf("\n");
			}
			cnt = 0;
			memset(vis,0,sizeof(vis));
			vis[0] = 1;
			if(!vis[belong[sa[i]]])
			{
				cnt ++;
				vis[belong[sa[i]]] = 1;
			}
		}
		if(!vis[belong[sa[i]]])
		{
			cnt++;
			vis[belong[sa[i]]] = 1;
		}
	}
	if(cnt >= num)
	{
		for(int i = sa[n-1];i < sa[n-1]+len;i++)
			printf("%c",str[i]-18);
		printf("\n");
	}
}

inline int max(int a,int b){return a > b ? a : b;}

char s[1008];
int main()
{
	int n;
	bool flag = true;
	while(scanf("%d",&n) != EOF && n)
	{
		if(!flag) printf("\n");
		else flag = false;
		memset(belong,0,sizeof(belong));
		int pos = 0,cat = 1;
		int l = 0,r = 0;
		for(int i = 1;i <= n;i++)
		{
			scanf("%s",s);
			int len = strlen(s);
			r = max(r,len);
			for(int j = 0;j < len;j++)
			{
				str[pos+j] = (int)(s[j]) + 18;
				belong[pos+j] = i;
			}
			str[pos+len] = cat++;
			pos = pos + len + 1;
		}
		str[pos] = 0;
		build(str,pos+1,150);
		getHeight(str,pos+1);
		int maxlen = 0;
		while(l <= r)
		{
			int mid = (l+r) >> 1;
			if(judge(pos+1,mid,n/2+1))
			{
				maxlen = mid;
				l = mid + 1;
			}
			else r = mid - 1;
		}
		if(maxlen == 0) printf("?\n");
		else print(pos+1,maxlen,n/2+1);
	}
	return 0;
}


 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值