DP(状压专题Ⅶ)

题意: 给定一个数字串s和正整数d,统计有多少不同的排列能被d整除(可以有前导零)

>> face <<

Strategy:状压DP, i是一个表示二进制的十进制数,'1’该位置对应的数被选上

状态: dp[i][k], 代表状态i, 被除余数为k的方案数;

目标: dp[(1 << n) - 1][0] -> 全选且被除为0

边界: dp[0][0] = 1 ;

合法判断: 每加入一个新数的时候, 判断这数是不是已经在状态i里面了, 并且判断该数有没有被用过, 如果被用过记得去重

转移方程:

$$dp[new][k\ast 10+新数]+=dp[old][k],k\in [0,mod-1]$$

attention: 去重

双倍经验: 状态很新颖,

#include <bits/stdc++.h>
#include <bits/extc++.h>
#define oo 0x3f3f3f3f
#define ll long long
#define db double
#define all(a) a.begin(), a.end()
#define met(a, b) memset(a, b, sizeof(a))
#define what_is(x) cerr << #x << " is " << x << endl;
#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 lowbit(x) x &(-x)
#define pi acos(-1.0)
using namespace std;
using namespace __gnu_pbds;
const int maxn = 20;
int mod;
const db eps = 1e-8;
int dp[1 << 10][1000], a[11];
bool used[10];
int main()
{
	int t;
	cin >> t;
	while (t--)
	{
		string s;
		cin >> s >> mod;
		int n = s.size();
		_rep(i, 1, n){
			a[i] = s[i-1] - '0';
		}
		met(dp, 0);
		dp[0][0] = 1;
		_rep(i, 0,(1 << n) - 1){

			met(used, 0);
			_rep(j, 1, n){

				if(i >> j - 1 & 1 || used[a[j]])
					continue;
				used[a[j]] = 1;
				_for(k, 0, mod){
					dp[i | 1 << j - 1][(k * 10 + a[j]) % mod] += dp[i][k];
				}
			}
		}
		cout << dp[(1 << n) - 1][0] << endl; 
	}
}

第一次回顾

感悟: 还是枚举新状态让旧状态推过来比较好, 思路和上面一样的,只不过(不要忘记清零dp数组, 然后就是注意_rep 和 _for, 不要拿着就敲, 还有就是used数组清零的位置每种状态对应一个used数组,所以在枚举新状态的时候就要清零)

#include <bits/stdc++.h>
#include <bits/extc++.h>
#define oo 0x3f3f3f3f
#define ll long long
#define db double
#define all(a) a.begin(), a.end()
#define met(a, b) memset(a, b, sizeof(a))
#define what_is(x) cout << #x << " is " << x << endl
#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 lowbit(x) x &(-x)
#define bin(x) cout << #x << " is " << bitset<sizeof(int) * 2>(x) << endl
#define pi acos(-1.0)
using namespace std;
const int maxn = 1002;
const int mod = 1e8;
int dp[1 << 10][maxn], a[maxn];
bool used[10];
int main()
{
	ios::sync_with_stdio(0);
	int t;
	cin >> t;
	while (t--)
	{
		string tmp;
		int d;
		cin >> tmp >> d;
		met(a, 0);
		
		met(dp, 0);
		dp[0][0] = 1;
		int n = tmp.size();
		_for(i, 0, n) a[i + 1] = tmp[i] - '0';
		_rep(i, 1, (1 << n) - 1)
		{
			met(used, 0);
			_rep(j, 1, n)
			{
				if (!(i >> j - 1 & 1) || used[a[j]])
					continue;
				_for(k, 0, d)
				{
					dp[i][(k*10+a[j])%d] += dp[i ^ 1 << j - 1][k];
				}
				used[a[j]] = 1;
			}
		}
		cout << dp[(1 << n) - 1][0] << endl;
	}
}