快速幂||取余运算 （分治算法）

#include<iostream>
using namespace std;

long b,p,k;
long skt=1;
int we,tsm;

int ksm(long b,long p,long k)
{
    while(p)
    {
    if(p%2!=0)
    {
    skt=skt*b%k;//之所以不用p=p-1是因为后面的位运算中p=p>>1，相当于p=p/2，多出的那个1被自动消去了，所以p=p-1可写可不写
    }
    b=b*b%k;//重初始化，新的式子是（b*b mod k）的p/2次方 mod k
    p=p>>1;//相当于p=p/2，当p=1时p位移为0
    }
    return skt%k;
}

int main()
{
    cin>>b>>p>>k;
    cout<<b<<"^"<<p<<" mod "<<k<<"="<<ksm(b,p,k);
}