BZOJ1195 [HNOI2006]最短母串 【状压dp】

题目

给定n个字符串（S1,S2,„,Sn），要求找到一个最短的字符串T，使得这n个字符串（S1,S2,„,Sn）都是T的子串。

输入格式

第一行是一个正整数n（n<=12），表示给定的字符串的个数。
以下的n行，每行有一个全由大写字母组成的字符串。每个字符串的长度不超过50.

输出格式

只有一行，为找到的最短的字符串T。在保证最短的前提下，
如果有多个字符串都满足要求，那么必须输出按字典序排列的第一个。

输入样例

2

ABCD

BCDABC

输出样例

ABCDABC

题解

写完状压dp后才知道可以用AC自动机水，而且十分简洁。。。

我还是说状压dp吧。。
思想简单却十分难写，，，

我们先将其它串的子串去掉，保证两两完全不包含
设\(f[s][i]\)表示当前串集合为\(s\)，最后一个为\(i\)号串的最小长度，同时记录一个\(pre\)数组记录状态转移的方向
再预处理一个\(at[i][j]\)表示在\(i\)串后接\(j\)串增加的长度

转移十分显然，枚举下一个不在集合中的串，如果方案更优则转移，如果同样优，利用pre数组还原串比较字典序

所有的字符串比较之类的都可以直接暴力

虽然代码长了一点丑了一点，但是跑得飞快

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define ULL unsigned long long int
#define LL long long int
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
#define mp(a,b) make_pair<int,int>(a,b)
#define cp pair<int,int>
using namespace std;
const int maxn = 15,maxm = (1 << 13),maxl = 55,INF = 1000000000;
struct String{
	char s[maxl];
	int len;
	ULL h[maxl];
}S[maxn];
cp pre[maxm][maxn];
int n,at[maxn][maxn],f[maxm][maxn];
ULL P[1000];
void exclude(){
	for (int i = 1; i <= n; i++){
		for (int j = 1; j <= n; j++){
			if (j != i && S[j].len >= S[i].len){
				int flag = false,L = S[j].len - S[i].len + 1;
				for (int k = 1; k <= L; k++){
					if (S[i].h[S[i].len] == S[j].h[k + S[i].len - 1] - S[j].h[k - 1] * P[S[i].len]){
						flag = true;
						break;
					}
				}
				if (flag){
					swap(S[i--],S[n--]);
					break;
				}
			}
		}
	}
	//for (int i = 1; i <= n; i++)
	//	printf("%s\n",S[i].s + 1);
}
void init(){
	REP(i,n) REP(j,n) if (i != j){
		for (int k = max(1,S[i].len - S[j].len + 1); k <= S[i].len; k++){
			int L = S[i].len - k + 1;
			if (S[j].h[L] == S[i].h[S[i].len] - S[i].h[k - 1] * P[L]){
				at[i][j] = L;
				break;
			}
		}
	}
	//REP(i,n) printf("%d ",at[1][i]); puts("");
}
char t1[1000],t2[1000];
int s1[55],s2[55],top1,top2,n1,n2;
int cmp(int a,int b,int c,int d,int x){
	top1 = 0;
	while (a && b){
		s1[++top1] = b;
		cp tmp = pre[a][b];
		a = tmp.first;
		b = tmp.second;
	}
	n1 = 0;
	for (int i = top1,last = 1; i; i--){
		int u = s1[i];
		for (int j = last; j <= S[u].len; j++){
			t1[++n1] = S[u].s[j];
		}
		last = at[u][s1[i - 1]] + 1;
	}
	if (x){
		int l = at[s1[1]][x] + 1;
		for (int i = l; i <= S[x].len; i++)
			t1[++n1] = S[x].s[i];
	}
	top2 = 0;
	while (c && d){
		s2[++top2] = d;
		cp tmp = pre[c][d];
		c = tmp.first;
		d = tmp.second;
	}
	n2 = 0;
	for (int i = top2,last = 1; i; i--){
		int u = s2[i];
		for (int j = last; j <= S[u].len; j++){
			t2[++n2] = S[u].s[j];
		}
		last = at[u][s2[i - 1]] + 1;
	}
	if (x){
		int l = at[s2[1]][x] + 1;
		for (int i = l; i <= S[x].len; i++)
			t1[++n2] = S[x].s[i];
	}
	for (int i = 1; i <= n1; i++)
		if (t1[i] != t2[i]) return t2[i] - t1[i];
	return 0;
}
void print(int x){
	int a = (1 << n) - 1,b = x;
	top1 = 0;
	//puts("");
	while (a && b){
		//puts(S[b].s + 1);
		s1[++top1] = b;
		cp tmp = pre[a][b];
		a = tmp.first;
		b = tmp.second;
	}
	n1 = 0;
	for (int i = top1,last = 1; i; i--){
		int u = s1[i];
		for (int j = last; j <= S[u].len; j++){
			t1[++n1] = S[u].s[j];
		}
		last = at[u][s1[i - 1]] + 1;
	}
	t1[n1 + 1] = '\0';
	printf("%s\n",t1 + 1);
}
void solve(){
	fill(f[0],f[0] + maxm * maxn,INF);
	int maxv = (1 << n) - 1;
	for (int i = 1; i <= n; i++) f[1 << i - 1][i] = S[i].len;
	for (int s = 0; s <= maxv; s++){
		for (int i = 1; i <= n; i++){
			if (!(s & (1 << i - 1))) continue;
			for (int j = 1; j <= n; j++){
				if ((s | (1 << j - 1)) != s){
					int e = (s | (1 << j - 1)),len = S[j].len - at[i][j];
					if (f[e][j] > f[s][i] + len){
						f[e][j] = f[s][i] + len;
						pre[e][j] = mp(s,i);
					}
					else if (f[e][j] == f[s][i] + len){
						int tmp = cmp(s,i,pre[e][j].first,pre[e][j].second,j);
						if (tmp > 0) pre[e][j] = mp(s,i);
					}
				}
			}
		}
	}
	int ans = 1;
	for (int i = 2; i <= n; i++){
		if (f[maxv][i] < f[maxv][ans] || (f[maxv][i] == f[maxv][ans] && cmp(maxv,i,maxv,ans,0) > 0)){
			ans = i;
		}
	}
	print(ans);
}
int main(){
	P[0] = 1;
	for (int i = 1; i < 1000; i++) P[i] = P[i - 1] * 107;
	scanf("%d",&n);
	for (int i = 1; i <= n; i++){
		scanf("%s",S[i].s + 1);
		S[i].len = strlen(S[i].s + 1);
		for (int j = 1; j <= S[i].len; j++){
			S[i].h[j] = S[i].h[j - 1] * 107 + S[i].s[j] - 'A';
		}
	}
	exclude();
	init();
	solve();
	return 0;
}