hdu 5833 Zhu and 772002

刘汝佳白皮书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;    }