結果

問題 No.1125 Without Parallelogram
ユーザー maspymaspy
提出日時 2020-07-22 22:09:24
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 138 ms / 1,000 ms
コード長 619 bytes
コンパイル時間 154 ms
コンパイル使用メモリ 10,976 KB
実行使用メモリ 30,168 KB
最終ジャッジ日時 2023-09-04 20:50:32
合計ジャッジ時間 10,675 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 134 ms
29,744 KB
testcase_01 AC 128 ms
29,876 KB
testcase_02 AC 128 ms
29,876 KB
testcase_03 AC 128 ms
29,888 KB
testcase_04 AC 132 ms
29,808 KB
testcase_05 AC 132 ms
29,916 KB
testcase_06 AC 133 ms
29,780 KB
testcase_07 AC 133 ms
29,892 KB
testcase_08 AC 129 ms
29,912 KB
testcase_09 AC 126 ms
29,856 KB
testcase_10 AC 124 ms
29,884 KB
testcase_11 AC 129 ms
29,956 KB
testcase_12 AC 129 ms
29,992 KB
testcase_13 AC 129 ms
29,872 KB
testcase_14 AC 126 ms
29,884 KB
testcase_15 AC 129 ms
29,844 KB
testcase_16 AC 129 ms
29,840 KB
testcase_17 AC 132 ms
30,052 KB
testcase_18 AC 129 ms
30,092 KB
testcase_19 AC 129 ms
30,124 KB
testcase_20 AC 131 ms
30,088 KB
testcase_21 AC 129 ms
30,100 KB
testcase_22 AC 132 ms
30,108 KB
testcase_23 AC 129 ms
30,096 KB
testcase_24 AC 128 ms
30,088 KB
testcase_25 AC 131 ms
30,104 KB
testcase_26 AC 133 ms
30,104 KB
testcase_27 AC 138 ms
30,168 KB
testcase_28 AC 133 ms
29,860 KB
testcase_29 AC 132 ms
29,800 KB
testcase_30 AC 129 ms
29,924 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import numpy as np

read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

def prime_table(N):
    is_prime = np.zeros(N, np.int64)
    is_prime[2:3] = 1
    is_prime[3::2] = 1
    for p in range(3, N, 2):
        if p * p >= N:
            break
        if is_prime[p]:
            is_prime[p * p::p + p] = 0
    return is_prime, np.where(is_prime)[0]

_, primes = prime_table(10**4)

N = int(read())

def main(N):
    p = primes[primes > N][0]
    x = np.arange(N)
    y = (x * x) % p
    return x, y

x, y = main(N)
for a, b in zip(x, y):
    print(a, b)
0