結果
問題 |
No.1809 Divide NCK
|
ユーザー |
|
提出日時 | 2022-01-14 23:06:10 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 68 ms / 2,000 ms |
コード長 | 1,860 bytes |
コンパイル時間 | 228 ms |
コンパイル使用メモリ | 82,076 KB |
実行使用メモリ | 63,980 KB |
最終ジャッジ日時 | 2024-11-20 13:51:24 |
合計ジャッジ時間 | 3,850 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 39 |
ソースコード
# import sys # input = sys.stdin.readline def mp():return map(int,input().split()) def lmp():return list(map(int,input().split())) def mps(A):return [tuple(map(int, input().split())) for _ in range(A)] def stoi(LIST):return list(map(int,LIST)) def itos(LIST):return list(map(str,LIST)) def bitA(X,A):return X & 1<<A == 1<<A import math import bisect import heapq import time from copy import copy as cc from copy import deepcopy as dc from itertools import accumulate from collections import Counter, defaultdict, deque def ceil(U,V):return (U+V-1)//V def modf1(N,MOD):return (N-1)%MOD+1 inf = int(1e18+20) mod = int(1e9+7) def seg_insu(x,M): return insu(x,M)-insu(x-k,M) def insu(x,M): ans = 0 now = M while now <= x: ans += x//now now *= M return ans def factorization(n): arr = [] temp = n for i in range(2, int(-(-n**0.5//1))+1): if temp%i==0: cnt=0 while temp%i==0: cnt+=1 temp //= i arr.append([i, cnt]) if temp!=1: arr.append([temp, 1]) if arr==[]: arr.append([n, 1]) return arr n,k,m = mp() mfact = factorization(m) #print(mfact) ans = 0 bl = len(mfact) u = [] for i,j in mfact: u.append((seg_insu(n,i)-seg_insu(k,i))//j) print(min(u)) # for i in range(1,1<<bl): # tar = [] # cnt = 1 # for j in range(bl): # if bitA(i,j): # cnt += mfact[j][0]**mfact[j][1] # tar.append(mfact[j]) # bunsi = [] # bunbo = [] # y = [] # y.append(cnt) # u = inf # v = inf # for num,f in tar: # u = min(u,seg_insu()) # now = inf # for num, f in tar: # now = min(now, (seg_insu(n,num)-seg_insu(k,num))//f) # if len(tar) % 2 == 1: # ans += now # else: # ans -= now # print(ans,tar) #print(ans)