CF251D Two Sets

Description

---

Philips非常喜欢数字。最近pwecar送给他了一个由n个非负整数组成的集合。Philips非常喜欢和Bing玩。他立即决定把他n个数字中的一部分送给Bing。为了让游戏更加有趣，Philips决定使得给她的数字集合满足如下条件：
我们用x1表示Philips的数字集合的xor值，用x2表示Bing的数字集合的xor值。要使得x1+x2尽可能地大。假如有多种划分集合的方法使得集合满足上述条件，Philips就要让x1尽可能地小。
Xor运算是“异或”，在Pascal语言中用”xor”表示，在C/C++/Java语言中用”^”表示。
帮助Philips按照上述方法划分集合。如果有多种合适的方法，就找出其中任意一种方法。请注意，Philips将一部分数字给了Bing之后，他可能就没有任何剩余的数字了。反之亦然，Philips也可以不给Bing任何数字。在这两种情况下，我们都假定空集的xor值为0。

对于30%的数据n<=10。
对于60%的数据n<=1000。
对于100%的数据n<=100000。

Solution

---

看到位运算的题目就要想到按位做了

从高位到低位贪心，设w为第i位上为1的数的数量，对m分奇偶性讨论
若m为奇数，最优肯定是第一集合分偶数个，第二集合分奇数个，贡献就是这一位上的1
若m为偶数，就是两边都分奇数个
讨论完奇偶性可以根据属于集合的情况列异或方程组，求一个任意解就行了

求解可以每次只加入一个式子，用之前的式子消去这个式子上的1，再用这个式子上剩余的1消去之前的式子（绕

Code

---

#include <stdio.h>
#include <string.h>
#include <bitset>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)

typedef long long LL;
const int N=200005;
const int L=62;

std:: bitset <N> a[L+1],r[L+1];
// int a[L+1][N],r[L+1][N];
int rec[N],ans[N];
int tot,n;

int main(void) {
    int n; LL sum=0; scanf("%d",&n);
    rep(i,1,n) {
        LL x; scanf("%lld",&x);
        sum^=x;
        rep(j,1,L) r[j][i]=(x&(1LL<<(j-1)))>0;
    }
    drp(j,L,1) {
        if (!((sum>>(j-1))&1)) {
            a[++tot]=r[j]; a[tot][n+1]=1;
            rec[tot]=n+2;
            rep(i,1,tot-1) if (a[tot][rec[i]]) a[tot]^=a[i];
            rep(i,1,n) if (a[tot][i]) {rec[tot]=i; break;}
            if (rec[tot]>n) {
                tot--;
                continue;
            }
            rep(i,1,tot-1) if (a[i][rec[tot]]) a[i]^=a[tot];
        }
    }
    drp(j,L,1) {
        if ((sum>>(j-1))&1) {
            a[++tot]=r[j]; a[tot][n+1]=0;
            rec[tot]=n+2;
            rep(i,1,tot-1) if (a[tot][rec[i]]) a[tot]^=a[i];
            rep(i,1,n) if (a[tot][i]) {rec[tot]=i; break;}
            if (rec[tot]>n) {
                tot--;
                continue;
            }
            rep(i,1,tot-1) if (a[i][rec[tot]]) a[i]^=a[tot];
        }
    }
    rep(i,1,tot) ans[rec[i]]=a[i][n+1];
    rep(i,1,n) printf("%d ", 2-ans[i]);
    return 0;
}