HDU 3001 【三进制状压DP】

最新推荐文章于 2024-08-05 20:06:07 发布

原创最新推荐文章于 2024-08-05 20:06:07 发布 · 517 阅读

0 ·

CC 4.0 BY-SA版权

DP。专栏收录该内容

49 篇文章

订阅专栏

本文介绍了一种结合了状压DP技术与三进制的应用案例，针对特定问题提出了有效的解决方法，并详细解释了如何通过三进制来处理可以重复访问节点的问题。

= =。。。我果真是数学弱渣。。。

普通状压，但是由于可以走两回所以变成了三进制的。。然后注意三进制的处理方法。。

关于状压的循环顺序，如果是那种给你格子的，就是一行一行遍历的，每一行都遍历所有可能符合情况的数，通常是用滚动数组加上部分维度或者是变量进行限制条件的记录。如果是需要顺序的，只要求最后遍历完的，就是在1<<n范围内，每个数都遍历一次n来找寻合适的下一状态。

#include <stdio.h>
#include <iostream>
#include <string.h>
using namespace std;
#define maxn 60000
int n,m;
int a[11][11];
int thr[11],ll[11],rr[11];
int dp[11][maxn];
int inner[maxn];
void init()
{
    thr[0]=1;ll[0]=0;rr[0]=0;
    memset(inner,1,sizeof(inner));
    for(int i=1;i<11;i++)
    {
        thr[i]=thr[i-1]*3;
        rr[i]=thr[i]-1;
        ll[i]=ll[i-1]+thr[i-1];
    }
    for(int i=0;i<maxn;i++)
    {
        int j=i;
        while(j)
        {
            if(j%3==0) {inner[i]=0;break;}
            j/=3;
        }
    }
}
int min(int x,int y){return x<y?x:y;}
int getmid(int pre,int num)
{
    int t=pre/thr[num];
    return t%3;
}
void search()
{
    int maxnum=rr[n];
    int ans=-1;
    memset(dp,-1,sizeof(dp));
    for(int i=0;i<=maxnum;i++)
    {
        for(int j=0;j<n;j++)
        {
            if(getmid(i,j)==2) continue;
            int now=i+thr[j];
            if(i==0)
            {
                dp[j][now]=0;
            }
            else
            {
                for(int k=0;k<n;k++) // pre node
                {
                    //if(k==j) continue;
                    if(dp[k][i]==-1||a[k][j]==-1) continue;
                    //printf("dp[%d][%d]=%d d=%d ",k,i,dp[k][i],a[k][j]);
                    if(dp[j][now]==-1) dp[j][now]=dp[k][i]+a[k][j];
                    else dp[j][now]=min(dp[j][now],dp[k][i]+a[k][j]);
                    
                }
            }
            if(now<=rr[n]&&now>=ll[n]&&inner[now]&&dp[j][now]!=-1)
            {
                //printf("dp[%d][%d]=%d\n",j,now,dp[j][now]);
                //for(int k=0;k<n;k++) printf("dp[%d][%d]=%d\n",k,i,dp[k][i]);
                if(ans==-1) ans=dp[j][now];
                else ans=min(ans,dp[j][now]);
                //cout<<endl;
            }
        }
    }
    //cout<<ll[n]<<" "<<rr[n]<<endl;
    cout<<ans<<endl;
}
int main()
{
    init();
    while(cin>>n>>m)
    {
        int x,y,z;
        memset(a,-1,sizeof(a));
        for(int i=0;i<n;i++) a[i][i]=0;
        for(int i=0;i<m;i++)
        {
            cin>>x>>y>>z;
            x--,y--;
            if(a[x][y]!=-1) a[x][y]=a[y][x]=min(a[x][y],z);
            else a[x][y]=a[y][x]=z; // consider the repeat road
        }
        //for(int i=0;i<n;i++) {for(int j=0;j<n;j++) cout<<a[i][j]<<" ";cout<<endl;}
        search();
    }
}

今天晚上去看SVG~~~~~~~ = w = ~~~~~