poj3229 哈密顿回路

最新推荐文章于 2019-02-20 20:49:43 发布

丿Smile灬晨星

最新推荐文章于 2019-02-20 20:49:43 发布

阅读量1k

点赞数

CC 4.0 BY-SA版权

分类专栏：动态规划

本文链接：https://blog.youkuaiyun.com/cxsmile/article/details/9327581

动态规划专栏收录该内容

76 篇文章

订阅专栏

本文介绍了一种通过路径 Floyd 处理结合状态压缩 DP 方法解决特定问题的技术。详细阐述了算法步骤，包括初始化、路径处理、状态转移等关键环节，并通过代码实现展示了实际应用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

很简单的题目，先将路径floyed处理，然后状态压缩DP即可。

ACcode：

#include<cstdio>
#include<cstring>

const int NS=15;
const int MS=1<<NS;
const double INF=1e12;

int n,m,x,y,f;
int ms,lim,res;
int dt[2]={80,120};
double tmp,k,z;
double time[NS];
double g[NS][NS];
double dp[MS][NS];

double Min(double a1,double b1)
{
    return a1<b1?a1:b1;
}

double Max(double a1,double b1)
{
    return a1>b1?a1:b1;
}

int getone(int st)
{
    int cnt=0;
    for (int i=st;i;i=i&(i-1)) cnt++;
    return cnt;
}

int main()
{
    while (~scanf("%d%d%lf",&n,&m,&k)&&(n+m+k))
    {
        ms=0,k*=12,res=-1,lim=1<<n;
        for (int i=0;i<m;i++)
        scanf("%d",&x),x--,ms|=1<<x;
        for (int i=0;i<n;i++)
        scanf("%lf",&time[i]);
        for (int i=0;i<n;i++)
        for (int j=0;j<n;j++)
        g[i][j]=i==j?0.0:INF;
        while (~scanf("%d%d%lf%d",&x,&y,&z,&f)&&(x+y+z+f))
        {
            x--,y--,z/=dt[f];
            g[y][x]=g[x][y]=Min(g[x][y],z);
        }

        for (int c=0;c<n;c++)
        for (int a=0;a<n;a++)
        for (int b=0;b<n;b++)
        g[a][b]=Min(g[a][b],g[a][c]+g[c][b]);

        for (int i=0;i<lim;i++)
        for (int j=0;j<n;j++)
        dp[i][j]=INF;
        for (int i=0;i<n;i++)
        dp[1<<i][i]=g[0][i]+time[i];

        for (int i=1;i<lim;i++)
        {
            for (int j=0;j<n;j++)
            {
                if (i&(1<<j))
                {
                    y=i^(1<<j);
                    for (int u=0;u<n;u++)
                    if (y&(1<<u))
                    dp[i][j]=Min(dp[i][j],dp[y][u]+g[u][j]+time[j]);
                }
//                printf("dp[%d][%d]=%.2lf\n",i,j,dp[i][j]);
            }
            if (((i&ms)==ms)&&(dp[i][0]<k*1.0)&&getone(i)>res)
            res=getone(i);
        }
        if (res<0) printf("No Solution\n");
        else printf("%d\n",res);
    }
    return 0;
}