#pragma GCC optimize("Ofast")
#pragma GCC optimize(1)
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC target("popcnt")
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef string str;
typedef pair<ll,ll> pii;
#define F first
#define S second
#define pb push_back
#define pq priority_queue
#define all(x) (x).begin(),(x).end()
#define bug(x) cerr << (x) << '\n'
#define yn(x) cout << (x==1?"Yes\n":"No\n")
#define PI 3.14159265358979323
#define Ststone ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0)
const ll max_n=2e5+20,max_int=2147483647,mod1=1e9+7,mod2=998244353;
const bool online = 0;
ll dx[8] = {1,0,0,-1,1,1,-1,-1};
ll dy[8] = {0,1,-1,0,1,-1,1,-1};
ll t,n,q,m,a,b,c,d,e;
vector<ll> v,v1,v2,v3,v4,vv;
char cc;
str s;

ll f(ll x,ll y){
    ll cnt = x%y;
    ll sum = 0;
    for(ll i=0;i<y;i++){
        if(cnt){
            sum += ((x/y)+1)*((x/y)+1);
            cnt--;
        }
        else
            sum += (x/y)*(x/y);
    }
    return sum;
}

void solve(){
    ll x;
    cin >> x >> a;
    ll ans = x+a*a;
    ll now = 2;
    while(f(a,now)+now*x<ans){
        ans = f(a,now)+now*x;
        now++;
    }
    cout << ans << '\n';

}

int main(){
    Ststone;
    if(online){
        #ifndef ONLINE_JUDGE
        freopen("input.txt","r",stdin);
        freopen("output.txt","w",stdout);
        #else
        #endif
    }
    cin >> t;
    //t = 1;
    for(ll i=1;i<=t;i++){
        //cout << "Case #"<<i<<": ";
        solve();
    }
}


/*
 a^2
 2*(a/2)^2 = a^2/2+x

*/