l = int(input()) ans = 1 mod = 10 ** 9 + 7 for _ in range(l): p, e = map(int,input().split()) ans *= ((pow(p, e + 2, mod) - 1) * pow((p - 1) * (p - 1), mod - 2, mod) - (e + 2) * pow((p - 1), mod - 2, mod)) % mod ans %= mod print(ans)