#int(input())
#map(int, input().split())
#list(map(int, input().split()))

T = int(input())
x = [0] * T
for i in range(T):
    x[i] = list(map(int, input().split()))

def f(c):
    return c ** 2 + X

for k in range(T):
    X, A = x[k]
    ans = A ** 2 + X
    for i in range(1, A+1):
        a, b = divmod(A, i)
        t = (i - b) * f(a) + b * f(a+1)
        ans = min(ans, t)
    print(ans)