結果

問題 No.890 移調の限られた旋法
ユーザー neterukun
提出日時 2019-09-21 11:32:50
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 682 ms / 2,000 ms
コード長 2,050 bytes
コンパイル時間 166 ms
コンパイル使用メモリ 81,920 KB
実行使用メモリ 146,836 KB
最終ジャッジ日時 2024-09-17 14:56:19
合計ジャッジ時間 24,620 ms
ジャッジサーバーID
(参考情報)
judge2 / judge6
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

class Combination:
'''MOD
O(N)
nCkO(1)'''
def __init__(self, n, MOD):
self.fact = [1]
for i in range(1, n + 1):
self.fact.append(self.fact[-1] * i % MOD)
self.inv_fact = [pow(self.fact[i], MOD - 2, MOD) for i in range(n + 1)]
self.MOD = MOD
def factorial(self, k):
"""k! O(1)"""
return self.fact[k]
def inverse_factorial(self, k):
"""k! O(1)"""
return self.inv_fact[k]
def permutation(self, k, r):
"""kPr O(1)"""
if k < r:
return 0
return (self.fact[k] * self.inv_fact[r]) % self.MOD
def combination(self, k, r):
"""kCr O(1)"""
if k < r:
return 0
return (self.fact[k] * self.inv_fact[k - r] * self.inv_fact[r]) % self.MOD
def make_divisors(n):
"""n
: O(sqrt(N))
: 12 -> [1, 2, 3, 4, 6, 12]
"""
divisors = []
for k in range(1, int(n**0.5) + 1):
if n % k == 0:
divisors.append(k)
if k != n // k:
divisors.append(n // k)
divisors = sorted(divisors)
return divisors
def gcd(a, b):
"""a, b(greatest common divisor:GCD)"""
if b == 0:
return a
return gcd(b, a%b)
n, k = map(int, input().split())
MOD = 10**9 + 7
comb = Combination(10**6 + 10, MOD)
gcd_nk = gcd(n, k)
div_list = make_divisors(gcd_nk)
li = []
memo = {}
for num in div_list:
if num == 1:
continue
total = n // num
selected = k // num
tmp = comb.combination(total, selected)
li.append(num)
memo[num] = tmp
li = li[::-1]
for i in range(len(li)):
for j in range(i):
if li[j] % li[i] == 0:
memo[li[i]] -= memo[li[j]]
memo[li[i]] %= MOD
ans = 0
for i in memo:
ans += memo[i]
ans %= MOD
print(ans)
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0