結果

問題 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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

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