結果
問題 | No.1358 [Zelkova 2nd Tune *] 語るなら枚数を... |
ユーザー |
![]() |
提出日時 | 2021-01-22 23:18:55 |
言語 | Rust (1.83.0 + proconio) |
結果 |
AC
|
実行時間 | 389 ms / 2,000 ms |
コード長 | 2,441 bytes |
コンパイル時間 | 15,662 ms |
コンパイル使用メモリ | 402,368 KB |
実行使用メモリ | 15,600 KB |
最終ジャッジ日時 | 2024-12-28 07:24:14 |
合計ジャッジ時間 | 17,594 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 |
ソースコード
// ---------- begin input macro ----------// reference: https://qiita.com/tanakh/items/0ba42c7ca36cd29d0ac8macro_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 ----------use std::io::Write;fn gcd(mut a: usize, mut b: usize) -> usize {while b > 0 {let r = a % b;a = b;b = r;}a}fn run() {let out = std::io::stdout();let mut out = std::io::BufWriter::new(out.lock());input! {t: usize,p: [([usize; 3], usize); t],}for (mut p, y) in p {p.sort();let (a, b, c) = (p[0], p[1], p[2]);let inf = y + 1;let mut dp = vec![inf; b];for i in 0.. {let pos = i * a % b;if i > 0 && pos == 0 {break;}dp[pos] = i * a;}let lcm = a * b / gcd(a, b);let mut ans = 0;const MOD: usize = 1_000_000_007;for i in 0..=(y / c) {let y = y - i * c;let r = y % b;if dp[r] <= y {ans += (y - dp[r]) / lcm + 1;}}ans %= MOD;writeln!(out, "{}", ans).ok();}}fn main() {run();}