1633：【例 3】Sumdiv

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const ll M=;
ll a,b,ans=,cnt;
ll p[],c[];
inline void divide(int n)
{
    for(int i=;i*i<=n;i++)
    {
        if(n%i==)
        {
            p[++cnt]=i;
            while(n%i==) n/=i,c[cnt]++;
        }
    }
    if(n>) p[++cnt]=n,c[cnt]=;
    //for(int i=1;i<=cnt;i++) cout<<p[cnt]<<' '<<c[cnt]<<' ';
}
inline int power(ll x,ll y)
{
    ll res=;
    while(y)
    {
        if(y&) res=res*x%M;
        x=x*x%M;
        y>>=;
    }
    return res;
}
int main()
{
    scanf("%lld%lld",&a,&b);
    divide(a);
    for(int i=;i<=cnt;i++)
    {
        if((p[i]-)%M==)
        {
            ans=(b*c[i]+)%M*ans%M;
            continue;
        }
        ll xx=power(p[i],b*c[i]+);
        xx=(xx-+M)%M;
        ll yy=p[i]-;
        yy=power(yy,M-);
        ans=ans*xx%M*yy%M;
    }
    printf("%lld\n",ans);
    return ;
}