import sys
input = sys.stdin.readline
sys.setrecursionlimit(10**7)
mod = 10**9+7

X = [[1, 1], [1, 1]]
Y = [[1, 1], [1, 1]]
def f():
    n, m, x = map(int, input().split())
    X[0][0] = 1
    X[0][1] = m
    X[1][0] = m
    X[1][1] = 1
    Y[0][0] = 1
    Y[0][1] = 0
    Y[1][0] = 0
    Y[1][1] = 1
    
    while n:
        if n%2:
            a = X[0][0]*Y[0][0] + X[0][1]*Y[1][0]
            b = X[0][0]*Y[0][1] + X[0][1]*Y[1][1]
            c = X[1][0]*Y[0][0] + X[1][1]*Y[1][0]
            d = X[1][0]*Y[0][1] + X[1][1]*Y[1][1]
            Y[0][0] = a%mod
            Y[0][1] = b%mod
            Y[1][0] = c%mod
            Y[1][1] = d%mod
        a = X[0][0]*X[0][0] + X[0][1]*X[1][0]
        b = X[0][0]*X[0][1] + X[0][1]*X[1][1]
        c = X[1][0]*X[0][0] + X[1][1]*X[1][0]
        d = X[1][0]*X[0][1] + X[1][1]*X[1][1]
        X[0][0] = a%mod
        X[0][1] = b%mod
        X[1][0] = c%mod
        X[1][1] = d%mod
        n //= 2
    
    print(Y[0][x])

s = int(input())
for _ in range(s):
    f()