結果

問題 No.1596 Distance Sum in 2D Plane
コンテスト
ユーザー naut3
提出日時 2024-01-30 15:11:46
言語 Rust
(1.83.0 + proconio)
結果
CE  
(最新)
AC  
(最初)
実行時間 -
コード長 2,720 bytes
コンパイル時間 30,298 ms
コンパイル使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-28 10:15:30
合計ジャッジ時間 32,201 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
コンパイルエラー時のメッセージ・ソースコードは、提出者また管理者しか表示できないようにしております。（リジャッジ後のコンパイルエラーは公開されます）
ただし、clay言語の場合は開発者のデバッグのため、公開されます。

コンパイルメッセージ
コンパイルが30秒の制限時間を超えました

ソースコード

diff # 

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

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 N = input!(usize);
    let M = input!(usize);

    let mut ans = ((N + N) * C(2 * N, N)) % MOD;

    for _ in 0..M {
        let t = input!(usize);
        let x = input!(usize);
        let y = input!(usize);

        if t == 1 {
            let mut d = C(x + y, x) * C(N - (x + 1) + N - y, N - y);
            d %= MOD;
            ans = ans + MOD - d;
            ans %= MOD;
        } else {
            let mut d = C(x + y, x) * C(N - x + N - (y + 1), N - x);
            d %= MOD;
            ans = ans + MOD - d;
            ans %= MOD;
        }
    }

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

const M_: usize = 4 * 100_001;
const MOD: usize = 1_000_000_007;

const fn comb() -> ([usize; M_], [usize; M_], [usize; M_]) {
    let mut f = [0; M_];
    let mut f_inv = [0; M_];
    let mut inv = [0; M_];

    f[0] = 1;
    f[1] = 1;
    f_inv[0] = 1;
    f_inv[1] = 1;
    inv[1] = 1;

    let mut i = 2;

    while i < M_ {
        f[i] = (f[i - 1] * i) % MOD;
        inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
        f_inv[i] = f_inv[i - 1] * inv[i] % MOD;

        i += 1;
    }

    return (f, f_inv, inv);
}

const _F: [usize; M_] = comb().0;
const _F_INV: [usize; M_] = comb().1;
const _INV: [usize; M_] = comb().2;

const fn C(n: usize, k: usize) -> usize {
    if n < k {
        return 0;
    } else {
        return _F[n] * (_F_INV[k] * _F_INV[n - k] % MOD) % MOD;
    }
}

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