厚脸皮求助

费马小定理求逆元出了错

#include<bits/stdc++.h>
using namespace std;
const long long mod=9901;
typedef long long ll;
ll a,b,cnt,ans;
ll ksm(ll x,ll y)
{
    ll ans=1;
    while(y>0)
    {
        if(y&1)
        {
            y--;
            ans=(ans*x)%mod;
        }
        y/=2;
        x=(x*x)%mod;
    }
    return ans%mod;
}
int main()
{
    long long i,j;
    scanf("%lld%lld",&a,&b);
    ans=1;
    for(i=2;i<=a;i++)
    {
        ll c=0;
        while(a%i==0)
        {
            a/=i;
            c++;
        }
        if(c)
            ans=((ans*(ksm(i,c*b+1)-1+mod))%mod*(ksm(i,mod-2)%mod)%mod)%mod;
    }
    if(a!=1) ans=((ans*(ksm(a,b+1)-1+mod))%mod*(ksm(a,mod-2)%mod)%mod)%mod;
    printf("%lld",ans);
    return 0;
}