BZOJ3105-新Nim游戏-线性基-拟阵-Nim博弈

3105: [cqoi2013]新Nim游戏

Time Limit: 10 Sec   Memory Limit: 128 MB
Submit: 1631   Solved: 960
[ Submit][ Status][ Discuss]

Description

传统的Nim游戏是这样的：有一些火柴堆，每堆都有若干根火柴（不同堆的火柴数量可以不同）。两个游戏者轮流操作，每次可以选一个火柴堆拿走若干根火柴。可以只拿一根，也可以拿走整堆火柴，但不能同时从超过一堆火柴中拿。拿走最后一根火柴的游戏者胜利。

本题的游戏稍微有些不同：在第一个回合中，第一个游戏者可以直接拿走若干个整堆的火柴。可以一堆都不拿，但不可以全部拿走。第二回合也一样，第二个游戏者也有这样一次机会。从第三个回合（又轮到第一个游戏者）开始，规则和Nim游戏一样。

如果你先拿，怎样才能保证获胜？如果可以获胜的话，还要让第一回合拿的火柴总数尽量小。

 

Input

第一行为整数k。即火柴堆数。第二行包含k个不超过10⁹的正整数，即各堆的火柴个数。

 

Output

 

输出第一回合拿的火柴数目的最小值。如果不能保证取胜，输出-1。

Sample Input

6
5 5 6 6 5 5

Sample Output

21

HINT

k<=100

Source

题解：

我们第一次拿完后，要使得剩下的火柴中不存在异或和为0的子集，否则对方会将先手必败的状态留给我们。

现在我们考虑剩下的火柴集合T。集合T的任意子集都不存在异或为0的子集 

=> 线性无关组 => 线性基是异或运算中的一个线性无关组

于是我们只需要构造一个线性基即可。

可是需要保留的元素和最大 => 需要构造一个极大线性无关组 => 极大线性基

根据拟阵证明(我也不会证!QAQ)

只需要贪心从大到小选择插入线性基中即可。

最后的答案就是所有元素总和 - 极大线性无关组的元素总和。

拟阵证明:

转自:传送门

那么现在设M=S,LM=S,L，S是所有数组成的矩阵，L是我们拿出来的极大线性无关组。

那么现在如果我们证得这个东西有遗传性并且也有交换性的话，那么他就是一个拟阵辣。其实主要就是证交换性的成立。

然后wyf大爷有这么一段证明。

设A,B∈I,且|A|<|B|,我们要证明存在x∈B,使得A∪{x}∈I.利用反证法，假设对于任意x∈B-A,均有A∪{x}不属于I,则B-A中的元素均在A的异或空间中，可由A的子集异或和表示。 
因此B中的元素都在A的异或空间中。那么必然有B的异或空间包含于A的异或空间。由|A|<|B|且A,B线性无关，显然矛盾。因此交换性存在，证毕

向量组的任意两个极大线性无关组等价。

代码：

/**
    我们第一次拿完后，要使得剩下的火柴中不存在异或和为0的子集，否则对方会将先手必败的状态留给我们。

    怎么使得剩下的火柴中不存在异或和为0的子集。
    这个时候我们想到了线性基。线性基就是一个异或运算的无关线性组。
    但是我们要求一个回合后的线性无关组最大。
    于是我们只需要求一个极大线性无关组即可。

    我们通过这n个元素从大到小构造出的线性基就是极大线性无关组。
*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 110;
const int MaxBasis = 62;
ll arr[maxn];
ll base[65];
struct LB {
    ll base[MaxBasis+10]; bool rel; int sz;
    vector<ll> Basis;
    LB() { memset(base, 0, sizeof(base)); rel = false; sz = 0; Basis.clear();}
    void init() {
        memset(base,0,sizeof(base));rel = false;sz = 0;Basis.clear();
    }
    bool add(ll x) {				///加入线性基中
        for (int i = MaxBasis; i >= 0; --i) {
            if (!(x >> i & 1)) continue;
            if (base[i]) x ^= base[i];
            else {
                for (int j = 0; j < i; ++j) if (x >> j & 1) x ^= base[j];
                for (int j = i+1; j <= MaxBasis; ++j) if (base[j] >> i & 1) base[j] ^= x;
                base[i] = x, ++sz;
                return true;
            }
        }
        rel = true; /// 表示有一个没插入进去,则说明当前线性基能表示出x，则如果放入x存在异或为0的结果。
        return false;
    }
    ll Max(ll ans = 0) {			///取最大
        for(int i=0;i<=MaxBasis;i++) ans = max(ans,ans^base[i]);
        return ans;
    }
    ll Min(ll ans = 0) {			///取最小
        for(int i=0;i<=MaxBasis;i++) ans = min(ans,ans^base[i]);
        return ans;
    }
    void GetBasis() {				///构造向量，用于之后求第k小
        for (int i = 0; i <= MaxBasis; ++i) if (base[i]) Basis.push_back(base[i]);
    }
    void umerge(struct LB &b) {			///合并线性基
        for(int i=0;i<=MaxBasis;i++) if(b.base[i])
            add(b.base[i]);
    }
    ll Min_Kth(ll k) {				///线性基中第k小
        if(rel) k--;				///线性基未满存在0
        if(k >= ((ll)1<<sz)) return -1;
        ll ans = 0;
        for(int i=0;i<sz;i++) if(k & ((ll)1<<i)){
            ans ^= Basis[i];
        }
        return ans;
    }
}lb;
int main()
{
    int n;scanf("%d",&n);
    ll ans = 0,sum = 0;
    for(int i=0;i<n;i++) scanf("%lld",&arr[i]),ans+=arr[i];
    sort(arr,arr+n);
    lb.init();
    for(int i=n-1;i>=0;i--) if(lb.add(arr[i])) sum += arr[i];
    printf("%lld\n",ans-sum);
    return 0;
}