刘汝佳白皮书160页原题,差不多,今天应该反思下一个月来的成绩和不足,应该朝那个方向努力。
#include<iostream>
#include<stdlib.h>
#include<queue>
#include<string.h>
#include<stdio.h>
using namespace std;
const int maxn=2001;
int prime[maxn],visit[maxn];
int N=0;
int b[2001][310];
long long a[310];
const int mod=1000000007;
void init()
{
for(long long i=2; i<maxn; i++)
{
if(!visit[i])
prime[N++]=i;
for(long long j=0; j<N; j++)
{
long long k=i*prime[j];
if(k>maxn)break;
visit[k]=1;
}
}
}
int mod_pow(int a,int b)
{
int res=1;
a%=mod;
while(b>0)
{
if(b&1)res=(res*a)%mod;
a=(a*a)%mod;
b/=2;
}
return res;
}
int rank1(int m,int n)
{
int k,r,u,i=0,j=0;
while(i<m&&j<n)
{
r=i;
for(k=i;k<m;k++)
if(b[k][j]){r=k;break;}
if(b[r][j])
{
if(r!=i)for(k=0;k<=n;k++)swap(b[r][k],b[i][k]);
for(u=i+1;u<m;u++)
if((b[u][j]))
for(k=i;k<=n;k++)b[u][k]^=b[i][k];
i++;
}
j++;
}
return i;
}
int main()
{
init();
int t,cnt=0;
scanf("%d",&t);
while(t--)
{
printf("Case #%d:\n",++cnt);
memset(b,0,sizeof(b));
int n;
scanf("%d",&n);
int max1=0;
for(int i=1;i<=n;i++)
{scanf("%I64d",&a[i]);
for(int j=0;j<N&&a[i]>=prime[j];j++)
{
if(a[i]%prime[j]==0)
{
long long tmp=0;
while(a[i]%prime[j]==0)
{
tmp++;
a[i]/=prime[j];
}
b[j][i]=tmp%2;
max1=max(max1,j);
}
}
}
int r=rank1(max1+1,n);
r=n-r;
cout<<mod_pow(2,r)-1<<endl;
}
return 0;
}