状态压缩DP题单

已于 2024-04-19 22:44:15 修改
阅读量408 收藏 8
点赞数 3
文章标签： 算法 c++ 数据结构
于 2024-04-18 22:59:30 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_62701031/article/details/137937407
版权
文章主要介绍了几种与IT技术相关的编程题目，涉及动态规划（如dp[i][s]表示状态转移）、最短路径算法（如Floyd-Warshall）、以及使用位操作优化的子集遍历。这些题目集中在计算最短路线、子集计数和状态转移等方面。
摘要由CSDN通过智能技术生成

技巧

for(int s = mask; s > 0; s = (s - 1) & mask) //遍历子集（mask将若干位上1->0的所有数
__builtin_popcount(s);//c++统计数位为1的和的函数

 

​​​​​​​P1433 吃奶酪（最短路）

dp(i, s) 表示从 i 出发经过的点的记录为 s 的路线距离最小值

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 20;
signed main()
{	
	int n; cin >> n;
	vector<double>x(n + 1), y(n + 1);
	vector<vector<double>> dp(n + 1, vector<double>(1 << (n + 1), 2e9));
	x[0] = 0, y[0] = 0;

	auto dis = [&](int i, int j) -> double
	{
		return sqrt((x[i] - x[j]) * (x[i] - x[j]) + (y[i] - y[j]) * (y[i] - y[j]));
	};
	for(int i = 1; i <= n; i ++)
	{
		cin >> x[i] >> y[i];
//		dp[i][1 << (i - 1)] = 0;
	}
	dp[0][1] = 0;
	for(int s = 1; s < (1 << (n + 1)); s ++)
	{
		for(int i = 0; i <= n; i ++)
		{
			if((s & (1 << (i))) == 0) continue;
			for(int j = 0; j <= n; j ++)
			{
				if(i == j) continue;
				if((s & (1 << (j))) == 0) {continue;}
//				cout << s << " " << i << " " << j << '\n';
				dp[i][s] = min(dp[i][s], dp[j][s - (1 << (i))] + dis(i, j)); 
			}
		}
	}
    double ans = -1;
    for(int i = 1; i <= n; i ++)
    {
        double s = dp[i][(1 << (n + 1)) -1];
        if(ans==-1||ans>s) ans=s;
    }
	cout << fixed << setprecision(2) << ans << '\n';
}

P8733 [蓝桥杯 2020 国 C] 补给（floyd）

 dp(i, s) 表示到达 i ，经过的点的记录为 s 的路线距离最小值

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 20;
signed main()
{	
	int n, D; cin >> n >> D;
	vector<double>x(n + 1), y(n + 1);
	vector<vector<double>> dp(n + 1, vector<double>((1 << n) + 10, 2e9)), f(n + 1, vector<double>(n + 1, 2e9));
	auto dis = [&](int i, int j) -> double
	{
		return sqrt((x[i] - x[j]) * (x[i] - x[j]) + (y[i] - y[j]) * (y[i] - y[j]));
	};
	for(int i = 0; i < n; i ++) cin >> x[i] >> y[i];
	
	for(int i = 0; i < n; i ++)
		for(int j = 0; j < n; j ++)
		{
			double now = dis(i, j);
			if(now <= D) f[i][j] = now;
		}
	
	for(int k = 0; k < n; k ++)
		for(int i = 0; i < n; i ++)
			for(int j = 0; j < n; j ++)
				f[i][j] = min(f[i][k] + f[k][j], f[i][j]);

	dp[0][1] = 0;
	for(int s = 1; s < (1 << n); s ++)
	{
		for(int i = 0; i < n; i ++)
		{
			if((s & (1 << i)) == 0) continue;
			for(int j = 0; j < n; j ++)
			{	
				if(i == j) continue;
				if((s & (1 << j)) == 0) continue;
				dp[i][s] = min(dp[i][s], dp[j][s - (1 << i)] + f[i][j]);
//				cout << i << " " << j << " " << f[i][j] << '\n';
			}
		}
	}
	double ans = 1e9;
	for(int i = 0; i < n; i ++)
	{
//		cout << i << " " << (1 << (n - 1)) << " " << dp[i][1 << (n - 1)] << '\n';
		double now = f[0][i] + dp[i][(1 << n) - 1];
		ans = min(ans, now);
	}
	cout << fixed << setprecision(2) << ans << '\n';
}

