t = int(input())
mod = 10 ** 9 + 7
for i in range(t):
    n,m = map(int,input().split())
    ans = pow(n,m,mod) * (pow(n,m,mod) + 1) * pow(2,mod-2,mod)
    val = pow(n,m//2,mod)
    ans -= val * (val+1)*(2*val+1) * pow(6,mod-2,mod)
    
    print(ans%mod)