uva11542 Square（异或方程组）

蒟蒻的第一篇blog，写得认真一点

好吧我承认几乎是抄LRJ书上的，但我把它详(fu)细(za)化了. . . 

【题目】

http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&category=27&page=show_problem&problem=2537

【题解】

设这n个数分别为a1,a2,…,an
∵ ai不含大于500的素因子 
∴ 将ai分解质因数 
例如：
n==4，这4个数分别为4,6,10,15
a1 == 2^2
a2 == 2^1 + 3^1
a3 == 2^1       + 5^1
a4 ==       3^1 + 5^1
若挑出几个，它们的乘积A为完全平方数，A == 2^p1 + 3^p2 + 5^p3 +…，则：p1,p2,…均为偶数 


因此，对于上面的例子，用x1,x2,…,xn表示a1,,a2,…,an取不取：xi==0表示不取ai，xi==1表示取ai，
则有：2*x1 + x2 + x3 ≡0(mod2)
             x2 + x4 ≡0(mod2)
             x3 + x4 ≡0(mod2)
其中，由于x*a ≡(x%2)*a (mod2)，因此可以把所有的xi的系数化为0或1

转化为下面的异或方程组：
          x2 xor x3 == 0 ① 
          x2 xor x4 == 0 ② 
          x3 xor x4 == 0 ③ 

对于一般的xor方程组：(n个未知数，m组方程，所有系数ai,bi,…和 ei 均为0或1)
a1*x1 xor b1*x2 xor c1*x3 xor d1*x4 xor …== e1
a2*x1 xor b2*x2 xor c2*x3 xor d2*x4 xor …== e2
……
am*x1 xor bm*x2 xor cm*x3 xor dm*x4 xor …== em
其解法类似于高斯消元：
对于x1,找到任意一个ai==1，用它把其他aj==1的方程中的这一项都消掉 
同样方法处理其他方程中的含x2,x3,…的项 

比如上例中：
1. 所有方程中x1系数都为0
2. ①,②中x2系数都为1，①消②(两式左右分别xor)得：x3 xor x4 == 0 ②'
    ③中x2系数为0，不管它 
3. ②',③中x3系数都为1，②'消③得0==0
消元后仅有2个方程：
          1*x2 xor 1*x3 xor 0*x4 == 0
                   1*x3 xor 1*x4 == 0
x1,x4可以任取，称为"自由变量"; x2,x3可以被推出 
对于有f个自由变量的xor方程组，共有2^f组解 


对于本题，不允许所有xi均为0，因此：ans=2^f-1

代码如下

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
typedef long long LL;
LL pri[105]={0,2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97,101,103,107,109,113,127,131,137,139,149,151,157,163,167,173,179,181,191,193,197,199,211,223,227,229,233,239,241,251,257,263,269,271,277,281,283,293,307,311,313,317,331,337,347,349,353,359,367,373,379,383,389,397,401,409,419,421,431,433,439,443,449,457,461,463,467,479,487,491,499};
int A[105][105]={0};
void jh(int* a,int* b)
{
	int t=*a;
	*a=*b;
	*b=t;
}
int gauss(int n,int m)//n个未知数，m组方程，返回非自由变量的个数 
{
	int i=1,j=1,k,l,r;
	while(i<=m&&j<=n)//正在处理第i个方程，第j个未知数 
	{
		for(r=i;r<=m;r++)
			if(A[r][j]==1) break;
		if(r<=m)
		{
			if(r!=i)
				for(k=j;k<=n;k++)
					jh(&A[i][k],&A[r][k]);
			for(k=i+1;k<=m;k++)
				if(A[k][j]==1)
					for(l=j;l<=n;l++)
						A[k][l]^=A[i][l];
			i++;
		}
		j++;
	}
	return i-1;
}
int main()
{
	LL x;
	int T,n,i,j,num=0,f;//num:构造方程的组数 
	scanf("%d",&T);
	for(;T>0;T--)
	{
		scanf("%d",&n);
		memset(A,0,sizeof(A));
		num=0;
		for(i=1;i<=n;i++)
		{
			scanf("%lld",&x);
			for(j=1;x>1;j++)
				while(x%pri[j]==0)
				{
					x/=pri[j];
					if(num<j) num=j;
					A[j][i]^=1;
				}
		}
		f=gauss(n,num);
		printf("%lld\n",(1LL<<n-f)-1);
	}
	return 0;
}