import sys
input=lambda:sys.stdin.readline().rstrip()
N=int(input())
for i in range(N):
	T,A=map(int,input().split())
	ans=float('inf')
	for j in range(1,A+1):
		ans=min(ans,j*T+(A%j)*(A//j+1)**2+(j-A%j)*(A//j)**2)
	print(ans)