結果

問題 No.187 中華風 (Hard)
ユーザー Kohei AsanoKohei Asano
提出日時 2020-01-01 18:09:05
言語 Rust
(1.72.1)
結果
WA  
実行時間 -
コード長 4,387 bytes
コンパイル時間 866 ms
コンパイル使用メモリ 149,008 KB
実行使用メモリ 4,504 KB
最終ジャッジ日時 2023-08-14 14:41:26
合計ジャッジ時間 4,214 ms
ジャッジサーバーID
(参考情報)
judge12 / judge11
このコードへのチャレンジ(β)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 111 ms
4,376 KB
testcase_03 AC 107 ms
4,376 KB
testcase_04 AC 124 ms
4,376 KB
testcase_05 AC 128 ms
4,380 KB
testcase_06 AC 129 ms
4,376 KB
testcase_07 AC 128 ms
4,380 KB
testcase_08 AC 83 ms
4,376 KB
testcase_09 AC 87 ms
4,380 KB
testcase_10 AC 84 ms
4,376 KB
testcase_11 WA -
testcase_12 AC 125 ms
4,380 KB
testcase_13 WA -
testcase_14 AC 27 ms
4,376 KB
testcase_15 WA -
testcase_16 WA -
testcase_17 AC 1 ms
4,376 KB
testcase_18 AC 1 ms
4,380 KB
testcase_19 AC 1 ms
4,376 KB
testcase_20 WA -
testcase_21 AC 1 ms
4,376 KB
testcase_22 AC 132 ms
4,380 KB
testcase_23 WA -
testcase_24 AC 1 ms
4,380 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
warning: unused imports: `max`, `min`
 --> Main.rs:2:16
  |
2 | use std::cmp::{max, min};
  |                ^^^  ^^^
  |
  = note: `#[warn(unused_imports)]` on by default

warning: unused imports: `HashMap`, `HashSet`
 --> Main.rs:3:24
  |
3 | use std::collections::{HashMap, HashSet};
  |                        ^^^^^^^  ^^^^^^^

warning: 2 warnings emitted

ソースコード

diff #

// =========
use std::cmp::{max, min};
use std::collections::{HashMap, HashSet};
use std::process::exit;

const MOD: i64 = 1000000007;

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, ) => {};
    // var... 変数の識別子, $t...型を一つよむ
    ($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, usize1) => {
        read_value!($iter, usize) - 1
    };
    // 配列の最後のNestではここで型が指定されてparseされる
    ($iter:expr, $t:ty) => {
        $iter.next().unwrap().parse::<$t>().expect("Parse error")
    };
}
// =========

// mod mの体におけるaの逆元
macro_rules! impl_modinv {
    ($U:ty) => {
        fn mod_inv(a: $U, m: $U) -> $U {
            let mut ab = (a as i64, m as i64);
            let mut uv = (1, 0);
            let mut t: i64;
            while ab.1 != 0 {
                t = ab.0 / ab.1;
                ab = (ab.1, ab.0 - t * ab.1);
                uv = (uv.1, uv.0 - t * uv.1);
            }
            // if ab.0 != 1 {
            //     // panic!("{} and {} are not coprime g={}", a, m, ab.0);
            //     println!("{:?}", -1);
            //     exit(0);
            // }
            let inv = uv.0 % m as i64;
            if inv < 0 {
                (inv + m as i64) as $U
            } else {
                inv as $U
            }
        }
    };
}
fn gcd(a: i64, b: i64) -> i64 {
    let (a, b) = if a < b { (b, a) } else { (a, b) };
    if b == 0 {
        return a;
    } else {
        return gcd(b, a % b);
    }
}

impl_modinv!(i64);
// mr[i].0 ... 互いに素
fn garner(mr: &mut Vec<(i64, i64)>, m: i64) -> i64 {
    mr.push((m, 0));
    // coef... mixed radixの係数, constants... 前まで求めた係数
    let mut coef: Vec<i64> = vec![1; mr.len()];
    let mut constants: Vec<i64> = vec![0; mr.len()];
    for i in 0..mr.len() - 1 {
        let v: i64 = (mr[i].1 - constants[i]) * mod_inv(coef[i], mr[i].0) % mr[i].0;
        for j in i + 1..mr.len() {
            constants[j] += coef[j] * v;
            constants[j] %= mr[j].0;
            coef[j] *= mr[i].0;
            coef[j] %= mr[j].0;
        }
    }
    constants[mr.len() - 1]
}

fn main() {
    input! {
        n: usize,
        mr: [(i64,i64);n]
    }
    let mut mr = mr.iter().map(|e| (e.1, e.0)).collect::<Vec<(i64, i64)>>();
    // 前処理 m を互いに素にする
    for i in 0..mr.len() {
        for j in i + 1..mr.len() {
            let mut g = gcd(mr[i].0, mr[j].0);
            // 解の条件チェック
            if (mr[i].1 - mr[j].1) % g != 0 {
                println!("{:?}", -1);
                exit(0);
            }
            //ひとまず互いに素にする
            mr[i].0 /= g;
            mr[j].0 /= g;

            //gi...mr[j].0に残らないやつ(とりきってない)
            let mut gi = gcd(mr[i].0, g);
            //mr[j]に残るのは必ず含まれるけど不純
            let mut gj = g / gi;
            // giが取り切れなかったのをとっていく
            while g != 1 {
                g = gcd(gi, gj);
                gi *= g;
                gj /= g;
            }
            mr[i].0 *= gi;
            mr[j].0 *= gj;
            // あまりの更新
            mr[i].1 %= mr[i].0;
            mr[j].1 %= mr[j].0;
        }
    }
    // let m: i64 = mr.iter().fold(1, |res, e| res * e.0);
    let mut ans = garner(&mut mr, MOD);
    if ans <= 0 {
        ans += MOD;
    }
    println!("{:?}", ans);
}
0