import sys, random
input = sys.stdin.readline
mod = 10 ** 9 + 7
rev2 = pow(2, mod - 2, mod)
rev6 = pow(6, mod - 2, mod)

for _ in [0] * int(input()):
	n, m = map(int, input().split());	n %= mod
	
	a = pow(n, m, mod)
	b = pow(n, m // 2, mod)

	s = rev2 * a % mod * (a + 1) % mod
	t = rev6 * b % mod * (b + 1) % mod * (b * 2 + 1) % mod

	ans = (s + mod - t) % mod
	print(ans)