洛谷 P1433 吃奶酪状态压缩DP

最新推荐文章于 2022-08-11 16:02:23 发布

自信的小螺丝钉

最新推荐文章于 2022-08-11 16:02:23 发布

阅读量358

点赞数

分类专栏： PAT甲练习文章标签：动态规划算法

本文链接：https://blog.csdn.net/qq_45791939/article/details/122881240

版权

PAT甲练习专栏收录该内容

27 篇文章 0 订阅

订阅专栏

该博客介绍了如何使用动态规划方法解决编程竞赛题目洛谷P1433——吃奶酪。博主首先给出了原题链接，并提供了C++代码实现，代码中计算了二维平面上各个点之间的距离，并通过动态规划更新点到原点的最短路径。最终求得所有点可达的最短总距离。

摘要由CSDN通过智能技术生成

练习动态规划中

本题参考博客：

P1433 吃奶酪（状压dp）

洛谷普及组 P1433 吃奶酪

原题链接：洛谷 P1433 吃奶酪

#include <bits/stdc++.h>
using namespace std;
double x[16],y[16],d[16][16],dp[16][40000];

double dis(int a,int b)
{
    return sqrt( (x[a]-x[b])*(x[a]-x[b])+(y[a]-y[b])*(y[a]-y[b]) );
}


int main()
{
    std::ios::sync_with_stdio(false);
    int n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>x[i]>>y[i];
    }
    memset(dp,127,sizeof(dp));//初始化点在各种状态下的距离
    for(int i=1;i<=n;i++)
    {
        for(int j=i+1;j<=n;j++)
        {
           d[i][j]=d[j][i]=dis(i,j);
        }
    }
    //第i点到(0,0)的距离，状态为1<<(i-1)
    for(int i=1;i<=n;i++)
    {
        dp[i][1<<(i-1)]=dis(0,i);
    }
    for(int r=1;r<=(1<<n)-1;r++)//遍历所有状态
    {
        for(int now=1;now<=n;now++)
        {
            if((r&(1<<(now-1)))==0)
                continue;//跳过未访问的点
            for(int pre=1;pre<=n;pre++)
            {
                if((r&(1<<(pre-1)))==0 ||now==pre)
                    continue;//如果pre不在当前状态已经走过的点当中或者now和pre是同一个点，就跳过
                dp[now][r]=min(dp[now][r],dp[pre][r-(1<<(now-1))]+d[now][pre]);
            }
        }
    }
    double res=1e8;
    for(int i=1;i<=n;i++)
    {
        res=min(dp[i][(1<<n)-1],res);
    }
    cout<<setiosflags(ios::fixed)<<setprecision(2)<<res<<endl;

    return 0;
}