MOD = 10**9+7
T = int(input())
for _ in range(T):
    N,M = map(int,input().split())
    a = pow(N,M,MOD)
    b = pow(N,M//2,MOD)
    ans = ((a+1)*a*pow(2,MOD-2,MOD)-b*(b+1)*(2*b+1)*pow(6,MOD-2,MOD))%MOD
    print(ans)