Cod sursa(job #1256855)

Utilizator sebinechitasebi nechita sebinechita Data 6 noiembrie 2014 22:30:48
Problema ADN Scor 30
Compilator cpp Status done
Runda Arhiva de probleme Marime 2.34 kb
#include <iostream>
#include <fstream>
#include <cstring>
#include <climits>
using namespace std;
ifstream fin("adn.in");
ofstream fout("adn.out");
#define MAX 30005
#define MAXN 18
#define INF INT_MAX/2;
int a[1<<MAXN][MAXN];
int p[1<<MAXN][MAXN];
char v[MAXN][MAX];
int pi[MAX], C[MAXN][MAXN];
int dr, st[MAXN+2];


int kmp(char a[], char b[])
{
    int n=strlen(a+1);
    int m=strlen(b+1);
    pi[1]=0;
    int i, k;
    for(i=2;i<=n;i++)
    {
        k=pi[i-1];
        while(k!=0 && a[k+1]!=a[i])
            k=pi[k];
        if(a[k+1]==a[i])
        {
            k++;
        }
        pi[i]=k;
    }
    k=0;
    for(i=1;i<=m;i++)
    {
        while(k!=0 && a[k+1]!=b[i])
        {
            k=pi[k];
        }
        if(a[k+1]==b[i])
        {
            k++;
        }
        if(k==n)
        {
            return 0;
        }
    }
    return n-k;
}

int main()
{
    int n, i, j, k, val, ind, nod, m;
    fin>>n;
    for(i=0;i<n;i++)
    {
        fin>>(v[i]+1);
    }
    for(i=0;i<n;i++)
    {
        for(j=0;j<n;j++)
        {
            if(i==j)
                continue;
            C[i][j]=kmp(v[i], v[j]);
        }
    }
    for(i=1;i<(1<<n);i++)
    {
        for(j=0;j<n;j++)
            p[i][j]=-1;
        if((i&(i-1))==0)
            continue;
        for(j=0;j<n;j++)
        {
            if((i&(1<<j))==0)
                continue;
            a[i][j]=INF;
            for(k=0;k<n;k++)
            {
                if((i&(1<<k))==0)
                    continue;
                if(k==j)
                    continue;
                if(C[k][j]+a[i-(1<<j)][k]<a[i][j])
                {
                    a[i][j]=C[k][j]+a[i-(1<<j)][k];
                    p[i][j]=k;
                }
            }
        }
    }

    val=INF;
    for(i=0;i<n;i++)
    {
        if(a[(1<<n)-1][i]<val)
        {
            val=a[(1<<n)-1][i];
            ind = i;
        }
    }

    st[++dr]=ind;
    nod=(1<<n)-1;
    while(ind!=-1)
    {
        st[++dr]=p[nod][ind];
        nod=(nod)-(1<<ind);
        ind=st[dr];
    }
    for(i=1;i<=n;i++)
        cout << st[i] << ' ';
    fout << v[st[1]]+1;
    for(i=2;i<=n;i++)
    {
        m=strlen(v[st[i]]+1);
        for(j=m-C[st[i]][st[i-1]]+1;j<=m;j++)
            fout << v[st[i]][j];
    }
}