CF 251D Two Sets

贪心+高斯消元

显然从高到低按位考虑。如果这一位有奇数个1，那么这一位最终对答案的贡献肯定是1。如果有偶数个1，对答案的贡献可能是0或2。要最大化答案，因此我们需要更多的2。列出方程，系数$a_{i,j}$表示第j个数的第i位的值。变量$x_i$表示第i个数是否归自己。要出现尽量多的2，只需让自己这边这一位的异或和为1即可。即a1x1^a2x2^…^anxn=1。判断如果无解说明不能做到，去掉该方程即可。没有任何限制作用的方程也可以去掉。

在值最大的前提下需要自己的异或和最小，同理列出另一套方程，令它们=0即可。

#include<cstdio>
#include<bitset>
#define N 100005
#define K 62
using namespace std;
namespace runzhe2000
{
    typedef long long ll;
    bitset<N> f[K], base[K];
    int pos[N], ans[N];
    void main()
    {
        int n, tot = 0; scanf("%d",&n);
        ll sum = 0, x;
        for(int i = 1; i <= n; i++)
        {
            scanf("%I64d",&x);
            sum ^= x; 
            for(int j = K-1; ~j; j--)
                if((x>>j) & 1) base[j][i] = 1;
        }
        for(int bit = 0; bit <= 1; bit++)
        {
            for(int i = K-1; ~i; i--)
            {
                if( ((sum>>i)&1) == bit)
                {
                    f[++tot] = base[i];
                    f[tot][n+1] = (bit^1);
                    for(int j = 1; j <= n; j++)
                        if(pos[j] && f[tot][j])
                            f[tot] ^= f[pos[j]];
                    int p = 0;
                    for(int i = 1; i <= n; i++)
                        if(f[tot][i])
                        {
                            p = i;
                            break;
                        }
                    if(!p) --tot;
                    else pos[p] = tot;
                }
            }
        }
        for(int i = n; i; i--)
            if(pos[i])
            {
                ans[i] = f[pos[i]][n+1];
                for(int j = i-1; j; j--)
                    if(f[pos[j]][i])
                        f[pos[j]][n+1] = f[pos[j]][n+1] ^ ans[i];
            } 
        for(int i = 1; i <= n; i++)
            printf("%d ",ans[i]?1:2);
    }
}
int main()
{
    runzhe2000::main();
}