#include <bits/stdc++.h>
using namespace std;

void solve() {
    int64_t X, A;
    cin >> X >> A;

    int64_t ans = LLONG_MAX;
    for (int64_t n = 1; n <= A; n++) {
        const int64_t a = A / n;
        const int64_t m = A % n;

        const int64_t curr_ans = n * X + (a + 1) * (a + 1) * m + a * a * (n - m);
        ans = min(ans, curr_ans);
    }
    cout << ans << endl;
}

int main() {
    int64_t T;
    cin >> T;
    while (T--) {
        solve();
    }
}