T=int(input())
for _ in range(T):
  x,a=map(int,input().split())
  ans=(1+x)*a
  for i in range(1,a):
    b=a//i
    r=a%i
    ans=min(ans,((b)**2)*(i-r)+((b+1)**2)*r+x*i)
  print(ans)