import sys
input = sys.stdin.readline

mod = 10 ** 9 + 7
inv2 = pow(2, mod - 2, mod)
inv6 = pow(6, mod - 2, mod)

for _ in [0] * int(input()):
    n, m = map(int, input().split())
    h = m // 2
    a = pow(n, m, mod)
    b = pow(n, h, mod)
    ans = a * (a + 1) % mod * inv2 % mod + mod - b * (b + 1) % mod * (2 * b + 1) % mod * inv6 % mod
    ans %= mod
    print(ans)