hdu 3001 Travelling 三进制状压Dp

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3001

附上一个写的 比较容易理解的链接：http://blog.sina.com.cn/s/blog_6ec19c780100ydws.html

题意：给定几个城市间的旅游费用，求把每个城市至少遍历一遍，且至多两遍的情况下的至少花费。否则输出-1。

代码写出来之后，会发现和二进制状压Dp差不多，只是多了一个附加条件，每个点可以重复一次，所以得用三进制状态表示，由于计算机中存储方式不是三进制，所以我们需要预处理出每种状态的三进制表示形式，便于枚举枚举每个位置的状态。这就是主要与二进制不同的地方，二进制是直接利用计算机位运算获取状态，相对来说，少了这一过程。

代码：

<span style="font-size:18px;">#include<iostream>
#include<sstream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<set>
#include<map>
#include<queue>
#include<vector>
#include<string>
#define LL long long
#define INF 0x7fffffff

using namespace std;

int n,m,Map[12][12];
int bit[11] = {1,3,9,27,81,243,729,2187,6561,19683,59049};
int dp[60000][10],num[60000][11];

void init(){//获取每种状态的三进制表示形式
    for(int i=0;i<bit[10];i++){
        int t=i;
        for(int j=0;j<10;j++)
            num[i][j]=t%3,t/=3;
    }
}

void Dp(){
    memset(dp,-1,sizeof(dp));//初始化为-1
    for(int i=0;i<n;i++)//把每个起点的费用初始化为0
        dp[bit[i]][i]=0;
    int ans=-1;
    for(int i=0;i<bit[n];i++){//枚举每一种状态
        int op=1;//所有点都至少遍历过一次的标志
        for(int j=0;j<n;j++){//和二进制一样，枚举每个可以作为起点（即跳转）的点
            if(num[i][j]==0) op=0;
            if(dp[i][j]==-1 || num[i][j]==0) continue;//不能跳转
            for(int k=0;k<n;k++){//枚举可以走的下个点（跳转目标点）
                if(Map[j+1][k+1]==-1 || k==j || num[i][k]>=2) continue;//下一个点不合法
                int s=i+bit[k];
                if(dp[s][k]==-1) dp[s][k]=dp[i][j]+Map[j+1][k+1];
                else dp[s][k]=min(dp[s][k],dp[i][j]+Map[j+1][k+1]);
            }
        }
        if(op){//如果每个点至少遍历过一次，更新最小费用
            for(int j=0;j<n;j++)
                if(dp[i][j]!=-1)
                    ans=(ans==-1)? dp[i][j]:min(ans,dp[i][j]);
        }
    }
    cout<<ans<<endl;
}


int main(){
    //freopen("D:\\in.txt","r",stdin);
    init();
    while(cin>>n){
        cin>>m;
        int a,b,c;
        memset(Map,-1,sizeof(Map));
        for(int i=0;i<m;i++){
            cin>>a>>b>>c;
            if(Map[a][b]==-1) Map[a][b]=Map[b][a]=c;
            else Map[a][b]=Map[b][a]=min(Map[a][b],c);
        }
        Dp();
    }
    return 0;
}
</span>