【jzoj4668】【腐败】【数论】【快速乘】

题目大意

给出一个序列，求两两gcd的积。

解题思路

可以发现每个质数对答案贡献是独立的，分解质因数后发现数是无序的，排序后可以用前缀和求答案，排序可以用计数排序，快速幂可以放在分解质因数的循环里，还有模数较大，要用快速乘，类似快速幂。

code

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
using namespace std;
int const maxn=30000,maxa=10000000;
LL const mod=100000000009;
int n,ss[maxn+10],a[maxn+10],b[maxn+10],cnt[50+10];
LL mult(LL x,LL y){
    if(x<y)swap(x,y);
    LL z=0;
    for(;y;){
        if(y%2==1)z=(z+x)%mod;;
        x=(x+x)%mod;
        y/=2;
    }
    return z;
}
LL Pow(LL x,LL y){
    LL z=1;
    for(;y;){
        if(y%2==1)z=mult(z,x);
        x=mult(x,x);
        y/=2;
    }
    return z;
}
int main(){
    freopen("d.in","r",stdin);
    freopen("d.out","w",stdout);
    scanf("%d",&n);
    fo(i,1,n)scanf("%d",&a[i]);
    fo(i,2,sqrt(maxa)){
        bool ok=1;
        fo(j,2,i-1)
            if(i%j==0){
                ok=0;
                break;
            }
        if(ok)ss[++ss[0]]=i;
    }
    LL ans=1;
    fo(i,1,ss[0]){
        memset(b,0,sizeof(b));
        fo(j,1,n)
            for(;a[j]%ss[i]==0;b[j]++,a[j]/=ss[i]);
        fo(j,1,n)cnt[b[j]]++;
        b[0]=0;
        fo(j,0,50)
            for(;cnt[j];b[++b[0]]=j,cnt[j]--);
        LL sum=0;
        fo(j,1,n)sum+=1ll*b[j]*(n-j+1);
        ans=mult(ans,Pow(ss[i],sum));
    }
    sort(a+1,a+n+1);
    int pre=1;
    fo(i,1,n){
        if(a[i]!=a[i-1])pre=i;
        ans=mult(ans,Pow(a[i],i-pre+1));
    }
    printf("%lld",ans);
    return 0;
}