結果
問題 | No.414 衝動 |
ユーザー | GrayCoder |
提出日時 | 2018-08-19 15:48:08 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 79 ms / 1,000 ms |
コード長 | 1,165 bytes |
コンパイル時間 | 100 ms |
コンパイル使用メモリ | 12,672 KB |
実行使用メモリ | 16,128 KB |
最終ジャッジ日時 | 2024-11-15 09:19:40 |
合計ジャッジ時間 | 1,642 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 57 ms
14,976 KB |
testcase_01 | AC | 33 ms
10,880 KB |
testcase_02 | AC | 33 ms
10,880 KB |
testcase_03 | AC | 78 ms
16,128 KB |
testcase_04 | AC | 63 ms
16,128 KB |
testcase_05 | AC | 79 ms
16,000 KB |
testcase_06 | AC | 78 ms
15,872 KB |
testcase_07 | AC | 30 ms
10,880 KB |
testcase_08 | AC | 59 ms
13,696 KB |
testcase_09 | AC | 79 ms
16,128 KB |
testcase_10 | AC | 31 ms
10,880 KB |
testcase_11 | AC | 34 ms
11,136 KB |
testcase_12 | AC | 32 ms
10,880 KB |
ソースコード
from sys import stdin, stdout input = lambda: stdin.readline() write = stdout.write def main(): M = int(input()) primes = [i for i in primes235(int(M ** 0.5))] for i in primes: d, mod = divmod(M, i) if not mod: print(d, i) break else: print(1, M) def primes235(limit): yield 2; yield 3; yield 5 if limit < 7: return modPrms = [7,11,13,17,19,23,29,31] gaps = [4,2,4,2,4,6,2,6,4,2,4,2,4,6,2,6] ndxs = [0,0,0,0,1,1,2,2,2,2,3,3,4,4,4,4,5,5,5,5,5,5,6,6,7,7,7,7,7,7] lmtbf = (limit + 23) // 30 * 8 - 1 lmtsqrt = (int(limit ** 0.5) - 7) lmtsqrt = lmtsqrt // 30 * 8 + ndxs[lmtsqrt % 30] buf = [True] * (lmtbf + 1) for i in range(lmtsqrt + 1): if buf[i]: ci = i & 7; p = 30 * (i >> 3) + modPrms[ci] s = p * p - 7; p8 = p << 3 for j in range(8): c = s // 30 * 8 + ndxs[s % 30] buf[c::p8] = [False] * ((lmtbf - c) // p8 + 1) s += p * gaps[ci]; ci += 1 for i in range(lmtbf - 6 + (ndxs[(limit - 7) % 30])): if buf[i]: yield (30 * (i >> 3) + modPrms[i & 7]) main()