def egcd(a, b): if a == 0: return b, 0, 1 else: g, y, x = egcd(b % a, a) return g, x - (b // a) * y, y mod=10**9+7 x,k=map(int,input().split()) gcd,a,b=egcd(mod-1,-k) b*=-1 b%=mod-1 n=pow(x,b,mod) print(n) assert 1<=n<=10**9 assert 1<=k<=5*10**8 assert k%2==1 assert x==pow(n,k,mod)