hdu 5833 Zhu and 772002 异或方程组高斯消元

ccpc网赛卡住的一道题

蓝书上的原题 但是当时没看过蓝书

今天又找出来看看 其实也不是特别懂 但比以前是了解了一点了

主要还是要想到构造异或方程组 异或方程组的消元只需要xor就好搞了

数学真的是硬伤啊……

（链接：蓝书161页详细讲解 我也在看……

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 #include<cstring>
 #include<string>
 #define cl(a,b) memset(a,b,sizeof(a))
 #define debug(x) cerr<<#x<<"=="<<(x)<<endl
 using namespace std;
 typedef long long ll;
 
 const int maxn=3e2+;
 const int maxx=2e3+;
 const int mod=1e9+;
 
 int n,maxp,tol;
 int prime[maxn];
 int A[maxn][maxn];
 
 ll powmod(ll a,ll x)
 {
     ll t=;
     while(x)
     {
         if(x&) t=t*a%mod;
         a=a*a%mod;
         x>>=;
     }
     return t;
 }
 
 void init()
 {
     bool notprime[maxx];
     cl(notprime,false);
     notprime[]=notprime[]=true;
     for(int i=; i<maxx; i++)
     {
         if(!notprime[i])
         {
             if(i*i>maxx) continue;
             for(int j=i*i; j<maxx; j+=i)
                 notprime[j]=true;
         }
     }
     tol=;
     for(int i=; i<maxx; i++)
         if(!notprime[i]) prime[tol++]=i;
 }
 
 ll Gauss(int m)
 {
     int i=,j=,r;
     while(j<m&&i<n)
     {
         r=i;
         for(int k=i; k<m; k++)
             if(A[k][j]){r=k;break;}
         if(A[r][j]){
             if(r!=i) for(int k=; k<=n; k++) swap(A[r][k],A[i][k]);
             for(int u=i+; u<m; u++) if(A[u][j])
                     for(int k=i; k<=n; k++) A[u][k]^=A[i][k];
             i++;
         }
         j++;
     }
     return i;
 }
 
 int main()
 {
     int T,cas=;
     scanf("%d",&T);
     init();
     while(T--)
     {
         int maxp=;
         ll x;
         scanf("%d",&n);
         memset(A,,sizeof(A));
         for(int i=; i<n; i++)
         {
             scanf("%lld",&x);
             for(int j=; j<tol; j++)
             {
                 while(x%prime[j]==)
                 {
                     maxp=max(maxp,j);
                     x/=prime[j];
                     A[j][i]^=;
                 }
             }
         }
         ll r=Gauss(maxp+);
         ll tmp=n-r;
         ll ans=powmod(,tmp)-;
         printf("Case #%d:\n",cas++);
         printf("%lld\n",ans);
     }
     return ;
 }
 
 /*
 
 3
 4
 4 6 10 15
 3
 3 3 4
 3
 2 2 2
 
 */