N=int(input())
mod=10**9+7

def sum1(x):
	return (x*(x+1)//2) % mod
def sum2(x):
	return (x*(x+1)*(2*x+1)//6)%mod
	
for _ in range(N):
	a,b=map(int,input().split())
	print((sum1(pow(a,b,mod))-sum2(pow(a,b//2,mod)))%mod)