結果
問題 | No.1253 雀見椪 |
ユーザー |
![]() |
提出日時 | 2020-10-09 23:30:09 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 68 ms / 2,000 ms |
コード長 | 1,952 bytes |
コンパイル時間 | 12,363 ms |
コンパイル使用メモリ | 376,716 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-20 14:22:27 |
合計ジャッジ時間 | 14,064 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 14 |
ソースコード
macro_rules! read_to_tuple { ( $( $t:ty ),* ) => {{ let mut input = String::new(); std::io::stdin().read_line(&mut input).unwrap(); let mut iter = input.split_whitespace(); ( $( iter.next().unwrap().parse::<$t>().unwrap() ),* ) }} } const MOD: i64 = 1_000_000_007; fn gcd(a: i64, b: i64) -> i64 { let r = a % b; if r != 0 { gcd(b, r) } else { b } } fn lcm(a: i64, b: i64) -> i64 { a / gcd(a, b) * b } fn pow(mut a: i64, mut n: i64) -> i64 { let mut res = 1; while n > 0 { if n & 1 > 0 { res = res * a % MOD; } a = a * a % MOD; n >>= 1; } res } fn solve() -> String { let (n, mut ag, mut bg, mut ac, mut bc, mut ap, mut bp) = { read_to_tuple!(i64, i64, i64, i64, i64, i64, i64) }; let f = |a: &mut i64, b: &mut i64| { let x = gcd(*a, *b); *a /= x; *b /= x; }; f(&mut ag, &mut bg); f(&mut ac, &mut bc); f(&mut ap, &mut bp); let g = |al: i64, bl: i64, ar: i64, br: i64| -> i64 { let mut res = 1; let x = lcm(bl, br); res = res * pow((al * (x / bl) + ar * (x / br)) % MOD, n) % MOD; res = res * pow(pow(x % MOD, n), MOD - 2) % MOD; res }; let h = |a: i64, b: i64| -> i64 { let mut res = 1; res = res * pow(a, n) % MOD; res = res * pow(pow(b, n), MOD - 2); res }; let mut ans = 0; ans = (ans + g(ag, bg, ac, bc)) % MOD; ans = (ans + g(ac, bc, ap, bp)) % MOD; ans = (ans + g(ap, bp, ag, bg)) % MOD; ans = (ans + MOD - h(ag, bg) * 2 % MOD) % MOD; ans = (ans + MOD - h(ac, bc) * 2 % MOD) % MOD; ans = (ans + MOD - h(ap, bp) * 2 % MOD) % MOD; ans = (1 + MOD - ans) % MOD; ans.to_string() } fn main() { let t = read_to_tuple!(usize); let answers = (0..t).map(|_| solve()).collect::<Vec<_>>(); println!("{}", answers.join("\n")); }