結果
問題 |
No.3307 Almost Equal
|
ユーザー |
![]() |
提出日時 | 2025-10-05 13:55:12 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 1 ms / 2,000 ms |
コード長 | 2,868 bytes |
コンパイル時間 | 10,764 ms |
コンパイル使用メモリ | 400,476 KB |
実行使用メモリ | 7,716 KB |
最終ジャッジ日時 | 2025-10-05 13:55:26 |
合計ジャッジ時間 | 12,664 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 46 |
ソースコード
// round(x) = floor(x + 0.5) // round(i*A/B) = floor(Ai/B + 0.5) = floor((2Ai + B) / 2B) fn main() { input! { a: u64, b: u64, c: u64, d: u64, } let mut a = a; let mut b = b; let mut c = c; let mut d = d; if a * d < b * c { std::mem::swap(&mut a, &mut c); std::mem::swap(&mut b, &mut d); } let g = gcd(a, b); a /= g; b /= g; let g = gcd(c, d); c /= g; d /= g; if (a, b) == (c, d) { println!("-1"); return; } // (a/b - c/d)x >= 1 // (ad - bc) / bd * x >= 1 let x = b * d / (a * d - b * c); let ans = x - (floor_sum(x, 2 * b, 2 * a, 2 * a + b) - floor_sum(x, 2 * d, 2 * c, 2 * c + d)); println!("{}", ans); } fn gcd(a: u64, b: u64) -> u64 { if b == 0 { a } else { gcd(b, a % b) } } // ---------- begin input macro ---------- // reference: https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8 #[macro_export] macro_rules! input { (source = $s:expr, $($r:tt)*) => { let mut iter = $s.split_whitespace(); input_inner!{iter, $($r)*} }; ($($r:tt)*) => { let s = { use std::io::Read; let mut s = String::new(); std::io::stdin().read_to_string(&mut s).unwrap(); s }; let mut iter = s.split_whitespace(); input_inner!{iter, $($r)*} }; } #[macro_export] macro_rules! input_inner { ($iter:expr) => {}; ($iter:expr, ) => {}; ($iter:expr, $var:ident : $t:tt $($r:tt)*) => { let $var = read_value!($iter, $t); input_inner!{$iter $($r)*} }; } #[macro_export] macro_rules! read_value { ($iter:expr, ( $($t:tt),* )) => { ( $(read_value!($iter, $t)),* ) }; ($iter:expr, [ $t:tt ; $len:expr ]) => { (0..$len).map(|_| read_value!($iter, $t)).collect::<Vec<_>>() }; ($iter:expr, chars) => { read_value!($iter, String).chars().collect::<Vec<char>>() }; ($iter:expr, bytes) => { read_value!($iter, String).bytes().collect::<Vec<u8>>() }; ($iter:expr, usize1) => { read_value!($iter, usize) - 1 }; ($iter:expr, $t:ty) => { $iter.next().unwrap().parse::<$t>().expect("Parse error") }; } // ---------- end input macro ---------- // ---------- begin floor sum ---------- // sum_{i = 0}^{n - 1} floor((ai + b) / m) pub fn floor_sum(n: u64, m: u64, mut a: u64, mut b: u64) -> u64 { let mut ans = 0; ans += b / m * n; if n % 2 == 0 { ans += a / m * (n / 2) * (n - 1); } else { ans += a / m * n * ((n - 1) / 2); } a %= m; b %= m; let p = a as u128 * n as u128 + b as u128; if p >= (m as u128){ ans += floor_sum((p / m as u128) as u64, a, m, (p % m as u128) as u64); } ans } // ---------- end floor sum ----------