HDU 3001 Travelling 三进制状压DP

题意：一个图，每个点都要遍历到，且每个点最多经过两次。

思路：用三进制表示状态，第i位上的数（0，1，2）代表第i+1个点经过的次数。

设dp[S][v]为状态S下，到达S状态中的点v的最小花费。

S0 为 状态S中第v-1位的值-1后的结果 ，其他位上的值不变。

那么状态转移方程：dp[S][v] = min(dp[S0][u] + maps[u][v],dp[S][v]).(u在状态S0中)

我的代码：

#include<cstdio>
#include<cstring>
#include<iostream>

#define min(a,b) (a < b ? a : b)
using namespace std;
typedef int LL;
const LL inf = 0x3f3f3f3f;
const LL maxn = 11;

LL n,m;
LL maps[maxn][maxn];
LL dp[177147][11];

LL mod_pow(LL n){
    LL res =  1,a = 3;
    while(n){
        if(n & 1) res *= a;
        a *= a;
        n >>= 1;
    }
    return res;
}

void init(){
    for(int i=0;i<=n;i++){
        for(int j=0;j<=n;j++){
            if(i == 0 || j == 0 || i == j) maps[i][j] = 0;
            else maps[i][j] = inf;
        }
    }
}

bool check(LL S){
    LL tmp;
    for(LL i = 0;i < n; i++){
        tmp = mod_pow(i);
        if((S / tmp) % 3 == 0) return false;
    }
    return true;
}

void solve(){
    LL S,S0,v,u,tmp;
    LL Ed = mod_pow(n);
    for(S = 0;S < Ed; S++){
        for(v = 1; v <= n; v++){
            dp[S][v] = inf;
            tmp = mod_pow(v - 1);
            if((S / tmp) % 3 == 0) continue;
            S0 = S - tmp;
            if(S0 == 0) dp[S][v] = 0;
            else{
                for(u = 1;u <= n; u++){
                    tmp = mod_pow(u - 1);
                    if(v == u || (S0 / tmp) % 3 == 0) continue;
                    dp[S][v] = min(dp[S][v],dp[S0][u] + maps[u][v]);
                }
            }
        }
    }
    LL res = inf;
    for(S = 0;S < Ed; S++){
        if(check(S)){
            for(v = 1; v <= n ; v++)
                res = min(res,dp[S][v]);
        }
    }
    if(res == inf) printf("-1\n");
    else  printf("%d\n",res);
}

int main(){
    while(~scanf("%d%d",&n,&m)){
        init();
        LL a,b,c;
        while(m--){
            scanf("%d%d%d",&a,&b,&c);
            maps[a][b] = maps[b][a] = min(maps[a][b],c);
        }
        solve();
    }
    return 0;
}