結果
問題 | No.890 移調の限られた旋法 |
ユーザー | brthyyjp |
提出日時 | 2021-03-26 07:07:19 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 139 ms / 2,000 ms |
コード長 | 715 bytes |
コンパイル時間 | 332 ms |
コンパイル使用メモリ | 82,176 KB |
実行使用メモリ | 98,588 KB |
最終ジャッジ日時 | 2024-05-05 17:34:52 |
合計ジャッジ時間 | 4,449 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 77 ms
90,368 KB |
testcase_01 | AC | 71 ms
90,496 KB |
testcase_02 | AC | 72 ms
90,624 KB |
testcase_03 | AC | 73 ms
90,368 KB |
testcase_04 | AC | 71 ms
90,752 KB |
testcase_05 | AC | 70 ms
90,624 KB |
testcase_06 | AC | 71 ms
90,752 KB |
testcase_07 | AC | 72 ms
90,240 KB |
testcase_08 | AC | 72 ms
90,368 KB |
testcase_09 | AC | 71 ms
90,240 KB |
testcase_10 | AC | 71 ms
90,752 KB |
testcase_11 | AC | 73 ms
90,368 KB |
testcase_12 | AC | 72 ms
90,880 KB |
testcase_13 | AC | 132 ms
98,304 KB |
testcase_14 | AC | 89 ms
98,588 KB |
testcase_15 | AC | 90 ms
98,432 KB |
testcase_16 | AC | 112 ms
98,432 KB |
testcase_17 | AC | 94 ms
97,792 KB |
testcase_18 | AC | 107 ms
97,888 KB |
testcase_19 | AC | 91 ms
94,720 KB |
testcase_20 | AC | 78 ms
93,056 KB |
testcase_21 | AC | 74 ms
90,752 KB |
testcase_22 | AC | 88 ms
96,512 KB |
testcase_23 | AC | 88 ms
98,152 KB |
testcase_24 | AC | 82 ms
94,976 KB |
testcase_25 | AC | 75 ms
91,520 KB |
testcase_26 | AC | 89 ms
98,048 KB |
testcase_27 | AC | 89 ms
98,528 KB |
testcase_28 | AC | 83 ms
95,616 KB |
testcase_29 | AC | 80 ms
94,004 KB |
testcase_30 | AC | 92 ms
97,408 KB |
testcase_31 | AC | 94 ms
94,720 KB |
testcase_32 | AC | 112 ms
97,300 KB |
testcase_33 | AC | 139 ms
97,664 KB |
testcase_34 | AC | 106 ms
97,408 KB |
ソースコード
N = 10**6+50 mod = 10**9+7 fac = [1]*(N+1) finv = [1]*(N+1) for i in range(N): fac[i+1] = fac[i] * (i+1) % mod finv[-1] = pow(fac[-1], mod-2, mod) for i in reversed(range(N)): finv[i] = finv[i+1] * (i+1) % mod def cmb1(n, r, mod): if r <0 or r > n: return 0 r = min(r, n-r) return fac[n] * finv[r] * finv[n-r] % mod n, k = map(int, input().split()) X = [0]*(n+1) ans = 0 for i in range(1, n): if n%i == 0: p = n//i if k%p == 0: x = cmb1(i, k//p, mod) X[i] += x ans += X[i] j = 2*i while j <= n: if n%j == 0: X[j] -= X[i] j += i #print(X) print(ans%mod)