z=10**9+7
N,M=map(int,input().split())
g=[1,1]
h=[1,1]
v=[0,1]
for i in range(2,N+1):
  g.append(g[-1]*i%z)
  v.append(-v[z%i]*(z//i)%z)
  h.append(h[-1]*v[-1]%z)
L=N-M
print(pow(M,L,z)*g[N]*h[L]*pow(L+1,z-2,z)%z)