結果

問題 No.1809 Divide NCK
ユーザー 👑 KazunKazun
提出日時 2022-01-14 21:31:03
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 49 ms / 2,000 ms
コード長 1,374 bytes
コンパイル時間 209 ms
コンパイル使用メモリ 82,560 KB
実行使用メモリ 58,112 KB
最終ジャッジ日時 2024-11-20 07:57:22
合計ジャッジ時間 3,008 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 39
権限があれば一括ダウンロードができます

ソースコード

diff # 

#素因数分解
def Prime_Factorization(N):
    if N==0:
        return [[0,1]]

    if N<0:
        R=[[-1,1]]
    else:
        R=[]

    N=abs(N)

    if N&1==0:
        C=0
        while N&1==0:
            N>>=1
            C+=1
        R.append([2,C])

    if N%3==0:
        C=0
        while N%3==0:
            N//=3
            C+=1
        R.append([3,C])

    k=5
    Flag=0
    while k*k<=N:
        if N%k==0:
            C=0
            while N%k==0:
                C+=1
                N//=k
            R.append([k,C])
        k+=2+2*Flag
        Flag^=1

    if N!=1:
        R.append([N,1])

    return R

#根基
def Radical(N):
    """Nが素因数分解 N=p^a*q^b*r^c ...となるとき, pqr...を返す.

    N:非負整数
    """

    assert N>=0,"Nは非負整数ではない."
    a=1

    if N&1==0:
        a*=2
        while N&1==0:
            N>>=1

    if N%3==0:
        a*=3
        while N%3==0:
            N//=3

    k=5
    Flag=0
    while k*k<=N:
        if N%k==0:
            a*=k
            while N%k==0:
                N//=k
        k+=2+2*Flag
        Flag^=1

    if N>1:
        a*=N
    return a

def f(a,N):
    res=0
    while a:
        res+=a//N
        a//=N
    return res

N,K,M=map(int,input().split())

Ans=float("inf")
for p,e in Prime_Factorization(M):
    Ans=min(Ans,(f(N,p)-f(N-K,p)-f(K,p))//e)

print(Ans)
0