t = int(input())
INF = 10**18
for _ in range(t):
    x, a = map(int, input().split())
    ans = INF
    for cnt in range(1, a + 1):
        q, r = divmod(a, cnt)
        res = cnt * x + (cnt - r) * q**2 + r* (q + 1)**2
        ans = min(ans, res)
    print(ans)