結果
問題 | No.890 移調の限られた旋法 |
ユーザー | convexineq |
提出日時 | 2019-09-20 22:14:06 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 168 ms / 2,000 ms |
コード長 | 1,174 bytes |
コンパイル時間 | 678 ms |
コンパイル使用メモリ | 87,260 KB |
実行使用メモリ | 99,888 KB |
最終ジャッジ日時 | 2023-10-12 19:37:06 |
合計ジャッジ時間 | 7,377 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 143 ms
99,552 KB |
testcase_01 | AC | 150 ms
99,436 KB |
testcase_02 | AC | 154 ms
99,556 KB |
testcase_03 | AC | 150 ms
99,500 KB |
testcase_04 | AC | 152 ms
99,560 KB |
testcase_05 | AC | 153 ms
99,524 KB |
testcase_06 | AC | 153 ms
99,384 KB |
testcase_07 | AC | 143 ms
99,424 KB |
testcase_08 | AC | 157 ms
99,524 KB |
testcase_09 | AC | 154 ms
99,284 KB |
testcase_10 | AC | 154 ms
99,384 KB |
testcase_11 | AC | 160 ms
99,600 KB |
testcase_12 | AC | 154 ms
99,568 KB |
testcase_13 | AC | 165 ms
99,888 KB |
testcase_14 | AC | 153 ms
99,384 KB |
testcase_15 | AC | 156 ms
99,404 KB |
testcase_16 | AC | 161 ms
99,828 KB |
testcase_17 | AC | 162 ms
99,788 KB |
testcase_18 | AC | 165 ms
99,820 KB |
testcase_19 | AC | 156 ms
99,444 KB |
testcase_20 | AC | 158 ms
99,404 KB |
testcase_21 | AC | 156 ms
99,412 KB |
testcase_22 | AC | 153 ms
99,380 KB |
testcase_23 | AC | 163 ms
99,596 KB |
testcase_24 | AC | 168 ms
99,532 KB |
testcase_25 | AC | 157 ms
99,692 KB |
testcase_26 | AC | 166 ms
99,388 KB |
testcase_27 | AC | 161 ms
99,584 KB |
testcase_28 | AC | 161 ms
99,524 KB |
testcase_29 | AC | 159 ms
99,544 KB |
testcase_30 | AC | 162 ms
99,880 KB |
testcase_31 | AC | 158 ms
99,812 KB |
testcase_32 | AC | 160 ms
99,868 KB |
testcase_33 | AC | 160 ms
99,876 KB |
testcase_34 | AC | 153 ms
99,736 KB |
ソースコード
SIZE=1000000; MOD=10**9+7 #998244353 #ここを変更する SIZE += 1 inv = [0]*SIZE # inv[j] = j^{-1} mod MOD fac = [0]*SIZE # fac[j] = j! mod MOD finv = [0]*SIZE # finv[j] = (j!)^{-1} mod MOD inv[1] = 1 fac[0] = fac[1] = 1 finv[0] = finv[1] = 1 for i in range(2,SIZE): inv[i] = MOD - (MOD//i)*inv[MOD%i]%MOD fac[i] = fac[i-1]*i%MOD finv[i]= finv[i-1]*inv[i]%MOD def choose(n,r): # nCk mod MOD の計算 if 0 <= r <= n: return (fac[n]*finv[r]%MOD)*finv[n-r]%MOD else: return 0 from math import gcd def divisor_list(N): #約数のリスト if N == 1: return [1] res = [] for i in range(1,N): if i*i >= N: break if N%i == 0: res.append(i) res.append(N//i) if i*i == N: res.append(i) return sorted(res) n,k = [int(i) for i in input().split()] g = gcd(n,k) if g == 1: print(0) exit() div = divisor_list(n) div = [i for i in div if g%i==0] d = {i:choose(n//i,k//i) for i in div} d[1] = 0 #print(d) for i in reversed(div): for j in div: if j%i == 0 and i != j: d[i] -= d[j] d[i] %= MOD #print(d) print(-d[1]%MOD)