SGU200 Cracking RSA

题目大意：给定m(m<=100)个数，它们都是由前t(t<=100)个素因子构成的。求有多少个这m个数的子集，使得他们的乘积是完全平方数。

完全平方数的每个素因子出现偶数次。

可以列出t个方程，m个变量的方程组，即t*m的01矩阵，0代表素因子出现偶数次，1代表出现奇数次。

问题即转化为求使每个方程结果都为 0 的解的个数。

时间复杂度为O（t^2*m）

#include<iostream>
#include<cstdio>
using namespace std;
const int maxx=10000;
int n,m;
int ans;
int p[10001];
bool vst[10001];
int b[101][101];
void Eular() {
	int i,j,k;
	for (i=2; i<=maxx; i++) {
		if (!vst[i]) p[++p[0]]=i;
		for (j=1; j<=p[0]&&i*p[j]<=maxx; j++) {
			k=i*p[j];
			vst[k]=1;
			if (i%p[j]==0) break;
		}
	}
}
void fac(int u,int v) {
	int i,t;
	if (v==1) return;
	for (i=1; i<=100; i++)
		if (v%p[i]==0) {
			t=0;
			while (v%p[i]==0) {
				v/=p[i];
				t++;
			}
			b[i][u]=t&1;
			if (v==1) return;
		}
}
void Gauss() {
	int i,j,k,t;
	for (i=1,t=0; i<=m&&t<=n; t++) {
		for (j=i; j<=m; j++)
			if (b[j][t]) break;
		if (j==m+1) continue;
		if (j!=i)
			for (k=t; k<=n; k++)
				swap(b[i][k],b[j][k]);
		for (j=i+1; j<=m; j++) {
			if (b[j][t])
				for (k=t; k<=n; k++)
					b[j][k]^=b[i][k];
		}
		i++;
	}
	ans=n-i+1;
}
int c[101];
void pow() {
	int i,j;
	if (ans==0) {
		printf("0\n");
		return;
	}
	c[0]=1;
	c[1]=2;
	for (i=1; i<ans; i++) {
		for (j=1; j<=c[0]; j++)
			c[j]<<=1;
		for (j=1; j<=c[0]; j++) {
			c[j+1]+=c[j]/10;
			c[j]%=10;
		}
		while (c[c[0]+1]) c[0]++;
	}
	c[1]--;
	for (i=1; i<=c[0]; i++)
		if (c[i]<0) {
			c[i+1]--;
			c[i]+=10;
		}
	while (c[c[0]]==0&&c[0]) c[0]--;
	for (i=c[0]; i; i--) printf("%d",c[i]);
}
int main() {
	int i,x;
	Eular();
	scanf("%d%d",&m,&n);
	for (i=1; i<=n; i++) {
		scanf("%d",&x);
		fac(i,x);
	}
	Gauss();
	pow();
	return 0;
}