import sys input = sys.stdin.readline mod=10**9+7 ANS=1 L=int(input()) for t in range(L): p,e=map(int,input().split()) X=-(e+1) + (p * (1-pow(p,e,mod)))*pow(1-p,mod-2,mod) + pow(p,e+1,mod) #print(X%mod*pow(p-1,mod-2,mod)) ANS=ANS*X%mod*pow(p-1,mod-2,mod)%mod print(ANS)