P7859 [COCI2015-2016#2] GEPPETTO（简单枚举状态）

 

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 4e5 + 10;
typedef long long ll;
int a[25][25]; 
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	cout.tie(0);
	int n, m; cin >> n >> m;
	for(int i = 1; i <= m; i ++)
	{
		int x, y; cin >> x >> y;
		a[x][y] = a[y][x] = 1;
	}
	int ans = 1;
	for(int s = 1; s < (1 << n); s ++)
	{
		int f = 0;
		for(int i = 1; i <= n; i ++)
		{
			if((s & (1 << (i - 1))) == 0) continue; 
			for(int j = i + 1; j <= n; j ++)
			{
				if((s & (1 << (j - 1))) == 0) continue;
				if(a[i][j] || a[j][i]) 
				{
					f = 1; 
//					cout << s << " " << i << " " << (s & (1 << (i - 1))) << " " << j << " "<< (s & (1 << (j - 1))) << '\n'; 
					break;
				}
			}
			if(f) break;
		}
//		cout << s << ' ' << f << '\n';
		if(!f) ans ++;
	}
	cout << ans << '\n';
}

P8687 [蓝桥杯 2019 省 A] 糖果

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 20;
/*
1 3 3
1 2 3
*/
signed main()
{	
	int n, m, k; cin >> n >> m >> k;
	vector<int>dp(1 << 21, -1), a(105);
	for(int i = 1; i <= n; i ++) 
	{
		int s = 0;
		for(int j = 1; j <= k; j ++)
		{
			int x; cin >> x;
			x --;
			s |= (1 << x);
		}
//		s |= (1 << x), s |= (1 << y), s |= (1 << z);
//		int s = (1 << x) + (1 << y) + (1 << z);
		a[i] = s;
	}
	dp[0] = 0;
	for(int s = 0; s < (1 << m); s ++)
	{
		if(dp[s] == -1) continue;
		for(int i = 1; i <= n; i ++)
		{
			if(dp[s | a[i]] == -1 || dp[s] + 1 < dp[s | a[i]])
			{
				dp[s | a[i]] = dp[s] + 1;
//				cout << s << " " << i << " " << a[i] << ' ' << dp[s | a[i]] << '\n';
			}
		}
	}
	cout << dp[(1 << m) - 1] << '\n';
}

[ABC332E] Lucky bag（SOS DP）

dp(s, i) 表示状态为 s 装入 i 个背包的最小值。

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 25; 
signed main()
{	
	ios::sync_with_stdio(0);
	cin.tie(0); cout.tie(0);
	int n, d; cin >> n >> d;
	vector<double>a(n + 1);
	vector<vector<double>>dp(1 << n, vector<double>(d + 1, 1e18));
	double sum = 0;
	for(int i = 0; i < n; i ++)
	{
		cin >> a[i];
//		cout << a[i] << ' ';
		sum += a[i];
	}
	double ave = sum / d;
//	cout << ave << '\n';
	for(int s = 0; s < (1 << n); s ++)
	{
		double res = 0;
		for(int j = 0; j < n; j ++) if((1 << j) & s) res += a[j];
		dp[s][1] = (res - ave) * (res - ave);
	}
	for(int i = 2; i <= d; i ++)
	{
		for(int s = 0; s < (1 << n); s ++)
		{
			for(int t = s; t > 0; t = (t - 1) & s)
			{
				dp[s][i] = min(dp[s][i], dp[t][i - 1] + dp[s ^ t][1]);
//				cout << s << " " << i << '\n';
			}
		}		
	}
	cout << fixed << setprecision(15) << dp[(1 << n) - 1][d] / d << '\n';
}

彩色路径（最短路/边数限制+折半存储）

CCF-CSP认证考试 202312-5 彩色路径 20/50/100分题解_ccf彩色路径-CSDN博客​​​​​​​

​​​​​​​E - Compatible Numbers（SOS DP）

对于数 a[ i ] 只要在 （（1 << 22）- 1）^ a[ i ] 中的子集找到即可

dp[ x ] 表示 x 的 子集中满足条件的数

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 25; 
signed main()
{	
	ios::sync_with_stdio(0);
	cin.tie(0); cout.tie(0);
	int n; cin >> n;
	vector<int>a(n + 1), dp(1 << 22, -1);
	for(int i = 1; i <= n; i ++) cin >> a[i], dp[a[i]] = a[i];
	for(int s = 0; s < (1 << 22); s ++)
	{
		for(int i = 0; i < 22; i ++)
		{
			if((s & (1 << i)) && dp[s ^ (1 << i)] != -1) dp[s] = dp[s ^ (1 << i)];
		}
	}
	for(int i = 1; i <= n; i ++)
	{
		cout << dp[((1 << 22) - 1) ^ a[i]] << " ";
	}
	cout << '\n';
}

