結果
問題 | No.1581 Multiple Sequence |
ユーザー | ygd. |
提出日時 | 2021-07-02 22:38:56 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,123 ms / 2,000 ms |
コード長 | 1,507 bytes |
コンパイル時間 | 439 ms |
コンパイル使用メモリ | 87,300 KB |
実行使用メモリ | 378,660 KB |
最終ジャッジ日時 | 2023-09-11 22:59:47 |
合計ジャッジ時間 | 18,666 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 157 ms
93,396 KB |
testcase_01 | AC | 162 ms
93,016 KB |
testcase_02 | AC | 958 ms
328,816 KB |
testcase_03 | AC | 1,041 ms
346,056 KB |
testcase_04 | AC | 463 ms
177,332 KB |
testcase_05 | AC | 1,091 ms
359,876 KB |
testcase_06 | AC | 631 ms
228,244 KB |
testcase_07 | AC | 452 ms
174,692 KB |
testcase_08 | AC | 542 ms
201,220 KB |
testcase_09 | AC | 978 ms
336,708 KB |
testcase_10 | AC | 756 ms
265,072 KB |
testcase_11 | AC | 333 ms
135,652 KB |
testcase_12 | AC | 568 ms
209,180 KB |
testcase_13 | AC | 215 ms
109,128 KB |
testcase_14 | AC | 1,034 ms
354,000 KB |
testcase_15 | AC | 816 ms
286,224 KB |
testcase_16 | AC | 369 ms
150,400 KB |
testcase_17 | AC | 1,097 ms
370,440 KB |
testcase_18 | AC | 315 ms
133,052 KB |
testcase_19 | AC | 950 ms
324,664 KB |
testcase_20 | AC | 960 ms
329,828 KB |
testcase_21 | AC | 1,016 ms
345,324 KB |
testcase_22 | AC | 701 ms
252,216 KB |
testcase_23 | AC | 1,123 ms
378,660 KB |
ソースコード
from collections import defaultdict from functools import lru_cache import sys sys.setrecursionlimit(50000000) def prime_list(n): is_prime = [0] * (n + 1) #素数かどうか。1なら素数。 is_prime[2] = 1 for i in range(3, n + 1, 2): is_prime[i] = 1 for i in range(2, int(n ** 0.5) + 1): if is_prime[i]: for j in range(i * i, n + 1, 2 * i): is_prime[j] = 0 #最後に全部並べる。 d = [] for j in range(n + 1): if is_prime[j] == 1: d.append(j) return d, is_prime def main(): dic = defaultdict(list) MAX = pow(10,5) + 100 L,dummy = prime_list(MAX) for i in L: a = i while a <= MAX: dic[a].append(i) a += i #print(dic) M = int(input()); MOD = pow(10,9) + 7 @lru_cache(maxsize=None) def solve(p,dic=dic,MOD=MOD): if p == 1: return 1 if p == 2: return 2 ret = solve(p-1) #print(p,ret) n = len(dic[p]) #print(dic[p]) for i in range(1,1<<n): val = 1 ppcnt = 0 for j in range(n): if (i>>j)&1 == 1: val *= dic[p][j] ppcnt += 1 if ppcnt%2 == 1: ret += solve(p//val) else: ret -= solve(p//val) ret %= MOD return ret%MOD print(solve(M)%MOD) if __name__ == '__main__': main()