状压DP hdu 3001 Travelling

【题意】10个点的TSP问题，要求经过每个点不超过两次。

【解析】用三进制状态压缩（表示经过某个点几次），几乎与用二进制完全相同（参考poj3311Hie With The Pie），只是不能使用位运算。可以先预处理每个状态的第k位是什么，也可以直接根据( s / (3^k) )%3来判断第k位是什么（s表示状态）。

【状态表示】dp[s][i]表示经过状态s后到达i点需要的最小花费。

【状态转移】dp[s][i] = min{ dp[s`][j] + dis[j][i] }，dis[j][i]为j到i的花费。

【注意】去重边。记得处理dp[3^k][k]为0。注意在何处能取得最终答案。

#include "iostream"
#include "cstdio"
#include "cmath"
#include "cstring"
using namespace std;
const int sub = 0x3f3f3f3f;

int n, m;
int mp[11][11];
int dp[59049][11];
int po[11];

int main()
{
    int i, j, a, b, c;
    po[0]=1;
    for(i = 1; i < 11; i++) po[i] = po[i-1]*3;
    while(~scanf("%d%d",&n,&m))
    {
        memset(mp,-1,sizeof(mp));
        for(int i = 0; i < m; i++)
        {
            scanf("%d%d%d",&a,&b,&c);
            if(mp[a][b]==-1) mp[a][b]=mp[b][a]=c;
            else mp[a][b]=mp[b][a]=min(mp[a][b],c);
        }
        if(n==0) {printf("0\n"); continue;}
        int ans = sub;
        for(int s = 0; s < po[n]; s++)
        {
            int flag = 1;
            for(i = 1; i <= n; i++)
            {
                dp[s][i] = sub;
                if((s/po[i-1])%3 == 0) { flag=0; continue;}
                if(s == po[i-1]) { dp[s][i]=0; continue;}
                for(j = 1; j <= n; j++)
                {
                    if(i==j || mp[j][i]==-1 || (s/po[j-1])%3==0) continue;
                    dp[s][i] = min(dp[s][i], dp[s-po[i-1]][j]+mp[j][i]);
                }
            }
            if(flag)
            {
                for(int i = 1; i <= n; i++)
                    ans = min(ans, dp[s][i]);
            }
        }
        if(ans == sub) printf("-1\n");
        else printf("%d\n",ans);
    }
    return 0;
}