poj 1795 DNA Laboratory

题目：http://poj.org/problem?id=1795

求一长度最小，字典序最小的字符串，使他包含所有输入的字符串。

状态dp，首先可以把本来就包含在其他字符串里面的字符串去掉。预处理出dis[i][j],dis[i][j]表示在i字符串的后面最少要加多少个字符，才能覆盖j字符串。

然后用dp[state][i]表示以字符串i开头的状态为state的最小长度。求出最小长度后，再每次选字典序最小的来凑答案。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<string>
#include<algorithm>
#include<map>
using namespace std;

#define rep(i,s,t) for(int i=s;i<t;i++)
#define repp(i,s,t) for(int i=s-1;i>=t;i--)
#define clr(a) memset(a,0,sizeof(a));

const int INF=(int)(1e9);
int dis[15][15],dp[1<<15][15],fail[15][110],len[15];
string s[15],in[15];
int t,n,m,all;
string ans;
bool vis[15];

inline void getFail(string s,int *p,int l){
	int j=-1;
	p[0]=-1;
	rep(i,1,l){
		while(j>=0 && s[i]!=s[j+1]) j=p[j];
		if(s[i]==s[j+1]) j++;
		p[i]=j;
	}
}

inline int getDis(string a,string b,int *p,int l1,int l2){
	int j=-1;
	rep(i,0,l1){
		while(j>=0 && a[i]!=b[j+1]) j=p[j];
		if(a[i]==b[j+1]) j++;
		if(j==l2-1) return 0;
	}
	return l2-(j+1);
}

inline void dfs(int num,int now,string tmp,int st){
	//if(tmp>ans) return;
	if(num==n){
		if(ans>tmp) ans=tmp;
	}
	else{
		int id=-1;
		string ttmp="zzZZ";
		rep(i,0,n){
			if(!vis[i]){
				if(dp[st][now] == dp[st^(1<<now)][i]+len[now]-(len[i]-dis[now][i])){
					string a=tmp;
					rep(j,len[i]-dis[now][i],len[i]){
						a+=s[i][j];
					}
					if(id==-1 || ttmp>a) ttmp=a,id=i;
				}
			}
		}
		if(id!=-1){
			vis[id]=1;
			dfs(num+1,id,ttmp,st^(1<<now));
		}
	}
} 
int main(){
	cin>>t;
	rep(ca,1,t+1){
		cin>>n;
		rep(i,0,n){
			cin>>in[i];
		}
		m=n;n=0;
		clr(vis);
		rep(i,0,m){
			int tag=1;
			rep(j,0,m){
				if(i!=j && !vis[j]){
					if(in[j].find(in[i])!=in[j].npos){
						tag=0;break;
					}
				}
			}
			if(tag) s[n++]=in[i];
			else vis[i]=1;
		}
		all=(1<<n)-1;
		rep(i,0,all+1) rep(j,0,n) dp[i][j]=INF;
		rep(i,0,n){
			len[i]=s[i].size();
			getFail(s[i],fail[i],len[i]);
			dp[1<<i][i]=len[i];
		}
		rep(i,0,n) rep(j,0,n){
			if(i!=j) dis[i][j]=getDis(s[i],s[j],fail[j],len[i],len[j]);
		}
		rep(i,0,all+1){
			rep(j,0,n){
				if((i&(1<<j)) && dp[i][j]!=INF){
					rep(k,0,n){
						if(!(i&(1<<k))){
							int nxt=(i|(1<<k));
							dp[nxt][k]=min(dp[nxt][k],dp[i][j]+(len[k]-(len[j]-dis[k][j])));	
						}
					}
				}
			}
		}
		int Min=dp[all][0],id=0;
		rep(i,1,n){
			if(dp[all][i]<Min || (dp[all][i]==Min && s[id]>s[i])){
				Min=dp[all][i];id=i;
			}
		}
		//cout<<Min <<" " <<id <<" "<<n <<endl;
		ans="zzzZZZ";
		clr(vis);vis[id]=1;
		dfs(1,id,s[id],all);
		cout<<"Scenario #"<<ca<<":"<<endl;
		cout<<ans<<endl;
		cout<<endl;
	}
	return 0;
}