import sys
input = sys.stdin.readline

mod = 10**9+7
T = int(input())
for _ in range(T):
    a, N = map(int, input().split())
    M = int((a*N)**0.5)
    ans = 0
    if (M+1)**2 <= a*N:
        M += 1
    for r in range(a):
        m = (M-r)//a
        ans += a*m*(m+1)*(2*m+1)//6+m*(m+1)*r + (r**2-1)//a*(m+1)
    ans = (N*M - ans - 1) % mod
    print(ans)