結果

問題 No.1555 Constructed Balancing Sequence
コンテスト
ユーザー chineristAC
提出日時 2021-01-06 07:32:12
言語 PyPy3
(7.3.17)
コンパイル:
pypy3 -mpy_compile _filename_
実行:
pypy3 _filename_
結果
RE  
実行時間 -
コード長 1,552 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 151 ms
コンパイル使用メモリ 85,120 KB
実行使用メモリ 135,528 KB
最終ジャッジ日時 2026-03-07 00:46:33
合計ジャッジ時間 4,651 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 15 RE * 27
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

def perferct_solve(N,K,A):
    mod = 998244353

    diff = [A[i] for i in range(N)]
    cum = [A[i] for i in range(N)]
    S = A[0]
    for i in range(1,N):
        diff[i] = S - A[i]
        if diff[i] < 0:
            return 0
        S += A[i]
        cum[i] = S

    dp = [[0 for minus in range(3*K+1)] for i in range(N)]

    for minus in range(3*K+1):
        first = A[0] - minus
        if -K<=first<=K:
            dp[0][minus] = 1
        if minus:
            dp[0][minus] += dp[0][minus-1]
            dp[0][minus] %= mod

    for i in range(1,N):
        for minus in range(3*K+1):

            if diff[i] and -K<=A[i]-minus<=K:
                dp[i][minus] += dp[i-1][0]
                dp[i][minus] %= mod

            if minus%2==diff[i] and -K<=(cum[i]-minus)//2<=K:
                dp[i][minus] += dp[i-1][3*K] - dp[i-1][minus//2+diff[i]]
                dp[i][minus] %= mod

            if diff[i]<=1:
                pre_minus_L = max(diff[i],-K+minus-A[i])
                pre_minus_R = min((minus+diff[i])//2,K+minus-A[i])
                if pre_minus_L<=pre_minus_R:
                    dp[i][minus] += dp[i-1][pre_minus_R] - dp[i-1][pre_minus_L-1] * (pre_minus_L>0)
                    dp[i][minus] %= mod

            if minus:
                dp[i][minus] += dp[i][minus-1]
                dp[i][minus] %= mod

    return dp[N-1][0]

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

assert 2<=N<=1500
assert 1<=K<=1500
assert len(A)==N
for i in range(N):
    assert abs(A[i])<=K

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