結果
問題 |
No.2125 Inverse Sum
|
ユーザー |
![]() |
提出日時 | 2022-12-31 22:00:34 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 62 ms / 2,000 ms |
コード長 | 932 bytes |
コンパイル時間 | 12,531 ms |
コンパイル使用メモリ | 376,916 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-26 21:38:49 |
合計ジャッジ時間 | 14,892 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
fn main() { let mut s = String::new(); std::io::stdin().read_line(&mut s).ok(); let mut itr = s.trim().split_whitespace(); let p: u64 = itr.next().unwrap().parse().unwrap(); let q: u64 = itr.next().unwrap().parse().unwrap(); let mut v = Vec::new(); let mut n = q; let mut i = 2u64; loop { if i * i > n { break; } if n % i == 0 { v.push(i); while n % i == 0 { n /= i; } } i += 1; } if n > 1 { v.push(n); } let mut n = q * q; let mut d = Vec::new(); d.push(1u64); for x in v.iter() { let l = d.len(); while n % x == 0 { for i in (d.len() - l)..d.len() { let y = d[i] * x; d.push(y); } n /= x; } } let mut ans = Vec::new(); for &d in d.iter() { if (d + q) % p == 0 && (q * q / d + q) % p == 0{ ans.push(((d + q) / p, (q * q / d + q) / p)); } } ans.sort(); ans.dedup(); println!("{}", ans.len()); for a in ans { println!("{} {}", a.0, a.1); } }