DP(状压专题四)

题意 :给定奶牛的数目, 给定一个熵, 在给定奶牛数组的的号码牌, 求最大无序组合的个数(无序指的是每两个奶牛的号码牌的差要大于给定的熵)

>> face <<

Strategy:状压DP, i是一个表示二进制的十进制数,'1’是第i位的牛被选了

状态: dp[j][i]在状态i中, 以j结尾的奶牛的最大无序数

目标:

$$\sum _{i=1}^{i\le n}dp[i][(1<<n)-1]$$

边界:

$$dp[i][1<<(i-1)]=1,i\in [1,2,\dots ,n](一只牛的无序队)$$

合法判断: 队尾必须包含在状态i中, 新队尾不在状态i中

转移方程: += 在原状态数组后加了一个确定的数,理应+=原状态的情况

$$dp[new][i|1<<new-1]+=dp[old][i]$$

attention:

双倍经验: 注意边界, 目标以及限制条件

#include <bits/stdc++.h>
#include <bits/extc++.h>
#define _rep(i, a, b) for (int i = (a); i <= (b); ++i)
#define _rev(i, a, b) for (int i = (a); i >= (b); --i)
#define _for(i, a, b) for (int i = (a); i < (b); ++i)
#define _rof(i, a, b) for (int i = (a); i > (b); --i)
#define ll long long
#define db double
#define eps 0.00001
#define all(x) x.begin(), x.end()
#define met(a, b) memset(a, b, sizeof(a))
#define what_is(x) cerr << #x << " is " << x << endl;
#define lowbit(x) x &(-x)
using namespace std;
const int maxn = 1 << 16;
ll dp[17][maxn]; //dp[i][j]表示以第i头牛结尾,选取状态为j的最大方案数
int n, k;
int a[maxn];
int main()
{
	cin >> n >> k;
	_rep(i, 1, n)
	{
		cin >> a[i];
	}
	_rep(i, 1, n)
	{
		dp[i][1 << i-1] = 1;
	}
	_rep(i, 0, (1 << n) - 1)
	{
		_rep(tail, 1, n)
		{ //枚举队尾
			if (i >> tail - 1 & 1)
			{
				_rep(newer, 1, n)
				{
					if (i >> newer - 1 & 1)
						continue;
					if (fabs(a[newer] - a[tail]) <= k)
						continue;
					dp[newer][i | 1 << newer - 1] += dp[tail][i];
				}
			}
		}
	}
	ll ans = 0;
	_rep(i, 1, n)
	{
		ans += dp[i][(1 << n) - 1];
	}
	cout << ans << endl;
}

回顾:: 不要设 变量名为 i, j, k (循环因子) , 推导时默认原数组里面有数字, 所以边界情况不能为dp[0][0] = 1, 然后就是 不开 long long 见祖宗