T = int(input())
for _ in range(T):
    X, A = map(int, input().split())
    ans = A * A + X
    for i in range(1, max(1, A // 2) + 1):
        b = A % i
        a = (A - b * (i + 1)) // i
        if a >= 0 and b >= 0 and a * i + b * (i + 1) == A:
            ans = min(ans, a * (i * i + X) + b * ((i + 1) * (i + 1) + X))
            # print(i,a,b,a*i+b*(i+1),ans)
        a = A % (i + 1)
        b = (A - a * i) // (i + 1)
        if a >= 0 and b >= 0 and a * i + b * (i + 1) == A:
            ans = min(ans, a * (i * i + X) + b * ((i + 1) * (i + 1) + X))
            # print(i,a,b,a*i+b*(i+1),ans)
    print(ans)