#include<cstdio>
const int mod=9901;
int inv[9901],i,j,x,y,a,b,n,ciur[8000],prime[8000];
void invers(int a,int b,int &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
	}
	else
	{
		int x0,y0;
		invers(b,a%b,x0,y0);
		x=y0;
		y=x0-(a/b)*y0;
	}
}
int pow(int a,int b)
{
	if(b==1)
		return a%mod;
	int x=pow(a,b/2);
	x=(x*x)%mod;
	if(b&1) x=(x*a)%mod;
	return x;
}
int main()
{
	freopen("sumdiv.in","r",stdin);
	freopen("sumdiv.out","w",stdout);
	inv[1]=1;
	for(i=2;i<mod;i++)
		if(!inv[i])
		{
			invers(i,mod,x,y);
			if(x<0) x=x+mod;
			x=x%mod;
			inv[i]=x;
			inv[x]=i;
		}
	for(i=2;i<8000;i++)
		if(ciur[i]==0)
			for(j=2*i;j<8000;j+=i)
				ciur[j]=1;
	for(i=2;i<8000;i++)
		if(ciur[i]==0) prime[++prime[0]]=i;
	scanf("%d %d",&a,&b);
	if(a!=1 && b) 
	{
		i=1;
		n=1;
		while(prime[i]*prime[i]<=a && n)
		{
			x=0;
			while(a%prime[i]==0)
				a/=prime[i],x++;
			if(x)
			{
				n=(n*(pow(prime[i],(x*b+1))-1)*inv[(prime[i]-1)%mod])%mod;
				if(n<0) n+=mod;
			}
			i++;
		}
		if(a!=1 && n)
			n=(n*(pow(a,b+1)-1)*inv[(a-1)%mod])%mod;
		if(n<0) n+=mod;
	}
	else
		n=1;
	printf("%d\n",n);
	fclose(stdout);
	return 0;
}