[Sdoi2010]古代猪文

最新推荐文章于 2020-10-07 19:13:36 发布

weixin_30527323

最新推荐文章于 2020-10-07 19:13:36 发布

阅读量76

点赞数

原文链接：http://www.cnblogs.com/ZUTTER/p/10371467.html

版权

扩展欧拉定理

\[a^k=a^{k\mod \varphi(p)+\varphi(p)}\ \ \ (mod\ p)\]

让后发现\(\varphi(p)\)不是质数不能求逆元

\[999911659=2\times 3\times 4679\times35617\]

分别用lucas定理求完以后再用中国剩余定理合并一下就行

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define M 1000001
#define LL long long 
using namespace std;

LL n,m,g,A[M],B[M],inv[M],k,ans1,ans2,ans3,ans4,p=999911658;

LL C(LL x,LL y,LL p)
{
    if(y>x) return 0;
    if(x<p && y<p) return A[x]*B[y]%p*B[x-y]%p;
    return C(x/p,y/p,p)*C(x%p,y%p,p)%p;
}

LL solve(LL p)
{
    LL ans=0; A[0]=B[0]=inv[1]=1;
    for(int i=2;i<p;i++) inv[i]=(p-p/i)*inv[p%i]%p;
    for(int i=1;i<p;i++) A[i]=A[i-1]*i%p;
    for(int i=1;i<p;i++) B[i]=B[i-1]*inv[i]%p;
    for(int i=1;i<=m;i++) if(n%i==0) ans=(ans+C(n,i,p)+C(n,n/i,p))%p;
    if(m*m==n) ans=(ans-C(n,m,p)+p)%p;
    return ans;
}

LL ksm(LL x,LL y,LL p)
{
    LL z=1;
    for(;y>1;y>>=1, x=x*x%p) if(y&1) z=z*x%p;
    return z*x%p;
}

int main()
{
    scanf("%lld%lld",&n,&g); m=sqrt(n);
    ans1=solve(2)*(p/2)%p; ans4=solve(35617)*(p/35617)%p*ksm(p/35617,35615,35617)%p; 
    ans3=solve(4679)*(p/4679)%p*ksm(p/4679,4677,4679)%p; ans2=solve(3)*(p/3)%p*(p/3%3)%p; 
    m=(ans1+ans2+ans3+ans4)%p;
    printf("%lld",ksm(g,m+p,p+1));
}

转载于:https://www.cnblogs.com/ZUTTER/p/10371467.html