mod = 10**9 + 7
inv2 = pow(2, mod-2, mod)
inv6 = pow(6, mod-2, mod)
T = int(input())
for _ in range(T):
	n, m = map(int,input().split())
	# 1 -> n^m のうち, 平方数を除いたものの, 総和!!
	v = pow(n, m, mod)
	w = pow(n, m//2, mod)
	print((v*(v+1)*inv2 - w*(w+1)*(2*w+1)*inv6) % mod)