Cod sursa(job #235827)

Utilizator andrei-alphaAndrei-Bogdan Antonescu andrei-alpha Data 25 decembrie 2008 23:01:06
Problema ADN Scor 20
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.84 kb
using namespace std;

#include <bitset>
#define II inline
#define FOR(i,a,b) for(int i=a;i<=b;++i)

#define IN  "adn.in"
#define OUT "adn.out"
#define N_MAX 1<<5
#define L_MAX 1<<15
#define oo 1<<30

int N,K;
char V[N_MAX][L_MAX];
int P[N_MAX][N_MAX],L[N_MAX];
int A[1<<18][19];
int pi[L_MAX];

II void print(int x,int last);
II void scan()
{
	freopen(IN,"r",stdin);
	freopen(OUT,"w",stdout);
	scanf("%d\n",&N);
	FOR(i,1,N)
	{
		fgets(V[i],1<<15,stdin);
		for(int j=0;;L[i] = j,++j)
			if(V[i][j] == '\n' || V[i][j] == '\0')
			{
				V[i][j] = '\0';
				break;
			}	
	}	
	K = (1<<N)-1;
}

II void prefix(char N[],int n)
{
	int k = 0;
	pi[0] = 0;
	
	FOR(i,1,n)
	{
		while(k && N[k] != N[i])
			k = pi[k-1];
		if(N[k] == N[i])
			++k;
		pi[i] = k;
	}	
}

II int kmp(char N[],char M[],int n,int m)
{
	int q = 0;
	
	FOR(i,0,m)
	{
		while(q && N[q] != M[i])
			q = pi[q-1];
		if(N[q] == M[i])
			++q;
		if(q==n)
			return oo;
	}	
	return q;
}

II void solve()
{
	memset(A,100,sizeof(A));
	A[0][0] = 0;
	
	FOR(i,1,N)
	FOR(j,1,N)
	if(i!=j)
	{
		prefix(V[j],L[j]);
		P[i][j] = kmp(V[j],V[i],L[j],L[i]);
		if(P[i][j] == oo)
			K ^= (1<<j)-1;
	}	
	
	FOR(i,1,K)
	FOR(j,1,N)
	{
		if(!(K & 1<<(j-1)) || !(1<<(j-1) & i) )
			continue;
		FOR(k,0,N)
		{
			if(k==j)
				continue;
			A[i][j] = min(A[i][j],A[i ^ 1<<(j-1)][k] + L[j] - P[k][j]);
		}
	}	
	
	int poz(0);
	A[K][0] = oo;
	FOR(i,1,N)
		if(A[K][i] < A[K][poz])
			poz = i;
	print(K,poz);	
}

II void print(int x,int last)
{
	if(x == 1<<(last-1))
	{
		printf("%s",V[last]);
		return;
	}	
	for(int fiu = 1;fiu <= N;++fiu)
		if(A[x ^ 1<<(last-1) ][fiu] + L[last] - P[fiu][last] == A[x][last])
		{
			print(x ^ 1<<(last-1),fiu);
			printf("%s",V[last]+P[fiu][last]);
		}		
}

int main()
{
	scan();
	solve();
	return 0;
}