const MOD: isize = 1e9 as isize + 7; fn power_mat(target: &Vec>, times: usize) -> Vec> { let n = target.len(); if times == 0 { let mut ret = vec![vec![0isize; n]; n]; for i in 0..n { ret[i][i] = 1; } return ret; } if times == 1 { return target.clone(); } let temp = power_mat(target, times/2); if times % 2 == 1 { calc_mat(&calc_mat(&temp, &temp), target) } else { calc_mat(&temp, &temp) } } fn calc_mat(left: &Vec>, right: &Vec>) -> Vec> { let sizei = left.len(); let sizej = right[0].len(); let sizek = left[0].len(); let mut ret = vec![vec![0isize; sizej]; sizei]; for i in 0..sizei { for j in 0..sizej { ret[i][j] = (0..sizek).map(|k| left[i][k] * right[k][j] % MOD).sum::() % MOD; } } ret } fn main() { let mut n = String::new(); std::io::stdin().read_line(&mut n).ok(); let n: usize = n.trim().parse().unwrap(); let mut abc = String::new(); std::io::stdin().read_line(&mut abc).ok(); let abc: Vec = abc.trim().split_whitespace().map(|s| s.parse().unwrap()).collect(); let a = abc[0]; let b = abc[1]; let c = abc[2]; let times = n - 1; let mat = vec![vec![1, -1, 0], vec![0, 1, -1], vec![-1, 0, 1]]; let mat = power_mat(&mat, times); let vals = vec![vec![a], vec![b], vec![c]]; let result = calc_mat(&mat, &vals); let ra = (result[0][0] + MOD) % MOD; let rb = (result[1][0] + MOD) % MOD; let rc = (result[2][0] + MOD) % MOD; println!("{} {} {}", ra, rb, rc); }