#!/usr/bin/env PyPy3

from collections import Counter, defaultdict, deque
import itertools
import re
import math
from functools import reduce
import operator
import bisect
from heapq import *
import functools
mod=10**9+7

import sys
input=sys.stdin.readline

t = int(input())
for _ in range(t):
    n,m=map(int,input().split())
    nm = pow(n,m,mod)
    ans = (1 + nm) * nm // 2
    ans %= mod
    nm2 = pow(n,m//2,mod)
    ans -= nm2 * (nm2 + 1) * (1 + 2 * nm2) // 6
    ans %= mod
    print(ans)