CodeForces 21 D.Traveling Graph（状压DP+Floyd）

Description
给出一张n个点m条边的无向图，每条边有边权，要求从1开始经过所有边至少一次后回到1处，问路径边权和最小值
Input
第一行两个整数n和m表示点数和边数，之后m行每行三个整数u,v,w表示u和v之间有一条边权为w的边(1<=n<=15,0<=m<=2000,1<=w<=1e4)
Output
输出从1开始经过所有边至少一次后回到1处的路径边权和最小值
Sample Input
3 3
1 2 1
2 3 1
3 1 1
Sample Output
3
Solution
如果是每条边只经过一次问题就变成求原图的欧拉回路，欧拉回路存在的重要条件是图连通且每个点度数为偶数，所以该问题可以转化为在原图的奇度数点（奇数点数量必为偶数）之间连边把奇度点变为偶度点使得加边后图的欧拉回路边权和最小，而为使路径边权和最小，加边的边权也必须尽可能小，对原图跑一边Floyd得到任意两点之间最短路dis[i][j]，之后状压DP，以dp[i]表示把i二进制中1所处位置对应的点的编号全部变成偶读点所需的最小代价，那么每次找到两个对应i状态中0位置的j,k，相当于找到两个奇度点然后连边，有转移方程
dp[ii]=min(dp[ii],dp[i]+dis[j][k])，其中ii=i^(1 << j)^(1 << k)
时间复杂度O(n*2^n)
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 16
int n,m,dp[1<<maxn],du[maxn],dis[maxn][maxn];
int Solve(int ans)
{
    for(int i=1;i<n;i++)//图不连通 
        if(du[i]&&dis[0][i]==INF)return -1;
    for(int i=0;i<n;i++)du[i]%=2;
    int N=1<<n,sta=0;
    for(int i=0;i<N;i++)dp[i]=INF;
    for(int i=0;i<n;i++)
        if(du[i]==0)sta^=(1<<i);
    dp[sta]=0;
    for(int k=sta;k<N;k++)
    {
        for(int i=0;i<n;i++)
            if(du[i]&&!(k&(1<<i)))
                for(int j=0;j<n;j++)
                    if(du[j]&&dis[i][j]!=INF&&!(k&(1<<j)))
                    {
                        int sta=k^(1<<i)^(1<<j);
                        dp[sta]=min(dp[sta],dp[k]+dis[i][j]);
                    }   
    }
    if(dp[N-1]==INF)return -1;
    return dp[N-1]+ans;
}
int main()
{
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=0;i<n;i++)
            for(int j=0;j<n;j++)
                dis[i][j]=INF;
        memset(du,0,sizeof(du));
        int ans=0;
        while(m--)
        {
            int u,v,w;
            scanf("%d%d%d",&u,&v,&w);
            u--,v--;
            du[u]++,du[v]++;
            ans+=w,dis[u][v]=dis[v][u]=min(dis[u][v],w);
        }
        for(int k=0;k<n;k++)
            for(int i=0;i<n;i++)
                for(int j=0;j<n;j++)
                    dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
        for(int i=0;i<n;i++)dis[i][i]=0;
        printf("%d\n",Solve(ans));
    }
    return 0;
}