BZOJ4903: [Ctsc2017]吉夫特

传送门
可以发现, ( n m ) ≡ 1 ( m o d   2 ) \binom{n}{m}\equiv 1(mod~2) (mn)1(mod 2) 当且仅当 m   a n d   n   =   m m~and~n~=~m m and n = m
m m m 二进制下为 n n n 的子集
那么可以直接写一个 3 18 3^{18} 318 的枚举子集 D P DP DP
但是还有一个 6 9 6^9 69 的做法
把数字分成前 9 9 9 位和后 9 9 9
f ( s 1 , s 2 ) f(s_1,s_2) f(s1,s2) 表示前 9 9 9 位为 s 1 s_1 s1,后 9 9 9 位为 s 2 s_2 s2 的超集的答案
那么对于一个数 x x x,分成 x 1 , x 2 x_1,x_2 x1,x2,转移的时候枚举 x 1 x_1 x1 的超集,更新的时候枚举 x 2 x_2 x2 的子集即可

# include <bits/stdc++.h>
using namespace std;
typedef long long ll;

const int maxn(1 << 9);
const int mod(1e9 + 7);

inline void Inc(int &x, int y) {
	x = x + y >= mod ? x + y - mod : x + y;
}

int n, f[maxn][maxn], sz = maxn - 1;

int main() {
	register int i, j, v, f1, f2, t, g;
	scanf("%d", &n);
	for (i = 1; i <= n; ++i) {
		scanf("%d", &v), f1 = v >> 9, f2 = v & sz, g = 1;
		for (t = j = sz ^ f1; ; j = (j - 1) & t) {
			Inc(g, f[sz ^ j][f2]);
			if (!j) break;
		}
		for (j = f2; ; j = (j - 1) & f2) {
			Inc(f[f1][j], g);
			if (!j) break;
		}
	}
	for (g = mod - n, i = 0; i <= sz; ++i) Inc(g, f[i][0]);
	printf("%d\n", g);
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值