結果

問題 No.1809 Divide NCK
ユーザー customaddonecustomaddone
提出日時 2022-01-14 22:00:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 64 ms / 2,000 ms
コード長 1,981 bytes
コンパイル時間 374 ms
コンパイル使用メモリ 82,392 KB
実行使用メモリ 64,132 KB
最終ジャッジ日時 2024-11-20 10:09:01
合計ジャッジ時間 3,498 ms
ジャッジサーバーID
(参考情報)
judge3 / judge1
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 48 ms
61,460 KB
testcase_01 AC 63 ms
63,248 KB
testcase_02 AC 48 ms
62,824 KB
testcase_03 AC 47 ms
62,212 KB
testcase_04 AC 47 ms
61,596 KB
testcase_05 AC 47 ms
61,852 KB
testcase_06 AC 47 ms
62,072 KB
testcase_07 AC 46 ms
62,560 KB
testcase_08 AC 50 ms
62,972 KB
testcase_09 AC 47 ms
61,068 KB
testcase_10 AC 48 ms
61,196 KB
testcase_11 AC 60 ms
62,072 KB
testcase_12 AC 61 ms
63,712 KB
testcase_13 AC 50 ms
63,436 KB
testcase_14 AC 54 ms
62,696 KB
testcase_15 AC 48 ms
62,076 KB
testcase_16 AC 62 ms
63,500 KB
testcase_17 AC 49 ms
62,436 KB
testcase_18 AC 61 ms
63,604 KB
testcase_19 AC 64 ms
62,664 KB
testcase_20 AC 62 ms
63,156 KB
testcase_21 AC 58 ms
62,584 KB
testcase_22 AC 47 ms
62,736 KB
testcase_23 AC 47 ms
61,672 KB
testcase_24 AC 47 ms
62,948 KB
testcase_25 AC 48 ms
61,832 KB
testcase_26 AC 48 ms
61,684 KB
testcase_27 AC 47 ms
62,224 KB
testcase_28 AC 48 ms
61,924 KB
testcase_29 AC 58 ms
62,588 KB
testcase_30 AC 58 ms
63,196 KB
testcase_31 AC 49 ms
62,068 KB
testcase_32 AC 56 ms
61,780 KB
testcase_33 AC 53 ms
62,376 KB
testcase_34 AC 57 ms
62,512 KB
testcase_35 AC 59 ms
62,076 KB
testcase_36 AC 55 ms
63,064 KB
testcase_37 AC 60 ms
62,336 KB
testcase_38 AC 54 ms
62,200 KB
testcase_39 AC 59 ms
64,132 KB
testcase_40 AC 62 ms
62,824 KB
testcase_41 AC 55 ms
62,156 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

from collections import defaultdict, deque, Counter
from heapq import heapify, heappop, heappush
import math
from copy import deepcopy
from itertools import combinations, permutations, product, combinations_with_replacement
from bisect import bisect_left, bisect_right

import sys

def input():
    return sys.stdin.readline().rstrip()
def getN():
    return int(input())
def getNM():
    return map(int, input().split())
def getList():
    return list(map(int, input().split()))
def getListGraph():
    return list(map(lambda x:int(x) - 1, input().split()))
def getArray(intn):
    return [int(input()) for i in range(intn)]

mod = 10 ** 9 + 7
MOD = 998244353
# sys.setrecursionlimit(10000000)
inf = float('inf')
eps = 10 ** (-15)
dy = [0, 1, 0, -1]
dx = [1, 0, -1, 0]

#############
# Main Code #
#############

def prime_factorize(n):
    divisors = []
    # 27(2 * 2 * 7)の7を出すためにtemp使う
    temp = n
    for i in range(2, int(math.sqrt(n)) + 1):
        if temp % i == 0:
            cnt = 0
            while temp % i == 0:
                cnt += 1
                # 素因数を見つけるたびにtempを割っていく
                temp //= i
            divisors.append([i, cnt])
    if temp != 1:
        divisors.append([temp, 1])
    if divisors == []:
        divisors.append([n, 1])

    return divisors

"""
x = 0なら余りは0
nCkは求められない どんな数か?
Mの因数はいくつあるか?という問題 分子と分母でMの個数を数えてみる
nCk = n! / (n - r)!r! とすると簡単
8*7*6*5*4*3*2*1/3*2*1 * 5*4*3*2*1

まず素数で割る
"""

N, K, M = getNM()

# 素因数の数を数える
def calc(x, m):
    tar = m
    res = 0
    while tar <= x:
        res += x // tar
        tar *= m
    return res

def count(n, r, m):
    return calc(n, m) - calc(n - r, m) - calc(r, m)

ans = inf
# 各因数について見ていく
for k, v in prime_factorize(M):
    ans = min(ans, count(N, K, k) // v)
print(ans)
0