#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using P = pair<ll,ll>;
#define fix(x) fixed << setprecision(x)
#define asc(x) x, vector<x>, greater<x>
#define rep(i, n) for(ll i = 0; i < n; i++)
#define all(x) (x).begin(),(x).end()
template<class T>bool chmin(T&a, const T&b){if(a>b){a=b;return 1;}return 0;}
template<class T>bool chmax(T&a, const T&b){if(a<b){a=b;return 1;}return 0;}

int main(){
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    int t;
    cin >> t;
    rep(houren,t){
        ll x,a,ans = (1LL<<60);
        cin >> x >> a;
        for(ll i=1;i<=a;i++){
            ll upd = x*i + (a/i)*(a/i)*(i-a%i) + (a/i+1)*(a/i+1)*(a%i);
            chmin(ans, upd);
        }
        cout << ans << '\n';
    }
    return 0;
}