蒟蒻的第一篇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;
}