結果

問題 No.1809 Divide NCK
ユーザー fukafukatanifukafukatani
提出日時 2022-01-15 17:32:33
言語 Rust
(1.77.0)
結果
AC  
実行時間 13 ms / 2,000 ms
コード長 2,113 bytes
コンパイル時間 2,486 ms
コンパイル使用メモリ 139,804 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-14 02:12:06
合計ジャッジ時間 3,940 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 13 ms
4,380 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 1 ms
4,380 KB
testcase_09 AC 1 ms
4,380 KB
testcase_10 AC 1 ms
4,376 KB
testcase_11 AC 10 ms
4,376 KB
testcase_12 AC 13 ms
4,376 KB
testcase_13 AC 3 ms
4,376 KB
testcase_14 AC 5 ms
4,380 KB
testcase_15 AC 1 ms
4,380 KB
testcase_16 AC 12 ms
4,380 KB
testcase_17 AC 1 ms
4,376 KB
testcase_18 AC 11 ms
4,380 KB
testcase_19 AC 12 ms
4,376 KB
testcase_20 AC 10 ms
4,376 KB
testcase_21 AC 9 ms
4,380 KB
testcase_22 AC 1 ms
4,376 KB
testcase_23 AC 1 ms
4,380 KB
testcase_24 AC 1 ms
4,376 KB
testcase_25 AC 1 ms
4,380 KB
testcase_26 AC 1 ms
4,376 KB
testcase_27 AC 1 ms
4,380 KB
testcase_28 AC 1 ms
4,380 KB
testcase_29 AC 10 ms
4,376 KB
testcase_30 AC 10 ms
4,380 KB
testcase_31 AC 3 ms
4,380 KB
testcase_32 AC 8 ms
4,376 KB
testcase_33 AC 6 ms
4,376 KB
testcase_34 AC 9 ms
4,380 KB
testcase_35 AC 9 ms
4,376 KB
testcase_36 AC 7 ms
4,376 KB
testcase_37 AC 11 ms
4,376 KB
testcase_38 AC 7 ms
4,380 KB
testcase_39 AC 12 ms
4,380 KB
testcase_40 AC 11 ms
4,376 KB
testcase_41 AC 7 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#![allow(unused_imports)]
use std::cmp::*;
use std::collections::*;
use std::io::Write;
use std::ops::Bound::*;

#[allow(unused_macros)]
macro_rules! debug {
    ($($e:expr),*) => {
        #[cfg(debug_assertions)]
        $({
            let (e, mut err) = (stringify!($e), std::io::stderr());
            writeln!(err, "{} = {:?}", e, $e).unwrap()
        })*
    };
}

fn main() {
    let v = read_vec::<i128>();
    let (n, k, m) = (v[0], v[1], v[2]);
    let mut sq = 1;
    for i in 1.. {
        if i * i > m {
            sq = i - 1;
            break;
        }
    }
    let primes = get_primes(sq + 1);

    let mut mm = m;
    let mut dividors = vec![];
    let mut q = vec![];
    for p in primes {
        if mm % p == 0 {
            dividors.push(p);
            q.push(0);
            while mm % p == 0 {
                q[dividors.len() - 1] += 1;
                mm /= p;
            }
        }
    }
    if mm > 1 {
        dividors.push(mm);
        q.push(1);
    }
    debug!(dividors);
    let ans = dividors
        .iter()
        .zip(q.iter())
        .map(|(&x, &q)| (count(x, n) - count(x, n - k) - count(x, k)) / q)
        .min()
        .unwrap();
    println!("{}", ans);
}

fn count(p: i128, n: i128) -> i128 {
    let mut cur = p;
    let mut ret = 0;
    while cur <= n {
        ret += n / cur;
        cur *= p;
    }
    ret
}

fn get_primes(n: i128) -> Vec<i128> {
    let mut is_prime = vec![true; n as usize + 1];
    let mut primes = Vec::new();
    is_prime[0] = false;
    is_prime[1] = false;

    for i in 2..n + 1 {
        if is_prime[i as usize] {
            primes.push(i);
            let mut j = 2 * i;
            while j <= n {
                is_prime[j as usize] = false;
                j += i;
            }
        }
    }
    primes
}

fn read<T: std::str::FromStr>() -> T {
    let mut s = String::new();
    std::io::stdin().read_line(&mut s).ok();
    s.trim().parse().ok().unwrap()
}

fn read_vec<T: std::str::FromStr>() -> Vec<T> {
    read::<String>()
        .split_whitespace()
        .map(|e| e.parse().ok().unwrap())
        .collect()
}
0