結果
問題 | No.1276 3枚のカード |
ユーザー | akakimidori |
提出日時 | 2020-10-30 23:09:23 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 685 ms / 2,000 ms |
コード長 | 3,763 bytes |
コンパイル時間 | 11,372 ms |
コンパイル使用メモリ | 403,176 KB |
実行使用メモリ | 8,396 KB |
最終ジャッジ日時 | 2024-07-22 02:50:41 |
合計ジャッジ時間 | 32,677 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 61 |
コンパイルメッセージ
warning: function `naive` is never used --> src/main.rs:55:4 | 55 | fn naive() { | ^^^^^ | = note: `#[warn(dead_code)]` on by default
ソースコード
// ---------- begin input macro ---------- // reference: https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8 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_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_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 ---------- // A の条件無視したものからAの条件に違反するものを引く // Bの倍数で違反するものは容易に消せるがBの約数の条件はどうしよう // fn naive() { input! { n: u128, } let mut ans = 0; for b in 1..=n { let rem = n / b - 1; ans += rem * (n - 2);// A全部 ans -= rem * (rem - 1);// AはBの倍数 } println!("{}", ans); // a はb の約数 for a in 1..=n { let q = n / a; if q - 1 < 2 { break; } for b in 2..=q { let rem = q / b - 1; if rem == 0 { break; } ans -= rem; } } ans %= 1_000_000_007; println!("{}", ans); } fn run() { input! { n: u128, } let mut ans = 0; let mut q = 1; let mut r = n; loop { let l = n / (q + 1); if l * l <= n { break; } let m = r - l; ans += m * (q - 1) * (n - 2);// A全部 ans -= m * (q - 1) * (q - 2);// AはBの倍数 r = l; q += 1; } for b in 1..=r { let rem = n / b - 1; ans += rem * (n - 2);// A全部 ans -= rem * (rem - 1);// AはBの倍数 } let mut map = std::collections::HashMap::new(); for i in 1..=r { map.insert(n / i, 0); } for i in 1..=q { map.insert(i, 0); } for (q, v) in map.iter_mut() { let n = *q; let mut val = 0u128; let mut q = 1; let mut r = n; loop { let l = n / (q + 1); if l * l <= n { break; } let m = r - l; val += (q - 1) * m; q += 1; r = l; } for i in 2..=r { val += n / i - 1; } *v = val; } let mut r = n; let mut q = 1; loop { let l = n / (q + 1); if l * l <= n { break; } let m = r - l; if q >= 3 { ans -= m * *map.get(&q).unwrap(); } r = l; q += 1; } for i in 1..=r { let q = n / i; if q >= 3 { ans -= *map.get(&q).unwrap(); } } ans %= 1_000_000_007; println!("{}", ans); } fn main() { // naive(); run(); }