N=int(input()) ans=1 mod=10**9+7 def mult(A,B): C=[[0 for i in range(2)] for j in range(2)] for i in range(2): for j in range(2): for k in range(2): C[i][j]+=A[i][k]*B[k][j] C[i][j]%=mod return C AS=[[[0,1],[1,1]]] # A^(2^i) for i in range(60): AS.append(mult(AS[i],AS[i])) for _ in range(N): C,D=map(int,input().split()) D%=mod-1 X=[1,1] A=[[0,1],[1,1]] for i in range(60): if C%2==1: X[0],X[1]=X[0]*AS[i][0][0]+X[1]*AS[i][0][1],X[0]*AS[i][1][0]+X[1]*AS[i][1][1] X[0]%=mod X[1]%=mod C>>=1 ans*=pow(X[1],D,mod) ans%=mod print(ans)