mod = 1000000007


def main():
    import sys
    from math import floor, sqrt
    input = sys.stdin.readline

    def floor_sqrt(a):
        x0 = floor(sqrt(a))
        res = 0
        for x in range(x0 - 1, x0 + 2):
            if x <= 0:
                continue
            if x ** 2 <= a:
                res = max(res, x)
        return res

    for _ in range(int(input())):
        a, N = map(int, input().split())
        ans = prev = floor_sqrt(a)
        n = 1
        while n < N:
            n_new = ((prev + 1) ** 2 - 1) // a + 1
            val = floor_sqrt(a * n_new)
            assert val ** 2 <= a * n_new < (val + 1) ** 2
            if n_new <= N:
                ans += val + (n_new - n - 1) * prev
            else:
                ans += prev * (N - n)
            ans %= mod
            n = n_new
            prev = val
        print(ans)


if __name__ == '__main__':
    main()