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