[POJ3986]Math teacher's homework && 数位DP

首先假设xi足够大 就会有

ans = (m1 + 1) * (m2 + 1) * ... *(mi-1 + 1) * (mi+1 + 1) *....

因为对于此时xi = x1 ^ x2 ^ x3 ^....^ xn ^ k 它被其余n-1个元素唯一确定 我们可以针对这一点进行DP

从高位往低位枚举 记录K的当前位aim 

用状态f[i][j] 表示考虑前 i 个元素， 当前位(记作bit位)为j的方案数

对于每一个考虑过的元素 我们视为其取到最大值（和数位DP的思想差不多吧） 那么当我们对xi进行讨论的时候 前i-1个元素的方案数是确定的 只需要对i+1 ~ n的元素进行DP 

对于每一个xj 如果设他当前位的值为p[j]

若p[j] == 1 

    则xj当前位取0的方案数就有 2^bit 种， 取1的方案数就有 xj % （2^bit) + 1种 (因为视为考虑过的数位取到最大值)

若p[j] == 0

   则xj当前位取0的方案数为 xj % （2^bit) + 1， 取1的方案数为0

设取0的方案数记为d0， 取1的方案数记为d1

那么转移方程可以写作（按照异或的计算规则进行转移）

f[k][0] = f[k-1][0] * d0 + f[k-1][1] * d1

f[k][1] = f[k-1][0] * d1 + f[k-1][0] * d0

如果发现当前位每个xi都取到最大值时与K的当前位不符 则不能均取最大值 退出循环

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<queue>
#include<cmath>
#include<deque>
#define SF scanf
#define PF printf
using namespace std;
typedef long long LL;
const int MAXN = 60;
const int mod = 1000000003;
LL f[MAXN+10][2];
LL pow2[MAXN+10];
bool a[MAXN+10];
void init() {
	pow2[0] = 1;
	for(int i = 1; i <= MAXN; i++) pow2[i] = pow2[i-1] << 1;
}
int A[MAXN+10], n, K;
int main() {
	init();
	while(~SF("%d%d", &n, &K) && n+K) {
		for(int i = 1; i <= n; i++) SF("%d", &A[i]);
		LL ans = 0;
		bool ok = true;
		for(int i = 30; i >= 0; i--) {
			for(int j = 1; j <= n; j++)
				a[j] = (bool)(pow2[i] & A[j]);
			int aim = (bool)(pow2[i] & K);
			for(int j = 1; j <= n; j++) {
				if(a[j]) {
					LL t = 1;
					int cur = 0;
					for(int k = 1; k < j; k++) {
						t = t * (A[k] % pow2[i] + 1) % mod;
						cur ^= a[k];
					}
					f[j][0] = f[j][1] = 0;
					f[j][cur] = t;
					for(int k = j+1; k <= n; k++) {
						LL d0 = a[k] ? pow2[i] : A[k] % pow2[i] + 1;
						LL d1 = a[k] ? A[k] % pow2[i] + 1 : 0;
						f[k][0] = (f[k-1][0] * d0 + f[k-1][1] * d1) % mod;
						f[k][1] = (f[k-1][0] * d1 + f[k-1][1] * d0) % mod;
					}
					(ans += f[n][aim]) %= mod;
				}
			}
			int high = 0;
			for(int j = 1; j <= n; j++) high ^= a[j];
			if(high != aim) {
				ok = false;
				break;
			}
		}
		if(ok) ans = (ans+1) % mod;
		PF("%d\n", (int)ans);
	}
}