BZOJ4328 始祖鸟 [高斯消元][bitset?没听说过。。]

什么$N^3,N=2000$要T？没听说过，小爷我就是不服，强行数组矩阵1A,就是那么华丽。。
其实还是可以用bitset的，那么友善的东西，丫丫常数也是极好的。

高斯消元灵活运用题，难度不大。
引用SD_le大大的题解吧，言简意赅。

把每只鸟在上游或下游看成一个0/1变量。
对于偶数条出边的鸟，要求它指向的鸟选1的为偶数个。
对于奇数条出边的鸟，如果它选1，要求它指向的鸟选1的为偶数个，否则为奇数个。

#include<bits/stdc++.h>
using namespace std;
int f[2005],n,a[2005][2005],t,m;
inline void read(int &res){
    static char ch;register int flag=1;
    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;
    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;
}
int main(){
    read(n);
    for(register int i=1;i<=n;i++){
        read(m);
        if(m&1)a[i][i]=a[i][0]=1;
        for(register int j=1;j<=m;j++)
            read(t),a[i][t]=1;
    }
    for(register int i=1;i<n;i++){
        register int p=0;
        for(register int j=i;j<=n;j++)
            if(a[j][i]){p=j;break;}
        if(!p)continue;
        if(p!=i){
            swap(a[p],a[i]);
        }
        for(register int j=1;j<=n;j++)
            if(a[j][i]&&j!=i)
                for(register int k=0;k<=n;++k)a[j][k]^=a[i][k];//就是直接硬钢！
    }
    t=0;
    for(register int i=n;i;i--){
        if(a[i][i]){
            f[i]=a[i][0];
            if(f[i])t++;
        }else if(a[i][0]){
            puts("Impossible");
            return 0;
        }
    }
    printf("%d\n",t);
    for(register int i=1;i<=n;i++)
        if(f[i])printf("%d ",i);
    return 0;
}

---

注意，本人已经将数据更新，上诉方法已经无法卡过了！！

new code

#include<bits/stdc++.h>
using namespace std;
int f[2005],n,a[2005],t,m;
bitset<2005>a[2005];
inline void read(int &res){
    static char ch;register int flag=1;
    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;
    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;
}
int main(){
    read(n);
    for(register int i=1;i<=n;i++){
        read(m);
        if(m&1)a[i][i]=a[i][0]=1;
        for(register int j=1;j<=m;j++)
            read(t),a[i][t]=1;
    }
    for(register int i=1;i<n;i++){
        register int p=0;
        for(register int j=i;j<=n;j++)
            if(a[j][i]){p=j;break;}
        if(!p)continue;
        if(p!=i){
            swap(a[p],a[i]);
        }
        for(register int j=1;j<=n;j++)
            if(a[j][i]&&j!=i)a[j]^=a[i];
    }
    t=0;
    for(register int i=n;i;i--){
        if(a[i][i]){
            f[i]=a[i][0];
            if(f[i])t++;
        }else if(a[i][0]){
            puts("Impossible");
            return 0;
        }
    }
    printf("%d\n",t);
    for(register int i=1;i<=n;i++)
        if(f[i])printf("%d ",i);
    return 0;
}