题意: 给你一个双向连通图,求 获得权值最大 的 哈密顿通路的 权值 和 这个权值对应的数目;
其中权值计算方法是 列如 ABCD 权值是a+b+c+d+ab+bc+cd 如果 A,B,C 和B,C,D 可构成三角形分别加上abc,bcd;
这个题 和poj 3311 很相像: 那个需要记录一个最后到达的地方 这个需要记录俩个罢了
DP[i][a][b]其中 i 二进制 中1表示这个点走过了 最后走的的 的是b>>a
因为对于已经走过了{1,2,3,4,,5,6,..,N} 为了推出下一次走到X 我们要获得权值 只和最后走的俩个点有关?
给出代码::
1 #include <cstdio> 2 #include <algorithm> 3 #include <cmath> 4 #include <cstdlib> 5 #include <cstring> 6 #include <iostream> 7 using namespace std; 8 long long dp[1<<13][13][13]; 9 long long num[1<<13][13][13]; 10 int n; 11 int val[13]; 12 bool Map[13][13]; 13 void solve() 14 { 15 for(int i=0;i<(1<<n);i++)for(int a=0;a<n;a++)for(int b=0;b<n;b++) 16 { 17 if(!Map[a][b])continue; 18 if(a==b)continue; 19 if((i&(1<<a))==0||(i&(1<<b))==0)continue; 20 if(i==((1<<a)+(1<<b))) 21 { 22 dp[i][a][b]=val[a]+val[b]+val[a]*val[b]; 23 num[i][a][b]=1; 24 continue; 25 } 26 else 27 { 28 for(int c=0;c<n;c++) 29 { 30 if(c==a||c==b||(!Map[b][c])) continue; 31 if((i&(1<<c))==0) continue; 32 long long tmp; 33 if(dp[i^(1<<a)][b][c]==-1)continue; 34 if(Map[a][c]) tmp=val[a]+val[a]*val[b]+val[a]*val[b]*val[c]; 35 else tmp=val[a]+val[a]*val[b]; 36 if(dp[i][a][b]<dp[i^(1<<a)][b][c]+tmp) 37 { 38 num[i][a][b]=num[i^(1<<a)][b][c]; 39 dp[i][a][b]=dp[i^(1<<a)][b][c]+tmp; 40 } 41 else if(dp[i][a][b]==dp[i^(1<<a)][b][c]+tmp) 42 num[i][a][b]+=num[i^(1<<a)][b][c]; 43 } 44 } 45 } 46 long long ans=-1; 47 long long ans_num=0; 48 for(int b=0;b<n;b++) for(int c=0;c<n;c++) 49 { 50 if(ans<dp[(1<<n)-1][b][c]) 51 { 52 ans=dp[(1<<n)-1][b][c]; 53 ans_num=num[(1<<n)-1][b][c]; 54 } 55 else if(ans==dp[(1<<n)-1][b][c]) 56 ans_num+=num[(1<<n)-1][b][c]; 57 58 } 59 if(ans==-1) printf("0 0\n"); 60 else printf("%lld %lld\n",ans,ans_num/2); 61 } 62 int main() 63 { 64 int q; 65 scanf("%d",&q); 66 while(q--) 67 { 68 memset(Map,false,sizeof(Map)); 69 memset(dp,-1,sizeof(dp)); 70 memset(num,0,sizeof(num)); 71 int m; 72 scanf("%d%d",&n,&m); 73 for(int i=0;i<n;i++)scanf("%d",&val[i]); 74 for(int i=1;i<=m;i++) 75 { 76 int a,b; 77 scanf("%d%d",&a,&b); 78 a--,b--; 79 Map[a][b]=Map[b][a]=true; 80 } 81 if(n==1) printf("%d 1\n",val[0]); 82 else 83 solve(); 84 } 85 return 0; 86 }
poj 2288 Islands and Bridges,布布扣,bubuko.com
原文地址:http://www.cnblogs.com/shuly/p/3856248.html