快速幂只能对前三个 求大佬来看一看哪里错了

#include<iostream>
using namespace std;
long long  b,p,k,ans;
long long b1,p1;
int main()
{
    cin>>b>>p>>k;
    ans=1;
    b1=b;
    p1=p;
    while(p>0)
    {
    if(p&1)
    ans*=b%k;
    b=(b*b)%k;                  
    p>>=1;
    }
    ans%=k;
    cout<<b1<<"^"<<p1<<" mod "<<k<<"="<<ans;
    return 0;
 }