for _ in range(int(input())):
    n, m = map(int, input().split())
    p = 1000000007
    a = pow(n, m, p)
    b = pow(n, m//2, p)
    print(((1+a)*a%p*500000004%p-b*(b+1)%p*(2*b+1)%p*166666668%p)%p)