【HDU 5833】Zhu and 772002（异或方程组高斯消元讲解）

#include<cstdio>
#include<stdio.h>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<queue>
#define INF 0x3f3f3f
#define MAX 2105
#define MOD 1000000007

using namespace std;

long long c[MAX][505],p[MAX],v[MAX],a[MAX],cnt,n;//c存矩阵，p存素数表，cnt代表素数的个数

void MakeTab()//打素数表
{
    int i,j;
    memset(v,0,sizeof(v));
    memset(p,0,sizeof(p));
    cnt=0;
    for(i=2; i<=2000; i++)
    {
        if(!v[i])
        {
            p[++cnt]=i;
            for(j=i; j<=2000; j+=i)
            {
                v[j]=1;
            }
        }
    }
}

int Rank()//计算秩
{
    int i,j,k,r,u;
    i=1;
    j=1;
    while(i<=cnt && j<=n)
    {
        r=i;
        while(!c[r][j] && r<=cnt)
            r++;
        if(c[r][j])
        {
            swap(c[i],c[r]);//如果发现了第r行第j列为1，就讲r行和i行行互换（初等行变换）
            for(u=i+1; u<=cnt; u++)
            {
                if(c[u][j])
                {
                    for(k=i; k<=n; k++)
                    {
                        c[u][k]=c[u][k]^c[i][k];//每找到一个未知数就对其进行亦或处理，去掉系数c[i][k]
                    }

                }
            }
            i++;
        }
        j++;
    }
    return i;
}

int main()
{
    MakeTab();
    int i,j,cns=1,T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%lld",&n);
        memset(c,0,sizeof(c));
        for(i=1; i<=n; i++)
        {
            scanf("%lld",&a[i]);
        }

        for(i=1; i<=n; i++)
        {
            for(j=1; j<=cnt; j++)
            {
                long long num=a[i];
                if(num%p[j]==0)
                {
                    while(num%p[j]==0)
                    {
                        num/=p[j];
                        c[j][i]=c[j][i]^1;
                    }
                }
            }
        }

        long long k=(n-Rank());
        long long ans=1;
        for(i=1; i<=k; i++)
            ans=(ans*2)%MOD;
        printf("Case #%d:\n",cns++);
        printf("%lld\n",ans-1);//去掉全0的解
    }
    return 0;
}