T = int(input())
inf = 10 ** 18
for _ in range(T):
    X, A = map(int, input().split())
    ans = inf
    for i in range(1, A + 1):
        q, r = divmod(A, i)
        val = q ** 2 * (i - r) + (q + 1) ** 2 * r
        val += X * i
        ans = min(ans, val)
    print(ans)