for _ in range(int(input())):
    X, A = map(int, input().split())
    ans = 10**18
    
    for c in range(1, A+1):
        s, a = divmod(A, c)
        cost = (s+1)*(s+1)*a+s*s*(c-a)+c*X
        ans = min(ans, cost)
    
    print(ans)