import sys
def MI(): return map(int,sys.stdin.readline().rstrip().split())


N,K,M = MI()
mod = 10**9+7

fac = [1]
for i in range(1,N):
    fac.append((fac[-1]*i) % mod)
fac_inv = [pow(fac[i],mod-2,mod) for i in range(N)]

a = 0
for i in range(1,N+1):
    if K % i == 0:
        a += fac[N-1]*fac_inv[N-i]*pow(N,N-i,mod)
        a %= mod
        
if M == 1:
    print(a)
else:
    ans = (pow(N,N,mod)-a)*pow(N-1,mod-2,mod)
    ans %= mod
    print(ans)