pku 1699 Best Sequence 状态压缩dp

http://poj.org/problem?id=1699

DFS+剪枝解法http://www.cnblogs.com/E-star/archive/2012/08/10/2631584.html

题意：

现在给出几个基因片段，要求你将它们排列成一个最短的序列，序列中使用了所有的基因片段，而且不能翻转基因。，这些基因可以重叠，只要一个基因的后段和一个基因的前端一样，就可以将其重叠链接在一起。现问将这些 基因全部排列出来，最短的长度为多少。

思路：

将n个基因片段的所有状态压缩，dp[i][j]表示状态i以j基因片段结尾的最短长度，则有

dp[i][j] = min(dp[i][j],dp[tmp][k] + len[j] - share[k][j])

tmp表示i状态除去j状态后，以k结尾的最短长度，然后dp[tmp][k] + len[j] - share[k][j]表示连接上后的长度。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <vector>
#define N 12
#define maxn 22
using namespace std;

const int inf = 0x7fffffff;

char str[N][maxn];
int dp[1<<N][maxn];
int len[N],share[N][N],n;

int getR(int x,int y)
{
    int i,j,l;
    int L = 0;
    for (l = 0; l <= len[x] && l <= len[y]; ++l)
    {
        bool flag = false;
        for (i = len[x] - l,j = 0; i < len[x] && j < l; ++i,++j)
        {
            if (str[x][i] != str[y][j])
            {
                flag = true;
                break;
            }
        }
        if (!flag) L = l;
    }
    return L;
}
void init()
{
    int i,j;
    memset(share,0,sizeof(share));
    for (i = 0; i < n; ++i)
    {
        for (j = 0; j < n; ++j)
        {
          share[i][j] = getR(i,j);
        }
    }
}
int main()
{
   // freopen("d.txt","r",stdin);
   int i,j,t,k;
   scanf("%d",&t);
   while (t--)
   {
       scanf("%d",&n);
       for (i = 0; i < n; ++i)
       {
           scanf("%s",str[i]);
           len[i] = strlen(str[i]);
       }
       init();//计算共享长度
       int nn = (1<<n) - 1;//存储所有的状态
       for (i = 1; i <= nn; ++i)
       {
           for (j = 0; j < n; ++j)
           {
               dp[i][j] = inf;//初始化
               if (i&(1<<j))//如果j存在于i状态里面
               {
                   int tmp = i - (1<<j);
                   if (tmp == 0)//如果只存在j状态
                   {
                       dp[i][j] = len[j];
                       continue;
                   }
                   for (k = 0; k < n; ++k)
                   {
                       if (tmp&(1<<k))
                       {
                           dp[i][j] = min(dp[i][j],dp[tmp][k] + len[j] - share[k][j]);
                       }
                   }
               }
           }
       }
       int Min = inf;
       for (i = 0; i < n; ++i)
       {
           if (dp[nn][i] < Min)
           {
               Min = dp[nn][i];
           }
       }
       printf("%d\n",Min);
   }
    return 0;
}