結果

問題 No.1555 Constructed Balancing Sequence
コンテスト
ユーザー maspy
提出日時 2021-01-11 00:42:12
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
(最新)
AC  
(最初)
実行時間 -
コード長 1,588 bytes
コンパイル時間 523 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 30,680 KB
最終ジャッジ日時 2024-06-22 19:29:17
合計ジャッジ時間 4,279 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 10 TLE * 1 -- * 31
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
import itertools

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

MOD = 998_244_353

def is_balanced(A, Acum):
    for a, b in zip(A[1:], Acum[1:]):
        if a > b:
            return False
    return True

def main(N, K, A):
    S = [0] + list(itertools.accumulate(A))
    if not is_balanced(A, S):
        return 0
    MAX = K + 10
    # x 以下の数列は対策済
    dp = [0] * (MAX + MAX + 1)
    v = A[-1] - 1
    dp[v] = 1
    for n in range(N - 2, -1, -1):
        newdp = [0] * (MAX + MAX + 1)
        imos = [0] * (MAX + MAX + 1)
        for v in range(-MAX, MAX + 1):
            if not dp[v]:
                continue
            if v < S[n] + A[n] - 1:
                low = A[n]
            else:
                low = -K
            """
            for b in range(low, A[n] + 1):
                to = max(b - 1, v // 2)
                newdp[to] += dp[v]
            """
            # low <= b <= An and b <= v//2 + 1
            L, R = low, min(A[n], v // 2 + 1)
            if L <= R:
                newdp[v // 2] += dp[v] * (R - L + 1)
            # low <= b <= An and v//2+2 <= b
            L, R = max(low, v // 2 + 2), A[n]
            L, R = L - 1, R - 1
            if L <= R:
                imos[L] += dp[v]
                imos[R + 1] -= dp[v]
        for i in range(-MAX, MAX):
            imos[i + 1] += imos[i]
            newdp[i] += imos[i]
        dp = newdp
    return sum(dp) % MOD

N, K = map(int, readline().split())
A = list(map(int, readline().split()))

print(main(N, K, A))
0