标签：md ch 记录 long include WR ans 排列组合 猪文

这就是传说中的数论全家桶么……

用到了Lucas定理，中国剩余定理，欧拉定理

大体而言可以拆数发现999911658=2*3*4679*35617

然后跑四遍CRT即可

#include<cstdio>
#include<cstring>
#include<string>
#include<cmath>
#define WR WinterRain
using namespace std;
const long long WR=10010000,mod=999911659;
long long tot,bs,val;
long long g[WR],b[5]={0,2,3,4679,35617},a[5];
long long read(){
    long long s=0,w=1;
    char ch=getchar();
    while(ch>'9'||ch<'0'){
        if(ch=='-') w=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9'){
        s=(s<<3)+(s<<1)+ch-48;
        ch=getchar();
    }
    return s*w;
}
long long quick_pow(long long a,long long b,long long md){
    long long ans=1,base=a;
    while(b>0){
        if(b&1) ans=ans*base%md;
        base=base*base%md;
        b>>=1;
    }
    return ans;
}
long long comb(long long n,long long m,long long md){
    if(n<m) return 0;
    return g[n]*quick_pow(g[m],md-2,md)*quick_pow(g[n-m],md-2,md);
}
long long Lucas(long long n,long long m,long long md){
    if(!n) return 1;
    if(n<m) return 0;
    return comb(n%md,m%md,md)*Lucas(n/md,m/md,md)%md;
}
void CRT(long long n,long long md){
    for(long long i=1;i<=n;i++){
        val=(val+a[i]*(md/b[i])%md*quick_pow(md/b[i],b[i]-2,b[i])%md)%md;
    }
}
int main(){
    tot=read(),bs=read();
    g[0]=1;
    /*long long tmp=mod-1;
    for(long long i=2;i<=sqrt(tmp);i++){
        if(tmp%i==0){
            bool f=true;
            for(long long j=2;j<=sqrt(i);j++){
                if(i%j==0){
                    f=false;
                    break;
                }
            }
            if(f){
                prlong longf("%d\n",i);
                tmp/=i;
            }
        }
    }
    printf("%lld",tmp);*/
    for(long long i=1;i<=4;i++){
        for(long long j=1;j<=b[i];j++) g[j]=g[j-1]*j%b[i];
        for(long long j=1;j<=sqrt(tot);j++){
            if(tot%j==0){
                a[i]=(a[i]+Lucas(tot,j,b[i]))%b[i];
                if(j*j!=tot) a[i]=(a[i]+Lucas(tot,tot/j,b[i]))%b[i];
            }
        }
    }
    CRT(4,mod-1);
    printf("%lld",quick_pow(bs,val,mod));
    return 0;
}

标签：md,ch,记录,long,include,WR,ans,排列组合,猪文
来源： https://www.cnblogs.com/WintersRain/p/16247827.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。