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;
}