const INF: usize = 1usize << 60;

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 = INF;
    for times in 1..=a {
        let rest = a % times;
        let base = a / times;
        result = result.min(x * times + (times-rest) * base.pow(2) + rest * (base+1).pow(2));
    }
    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);
    }
}