結果

問題 No.1809 Divide NCK
コンテスト
ユーザー ygd.
提出日時 2022-01-14 21:56:48
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 48 ms / 2,000 ms
コード長 1,474 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 272 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 57,728 KB
最終ジャッジ日時 2024-11-20 09:57:22
合計ジャッジ時間 3,073 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys
#input = sys.stdin.readline #文字列につけてはダメ
#input = sys.stdin.buffer.readline #文字列につけてはダメ
#sys.setrecursionlimit(1000000)
#import bisect
#import itertools
#import random
#from heapq import heapify, heappop, heappush
#from collections import defaultdict 
#from collections import deque
#import copy
#import math
#from functools import lru_cache
#MOD = pow(10,9) + 7
#MOD = 998244353

def prime_factorize(n):
    ret = []
    #if n == 1:
    #    ret.append((1,1))
    #    return ret
    cnt = 0
    while n % 2 == 0:
        cnt += 1
        n //= 2
    if cnt > 0:
        ret.append((2,cnt))
    i = 3
    while i * i <= n:
        cnt = 0
        while n % i == 0:
            cnt += 1
            n //= i
        else:
            if cnt != 0: #cnt==0の時は追加しない
                ret.append((i,cnt))
            i += 2
    if n != 1:
        ret.append((n,1))
    return ret #(素因数、何乗)

def count(A,p): #A!に素因数pが何個含まれるか
    ret = 0
    while A > 0:
        ret += A//p
        A = A//p
    #print(ret)
    return ret

def main():
    N,K,M = map(int,input().split())
    L = prime_factorize(M)
    #print(L)
    LN = []
    for x,c in L:
        cnt = count(N,x) - count(K,x) - count(N-K,x)
        #print("x",x,"c",c,"cnt",cnt)
        LN.append(cnt//c)
    ans = min(LN)
    if ans < 0:
        print(0)
    else:
        print(ans)



if __name__ == '__main__':
    main()
0