ACM新手DAY 21 莫比乌斯反演&欧拉函数

题解

I - Super A^B mod C

给定A，B，C，您应该快速计算$A^{B}modC$的结果$(1<=A，C<=10000000000，1<=B<=10^{1000000})$。

• 欧拉降幂
• 先用欧拉降幂公式¹处理指数部分，然后快速幂求解。输入部分也要用字符串（或者数组）储存b的值。
• 对应代码

//i题，欧拉降幂+快速幂
#include<iostream>
#include<string.h>
using namespace std;
#define ll long long
int eular(int n)           //求单个数的欧拉函数
{
    int ans = n;
    for(int i = 2; i*i <= n; i++)
    {
        if(n % i == 0)
        {
            ans = ans/i*(i-1);       //ans=ans*(1-1/i)=ans*(i-1)/i
            //        =ans/i*(i-1) //先除防止溢出
            while(n % i == 0)      //消除i因子
                n /= i;
        }
    }
    if(n > 1)ans = ans/n*(n-1);    //最后可能还剩下一个质因数没有除
    return ans;
}
//快速幂部分
ll solve(ll a,ll b,ll c)
{
    ll ans=1;
    while(b)
    {
        if(b&1)
            ans=ans*a%c;
        a=a*a%c;
        b>>=1;
    }
    return ans % c;
}
ll a,c;
string b;
int main()
{
    while(cin >> a>> b >> c)
    {
        ll cc=eular(c);//首先获得m也就是c对应的值
        ll sum=0;
        for(int i=0; i<b.size(); i++)
        {
            sum=(sum*10+b[i]-'0') % cc;
        }
        sum += cc;//欧拉降幂部分
        cout<<solve(a,sum,c)<<endl;
    }
    return 0;
}

---

1. 欧拉降幂公式：$a^b\%m=a^{b\%\varphi (m)+\varphi (m)}\%m$ ↩︎