結果

問題 No.2527 H and W
コンテスト
ユーザー naut3
提出日時 2023-11-04 01:43:54
言語 Rust
(1.83.0 + proconio)
結果
AC  
実行時間 36 ms / 2,000 ms
コード長 3,193 bytes
コンパイル時間 13,420 ms
コンパイル使用メモリ 378,196 KB
実行使用メモリ 25,472 KB
最終ジャッジ日時 2024-09-25 21:49:53
合計ジャッジ時間 15,719 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 23
権限があれば一括ダウンロードができます
コンパイルメッセージ
warning: unused variable: `i`
  --> src/main.rs:47:17
   |
47 |     let (f, fi, i) = calc_combtable();
   |                 ^ help: if this is intentional, prefix it with an underscore: `_i`
   |
   = note: `#[warn(unused_variables)]` on by default

ソースコード

diff # 

#![allow(non_snake_case, unused_imports, unused_must_use)]
use std::io::{self, prelude::*};
use std::str;

const MOD: usize = 998_244_353;
const M: usize = 1_000_010;

fn calc_combtable() -> (Vec<usize>, Vec<usize>, Vec<usize>) {
    let mut f = vec![0; M];
    let mut fi = vec![0; M];
    let mut i = vec![0; M];

    f[0] = 1;
    f[1] = 1;
    fi[0] = 1;
    fi[1] = 1;
    i[1] = 1;

    for j in 2..M {
        f[j] = (f[j - 1] * j) % MOD;
        i[j] = MOD - i[MOD % j] * (MOD / j) % MOD;
        fi[j] = fi[j - 1] * i[j] % MOD;
    }

    return (f, fi, i);
}

fn main() {
    let (stdin, stdout) = (io::stdin(), io::stdout());
    let mut scan = Scanner::new(stdin.lock());
    let mut out = io::BufWriter::new(stdout.lock());

    macro_rules! input {
        ($T: ty) => {
            scan.token::<$T>()
        };
        ($T: ty, $N: expr) => {
            (0..$N).map(|_| scan.token::<$T>()).collect::<Vec<_>>()
        };
    }

    let H = input!(usize);
    let W = input!(usize);
    let K = input!(usize);

    let mut ans = 0;
    let (f, fi, i) = calc_combtable();

    macro_rules! comb {
        ($N: expr, $K: expr) => {
            if $N < $K {
                0
            } else {
                f[$N] * (fi[$K] * fi[$N - $K] % MOD) % MOD
            }
        };
    }

    // h \in 0..=H, w \in 0..=W, (H - h) * (W - w) == K
    let divs = calc_divisors(K);

    for d in divs {
        if d == K / d {
            if H >= d && W >= d {
                ans += 2 * comb!(H, H - d) * comb!(W, W - d);
            }
        } else {
            // H - h = d, W - w = K / d
            if H >= d && W >= K / d {
                ans += comb!(H, H - d) * comb!(W, W - K / d);
            }

            // // W - w = d, H - h = K / d
            if W >= d && H >= K / d {
                ans += comb!(W, W - d) * comb!(H, H - K / d);
            }
        }

        ans %= MOD;

        // println!("{} {}", d, ans);
    }

    ans *= 499122177;
    ans %= MOD;

    writeln!(out, "{}", ans);
}

fn calc_divisors(n: usize) -> Vec<usize> {
    let mut ret = vec![];

    for i in 1..=n {
        if i * i > n {
            break;
        }

        if n % i != 0 {
            continue;
        }

        ret.push(i);

        if n / i != i {
            ret.push(n / i);
        }
    }

    ret.sort();
    return ret;
}

struct Scanner<R> {
    reader: R,
    buf_str: Vec<u8>,
    buf_iter: str::SplitWhitespace<'static>,
}
impl<R: BufRead> Scanner<R> {
    fn new(reader: R) -> Self {
        Self {
            reader,
            buf_str: vec![],
            buf_iter: "".split_whitespace(),
        }
    }
    fn token<T: str::FromStr>(&mut self) -> T {
        loop {
            if let Some(token) = self.buf_iter.next() {
                return token.parse().ok().expect("Failed parse");
            }
            self.buf_str.clear();
            self.reader
                .read_until(b'\n', &mut self.buf_str)
                .expect("Failed read");
            self.buf_iter = unsafe {
                let slice = str::from_utf8_unchecked(&self.buf_str);
                std::mem::transmute(slice.split_whitespace())
            }
        }
    }
}
0