[acm/icpc2016ChinaFinal][CodeforcesGym101194] Mr. Panda and Fantastic Beasts 后缀自动机

这题解法好多耶

听到有后缀自动机的写法,就为此专门补习后缀自动机的知识

把第2-n个字符串两两间用符号隔开(这个符号可以相同,只要不是小写字母即可),拼成长串插入后缀自动机

扫一遍第一个字符串,和后缀自动机求最长公共字串,失配的时候更新答案


后缀自动机也不只有一种方法可以做,可以在第2-n个串任意两个字符串间插入同一个特殊字符以后构建后缀自动机求最长公共字串,可以建广义后缀自动机,还可以对第一个串建后缀自动机后dp。题目细节不多写起来不难,很快就写完拍完。


除此之外,后缀数组和hash+二分答案都可以解决这道题

#include <iostream>
#include <cstdio>
#include <cstring>

#define N 600050

using namespace std;
char T[N],S[N];
int ch[N][28],len[N],link[N],last,cnt,rt;
int ans,cur,l,al,ar;

void add(int pos) {
	int x = S[pos] - 'a' + 1 , p = last , np = last = ++cnt;
	len[np] = len[p] + 1;
	
	while (p && !ch[p][x]) ch[p][x] = np , p = link[p];
	if (!p) 
		link[np] = rt;
	else {
		int q = ch[p][x];
		if (len[q] == len[p] + 1)
			link[np] = q;
		else {
			int nq = ++cnt;
			len[nq] = len[p] + 1;
			memcpy(ch[nq],ch[q],sizeof(ch[q]));
			link[nq] = link[q];
			link[np] = link[q] = nq;
			while (p && ch[p][x] == q) ch[p][x] = nq , p = link[p];
		}
	}
	return ;
}

void ut(int &al,int &ar,int l,int r) {
	if (r-l+1 < ar-al+1) {
		al = l; ar = r; return ;
	} else {
		int len = ar - al + 1;
		for (int j=0;j<len;j++) 
			if (T[al+j] != T[l+j]) {
				if (T[al+j] > T[l+j]) al = l , ar = r;
				return ;
			}
	}
	return ;
}

void gio(int pos) {
	int x = T[pos] - 'a' + 1;
	if (ch[cur][x]) 
		cur = ch[cur][x] , l = len[ link[cur] ] + 1;
	else {
		while (cur && !ch[cur][x]) {
			if (l+1 <= ans) {
				ans = l+1;
				ut(al,ar,pos-l,pos);
			}
			cur = link[cur] , l = len[ link[cur] ] + 1;
		}
		if (!cur) 
			cur = rt , l = 0;
		else 
			cur = ch[cur][x] , l = len[ link[cur] ] + 1;
	}
}

int main() {
	int cas = 0; scanf("%d",&cas);
	for (int _=1;_<=cas;_++) {
		printf("Case #%d: ",_);
		len[0] = -1;
		ans = (1<<30); al = 0; ar = (1<<30);
		last = cur = rt = cnt = 1;
		l = 0;
		int num; scanf("%d",&num);
		scanf("%s",T+1);
		for (int i=2;i<=num;i++) {
			scanf("%s",S+1); 
			int n = strlen(S+1) + 1; 
			S[n] = '{';
			for (int j=1;j<=n;j++) 
				add(j);
		}
		
		int n = strlen(T+1);
		for (int i=1;i<=n;i++) 
			gio(i);
		
		if (ans > n) 
			puts("Impossible");
		else {
			for (int i=al;i<=ar;i++) printf("%c",T[i]);
			printf("\n");
		}
		for (int i=1;i<=cnt;i++) 
			for (int j=1;j<=27;j++) ch[i][j] = 0;
		for (int i=1;i<=cnt;i++)
			link[i] = len[i] = 0;
	}
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值