首页 > 代码库 > UVA 11542 Square ——线性基
UVA 11542 Square ——线性基
【题目分析】
每个数没有超过500的因子。很容易想到把每一个数表示成一个二进制的数。
(0代表该质数的次数为偶数,1代表是奇数)
然后问题转化成了选取一些二进制数,使他们的异或和为0。
高斯消元,2^(自由元)即为答案,需要把空集的情况减去,所以减一。
然而发现并不需要知道哪些是自由元,所以只需要用线性基去维护即可。
然后代码就呼之欲出了。
【代码】
#include <cstdio> #include <cstring> #include <cmath> #include <cstdlib> #include <map> #include <set> #include <queue> #include <string> #include <iostream> #include <algorithm> using namespace std; #define maxn 500005 #define ll long long #define inf 0x3f3f3f3f #define F(i,j,k) for (int i=j;i<=k;++i) #define D(i,j,k) for (int i=j;i>=k;--i) void Finout() { #ifndef ONLINE_JUDGE freopen("in.txt","r",stdin); // freopen("out.txt","w",stdout); #endif } int Getint() { int x=0,f=1; char ch=getchar(); while (ch<‘0‘||ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();} while (ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-‘0‘; return x*f; } ll Getll() { ll x=0,f=1; char ch=getchar(); while (ch<‘0‘||ch>‘9‘) {if (ch==‘-‘) f=-1; ch=getchar();} while (ch>=‘0‘&&ch<=‘9‘) x=x*10+ch-‘0‘; return x*f; } int pri[205],top; void init() { F(i,2,500) { int flag=1; F(j,2,sqrt(i)) if (i%j==0) flag=0; if (flag) pri[++top]=i; } // F(i,1,top) cout<<pri[i]<<" "; // cout<<"over"<<endl; } int t,a[205][205],n,cnt; int lb[205][205],hav[205]; int main() { Finout(); init(); scanf("%d",&t); // cout<<t<<endl; while (t--) { // cout<<"test "<<t<<endl; memset(hav,0,sizeof hav); memset(lb,0,sizeof lb); memset(a,0,sizeof a); cnt=0; scanf("%d",&n); F(i,1,n) { ll x; scanf("%lld",&x); // cout<<"now is "<<x<<endl; F(j,1,top) { while (x%pri[j]==0) { a[i][j]=!a[i][j]; x/=pri[j]; // cout<<j<<" ^ 1"<<endl; } } // cout<<i<<": "; // D(j,top,1) cout<<a[i][j]; cout<<endl; } F(i,1,n) { int flag=0; D(j,top,1) { if (a[i][j]) { // cout<<"have in "<<j<<endl; if (!hav[j]) { // cout<<"put a"<<endl; // D(k,top,1) cout<<a[i][k]; cout<<endl; D(k,j,1) lb[j][k]=a[i][k]; hav[j]=1; flag=1; break; } else { // cout<<"star to ^"<<endl; D(k,j,1) a[i][k]^=lb[j][k]; // D(k,top,1) cout<<a[i][k]; cout<<endl; } } } if (!flag) cnt++; } cout<<(1LL<<cnt)-1<<endl; } }
UVA 11542 Square ——线性基
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。