阶乘取模预处理

int fact[MAX_P];//预处理n! mod p的表 O(p)

//分解n!=ap^e,返回a mod p O(log_p n)
int mod_fact(int n,int p,int& e)
{
    e=0;
    if(n==0) return 1;

    //计算p的倍数的部分
    int res=mod_fact(n/p,p,e);
    e+=n/p;

    //由于(p-1)!=-1,因此(p-1)!^(n/p)只需要知道n/p的奇偶性就可以计算了。
    if(n/p%2!=0) return res*(p-fact[n%p])%p;
    return res*fact[n%p]%p;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值