2016夏季练习——dp

来源：HDU3001

三进制的状态压缩

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int st[12];
int vis[60000][12];
int dp[60000][12];
int mp[12][12];
int n,m;
int ans;
void ini()
{
	st[0]=1;
	for(int i=1;i<=10;i++)
		st[i]=st[i-1]*3;
		//3^n
	/*for(int i=0;i<10;i++)
		cout<<st[i]<<" ";*/
	for(int i=0;i<=st[10];i++){
		int tmp=i;
		for(int j=0;j<=10;j++){
			vis[i][j]=tmp%3;
			tmp/=3;
			/*cout<<i<<" "<<j<<" "<<vis[i][j]<<endl;
			getchar();*/
		}
		//求三进制
	}
}
bool flag;
int main(){
	ini();
	int d1,d2,s;
	while(scanf("%d%d",&n,&m)!=EOF){
		memset(dp,0x3f,sizeof(dp));
		memset(mp,0x3f,sizeof(mp));
		for(int i=0;i<n;i++)
			dp[st[i]][i]=0;
		//自己不能到自己
		for(int i=0;i<m;i++)
		{
			scanf("%d%d%d",&d1,&d2,&s);
			d1--;
			d2--;
			mp[d1][d2]=mp[d2][d1]=min(mp[d1][d2],s);
		}
		//记录路径
		ans=0x3f3f3f3f;
		for(int i=0;i<st[10];i++){
			flag=true;
			for(int j=0;j<n;j++){
				if(vis[i][j]==0) flag=false;
				if(dp[i][j]==0x3f) continue;
				for(int k=0;k<n;k++) {
					if(j==k) continue;
					if(vis[i][k]>=2) continue;
					if(mp[j][k]==0x3f) continue;
					dp[i+st[k]][k]=min(dp[i+st[k]][k],dp[i][j]+mp[j][k]);
				}
 			}
 			if(flag){
				for(int j=0;j<n;j++){
					ans=min(ans,dp[i][j]);
				}
 			}
		}
		if(ans==0x3f3f3f3f) ans=-1;
		cout<<ans<<endl;
	}
	return 0;
}