NKOI 1905 慢跑小路

【S2状态压缩】慢跑小路
Time Limit:10000MS  Memory Limit:65536K

Total Submit:31 Accepted:16 
Case Time Limit:1000MS
Description
Gord准备参加马拉松赛跑，他在家的后面的一个公园进行练习。公园里面有很多设置有座椅和饮用水的休息点，庞大的慢跑小路构成的道路网将这些休息点连接了起来。Gord想要找出最短的一条慢跑小路路径，要求这条路径将每条小路都至少经过了一次。
Input
输入有多组数据，对于每组数据： 

第一行，两个整数n和m，分别表示休息点的数量和小路的数量。 

接下来m行，每行包含三个整数x,y,z,代表一条小路的信息，在休息点x和休息点y之间有条长度为z的慢跑小路。 

注意：两个休息点之间有可能有多条慢跑小路直接相连。Gord可以选择任意休息点作为起点，但是路径的终点和起点必须相同。 

单独的一行一个数字0表示输入结束。 

Output
对于每组测试数据，输出一行，每行一个整数，表示最短路径的长度。
Sample Input

4 5
1 2 3
2 3 4
3 4 5
1 4 10
1 3 12
0

Sample Output

41

Hint
n <= 15 

m < 1000 

z<=10000
Source
Waterloo local 2002

首先我们很容易知道这个图是一个无向图，那么在所求的回路内所有点的度数之和一定为一个偶数

我们分析可以发现，图中的欧拉回路的最短长度为所有边的长度+回到出发点的多走的长度，而这多走的边的端点一定是度数为奇的点

那么我们先求出来所有点的度数，对于一个点，如果他的度数为偶数，我们就可以不用讨论

我们求出所有度数为奇的点的集合s，对于s的子集i，我们设f[i]表示将i集合中的点两两配对所需的最小代价，那么得到

方程：

f[i|(1<<j)|(1<<k)]=f[i]+map[j][k]    j,k∈s且j,k∉i

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
const int maxn=20,inf=1e9;
int n,m,ans,s;
int map[maxn][maxn],f[1<<16],d[maxn];
inline void _read(int &x){
char t=getchar();bool sign=true;
while(t<'0'||t>'9')
{if(t=='-')sign=false;t=getchar();}
for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
if(!sign)x=-x;
}
bool ok(int x){
for(int i=0;i<n;i++)if(!(s&(1<<i))&&(x&(1<<i)))return 0;
return 1;
}
void clear(){
ans=0;
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
map[i][j]= i==j?0:inf;
memset(d,0,sizeof(d));
}
int main(){
int i,j,a,b,c,k;
while(scanf("%d",&n)&&n){
_read(m);
clear();
for(i=1;i<=m;i++){
_read(a);_read(b);_read(c);
ans+=c;
a--,b--;//将点从0开始编号
d[a]++,d[b]++;
if(map[a][b]>c)map[a][b]=map[b][a]=c;
}
s=0;
for(i=0;i<n;i++)
if(d[i]&1)s|=(1<<i);
for(i=1;i<=s;i++)f[i]=inf;
for(k=0;k<n;k++)
for(i=0;i<n;i++)
for(j=0;j<n;j++)
map[i][j]=min(map[i][k]+map[k][j],map[i][j]);//Floyd求出最短路
for(i=0;i<=s;i++)
if(ok(i))
for(j=0;j<n;j++)
if((s&(1<<j))&&(!(i&(1<<j))))
for(k=j+1;k<n;k++)
if((s&(1<<k))&&(!(i&(1<<k))))
f[i|(1<<j)|(1<<k)]=min(f[i|(1<<j)|(1<<k)],f[i]+map[j][k]);
printf("%d\n",ans+f[s]);
}
}