hdu5728PowMod+欧拉函数和+k的无穷次方取膜

最新推荐文章于 2020-03-17 17:06:06 发布

xtulollipop

最新推荐文章于 2020-03-17 17:06:06 发布

阅读量878

点赞数

CC 4.0 BY-SA版权

分类专栏：数学问题

本文链接：https://blog.youkuaiyun.com/xtulollipop/article/details/52067221

数学问题专栏收录该内容

64 篇文章

订阅专栏

本文介绍了一种计算特定形式幂次表达式的模运算的方法，利用欧拉函数的性质来简化计算过程，并通过素数筛选和快速幂运算提高效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

Problem Description
Declare:
$k=\sum_{i=1}^{m}φ(i∗n) mod 1000000007$

n is a square-free number.

φ is the Euler’s totient function.

find:
ans= $k^{k^{k^{k^{...^{k}}}}}$ mod p

There are infinite number of k

Input
Multiple test cases(test cases ≤100), one line per case.

Each line contains three integers, n,m and p.

1≤n,m,p≤107

Output
For each case, output a single line with one integer, ans.

Sample Input

1 2 6
1 100 9

Sample Output

4
7

Author
HIT

Source
2016 Multi-University Training Contest 1

这里写图片描述

证明参考这位菊苣的吧http://blog.youkuaiyun.com/wust_zzwh/article/details/51966450
后边那个k的无穷个k的那个参考八中oj3884
http://blog.youkuaiyun.com/xtulollipop/article/details/52065860

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<iostream>
#include<map>
#include<queue>
using namespace std;

#define LL long long

const int mod=1000000007;
const int maxn=10000005;
bool check[maxn];          //用于打表记录的中间量
LL sumPhi[maxn];           //前i个的欧拉函数和
int cnt,phi[maxn],prime[maxn];  //素数个数，欧拉表，素数表
//素数表是第几个素数是什么，欧拉表是i的欧拉是phi[i];
void init(){               //素数+欧拉表
    phi[1]=1;
    cnt=0;
    for(int i=2;i<maxn;i++){
        if(!check[i]){
            phi[i]=i-1;
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt;j++){
            if(i*prime[j]>=maxn)break;
            check[i*prime[j]]=true;
            if(i%prime[j]==0){
                phi[i*prime[j]]=phi[i]*prime[j];
                break;
            }
            else{
                phi[i*prime[j]]=phi[i]*(prime[j]-1);
            }
        }
    }
    sumPhi[0]=0;
    for(int i=1;i<maxn;i++) sumPhi[i]=(sumPhi[i-1]+phi[i])%mod;
}
LL phiz(int n,int m){       //计算k，
    if(n==1) return sumPhi[m];
    if(m==1) return phi[n];
    if(m<1) return 0;

    for(int i=0;i<cnt;i++){
        if(n%prime[i]==0){  //n的最小质数prime[i]
            LL temp1=(prime[i]-1)*phiz(n/prime[i],m)%mod;
            LL temp2=phiz(n,m/prime[i])%mod;
            return (temp1+temp2)%mod;
        }
    }
    return 0;
}
LL pow(LL a,LL b,int p){  //快速幂
    LL ret=1;
    a%=p;
    while(b){
        if(b&1) ret=(ret*a)%p;
        a=(a*a)%p;
        b/=2;
    }
    return ret;
}
LL powe(LL k,int p){  //k^b%p
    if(p==1) return 0;
    LL temp=powe(k,phi[p]);
    return pow(k,temp+phi[p],p);
}
int main(){
    init();
    int n,m,p;
    while(scanf("%d %d %d",&n,&m,&p)!=EOF){
        LL k=phiz(n,m);
        printf("%I64d\n",powe(k,p));
    }
    return 0;
}