fn solve() -> usize {
    let mut xa = String::new();
    std::io::stdin().read_line(&mut xa).ok();
    let xa: Vec<usize> = xa.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();
    let x = xa[0];
    let a = xa[1];

    let mut result = 1usize << 60;
    for cnt in 1..=a {
        let base = a / cnt;
        let rest = a % cnt;
        result = result.min(cnt * x + (base+1) * (base+1) * rest + base * base * (cnt-rest));
    }
    result
}

fn main() {
    let mut t = String::new();
    std::io::stdin().read_line(&mut t).ok();
    let t: usize = t.trim().parse().unwrap();
    let mut result = Vec::with_capacity(t);
    for _ in 0..t {
        result.push(solve());
    }
    for &v in result.iter() {
        println!("{}", v);
    }
}