阶乘取模预处理

来源:互联网 发布:傲剑弓箭升级数据 编辑:程序博客网 时间:2024/06/09 19:00
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;}
原创粉丝点击