結果

問題 No.97 最大の値を求めるくえり
ユーザー maspymaspy
提出日時 2023-05-12 21:15:58
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
AC  
実行時間 2,947 ms / 5,000 ms
コード長 1,281 bytes
コンパイル時間 232 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 50,396 KB
最終ジャッジ日時 2024-05-06 10:43:49
合計ジャッジ時間 38,992 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,333 ms
45,792 KB
testcase_01 AC 1,828 ms
45,848 KB
testcase_02 AC 2,213 ms
50,276 KB
testcase_03 AC 2,084 ms
50,268 KB
testcase_04 AC 1,417 ms
49,892 KB
testcase_05 AC 1,297 ms
49,764 KB
testcase_06 AC 1,292 ms
49,884 KB
testcase_07 AC 1,349 ms
49,924 KB
testcase_08 AC 1,435 ms
49,888 KB
testcase_09 AC 1,555 ms
49,760 KB
testcase_10 AC 1,913 ms
49,884 KB
testcase_11 AC 2,289 ms
50,020 KB
testcase_12 AC 2,559 ms
49,764 KB
testcase_13 AC 2,947 ms
49,760 KB
testcase_14 AC 1,819 ms
50,396 KB
testcase_15 AC 1,898 ms
50,272 KB
testcase_16 AC 1,929 ms
50,280 KB
testcase_17 AC 1,863 ms
50,048 KB
testcase_18 AC 1,835 ms
50,324 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
import numpy as np

MOD = 100003

N, _, *Q = map(int, read().split())
mask = (1 << 32) - 1


def generate_A(N):
    xor128_x = 123456789
    xor128_y = 362436069
    xor128_z = 521288629
    xor128_w = 88675123

    for _ in range(N):
        t = xor128_x ^ (xor128_x << 11)
        t &= mask
        xor128_x = xor128_y
        xor128_y = xor128_z
        xor128_z = xor128_w
        xor128_w = xor128_w ^ (xor128_w >> 19) ^ (t ^ (t >> 8))
        xor128_w &= mask
        yield xor128_w % MOD


A = np.array(list(generate_A(N)), np.int32)


def strategy_1(A):
    answer = [0] * MOD
    S = np.zeros_like(A)
    for i in range(MOD):
        answer[i] = S.max()
        S += A
        S[S >= MOD] -= MOD
    return answer


def strategy_2(A):
    answer = [0] * MOD
    isin_A = np.zeros(MOD, np.bool_)
    isin_A[A] = 1
    x = np.arange(99400, MOD)
    S = np.zeros_like(x)
    for i in range(1, MOD):
        inv = pow(i, MOD - 2, MOD)
        S += x
        S[S >= MOD] -= MOD
        answer[inv] = 99400 + np.where(isin_A[S])[0][-1]
    return answer


if N < 2000:
    answer = strategy_1(A)
else:
    answer = strategy_2(A)
for q in Q:
    print(answer[q])
0