快速幂取模

快速幂取模的的核心代码（x^y % z）

int powMod(int x, int y, int k) {
    int result = 1 % k;
    x %= k;
    while (y != 0) {
        if ((y & 1) == 1) {
            result = (result * x) % k;
        }
        x = (x * x) % k;
        y >>= 1;
    }
    return result;
}

一道清华大学机考题：
描述

N<k时，root(N,k) = N，否则，root(N,k) = root(N’,k)。N’为N的k进制表示的各位数字之和。输入x,y,k，输出root(x^ y,k)的值 (这里 ^ 为乘方，不是异或)，2=<k<=16，0<x,y<2000000000，有一半的测试点里 x^y会溢出int的范围(>=2000000000)

输入描述：

每组测试数据包括一行，x(0<x<2000000000), y(0<y<2000000000), k(2<=k<=16)

输出描述：
输入可能有多组数据，对于每一组数据，root(x^y, k)的值

输入：
4 4 10
输出：
4

答案：

#include <iostream>
using namespace std;
long long root(long long x,long long y,int k){
    long long ans = 1 % k;
    while(y != 0){
        if(y % 2 == 1){
            ans = x * ans%k;
        }
        y = y / 2;
        x = x * x % k;
    }
    return ans;
}
int main(){
    int x , y, k;
    cin >> x >> y >> k;
    int ans = root(x,y,k-1);
    if(ans == 0){
        ans = k - 1;
    }
    cout << ans << endl;
    return 0;
}