hdu3929 容斥+扩展Lucas定理

最新推荐文章于 2023-04-17 18:53:06 发布

羁绊残阳

最新推荐文章于 2023-04-17 18:53:06 发布

阅读量465

点赞数

分类专栏： ACM_组合数学

本文链接：https://blog.csdn.net/jibancanyang/article/details/52629868

版权

ACM_组合数学专栏收录该内容

9 篇文章 0 订阅

订阅专栏

分析：

根据扩展lucas定理， $(1+x)^n$ 含有的奇数项的个数是 $n$ 的二进制表达中数字1的个数。
i & -i得到的是 $i$ 的二进制表达最右边的 1 <script type="math/tex" id="MathJax-Element-114">1</script>的位置代表的十进制值。我们每次都减去它一个1，这样就有个函数可以算1的个数。
函数实现如下：

int get1(int x) {
    return x ? get(x - (x & -x)) + 1 : 0;
}

接下来容斥。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define pr(x) cout << #x << ": " << x << "  " 
#define pl(x) cout << #x << ": " << x << endl;

struct jibancanyang
{
    long long ans, n, A[22];

    int get(long long x) {
        return x ? get(x - (x & -x)) + 1 : 0;
    }


    void dfs(int ith, long long x, int f) {
        ans += ((long long)1 << get(x)) * f;
        for (int i = ith + 1; i <= n; ++i) 
            dfs(i, x & A[i], -2 * f);
    }

    void fun() {
        int T;
        scanf("%d", &T);
        for (int cas = 1; cas <= T; ++cas) {
            scanf("%lld", &n);
            for (int i = 1; i <= n; ++i) scanf("%lld", &A[i]);
            ans = 0;
            for (int i = 1; i <= n; ++i) 
                dfs(i, A[i], 1);
            printf("Case #%d: %lld\n", cas, ans);
        }
    }
}ac;

int main()
{
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
    //freopen("out.txt", "w", stdout);
#endif
    ac.fun();
    return 0;
}