uva 11542 - Square（高斯异或）

题意：给你n个数，从中选1个或多个，使得选出的整数的乘积是完全平方数，问有多少种选法?

解析：

1. xi=1表示选择第i个数，否则不选 
2.     4,6,10,15的因子只有2,3,5;可分解为4=2^2*3^0*5^0->(2,0,0)... 
3.      
4.     那么选出来的数的乘积就是2^(2*x1+x2+x3)*3^(x2+x4)*5^(x3+x4)，要使其为完全平方数，则各个指数必须是2的倍数， 
5. 即各个指数中各个和数的异或值为0.由于相同数的异或值为0，所以可以将xi的系数换成0或1 
6.     上例可得方程： 
7.                 x2^x3=0 
8.                 x2^x4=0 
9.                 x3^x4=0 
10.     得到矩阵(横坐标为各素数因子，纵坐标为n个输入值) 
11.             |0 1 1 0 |0| 
12.             |0 1 0 1 |0| 
13.             |0 0 1 1 |0| 
14.     用高斯消元进行异或消元，最后得到自由变量个数ans，那么可能取法有2^ans种(0,1取法)，还要减去全不取的情况 

#include<iostream>
#include<cstdio>
#include<math.h>
#include<vector>
#include<string.h>
using namespace std;
typedef long long LL;
#define N 505
int A[N][N];

int prime[N],a[N];
int inint()
{
    int i,j,k,c=0;
    memset(a,0,sizeof(a));
    for(i=2; i<N; i++)
    {
        if(!a[i]) prime[c++]=i;
        for(j=0; j<c&&(k=i*prime[j])<N; j++)
        {
            a[k]=1;
            if(i%prime[j]==0) break;
        }
    }
    return c;
}

int gauss(int n,int m)
{
    int i=0,j=0,k,r,t;
    while(i<m&&j<n)
    {
        r=i;
        for(k=i;k<m;k++)
          if(A[k][j]) {  r=k; break;  }
        if(A[r][j])
        {
            if(r!=i)
            for(k=0;k<=n;k++) swap(A[r][k],A[i][k]);

            for(t=i+1;t<m;t++)
              if(A[t][j])
              for(k=i;k<=n;k++)
                   A[t][k]^=A[i][k];
            i++;
        }
        j++;
    }
    return i;
}

int main()
{
    int i,j,k,n,t,res,ge,numf;
    LL ans,x;
    int m=inint();
    scanf("%d",&t);
    while(t--)
    {
        memset(A,0,sizeof(A));
        scanf("%d",&n);
        numf=0;
        for(i=0; i<n; i++)
        {
            scanf("%lld",&x);
            for(j=0; j<m; j++)
            {
                if(x%prime[j]==0)
                {
                    numf=max(numf,j);
                    while(x%prime[j]==0)
                    {
                        x/=prime[j];
                        A[j][i]^=1;
                    }
                }
            }
        }
        numf++;
        k=gauss(n,numf);
        ans=1;
        printf("%lld\n",(ans<<(n-k))-1);
    }
    return 0;
}