from collections import *
from itertools import *
from functools import *
from heapq import *
import sys,math
input = sys.stdin.readline


def answer():
    X,A = map(int,input().split())
    L = math.floor(math.sqrt(X))
    INF =  (1<<30)
    dp = [INF]*(A+1)
    
    dp[0] = 0
    for i in range(1,A+1):
        if i<L-10:
            dp[i] = i**2 + X
            continue
        dp[i] = min(k**2 + X + dp[i-k] for k in range(max(1,L-10),min(i,L+10)+1))
    print(dp[A])
for _ in range(int(input())):
    answer()