POJ 3986 数位dp

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int mod = 1000000003;
const int maxn = 1 << 6;
int n, k, inp[maxn], dp[maxn][maxn][2], pw[maxn];
long long solve(int i, int pre, int len)
{
	long long cur = 0, res = 0;
	pre &= ~((1 << len) - 1);
	if (n + 1 == i) return !pre;
	int& k = dp[i][len][(pre & (1 << len)) ? 1 : 0];
	if (k != -1) return k;
	for (int j = 31; j >= 0; j--)
		if (inp[i] & (1 << j))
			res += solve(i + 1, pre ^ cur, max(j, len)) * pw[min(j, len)] % mod, res %= mod, cur |= 1 << j;
	return k = res;
}
void init()
{
	pw[0] = 1;
	for (int i = 1; i < maxn; i++)
		pw[i] = (pw[i - 1] * 2) % mod;
}
int main(int argc, char const *argv[])
{
	init();
	while (~scanf("%d%d", &n, &k) && n + k)
	{
		for (int i = 1; i <= n; i++)
			scanf("%d", &inp[i]), inp[i]++;
		memset(dp, -1, sizeof(dp));
		printf("%lld\n", solve(1, k, 0));
	}
	return 0;
}


dp[i][pre][len]表示Xi~Xn[最低位为第len位的前缀]为pre的方案数

dp[i][pre][len]= Sigma(dp[i+1][pre^cur][max(j,len)] * 2^min(j,len))

详细来自:http://blog.csdn.NET/snowdut/article/details/51581438


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值