状态dp——吃奶酪

题目

原题传送门

分析

题目大致就是现在给定了n个点的坐标，从原点出发，求遍历这些点所需的最短距离。是旅行商问题。
状态dp能够解决最大20的问题。这个题目最大是15，所以用状态dp。
其中状态k（001101）表示的是共6个点，其中点1，3，4被走过了，在这个路径中。

1. 判断点i有没有在路径中：k&(1<<(i-1))==1表示i在路径k中
2. 要从路径k中删除点i，回到点i之前的点j：
   dp[j][k-(1<<(i-1))]

大致框架

//一大堆初始化：
//把dp中的每一个值都设置为无穷大。因为是double类型的所以可以使用memset(dp,127,sizeof(dp));
//距离数组：确定点i和点j的距离（i-j和j-i一样）
//dp数组：从原点直接到i点的距离设置好d[0][i];
//进行枚举：三层循环
//所有的状态：也就是路径k，从1-(1<<n)-1
//所有的点：从1-n，需要判断当前的点i是否在路径k中了
//在路径k中能够到达当前点的点：从1-n
//找到dp[i][(1<<n)-1]中的最小值

代码

#include<bits/stdc++.h>
using namespace std;
double dp[18][(1<<18)];//dp数组，dp[i][k]表示的是经过路径k（状态集合k）到达点i的所用的最短距离
double d[20][20];//距离数组，d[i][j]表示的是点i到点j的距离
double x[20],y[20];//x[i],y[i]表示的是点i的横纵坐标
int n;//点的个数
double res;//走遍所有点所需要的最短距离，从dp[i][(1<<n)-1]中找。i从1-n，i是终点


int main()
{
    ios::sync_with_stdio(false);
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>x[i]>>y[i];
    }
    //先进行一大堆初始化
    //先把dp数组里面的每一个值都弄成无穷大
    memset(dp,127,sizeof(dp));
    //距离数组：确定点i和点j的距离
    res=dp[0][0];
    x[0]=0;
    y[0]=0;
    for(int i=0;i<n;i++){
        for(int j=i+1;j<=n;j++){
            d[i][j]=sqrt(pow((x[i]-x[j]),2)+pow((y[i]-y[j]),2));
            d[j][i]=d[i][j];
        }
    }
    //dp数组：从原点直接到i点的距离设置好d[0][i];

    for(int i=1;i<=n;i++){
        dp[i][(1<<(i-1))]=d[0][i];
    }
    //进行枚举：三层循环
    //所有的状态：也就是路径k，从1-(1<<n)-1
    //所有的点：从1-n，需要判断当前的点i是否在路径k中了
    //在路径k中能够到达当前点的点：从1-n
    for(int k=1;k<(1<<n);k++){
        for(int i=1;i<=n;i++){
            if((k&(1<<(i-1)))==0)
                continue;
            for(int j=1;j<=n;j++){
                if(j==i)
                    continue;
                if((k&(1<<(j-1)))==0)
                    continue;
                dp[i][k]=min(dp[i][k],dp[j][k-(1<<(i-1))]+d[i][j]);
            }
        }
    }
    //找到dp[i][(1<<n)-1]中的最小值
    
    for(int i=1;i<=n;i++){
        res=min(res,dp[i][(1<<n)-1]);
    }
    printf("%.2f\n",res);
    return 0;
}