这是一道入门的状态压缩dp,只要理解了dp是根据不包括这一节点前的集合最短长度进行选择,以及利用二进制码进行状态压缩,表示集合即可解出
状态压缩仅仅是集合,而F这个二维数组则蕴含着顺序,代码中的注释已经很详细了,有问题可以评论区问
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<map>
using namespace std;
int n;
double F[16][36000], l[16][16];
//F的行代表目前吃下的奶酪,列代表吃下这块奶酪时的可能集合,l是点与点之间的距离
struct coordinate
{
double x, y;
}arr[16];//用于保存各个奶酪的坐标
double length(coordinate a, coordinate b)//计算各奶酪之间的距离
{
return sqrt(pow( a.x - b.x,2) + pow(a.y - b.y,2));
}
int main()
{
cin >> n;
for (int i = 1; i <=n; i++)
cin >> arr[i].x >> arr[i].y;
memset(F, 127, sizeof(F));//将F的所有可能初始化为无限大
for(int i=0;i<=n;i++)//计算每个奶酪的距离
for (int j = i + 1; j <=n; j++)
{
l[i][j] = length(arr[i], arr[j]);
l[j][i] = l[i][j];//无向图没有方向之分
}
for (int i=1; i<=n; i++)//吃下的第一个奶酪是从初始点(0,0)到任一奶酪,借此初始化F
F[i][1 << (i-1)] = l[0][i];//这是吃下一块奶酪的情况
for (int i = 1; i <1 << n ; i++)//遍历每一种情况
{
for (int t = 1; t <=n; t++)//试图寻找刚被吃掉的奶酪
{
if (!(i &1 <<(t-1)))//1&x,只和x是0和1有关,0&x则都是0
continue;//没吃的奶酪不用考虑
for (int g = 1; g <=n; g++)//试图寻找上一个被吃掉的奶酪
{
if (t== g || !(i & 1 << (g-1)))
continue;
F[t][i] = min(F[t][i], F[g][i &~(1 << t-1)] + l[g][t]);//选择最短方案
}
}
}
double ans = 1 << 15 - 1;
for (int i = 1; i <= n; i++)
{
double temp = F[i][(1 << n) - 1];//吃下所有奶酪的情况是全1
if (ans >temp )
ans = F[i][(1 << n) - 1];
}
printf("%.2lf", ans);
return 0;
}