P1896 [SCOI2005] 互不侵犯

 dp[i][j][s]:我们已经放了i行,用了j个国王,第i行的状态为s的方案数​​​​​​​

#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef pair<int, int> PII;
const int N = ((1 << 10) + 10), mod = 998244353;
int dp[10][100][1 << 9];
signed main() 
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	int n, k; cin >> n >> k;
	vector<int>ok, cnt(1 << 9);
	for(int s = 0; s < (1 << n); s ++)
	{
		if(s & ((s << 1) | (s >> 1))) continue;
		ok.push_back(s);
		int temp = __builtin_popcount(s);
		cnt[s] = temp;
//		cout << s << " " << temp << '\n';
	}
	dp[0][0][0] = 1;
	for(int i = 1; i <= n; i ++)
	{
		for(auto s1 : ok)
		{
			for(auto s2 : ok)
			{
				if(s1 & (s2 | (s2 << 1) | (s2 >> 1))) continue;
				for(int j = 0; j <= k; j ++)
				{
					if(j - cnt[s1] >= 0) dp[i][j][s1] += dp[i - 1][j - cnt[s1]][s2];
				}	
			}
		}
	}
	int ans = 0;
	for(auto s : ok) ans += dp[n][k][s];
	cout << ans << '\n';
}

不要麻酱
关注
九、状态压缩 DP
yake1965
03-27 553
维护一个集合的状态，常常用 int 或 bool 类型的数组进行表达与描述。如果这个集合包含的元素个数很少，就把描述集合状态的数组压缩到一个整数中，这种 用一个整数来描述一个集合的方法就叫作“状态压缩”。而当状态的某个维度存储的是一个集合的状态，这类问题就称作：状态压缩动态规划。
算法：状态压缩dp
m0_63997099的博客
05-01 1572
状态压缩dp实际上比其他dp要容易想一点，因为状态很容易定义，并且使用状压那么数据就必须很小比如int32位，longlong64位，数据大小最多那么大。是一种在解决某些特定类型的问题时非常有效的技术，特别是在涉及到限制状态或配置数量较多，但可以用整数或位操作来紧凑表示的情况。这种方法通常用于解决优化问题，特别是在图论、搜索问题、排列组合问题中非常常见，如旅行商问题（TSP）、位掩码问题等。状态压缩DP的核心思想是使用整数的二进制表示来表示复杂的状态集合。每一个比特位可以代表一个元素的某种状态（通常是存在或
状态压缩DP
qq_54783066的博客
11-12 852
状压dp
蓝桥杯：国际象棋【状态压缩dp】
weixin_62953519的博客
04-05 457
写了有五六道状态压缩dp的题，本菜发现这种题目其实和线性dp差不多，只不过是多了一个状态压缩，即将某个状态压缩到一串二进制中； 先看题： 因为行数很小，所以不难想到是状态压缩dp； 这道题和洛谷的炮兵阵地很像，都是一行的状态可以牵扯到上两个状态，所以我们在对i行进行状态转移方程的时候，我们得考虑第i-1和i-2行的状态，然后这道题目从某种意义上来说，比炮兵阵地要难一些，因为这个还有一个马的个数限制，所以我们得开4个维度的dp数组； 我们先明确一下，状态集合的定义，即dp数组的定义： f[.
状态压缩dp模型整理
m0_74120503的博客
01-22 964
状态压缩dp的核心在于将每一层视为一个状态，这些状态用二进制来表示，每一层由上一层转移而来，那么这么来看的话，我们只用预处理本层的每一个状态可以由上一层的哪些状态转移而来，然后由此来更新本层每个状态即可。这种一般由于状态不能选满，所以看似时间复杂度很高，但实际循环跑不满。
状态压缩dp(状压dp)
热门推荐
唐宋元明清的博客
07-16 3万+
注：在涉及到位运算时，一定要注意位运算的优先级。该加的括号一定要加 状压dp是一类比较难理解的dp; 在讲状压dp之前，我们应该清楚所有的dp是解决多阶段决策最优化问题的一种思想方法； 请注意多阶段这三个字： 经过前面三种背包的学习，可以发现如何定义状态是解决动态规划最重要的一步； 状态的定义也就决定了相当于阶段的划分； 在背包问题中，我们通过物品的件数i和背包的容量j来定义状态或者说...
动态规划：状态压缩DP - 愤怒的小鸟
07-09
### 动态规划：状态压缩DP - 愤怒的小鸟 #### 题目介绍 本题源自AcWing平台编号为524的题目《愤怒的小鸟》。题目要求寻找最少数量的抛物线，使得每条抛物线都能覆盖给定的一组点中的若干个点。具体而言，给定一组...
【状压DP】易懂讲解状态压缩/状态压缩DP
qq_35684989的博客
05-31 3276
本文是状态压缩DP的总结
cf 1316E(贪心优化的状态压缩dp) 好题
01-03
题意： 给你n个人，从中选出p个球员和k个观众，第i个人作为观众产生价值ai，第i个人作为j号球员产生价值Ci,j ，求最大价值 (2≤n≤105,1≤p≤7,1≤k,p+k≤n,ai<109,ci,j<109) 输入 第1行输入n p k ...
【算法】状态压缩dp
KbS_gQop的博客
08-14 217
精通状态压缩dp
CodeForces - 1102F 状态压缩dp
m0_62167326的博客
08-29 130
状态转移
16 合并区间
weixin_47521269的博客
11-12 750
【代码】16 合并区间。
数据库范式、MySQL 架构、算法与树的深入解析
lrh35的博客
11-11 1150
在数据库设计中，范式是一系列规则，用于确保数据的组织和存储具有良好的结构、完整性以及最小化的数据冗余。如果不遵循范式设计，数据可能会以平铺式罗列，仅使用冒号、分号等简单分隔。这种方式存在诸多弊端，例如在数据搜索时极为不便，数据修改也会变得十分棘手。由此，范式的概念应运而生。
点云配准之点到点，点到面，点到线ICP，NDT算法介绍
qq_54375651的博客
11-10 1525
点云配准在人工智能，自动驾驶，SLAM等任务中具有十分关键的作用，本文将介绍点云配准的点到点，点到线，点到面ICP和NDT算法
1-Equity-Transformer：求解NP-Hard Min-Max路由问题的顺序生成算法（AAAI-24)(完)（code）
太极生两鱼
11-12 159
最小最大路由问题旨在通过智能体合作完成任务来最小化多个智能体中最长行程的长度。这些问题包括对现实世界有重大影响的应用场景，但已知属于NP-hard问题。现有方法在大规模问题上面临挑战，尤其是在需要协调大量智能体以覆盖数千个城市的情况下。本文提出了Equity-Transformer来解决大规模最小最大路由问题。首先，我们采用顺序规划方法来处理最小最大路由问题，使我们能够利用强大的序列生成器（例如，Transformer）。其次，我们提出了关键的归纳偏差，以确保智能体之间的公平工作量分配。
LeetCode题解：3.无重复字符的最长子串【python 题解超详细】
m0_74420622的博客
11-10 785
LeetCode题解：3.无重复字符的最长子串【python 题解超详细】，滑动窗口与暴力求解，知识拓展：什么是哈希函数
每日小题--买股票的最佳时机
最新发布
m0_74977981的博客
11-14 488
毋庸置疑，一定需要数组遍历结束后才能确定这么一个买入点和卖出点，所以此时想要提升算法的话就要找到一个只需要遍历一遍就能出结果的方法，即将时间复杂度控制在O(n)；遍历数组的过程中去找数组中是否存在一个低价点，并且在低价点的基础上去判断有无上升空间；即杜绝了{4，3，2，1}的情况出现，买入后无法卖出。在第 2 天（股票价格 = 1）的时候买入，在第 5 天（股票价格 = 6）的时候卖出，最大利润 = 6-1 = 5。后出价最高的一天卖出，简而言之，就是找到一个比起始小，并且结尾在一个有最大差距的位置；
以太坊系地址衍生算法分层确定性生成逻辑
J_Lee
11-12 249
以太坊系地址衍生算法分层确定性生成逻辑
动态规划---算法中的坚持仙蛊
2301_77174919的博客
11-14 432
然后再写代码，如果代码没通过就打印dp数组，看看是不是和自己预先推导的哪里不一样。如果打印出来和自己预先模拟推导是一样的，那么就是自己的递归公式、初始化或者遍历顺序有问题了。如果和自己预先模拟推导的不一样，那么就是代码实现细节有问题。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值