ZOJ 2647 欧拉函数

最新推荐文章于 2021-03-17 19:06:12 发布

Gipsyu

最新推荐文章于 2021-03-17 19:06:12 发布

阅读量138

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/Gipsy_Danger/article/details/79853847

版权

数论专栏收录该内容

49 篇文章 0 订阅

订阅专栏

另一个版本的BZOJ 3884

AC代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll quickmod(ll a,ll b,ll m){
    ll ans=1;
    while(b){
        if(b&1) ans=(ans*a)%m;
        b>>=1;
        a=a*a%m;
    }
    return ans;
}
ll phi(ll x){
    ll ans=x;
    for(ll i=2;i*i<=x;i++){
        if(x%i==0){
            ans=ans/i*(i-1);
            while(x%i==0) x/=i;
        }
    }
    if(x>1) ans=ans/x*(x-1);
    return ans;
}
ll dfs(ll p,ll x){
    if(x==1)
        return 0;
    if(x==p)
        return 0;
    ll y=0;
    while(x%p==0){
        x/=p;
        ++y;
    }
    return quickmod(p,(dfs(p,phi(x))-y%phi(x)+phi(x))%phi(x),x)*pow(p,y);
}
int main(){
    ll p,m;
    int flag=0;
    while(cin>>p>>m){
    	if(flag){
    		cout<<endl;	
		}
		flag=1;
        ll x=m;
        for(ll i=1;i<x;++i)
            m=m*i;
        cout<<dfs(p,m)%m<<endl;
    }
}