【深搜+状态压缩+动态规划(状压dp)】洛谷 P1433 吃奶酪详细解析及AC代码【三合一】-CSDN博客

本文链接：https://blog.csdn.net/flyunicorninsky/article/details/134585287

前言

最近了解了一些研究生的事情，很焦虑，三无产品的大学生怎么找到研究生读😭但是CPS认证还是要争取一下的，说不定以后转战leetcode

题目

题目描述

房间里放着 $n$ 块奶酪。一只小老鼠要把它们都吃掉，问至少要跑多少距离？老鼠一开始在 $(0, 0)$ 点处。

输入格式

第一行有一个整数，表示奶酪的数量 $n$ 。

第 $2$ 到第 $(n + 1)$ 行，每行两个实数，第 $(i + 1)$ 行的实数分别表示第 $i$ 块奶酪的横纵坐标 $x_i, y_i$ 。

输出格式

输出一行一个实数，表示要跑的最少距离，保留 $2$ 位小数。

样例 #1

样例输入 #1

样例输出 #1

7.41

提示

数据规模与约定

对于全部的测试点，保证 $1\leq n\leq 15$ ， $|x_i|, |y_i| \leq 200$ ，小数点后最多有 $3$ 位数字。

提示

对于两个点 $x_1,y_1)$ ， $x_2, y_2)$ ，两点之间的距离公式为 $\sqrt{(x_1-x_2)^2+(y_1-y_2)^2}$ 。

$2022.7.13$ ：新增加一组 $\text{Hack}$ 数据。

题目分析

一开始还是想用最朴素的dfs深搜来解决，谁让我是dfs的忠实粉丝，但是很遗憾，朴素的深搜会有一半的TLE。那么我就想了，该不会还要剪枝吧？
于是最朴素的我又使用了最朴素的剪枝，在dfs开头判断if(total>res)return;OK，再次提交，结果过了几个点，还是有四五个超时TLE。哎呦想破脑袋都想不到了还有什么地方可以剪枝了。
我觉得也许是不能使用dfs，但是俗话说得好，回溯搜索很简单，难在剪枝，特别是找哪里可以剪。
还是求助题解，第一次知道了状态压缩这种方法！
观察数据，15个奶酪！这数据并不大。适合使用状态压缩
根据题意，老鼠走的路径是无后效性的，只要经过的点一致，所在的点也一致，接下来所要走的路径就是等价的。换句话说，可以使用0和1来记录点是否走过，而当“前面经过的点固定”，那么并不会影响到后面的数值，剩下点的最短距离不会因为前面的顺序或者数值而改变，那么使用这种方法就是可行的。
相反如果说，老鼠还有一个体力值，体力值消耗有一个函数，你前面跑的距离虽然是前面的最佳的，但是你可能没有考虑全局，后面没有体力了，于是前面的决策选择影响力后面的结果，就是说是有后效性的，就不能使用这种方法，简单来说，只有局部最优不影响全局最优才能用，（局部最优是全局最优的一部分？）。
所以可以用一个二进制数来记录走过的点，另一个数记录老鼠所在的点，如果之后搜到的答案比这个点要大，就不继续搜即可，大致思想就是这样。
于是定义一个数组dp来记录走过固定点的最优距离。当同样的固定点和出发点不优于最短距离的将会被剪枝，OK！

知识点传送

状压dp，全称状态压缩dp，是一种在数据范围不大（n<=20左右）的情况下，代替暴搜的常用方法
状态常以二进制形式表示，从左往右第k位表示第k个元素的选择情况，0表示未被选择，1表示已选
例如二进制下110101，表示第1、3、5、6个元素已被选择
因为涉及二进制，所以简单的位运算必须掌握

注意事项

1.朴素深搜过不了
2.朴素深搜+朴素剪枝过不了
3.状态压缩的dp数组开大一点
4.位运算记得是i-1，因为第一个点是0初始点，记得加括号。

代码

失败了三次终于拿下

#include<iostream>
#include<cmath>
using namespace std;

int n;
double total=0,a,b,res=90000;
int visited [20]= {0};
double dp[2097153][20]={0}; 
struct point {
	double x=0,y=0;
} cheese[20];
double distance(point A,point B)
{
	return sqrt((A.x-B.x)*(A.x-B.x)+(A.y-B.y)*(A.y-B.y));
}
void dfs(int now,int left,int zhuangtai)
{
	if(total>res)//剪枝 
		return;
	if(left==0) { //没有奶酪了
		if(total<res)
			res=total;
		return;
	}
	for(int i = 1; i <= n ; i ++) {
		if(!visited[i]) {
			int zt=zhuangtai+(1<<(i-1));
			double dis=distance(cheese[now],cheese[i]) ;
			if(dp[zt][i]!=0&&dp[zt][i]<=dis+total)continue;
			visited[i]=1;//标记已访问
			total+=dis;//距离加
			dp[zt][i]=total;//保存状态最优解
			dfs(i,left-1,zt);//更新最新点，奶酪减一
			total-=dis;//恢复现场
			visited[i]=0;//恢复为未访问
		}
	}


}
int main()
{
	cin>>n;
	for(int i=1; i<=n; i++) {
		cin>>a>>b;
		cheese[i].x=a;
		cheese[i].y=b;
	}
	//visited[0]=1;
	dfs(0,n,0);
	printf("%.2lf",res);
	return 0;
}

后话

额外测试用例

因为WA了两次获得了两个个用例

样例输入 #2

样例输出 #2

17.15

样例输入 #3

样例输出 #3

21.73

纯状压动态规划解法

因为本题还有纯动态规划的方法，但是专题侧重讲搜索，这里就不多进行扩展，可以看一下状态压缩经典的NOIP2016愤怒的小鸟

王婆卖瓜

感觉有收获或者想跟上我的进度刷题的，可以点个关注，或者点赞收藏评论都可以！

题目来源

洛谷链接