import sys input = sys.stdin.buffer.readline sys.setrecursionlimit(10 ** 7) mod = 10 ** 9 + 7 def F(p, e): n = e + 1 res = (pow(p, e + 1, mod) - 1) * pow(p - 1, mod - 2, mod) % mod res = p * res - e - 1 res = res * pow(p - 1, mod - 2, mod) % mod return res ans = 1 L = int(input()) for _ in range(L): p, e = map(int, input().split()) ans = ans * F(p, e) % mod print(ans)