結果

問題 No.97 最大の値を求めるくえり
ユーザー maspy
提出日時 2023-05-12 21:15:58
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
AC  
実行時間 3,868 ms / 5,000 ms
コード長 1,281 bytes
コンパイル時間 1,128 ms
コンパイル使用メモリ 12,416 KB
実行使用メモリ 50,408 KB
最終ジャッジ日時 2024-12-26 11:58:49
合計ジャッジ時間 53,721 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

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