2020牛客多校第二场E题 Exclusive OR FWT

Exclusive OR

题意

n n n个数 a 1 , a 2 , . . . , a n a_1,a_2,...,a_n a1,a2,...,an, 求从中取出 i i i个数的最大异或和,输出 i = 1 , 2 , . . . , n i=1,2,...,n i=1,2,...,n的答案

题解

FWT: C x = ∑ i ⊕ j = x A i B j C_x=\displaystyle\sum_{i \oplus j=x} A_iB_j Cx=ij=xAiBj
这是一个卷积的形式,只不过平常的卷积运算符是 + + +
我们用 f 1 [ x ] f_1[x] f1[x]记录 a 1 , a 2 , . . . , a n a_1,a_2,...,a_n a1,a2,...,an x x x是否出现
那么 f 2 [ x ] = ∑ i ⊕ j = x f 1 [ i ] f 1 [ j ] f_2[x]=\displaystyle\sum_{i \oplus j=x} f_1[i]f_1[j] f2[x]=ij=xf1[i]f1[j]代表的含义就是

选取两个数 i , j i,j i,j的异或值为 x x x的种数

同理, f k [ x ] = ∑ i ⊕ j = x f k − 1 [ i ] f 1 [ j ] f_k[x]=\displaystyle\sum_{i \oplus j=x} f_{k-1}[i]f_1[j] fk[x]=ij=xfk1[i]f1[j]代表的含义为

在选取了 k − 1 k-1 k1个数能得到的异或和数中选取一个 i i i,在初始数字中选取一个 j j j, 使得 i ⊕ j = x i \oplus j = x ij=x的种数

所以这就是一个FWT的题
但是我们不需要知道取得 x x x的方案数有多少种,我们只关心能不能取到,所以每次卷积完之后,都将能取到的 x x x设为 1 1 1

另外当 i > 19 i>19 i>19的情况,题目说的比较清楚了
题解

代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAX = 2e5 + 10;

int N;
int a[1 << 18], b[1 << 18], n;
int ans[MAX];

void get() { for (int i = 0; i < n; i++) a[i] *= b[i]; }

void XOR(int *f, int x = 1) {
    for (int o = 2, k = 1; o <= n; o <<= 1, k <<= 1)
        for (int i = 0; i < n; i += o)
            for (int j = 0; j < k; j++) {
                int a0 = f[i + j], a1 = f[i + j + k];
                f[i + j] = a0 + a1;
                f[i + j + k] = a0 - a1;
                if (x == -1) {
                    f[i + j] /= 2;
                    f[i + j + k] /= 2;
                }
            }
}


int main() {

    scanf("%d", &N); n = 1 << 18;
    for (int i = 1; i <= N; i++) {
        int x; scanf("%d", &x);
        a[x] = b[x] = 1;
        ans[1] = max(ans[1], x);
    }
    XOR(b,1);
    for (int i = 2; i <= 20; i++) {
        XOR(a, 1); get(); XOR(a, -1);
        for (int j = 0; j < n; j++)
            if (a[j]) {
                a[j] = 1;
                ans[i] = max(ans[i], j);
            }
    }
    for (int i = 20; i <= N; i++)
        ans[i] = ans[i - 2];
    for (int i = 1; i <= N; i++)
        printf("%d ", ans[i]);

    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值