HDU 5728 PowMod

具体推导过程可以看这篇博客..我是看这个看懂的，其余的人都是直接给的结果..

http://blog.csdn.net/wust_zzwh/article/details/51966450

分两步：1.求k，2.求答案。

对于第二步：即bzoj 3884，http://blog.csdn.net/skywalkert/article/details/43955611 看了这个就会了。

对于第一步：请看http://blog.csdn.net/wust_zzwh/article/details/51966450 详细推导过程。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<map>
#include<set>
#include<queue>
#include<stack>
#include<iostream>
using namespace std;
typedef long long LL;

const LL mod=1000000007;
const int maxn=10000010;
int phi[maxn],prime[maxn],sz;
bool check[maxn];
LL sumPhi[maxn],a[maxn],cnt;
LL n,m,P,k,ans;

LL pow(LL a,LL b,LL p)
{
    LL res=1; a%=p;
    while(b)
    {
        if(b&1) res=(res*a)%p;
        a=(a*a)%p, b/=2;
    }
    return res;
}

LL f(LL p)
{
    if(p==2) return k%2;
    return pow(k,f(phi[p])+(LL)phi[p],p);
}

LL get(int pos,LL x,LL y)
{
    if(x==1) return sumPhi[y];
    if(y==1) return phi[x];
    if(y<1) return 0;
    return (phi[a[pos]]*get(pos-1,x/a[pos],y)%mod+get(pos,x,y/a[pos])%mod)%mod;
}

void Eorue()
{
    phi[1]=1;
    sz=0;
    for(int i=2;i<maxn;i++)
    {
        if(!check[i]) { phi[i]=i-1; prime[sz++]=i; }
        for(int j=0;j<sz;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]+(LL)phi[i])%mod;
}

void fun(LL x)
{
    cnt=0;
    for(int i=0;i<sz;i++)
    {
        if(x==1) break;
        if(x%(LL)prime[i]==0)
        {
            a[cnt++]=(LL)prime[i];
            x=x/(LL)prime[i];
        }
    }
}

int main()
{
    Eorue();
    while(~scanf("%lld%lld%lld",&n,&m,&P))
    {
        fun(n);//质因数分解
        k=get(cnt-1,n,m);//计算k
        ans=f(P);
        printf("%lld
",ans);
    }
    return 0 ;
}