結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1,334 ms
45,928 KB
testcase_01 AC 1,756 ms
45,716 KB
testcase_02 AC 2,053 ms
50,400 KB
testcase_03 AC 2,018 ms
50,404 KB
testcase_04 AC 1,438 ms
49,888 KB
testcase_05 AC 1,288 ms
50,284 KB
testcase_06 AC 1,320 ms
50,016 KB
testcase_07 AC 1,365 ms
50,404 KB
testcase_08 AC 1,456 ms
50,396 KB
testcase_09 AC 1,606 ms
50,144 KB
testcase_10 AC 1,876 ms
50,144 KB
testcase_11 AC 2,270 ms
50,488 KB
testcase_12 AC 2,632 ms
50,272 KB
testcase_13 AC 2,944 ms
50,400 KB
testcase_14 AC 1,811 ms
50,656 KB
testcase_15 AC 1,839 ms
50,276 KB
testcase_16 AC 1,866 ms
50,652 KB
testcase_17 AC 1,804 ms
50,540 KB
testcase_18 AC 1,869 ms
50,656 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