uva11542(高斯消元)

题意:
给定一些数字,保证这些数字质因子不会超过500,求这些数字中选出几个,乘积为完全平方数,问有几种选法


思路:

对每个数字分解成质因子后,发现如果要是完全平方数,选出来的数字的每个质因子个数都必然要是偶数,这样每个质因子可以列出一个异或的方程,如果数字包含质因子,就是有这个未知数,然后进行高斯消元,求出自由变量的个数,每个自由变量可以选或不选,这样的情况就是(2^个数),然后在扣掉什么都不选的1种就是答案了


#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
using namespace std;

const ll N = 510;

int t, n, a[N][N], Max, vis[N], pn = 0;
ll prime[N];

void get_prime() {
	for (ll i = 2; i < N; i++) {
		if (vis[i]) 
			continue;
		prime[pn++] = i;
		for (ll j = i * i; j < N; j += i)
			vis[j] = 1;
	}
}

int gauss() {
	int i = 0, j = 0;
	while (i <= Max && j < n) {
		int k = i;
		for (; k <= Max; k++)
			if (a[k][j]) break;
		if (k != Max + 1) {
			for (int l = 0; l <= n; l++)
				swap(a[i][l], a[k][l]);
			for (int k = i + 1; k <= Max; k++) {
				if (a[k][j]) {
					for (int l = j; l <= n; l++)
						a[k][l] ^= a[i][l];
				}
			}
			i++;
		}
		j++;
	}
	return n - i;
}

int main() {
	get_prime();
	scanf("%d", &t);
	while (t--) {
		scanf("%d", &n);
		ll x;
		Max = 0;
		memset(a, 0, sizeof(a));
		for (int i = 0; i < n; i++) {
			scanf("%lld", &x);
			for (int j = 0; j < pn && prime[j] <= x; j++) {
				while (x % prime[j] == 0) {
					a[j][i] ^= 1;
					Max = max(Max, j);
					x /= prime[j];
				}
			}
		}
		printf("%lld\n", (1LL << (gauss())) - 1);
	}
	return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值