HDU-5833-Zhu and 772002(高斯消元)

32 篇文章 0 订阅

链接:http://acm.hdu.edu.cn/showproblem.php?pid=5833

题意:

给出n个数,每次可以从中取任何数,问有多少种取法使得取出来的数的乘积是完全平方数。

题解:

首先显然可以去掉偶数幂的质因子,每个数就变成了不同质数相乘的形式,已知题意中只有大约300多个质因子,直接列成质因子的01矩阵形式,高斯消元,求一下自由变元个数temp,显然答案就是temp^1-1。

CODE:

#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
#define bug cout<<"bug\n"
const int maxn = 1e5+7;
#define MAXN 2030
const long long mod = 1e9+7;
bool not_prime[MAXN*30];
long long prime[MAXN];
int num_prime;
int a[350][350];
void get_prime()
{
    long long i,j;
    memset(not_prime,0,sizeof(not_prime));
    num_prime=0;
    not_prime[1]=1;
    for(i=2; i<MAXN; ++i)
        if(!not_prime[i])
        {
            num_prime++;
            prime[num_prime]=i;
            for(j=i*i; j<MAXN; j+=i)
                not_prime[j]=1;
        }
}
int gauss(int m, int n)
{
    int i,j,k,r;
    for(i=j=0; i<m&&j<n ; ++i,++j)
    {
        r=i;
        for(int c=i; c<m; ++c)
        if(a[c][j]){r=c;break;}
        if(a[r][j])
        {
            if(r!=i)for(int c=0; c<=n; ++c)swap(a[r][c],a[i][c]);
            for(k=i+1; k<m; ++k)if(a[k][j])
                for(int c=i; c<=n; ++c)
                    a[k][c]^=a[i][c];
        }
        else --i;
    }
    return i;
}
int main()
{
    int T;
    int n,m;
    long long p;
    scanf("%d",&T);
    int cas=1;
    get_prime();
    while(T--)
    {
        memset(a,0,sizeof(a));
        scanf("%d",&n);
        int M=0;
        for(int i=0; i<n; ++i)
        {
            scanf("%I64d",&p);
            for(int j=1; p>=prime[j] && j<=num_prime; ++j)
            {
                while(p%prime[j]==0)
                {
                    M=max(M,j);
                    p/=prime[j];
                    ++a[i][j-1];
                }
                a[i][j-1]%=2;
            }
        }
        long long temp=n-gauss(n,M);
        long long ans=1;
        for(int i=0; i<temp; ++i){ans<<=1LL;ans%=mod;}
        printf("Case #%d:\n%I64d\n",cas++,ans-1);
    }
    return 0;
}
/*
100
3
3 3 4
3
2 2 2
5
4 4 4 4 4
9
4 4 4 4 4 5 5 5 5
6
4 4 80 20 5 5
*/

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值