魔法弹(线性基)

魔法弹(线性基)

给出n个m位数,问它们张成的异或集合的大小。

关于线性基,这篇博文讲的太好了:

所谓线性基,就是线性代数里面的概念。一组线性无关的向量便可以作为一组基底,张起一个线性的向量空间,这个基地又称之为线性基。这个线性基的基底进行线性运算,可以表示向量空间内的所有向量,也即所有向量可以拆成基底的线性组合。
在ACM领域,线性基主要用来处理有关异或和的极值问题。根据异或按照二进制数位进行的方式,我们可以把一个数字拆成它的二进制表示形式,而这个二进制形式拆成一位一位的,可以用向量来表示。显然,一组线性无关的向量可以张起一个向量空间,我们同样可以考虑构造这样一组数字的二进制形式组成的线性基。在这个线性基里面,我可以通过基底的线性组合、异或运算,表示所有的异或结果。
要构造这样的线性基,首先要满足一个性质:任取向量组中两个向量a、b,把其中一个替换成a xor b,这组向量线性组合得到的线性空间相同。这个性质对于xor来说是显然的,因为xor偶数具有抵消性。有了这个性质,就可以利用它对加入的向量进行修改。基本思路是:对于一个新加入的数字,从最高为考试往后扫,如果某一位该向量是1,但是之前的向量组中没有一个这一位是1,那么这一位就放如这个数字并且break;如果之前已经有该位为1的,那么该数字异或该位之前对应的数字。最后的结果是,新加入的数字x,要么被添加到某一个二进制位,要么变成0,说明这个数字对应的二进制向量不属于极大线性无关组。
线性基构造完毕,我们就要求其极值。很多时候我们需要求一堆数字最大的异或和,这个时候就要用到线性基。线性基构造完毕之后,根据它的定义,任何线性基所张的向量空间的所有向量都能够表示为其线性组合。所以我们只需要找,这个线性基中最大的向量即可。做法就是,从高位开始取线性基的每一个基底,如果异或能够使得结果变大,那么就异或。同时,这个还支持有初始值的最大值,因为线性基可以表示所有的向量,所以即使有一个初始值也是能够一起计算最大值的。

补充一下:n个的线性基能张成\(2^n\)个不同的元素,同时,有\(2^{n-1}\)个元素满足第i位为1(当然前提是线性基中有第i位为1的数)。可以这样考虑:找到一个第i位为1的线性基x,那么剩下n-1个元素随便异或,能张成\(2^{n-1}\)个数。为了让第i位为1,这些数中第i位为0的数需要异或x,否则就不用异或x。因此,这\(2^{n-1}\)个数都一一对应一个第i位为1的数。

对于这道题,尝试按位考虑。对于一个位,由于在线性基张成的空间中有\(2^{n-1}\)个数这个位为1,因此它对总和的贡献是\(x2^{n-1}\),x表示这个位的权值。

#include <cstdio>
#include <cstring>
using namespace std;

typedef long long LL;
const int maxn=405, mod=1e9+7;
int n, m, A[maxn];
struct Bigint{
    int a[maxn];  //从低到高
    bool empty(){
        int flag=1;
        for (int i=0; i<=m; ++i) if (a[i]) flag=0;
        return flag;
    }
    Bigint(){ memset(a, 0, sizeof(a)); }
    Bigint(char *s):Bigint(){
        for (int i=0; i<m; ++i) a[m-i-1]=s[i]-48;
        for (int i=0; i<m; ++i) A[i]|=a[i]; }
    int& operator[](int x){ return a[x]; }
    int value(){
        int re=0;
        for (int i=m; ~i; --i)
            re=(re*2+a[i])%mod;
        return re;
    }
}p[maxn];
Bigint operator ^(const Bigint &a, const Bigint &b){
    Bigint ans;
    for (int i=0; i<=m; ++i) ans[i]=a.a[i]^b.a[i];
    return ans; }

int fpow(int a, int x){
    LL base=a, ans=1;
    for (; x; x>>=1, (base*=base)%=mod)
        if (x&1) (ans*=base)%=mod;
    return ans; }

int main(){
    scanf("%d%d", &n, &m);
    char s[maxn]; Bigint t;
    for (int i=1; i<=n; ++i){
        scanf("%s", s);
        t=Bigint(s);
        for (int j=m; ~j; --j){
            if (!t[j]) continue;
            if (p[j].empty()) p[j]=t;
            t=t^p[j];
        }
    }
    LL ans=0, cnt=0, mi=1;
    for (int i=0; i<=m; ++i) if (!p[i].empty()) ++cnt;
    for (int i=0; i<=m; ++i){
        if (i) mi<<=1; mi%=mod;
        if (!A[i]) continue;
        ans+=fpow(2, cnt-1)*mi; ans%=mod;
    }
    printf("%lld\n", ans);
    return 0;
}
posted @ 2018-10-05 15:46 pechpo 阅读(...) 评论(...) 编辑 收藏