HDU-5833-Zhu and 772002（高斯消元）

链接：http://acm.hdu.edu.cn/showproblem.php?pid=5833

题意：

给出n个数，每次可以从中取任何数，问有多少种取法使得取出来的数的乘积是完全平方数。

题解：

首先显然可以去掉偶数幂的质因子，每个数就变成了不同质数相乘的形式，已知题意中只有大约300多个质因子，直接列成质因子的01矩阵形式，高斯消元，求一下自由变元个数temp，显然答案就是temp^1-1。

CODE：

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define bug cout<<"bug\n"
const int maxn = 1e5+7;
#define MAXN 2030
const long long mod = 1e9+7;
bool not_prime[MAXN*30];
long long prime[MAXN];
int num_prime;
int a[350][350];
void get_prime()
{
    long long i,j;
    memset(not_prime,0,sizeof(not_prime));
    num_prime=0;
    not_prime[1]=1;
    for(i=2; i<MAXN; ++i)
        if(!not_prime[i])
        {
            num_prime++;
            prime[num_prime]=i;
            for(j=i*i; j<MAXN; j+=i)
                not_prime[j]=1;
        }
}
int gauss(int m, int n)
{
    int i,j,k,r;
    for(i=j=0; i<m&&j<n ; ++i,++j)
    {
        r=i;
        for(int c=i; c<m; ++c)
        if(a[c][j]){r=c;break;}
        if(a[r][j])
        {
            if(r!=i)for(int c=0; c<=n; ++c)swap(a[r][c],a[i][c]);
            for(k=i+1; k<m; ++k)if(a[k][j])
                for(int c=i; c<=n; ++c)
                    a[k][c]^=a[i][c];
        }
        else --i;
    }
    return i;
}
int main()
{
    int T;
    int n,m;
    long long p;
    scanf("%d",&T);
    int cas=1;
    get_prime();
    while(T--)
    {
        memset(a,0,sizeof(a));
        scanf("%d",&n);
        int M=0;
        for(int i=0; i<n; ++i)
        {
            scanf("%I64d",&p);
            for(int j=1; p>=prime[j] && j<=num_prime; ++j)
            {
                while(p%prime[j]==0)
                {
                    M=max(M,j);
                    p/=prime[j];
                    ++a[i][j-1];
                }
                a[i][j-1]%=2;
            }
        }
        long long temp=n-gauss(n,M);
        long long ans=1;
        for(int i=0; i<temp; ++i){ans<<=1LL;ans%=mod;}
        printf("Case #%d:\n%I64d\n",cas++,ans-1);
    }
    return 0;
}
/*
100
3
3 3 4
3
2 2 2
5
4 4 4 4 4
9
4 4 4 4 4 5 5 5 5
6
4 4 80 20 5 5
*/