0x35 高斯消元与线性空间

颓了十天回来做题果然……

感觉还是很有收获的，这两以前都没学过

bzoj1013: [JSOI2008]球形空间产生器sphere

poj1830(upd)

之前做得很烂还被 D飞*2 了。。重做一次

对于一个灯，把能把它点亮的其他灯设为1，然后高斯消元。

注意在这里的系数只是一个判断的手段，判断是否受到影响，并不是参与运算的（之前纠结了很久）

用二进制状压判断无解多解比较方便。（之前直接判当前位是不是0直接break错的一逼）

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;

int n,q[];
void guess()
{
    for(int j=;j<=n;j++)
    {
        for(int i=j+;i<=n;i++)
            if(q[i]>q[j])swap(q[i],q[j]);
        if(q[j]<=)
        {
                 if(q[j]==)printf("%d\n", (<<(n-j+)) );
            else if(q[j]==)printf("Oh,it's impossible~!!\n");
            return ;
        }

        int z;
        for(int i=n;i>=;i--)
            if(q[j]&(<<i)){z=i;break;}
        for(int i=;i<=n;i++)
        {
            if(i==j)continue;
            if(q[i]&(<<z))q[i]^=q[j];
        }
    }
    printf("1\n");
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int x,y;
        scanf("%d",&n);
        memset(q,,sizeof(q));
        for(int i=;i<=n;i++)scanf("%d",&x),q[i]^=x;
        for(int i=;i<=n;i++)scanf("%d",&x),q[i]^=x;
        for(int i=;i<=n;i++)q[i]|=(<<i);
        while(scanf("%d%d",&x,&y)!=EOF)
        {
            if(x==&&y==)break;
            q[y]|=(<<x);
        }
        guess();
    }
    return ;
}

poj1830

bzoj4004: [JLOI2015]装备购买

hdu3949(upd):没开LL见祖宗系列（记得是1LL<<i啊）
裸的线性基，但是求第k大的时候稍微注意下，求第k大的时候实际上求的是第k+1大，因为0会被看作第0大的……而且0需要特判存不存在，假如存在那么k--，否则这里就和前面的多算抵消了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef unsigned long long uLL;

uLL lt[],zero;
void insert(uLL x)
{
    for(int i=;i>=;i--)
        if(x&(1LL<<i))
        {
            if(lt[i]==){lt[i]=x;return ;}
            else x^=lt[i];
        }
    zero=;
}
int plen;uLL p[];
void rebuild()
{
    for(int i=;i>=;i--)
        if(lt[i]>)
        {
            for(int j=i-;j>=;j--)
                if(lt[i]&(1LL<<j))lt[i]^=lt[j];
        }
    plen=;
    for(int i=;i<=;i++)
        if(lt[i]>)p[++plen]=lt[i];
}
void findkth(uLL x)
{
    if(x>=(1LL<<plen)){printf("-1\n");return ;}
    uLL ans=;
    for(int i=plen;i>=;i--)
        if(x&(1LL<<i-))ans^=p[i];
    printf("%I64u\n",ans);
}

int main()
{
    int T,T_T=;
    scanf("%d",&T);
    while(T--)
    {
        int n;uLL x;
        scanf("%d",&n); zero=;
        memset(lt,,sizeof(lt));
        for(int i=;i<=n;i++)
            scanf("%I64u",&x), insert(x);
        rebuild();

        int Q;
        scanf("%d",&Q);
        printf("Case #%d:\n",++T_T);
        while(Q--)
        {
            scanf("%I64u",&x);x-=zero;
            findkth(x);
        }
    }
    return ;
}

hdu3949