hdu 3001 Travelling 状压dp

题意:给一张城市地图,可以以任何城市为起点,求走完所有城市的最小费用,所有城市最多走过两遍。
思路:因为我们要控制经过城市次数,所以用三进制表示,注意有重边。
dp[i][j]代表i状态时我们在j城市的最小费用
当我们要从j城市走到k城市的时候,要判断k城市的访问次数,然后就很容易转移
dp[i+ck[k]][k]=min(dp[i][j]+Map[j][k],dp[i+ck[k]][k]);
就最后求答案的时候要找到满足走完城市的状态的时候比较麻烦

#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include<stack>
#include<string>
#include<vector>
#include<map>
#include<set>
using namespace std;
typedef long long LL;
#define maxn 200000
#define f(x) (x*1.0)
#define inf 0x3f3f3f3f
#define maxm maxn*maxn
#define lowbit(x) (x&(-x))
#define cheak(i) printf("%d ",i)
#define lson(x) (splay[x].son[0])
#define rson(x) (splay[x].son[1])
#define rfor(i,a,b) for(i=a;i<=b;++i)
#define lfor(i,a,b) for(i=a;i>=b;--i)
#define mem(a,b) memset(a,b,sizeof(a))
#define mec(a,b) memcpy(a,b,sizeof(b))
int dp[maxn][15];//i状态 当前在j位置 第k次访问 
int Map[15][15],num[maxn][11],ck[15];
void init()
{
    int i;
    rfor(i,1,maxn-1)
    {
        int x=i,len=-1;
        while(x)
        {
            num[i][++len]=x%3;
            x/=3;
        }
    }
    ck[1]=1;
    rfor(i,2,12) ck[i]=ck[i-1]*3;
}
int main()
{
    int i,j,k,n,m,u,v,w;
    init();
    while(~scanf("%d%d",&n,&m))
    {
        rfor(i,1,(ck[n+1])-1)
         rfor(j,1,n) dp[i][j]=inf;
        rfor(i,1,n)
         rfor(j,1,n) Map[i][j]=inf;
        //rfor(i,1,n) printf("*%d ",ck[i]);printf("\n");
        rfor(i,1,n) dp[ck[i]][i]=0;
        rfor(i,1,m)
        {
            scanf("%d%d%d",&u,&v,&w);
            if(w<Map[u][v])
            Map[u][v]=Map[v][u]=w;
        }
        int ans=inf;
        rfor(i,0,(ck[n+1])-1)
        {
            int x=i,t=0,flag=0;
            while(x)
            {
                if(x%3) t++;
                else break;
                x/=3;
            }
            if(t==n) flag=1;
            rfor(j,1,n)
            {
                if((i/ck[j])%3)
                {
                    rfor(k,1,n)
                    {
                        if((i/ck[k])%3>=2) continue;
                        dp[i+ck[k]][k]=min(dp[i][j]+Map[j][k],dp[i+ck[k]][k]);
                    }
                }
                if(flag) ans=min(dp[i][j],ans);
            }
        }
        if(ans<inf)
        printf("%d\n",ans);
        else printf("-1\n");
    }
    return 0;
}

你可能感兴趣的:(poj,状压dp)