def product_mod(*X): y=1 for x in X: y*=x%Mod; y%=Mod return y #================================================== import sys input=sys.stdin.readline Ans=1; Mod=10**9+7 L=int(input()) for _ in range(L): P,e=map(int,input().split()) alpha=pow(P-1,Mod-2,Mod) Ans*=product_mod(pow(P,e+2,Mod)-(P-1)*(e+2)-1, alpha**2) Ans%=Mod print(Ans%Mod)