B - Binomial Gym - 102576B(Lucas定理,SOS DP)

5 篇文章 0 订阅
2 篇文章 0 订阅

在这里插入图片描述

题意:
n个数,求多少个数满足 ( a [ i ] a [ j ] ) \tbinom{a[i]}{a[j]} (a[j]a[i])为奇数。

思路:
在这里插入图片描述
由卢卡斯定理可以得到,如果存在n的某个二进制位为0,且m的该二进制位为1,那么就有
( n m ) m o d    2 = 0 \tbinom{n}{m}\mod2=0 (mn)mod2=0

所以只有m是 n n n的子集的时候 ( n m ) \tbinom{n}{m} (mn)才会是奇数。

那么实际上就是求子集和,这个直接用SOS DP就好了。

#pragma optimize("-O3")

#include <cstdio>
#include <cassert>
#include <cmath>
#include <vector>
#include <cstring>

using namespace std;

typedef long long ll;
const int maxn = (1 << 21) + 5;
int dp[maxn],cnt[maxn];

int main() {
    int T;scanf("%d",&T);
    while(T--) {
        memset(dp,0,sizeof(dp));
        memset(cnt,0,sizeof(cnt));
        int n;scanf("%d",&n);
        for(int i = 1;i <= n;i++) {
            int x;scanf("%d",&x);
            dp[x]++;
            cnt[x]++;
        }
        int all = 1 << 21;
        for(int j = 0;j <= 21;j++) {
            for(int i = all;i >= 0;i--) {
                if(i & (1 << j)) {
                    dp[i] += dp[i ^ (1 << j)];
                }
            }
        }
        ll ans = 0;
        for(int i = 0;i <= all;i++) {
            if(cnt[i]) {
                ans += 1ll * dp[i] * cnt[i];
            }
        }
        printf("%lld\n",ans);
    }
    return 0;
}

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值