题意:
有n个数,选任意个数,使得乘积为一个完全平方数。
解题思路:
题目保证这n个数的素因数不超过2000,也就是其素数的个数不超过
303
,构造一个
303∗n
的
01
矩阵,进行高斯消元
ans=2
^
(n−秩)−1
。
矩阵的构造:分解每一个数,统计其素因数出现的次数,奇数个在相应位置记为
1
,否则记为
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<queue>
#include<map>
#include<vector>
using namespace std;
typedef long long ll;
const ll mod=1000000007;
int prime[10000];//存储了小于n的素数
bool visit[10000];//数组标记了是否为素数
//返回值num表示小于n的素数的个数
ll bigpow(int x, int y,ll k)
{
ll ret = 1;
ll tmp = x%k;
while (y > 0)
{
if (y & 1)
ret =ret* tmp%k;
y >>= 1;
tmp =tmp* tmp%k;
}
return ret%k;
}
int init_prim(int n)
{
memset(visit, true, sizeof(visit));
int num = 0;
for (int i = 2; i <= n; ++i)
{
if (visit[i] == true)
{
num++;
prime[num] = i;
}
for (int j = 1; ((j <= num) && (i * prime[j] <= n)); ++j)
{
visit[i * prime[j]] = false;
if (i % prime[j] == 0) break; //点睛之笔
}
}
return num;//返回质数个数
}
const int maxn=350;
int equ,var;//等式数,变量数
int a[maxn][maxn],x[maxn];//矩阵(多一列最终状态),答案
int free_x[maxn],free_num;//自由变元
int Gauss()
{
int max_r,col,k;
free_num=0;
for(k=0,col=0;k<equ&&col<var;k++,col++)
{
max_r=k;
for(int i=k+1;i<equ;i++)
{
if(abs(a[i][col])>abs(a[max_r][col]))
max_r=i;
}
if(a[max_r][col]==0)
{
k--;
free_x[free_num++]=col;
continue;
}
if(max_r!=k)
{
for(int j=col;j<var+1;j++)
swap(a[k][j],a[max_r][j]);
}
for(int i=k+1;i<equ;i++)
{
if(a[i][col]!=0)
{
for(int j=col;j<var+1;j++)
a[i][j]^=a[k][j];
}
}
}
for(int i=k;i<equ;i++)
if(a[i][col]!=0)
return -1;
if(k<var) return var-k;
for(int i=var-1;i>=0;i--)
{
x[i]=a[i][var];
for(int j=i+1;j<var;j++)
x[i]^=(a[i][j]&&x[j]);
}
return 0;
}
int main()
{
int cas=0;
int T;
scanf("%d",&T);
int num=init_prim(2000);
while(T--)
{
int n;
ll m,k;
scanf("%d",&n);
equ=303,var=n;//等式数,变量数
memset(a,0,sizeof(a));
for(int i=0;i<n;i++)
{
scanf("%I64d",&m);
k=m;
for(int j=1;prime[j]<=k&&j<=303;j++)/// 构造01矩阵
{
int nn=0;
while(k%prime[j]==0)
{
k/=prime[j];
nn++;
}
if(nn&1)
a[j-1][i]=1;
}
}
int zhi=Gauss(); ///返回值为n-秩
ll ans=(bigpow(2,zhi,mod)-1+mod)%mod;
printf("Case #%d:\n%I64d\n",++cas,ans);
}
return 0;
}