T = int(input())
for t in range(T):
    X , A = map(int,input().split())
    dp = [0] + [10**20] * A
    XB = int(X**0.5)+2
    XA = min(max(1,A-2),int(X**0.5)-2)
    for j in range(XA,XB):
        flag = 0
        for i in range(A):
            if i + j > A:
                continue
            tmp = dp[i] + j * j + X
            if dp[i + j] > tmp:
                dp[i + j] = tmp
                flag = 1
        #print(dp)
        if flag == 0:
            break
print(dp[-1])