結果

問題 No.1809 Divide NCK
ユーザー customaddonecustomaddone
提出日時 2022-01-14 22:00:19
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 131 ms / 2,000 ms
コード長 1,981 bytes
コンパイル時間 371 ms
コンパイル使用メモリ 87,144 KB
実行使用メモリ 77,520 KB
最終ジャッジ日時 2023-08-12 19:35:05
合計ジャッジ時間 7,400 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 116 ms
72,876 KB
testcase_01 AC 131 ms
77,348 KB
testcase_02 AC 115 ms
73,308 KB
testcase_03 AC 113 ms
73,064 KB
testcase_04 AC 114 ms
73,248 KB
testcase_05 AC 116 ms
73,440 KB
testcase_06 AC 115 ms
72,904 KB
testcase_07 AC 112 ms
73,444 KB
testcase_08 AC 115 ms
73,312 KB
testcase_09 AC 114 ms
73,128 KB
testcase_10 AC 115 ms
73,352 KB
testcase_11 AC 130 ms
77,340 KB
testcase_12 AC 131 ms
77,396 KB
testcase_13 AC 116 ms
77,424 KB
testcase_14 AC 120 ms
77,404 KB
testcase_15 AC 117 ms
77,428 KB
testcase_16 AC 131 ms
77,376 KB
testcase_17 AC 118 ms
77,520 KB
testcase_18 AC 131 ms
76,996 KB
testcase_19 AC 130 ms
77,412 KB
testcase_20 AC 127 ms
77,384 KB
testcase_21 AC 129 ms
77,000 KB
testcase_22 AC 116 ms
73,340 KB
testcase_23 AC 112 ms
72,860 KB
testcase_24 AC 111 ms
73,092 KB
testcase_25 AC 114 ms
73,148 KB
testcase_26 AC 112 ms
73,516 KB
testcase_27 AC 113 ms
73,252 KB
testcase_28 AC 113 ms
73,356 KB
testcase_29 AC 128 ms
77,244 KB
testcase_30 AC 131 ms
77,260 KB
testcase_31 AC 118 ms
77,232 KB
testcase_32 AC 126 ms
77,196 KB
testcase_33 AC 128 ms
77,260 KB
testcase_34 AC 128 ms
76,888 KB
testcase_35 AC 127 ms
77,320 KB
testcase_36 AC 122 ms
77,304 KB
testcase_37 AC 128 ms
76,892 KB
testcase_38 AC 123 ms
77,128 KB
testcase_39 AC 131 ms
77,324 KB
testcase_40 AC 129 ms
77,324 KB
testcase_41 AC 127 ms
77,160 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