T = int(input())
for _ in range(T):
    X,A = map(int,input().split())
    ans = A**2+X
    for n in range(2,A+1):
        q = A//n
        r = A%n
        ans = min(ans,(q+1)**2*r+q**2*(n-r)+n*X)
    print(ans)