ZOJ 2674 Strange Limit

欧拉函数。

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define MAXSIZE 100003
#define N 100003
using namespace std;

long long prime[];
long long flag[];
long long cn;

void sushu()
{
    memset(flag,,sizeof(flag));
    cn=;
    long long i,j;
    for (i=;i<MAXSIZE+;i++)
    {
        if(!flag[i])
             prime[cn++]=i;
        for(j=;(j<cn)&&(prime[j]*i<MAXSIZE+);j++)
        {
            flag[i*prime[j]]=;
        }

    }
}
long long fast(long long a,long long p,long long m)
{
    long long r,k;
    if(p==)
        return %m;
    else if(p==)
        return a%m;
    r=a%m;
    k=;
    while(p>)
    {
        if((p&)!=)
            k=(k*r)%m;
        r=(r*r)%m;
        p>>=;
    }
       return (r*k)%m;
}
long long fac(long long m)
{
    long long i,res=;
    for(i=m;i>=;i--)
        res*=i;
    return res;
}
long long ola(long long n)
{
    long long i,n1=n;
    long long res=n;
    for(i=;(i<cn)&&prime[i]*prime[i]<=n;i++)
    {
        if(n%prime[i]==)
        {
            res=res/prime[i]*(prime[i]-);
            while(n1%prime[i]==)
                n1/=prime[i];
        }
    }
    if(n1!=)
        res=res/n1*(n1-);
    return res;
}
long long solve(long long a,long long m)
{
    if(m==)
        return ;
    else
    {
        long long tmp=ola(m);           //phi[m];     //euler(m);
        return fast(a,tmp,m)*fast(a,solve(a,tmp),m);
    }
}

int main()
{
    long long m,m1;
    long long p;
    sushu();
    bool tag=false;
    while(scanf("%lld%lld",&p,&m1)!=EOF)
    {
        if(tag)
        {
            puts("");
        }
        else tag=true;
        m=fac(m1);
        printf("%lld\n",solve(p,m)%m);
    }
    return ;
}