5418
状压DP+Floyed算法求多源最短路径
参考:Floyed算法、C++中的位运算符
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;//不会溢出,大概大小是10^9多一点
const int MXS=(1<<16)+5;
const int MXN=20;
int n,m,g[MXN][MXN],dp[MXS][MXN],c[MXS];
int main(){
int t,flag=0,num,x,y,z;
for(int i=0;i<16;i++)c[1<<i]=i+1;//初始化
scanf("%d",&t);
while(t--){
memset(g,inf,sizeof(g));
memset(dp,inf,sizeof(dp));
scanf("%d%d",&n,&m);
num=(1<<n)-1;
for(int i=1;i<=n;i++)g[i][i]=0;
while(m--){
scanf("%d%d%d",&x,&y,&z);
g[x][y]=g[y][x]=min(g[x][y],z);
}
//Floyed算法
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
for(int k=1;k<=n;k++){
g[j][k]=min(g[j][k],g[j][i]+g[i][k]);
}
}
}
dp[1][1]=0;
for(int i=1;i<num;i++){
for(int last=c[-i&i],j=i;j>0;j&=j-1,last=c[-j&j]){//枚举集合元素
if(dp[i][last]==inf)continue;//未出现
int item,next;
for(int k=num^i;k>0;k&=k-1){//枚举集合外的元素求最小
next=c[-k&k],item=i|(-k&k);//|是按位或
dp[item][next]=min(dp[item][next],dp[i][last]+g[last][next]);//g[][]是两点间最短路径
}
}
}
for(int i=2;i<=n;i++){
if(dp[num][1]>dp[num][i]+g[1][i])dp[num][1]=dp[num][i]+g[1][i];
}
printf("%d\n",dp[num][1]);
}
return 0;
}