POJ-1699 Best Sequence 状态压缩DP
题目链接:http://poj.org/problem?id=1699
太爽了这题,1AC。容易想到用状态压缩DP来做,f[k][i][j]表示当前 i 状态有 k 个串并且串以 j 结尾的最短串。则 f[k][i][rt]=Min{ f[k][i][rt] , f[k][st][r]+len[j]+g[r][j] }。其中状态 i 用位运算表示所包含的具体的串,1010表示包含串2和4。这里要注意转移方程中的rt,要考虑串包含的情况。
1 //STATUS:C++_AC_0MS_324KB 2 #include<stdio.h> 3 #include<stdlib.h> 4 #include<string.h> 5 #include<math.h> 6 #include<iostream> 7 #include<string> 8 #include<algorithm> 9 #include<vector> 10 #include<queue> 11 #include<stack> 12 #include<map> 13 using namespace std; 14 #define LL __int64 15 #define pii pair<int,int> 16 #define Max(a,b) ((a)>(b)?(a):(b)) 17 #define Min(a,b) ((a)<(b)?(a):(b)) 18 #define mem(a,b) memset(a,b,sizeof(a)) 19 #define lson l,mid,rt<<1 20 #define rson mid+1,r,rt<<1|1 21 const int N=11,INF=0x3f3f3f3f,MOD=10000,STA=8000010; 22 const double DNF=1e13; 23 24 char s[N][22]; 25 int f[2][1<<N][11],len[N],g[N][N],next[2][1<<N],cnt[2],vis[1<<N]; 26 int T,n; 27 28 void getg() 29 { 30 int i,j,k,p; 31 for(i=0;i<n;i++){ 32 for(j=0;j<n;j++){ 33 for(k=0;k<=len[i];k++){ 34 for(p=k;p<len[i] && p-k<len[j];p++) 35 if(s[i][p]!=s[j][p-k])break; 36 if(p==len[i] || p-k==len[j]){ 37 g[i][j]=Min(len[i]-k,len[j]); 38 break; 39 } 40 } 41 } 42 } 43 } 44 45 int main() 46 { 47 // freopen("in.txt","r",stdin); 48 int i,j,p,st1,k,r,st2,ans,rt; 49 scanf("%d",&T); 50 while(T--) 51 { 52 scanf("%d",&n); 53 for(i=0;i<n;i++){ 54 scanf("%s",s[i]); 55 len[i]=strlen(s[i]); 56 } 57 getg(); 58 59 mem(cnt,0); 60 mem(vis,0); 61 for(i=0;i<n;i++){ 62 next[0][cnt[0]]=1<<i; 63 f[0][1<<i][i]=len[i]; 64 cnt[0]++; 65 } 66 mem(f[p=1],INF); 67 for(i=1;i<n;i++){ 68 for(k=0;k<cnt[!p];k++){ 69 st1=next[!p][k]; 70 for(r=0;r<n;r++){ 71 if(st1&(1<<r)){ 72 for(j=0;j<n;j++){ 73 if(st1&(1<<j))continue; 74 st2=st1|(1<<j); 75 if(!vis[st2]){ 76 vis[st2]=1; 77 next[p][cnt[p]++]=st2; 78 } 79 if(g[r][j]==len[j])rt=r; 80 else rt=j; 81 f[p][st2][rt]=Min(f[p][st2][rt],f[!p][st1][r]+len[j]-g[r][j]); 82 } 83 } 84 } 85 } 86 mem(vis,0); 87 cnt[p=!p]=0; 88 mem(f[p],INF); 89 } 90 91 ans=INF; 92 k=(1<<n)-1; 93 for(i=0;i<n;i++) 94 ans=Min(ans,f[!p][k][i]); 95 printf("%d\n",ans); 96 } 97 return 0; 98 }