标签：return int mo x% SDOI2010 base 1951 猪文 jc

数论大杂烩，留坑代填，先放代码：

/**************************************************************
    Problem: 1951
    User: JBLee
    Language: C++
    Result: Accepted
    Time:152 ms
    Memory:2856 kb
****************************************************************/
 
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int maxn=1e5+7;
const int mod=999911658;
int n,g;
int sum;
int m[5]={0,2,3,4679,35617};
int jc[maxn];
int x,y;
int china[maxn];
int ksm(int x,int n,int mo){
    int base=1;
    while(n){
        if(n&1) base=base*x%mo;
        x=x*x%mo;
        n>>=1; 
    }
    return base;
}
int exgcd(int a,int &x,int b,int &y){
    if(!b){
        x=1;y=0;
        return a;
    }
    int gcd=exgcd(b,x,a%b,y);
    int tmp=x;x=y;y=tmp-(a/b)*y;
    return gcd;
}
int inv(int a,int mo){
    exgcd(a,x,mo,y);
    return (x%mo+mo)%mo;
}
void work(int p){
    jc[0]=1,jc[1]=1;
    for(int i=2;i<=p;i++) jc[i]=jc[i-1]*i%p;
}
int comb(int n,int m,int p){
    if(m>n) return 0;
    return jc[n]*inv(jc[n-m],p)%p*inv(jc[m],p)%p;
}
int lucas(int n,int m,int p){
    if(!m) return 1;
    else return comb(n%p,m%p,p)*lucas(n/p,m/p,p)%p;
}
int solve(int p){
    work(p);
    int ans=0;
    for(int i=1;i*i<=n;i++){
        if(n%i==0){
            ans=(ans+lucas(n,i,p))%p;
            if(n/i!=i) ans=(ans+lucas(n,n/i,p))%p;
        }
    }
    return ans;
}
int china_mod(){
    for(int i=1;i<=4;i++) china[i]=solve(m[i]);
    int lcm=1,ans=0,mi;
    for(int i=1;i<=4;i++) lcm*=m[i];
    for(int i=1;i<=4;i++){
        mi=lcm/m[i];
        exgcd(mi,x,m[i],y);
        ans=((ans+mi*x*china[i])%lcm+lcm)%lcm;
    }
    return ans;
} 
signed main(){
    scanf("%lld%lld",&n,&g);
    if(g%(mod+1)==0){
        printf("0\n");
        return 0;
    }
    printf("%lld\n",ksm(g,china_mod(),mod+1));
    return 0;
} 

标签：return,int,mo,x%,SDOI2010,base,1951,猪文,jc
来源： https://www.cnblogs.com/LJB666/p/11779219.html

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