poj2288(状态压缩dp)

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <iostream>
#include <algorithm>
using namespace std;
long long v[15],G[15][15];
long long dp[(1<<14)][15][15],path[(1<<14)][15][15];
int main(){                                     //dp[s][i][j]表示状态是s，最后两个点是i和j
long long t,n,m,i,j,k,x,y,s,s_tmp,ans,sum,tmp;
scanf("%I64d",&t);
while(t--){
scanf("%I64d%I64d",&n,&m);
for(i=0;i<n;i++)
scanf("%I64d",&v[i]);
memset(G,0,sizeof(G));
memset(dp,-1,sizeof(dp));
memset(path,0,sizeof(path));
for(i=0;i<m;i++){
scanf("%I64d%I64d",&x,&y);
x--,y--;
G[x][y]=G[y][x]=1;
}
if(n==1){                               //一个点的时候要进行特判
printf("%d 1\n",v[0]);
continue;
}
for(i=0;i<n;i++){
for(j=0;j<n;j++){
if(i!=j&&G[i][j]){
dp[(1<<i)+(1<<j)][i][j]=v[i]+v[j]+v[i]*v[j];
path[(1<<i)+(1<<j)][i][j]=1;
}                               //先预处理出所有相邻的情况
}
}
for(s=0;s<(1<<n);s++){
for(i=0;i<n;i++){
if(s&(1<<i))
for(j=0;j<n;j++){
if((s&(1<<j))&&i!=j&&G[i][j])
for(k=0;k<n;k++){           //找出一个三个点都在的状态进行转移
if((s&(1<<k))&&i!=k&&j!=k&&G[j][k]){
s_tmp=s-(1<<i);
if(dp[s_tmp][j][k]!=-1){
tmp=v[i]+v[i]*v[j]+dp[s_tmp][j][k];
if(G[i][k])     //判断是否成环
tmp+=(v[i]*v[j]*v[k]);
if(tmp>dp[s][i][j]){
dp[s][i][j]=tmp;
path[s][i][j]=path[s_tmp][j][k];
}
else if(tmp==dp[s][i][j])
path[s][i][j]+=path[s_tmp][j][k];
}
}
}
}
}
}
ans=-1,sum=0;
for(i=0;i<n;i++){
for(j=0;j<n;j++){
if(i==j)
continue;                       //遍历最后两个点判断最后的最大值
if(dp[(1<<n)-1][i][j]>ans){
ans=dp[(1<<n)-1][i][j];
sum=path[(1<<n)-1][i][j];
}
else if(dp[(1<<n)-1][i][j]==ans)
sum+=path[(1<<n)-1][i][j];
}
}
if(ans==-1)
puts("0 0");
else
printf("%I64d %I64d\n",ans,sum/2);
}
return 0;
}