codeforce338 D. Multipliers 费马小定理降幂

codeforces338 D. Multipliers

给定组成一个数的质因子序列,求其所有因子的乘积

一个质因子p的贡献为(p^1p^2p^3...p^c)*(c1+1)(c2+1)...(cn+1)/(c+1)

化简得                    \prod_{i=1}^{m}(p^{(1+c)*c/2})^{mul/(c+1)}=\prod_{i=1}^{m}p^{c*mul/2}

由于c*mul/2可能为一个大数,根据费马小定理 p^{x\%(mod-1)} \equiv 1(mod)进行降所以对

所以对c*mul/2mod-1但是mul可能很大,需要在计算过程中不断取模,

mul\%(mod-1)*c/2\ =?=\ c*mul/2\%(mod-1)

因为除法不满足取模性质,所以不可以直接外提,且mod-1不是质数,2不存在逆元,考虑将模数扩大2倍消除其影响mod*2-2

#include<bits/stdc++.h>
using namespace std;
const int MAX=2e5+5;
const long long MOD=1e9+7;
int n,a[MAX];
map<int,long long>mp;
map<int,long long>::iterator p;
long long pow(long long a,int k,long long mod)
{
    long long b=1;
    while(k)
    {
        if(k&1) b=b*a%mod;
        a=a*a%mod;
        k>>=1;
    }
    return b;
}
int main()
{
    long long x;
    while(~scanf("%d",&n))
    {
        mp.clear();
        for(int i=0;i<n;++i)
        {
            scanf("%d",&x);
            if(mp.find(x)!=mp.end()) ++mp[x];
            else mp[x]=1;
        }
        long long x=1;
        for(p=mp.begin();p!=mp.end();++p)
            x=x*(p->second+1)%(MOD*2-2);
        long long mul=1;
        for(p=mp.begin();p!=mp.end();++p)
        {
            mul=mul*pow(p->first,(x*p->second/2)%(MOD-1),MOD)%MOD;
        }
        printf("%lld\n",mul);
    }
    return 0;
}

 

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

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值