跳到主要內容

uva 11542

壓根兒想不到是高斯,首次學習高斯竟然可以這樣用..

//
// GGGGGGGGGGGGG CCCCCCCCCCCCC AAA
// GGG::::::::::::G CCC::::::::::::C A:::A
// GG:::::::::::::::G CC:::::::::::::::C A:::::A
// G:::::GGGGGGGG::::G C:::::CCCCCCCC::::C A:::::::A
// G:::::G GGGGGG C:::::C CCCCCC A:::::::::A
//G:::::G C:::::C A:::::A:::::A
//G:::::G C:::::C A:::::A A:::::A
//G:::::G GGGGGGGGGGC:::::C A:::::A A:::::A
//G:::::G G::::::::GC:::::C A:::::A A:::::A
//G:::::G GGGGG::::GC:::::C A:::::AAAAAAAAA:::::A
//G:::::G G::::GC:::::C A:::::::::::::::::::::A
// G:::::G G::::G C:::::C CCCCCC A:::::AAAAAAAAAAAAA:::::A
// G:::::GGGGGGGG::::G C:::::CCCCCCCC::::C A:::::A A:::::A
// GG:::::::::::::::G CC:::::::::::::::C A:::::A A:::::A
// GGG::::::GGG:::G CCC::::::::::::C A:::::A A:::::A
// GGGGGG GGGG CCCCCCCCCCCCCAAAAAAA AAAAAAA
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <climits>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <cctype>
#include <utility>
#include <ctime>
using namespace std;
#ifdef DEBUG
#define VAR(a,b) decltype(b) a=(b)
#define debug(...) printf("DEBUG: "),printf(__VA_ARGS__)
#define gettime() end_time=clock();printf("now running time is %.7f\n",(float)(end_time - start_time)/CLOCKS_PER_SEC);
#else
#define VAR(a,b) __typeof(b) a=(b)
#define debug(...)
#define gettime()
#endif
typedef unsigned int uint;
typedef long long int Int;
typedef unsigned long long int UInt;
#define Set(a,s) memset(a,s,sizeof(a))
#define Write(w) freopen(w,"w",stdout)
#define Read(r) freopen(r,"r",stdin)
#define Pln() printf("\n")
#define I_de(x,n)for(int i=0;i<n;i++)printf("%d ",x[i]);Pln()
#define De(x)printf(#x"%d\n",x)
#define For(i,x)for(int i=0;i<x;i++)
#define CON(x,y) x##y
#define Pmz(dp,nx,ny)for(int hty=0;hty<ny;hty++){for(int htx=0;htx<nx;htx++){\
printf("%d ",dp[htx][hty]);}Pln();}
#define M 505
#define PII pair<int,int>
#define PB push_back
#define oo INT_MAX
#define Set_oo 0x3f
#define FOR(a,b) for(VAR(a,(b).begin());a!=(b).end();++a)
#define eps 1e-6
#define X first
#define Y second
clock_t start_time=clock(), end_time;
bool xdy(double x,double y){return x>y+eps;}
bool xddy(double x,double y){return x>y-eps;}
bool xcy(double x,double y){return x<y-eps;}
bool xcdy(double x,double y){return x<y+eps;}
int min3(int x,int y,int z){
int tmp=min(x,y);
return min(tmp,z);
}
int max3(int x,int y,int z){
int tmp=max(x,y);
return max(tmp,z);
}
int mz[105][300];
int seize[M];
int prime[M],pnum;
void pre(){
pnum=0;
for(int i=2;i<M;i++){
if(!seize[i]){
prime[pnum++]=i;
for(int j=i+i;j<M;j+=i){
seize[j]=1;
}
}
}
}
int n,m;
Int x[105];
int gaus_xor(){
int i,j;
i=j=0;
while(i<n&&j<m){
// for(int jj=0;jj<n;jj++){
// for(int kk=0;kk<m;kk++)printf("%d ",mz[kk][jj]);
// Pln();
// }Pln();
int r=i;
for(int k=i;k<n;k++){
if(mz[j][k]){
r=k;
break;
}
}
if(mz[j][r]){
if(r!=i){//change
for(int k=0;k<m;k++)swap(mz[k][r],mz[k][i]);
}
for(int k=i+1;k<n;k++){
if(mz[j][k]){
for(int a=j;a<m;a++){
mz[a][k]^=mz[a][i];
}
}
}
i++;
}
j++;
}
return n-i;
}
int main() {
ios_base::sync_with_stdio(0);
int test;
pre();
scanf("%d",&test);
while(test--){
m=0;
Set(mz,0);
scanf("%d",&n);
for(int i=0,j;i<n;i++){
scanf("%lld",&x[i]);
Int t=x[i];
for(j=0;j<pnum&&t!=1;j++){

int exp=0;
while(t%prime[j]==0){
t/=prime[j];
m=max(m,j);
exp++;
// debug("%d %lld\n",prime[j],x[i]);
}
mz[j][i]=exp&1;
// debug("%d %lld %d %d\n",i,x[i],j,mz[i][j]);
}
}
// debug("%d\n",m);
m++;
printf("%lld\n",(1LL<<gaus_xor())-1);
}















}

留言