結果

問題 No.913 木の燃やし方
ユーザー roarisroaris
提出日時 2023-03-16 05:40:29
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 4,084 bytes
コンパイル時間 170 ms
コンパイル使用メモリ 81,768 KB
実行使用メモリ 239,844 KB
最終ジャッジ日時 2023-10-18 12:50:14
合計ジャッジ時間 46,988 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 44 ms
55,524 KB
testcase_01 AC 42 ms
55,524 KB
testcase_02 AC 41 ms
55,524 KB
testcase_03 AC 197 ms
78,964 KB
testcase_04 AC 200 ms
78,888 KB
testcase_05 WA -
testcase_06 AC 186 ms
78,920 KB
testcase_07 AC 158 ms
77,392 KB
testcase_08 AC 1,499 ms
189,948 KB
testcase_09 AC 1,429 ms
187,312 KB
testcase_10 AC 1,416 ms
187,176 KB
testcase_11 AC 1,446 ms
190,248 KB
testcase_12 AC 1,518 ms
198,184 KB
testcase_13 AC 1,454 ms
191,496 KB
testcase_14 AC 1,383 ms
185,876 KB
testcase_15 AC 1,363 ms
186,772 KB
testcase_16 AC 1,472 ms
213,612 KB
testcase_17 AC 1,432 ms
210,920 KB
testcase_18 AC 1,401 ms
210,040 KB
testcase_19 AC 1,220 ms
239,844 KB
testcase_20 AC 1,456 ms
220,296 KB
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 AC 1,228 ms
225,220 KB
testcase_26 AC 1,420 ms
197,236 KB
testcase_27 AC 1,192 ms
217,020 KB
testcase_28 AC 1,502 ms
183,252 KB
testcase_29 AC 1,541 ms
188,356 KB
testcase_30 WA -
testcase_31 AC 1,313 ms
211,124 KB
testcase_32 AC 1,326 ms
215,252 KB
testcase_33 AC 1,299 ms
218,004 KB
testcase_34 AC 1,612 ms
201,196 KB
testcase_35 AC 1,572 ms
202,748 KB
testcase_36 AC 1,574 ms
209,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
input = sys.stdin.readline
from collections import *

class ConvexHullTrick:
    def __init__(self, getmin=1, a_inc=0, x_inc=1):
        self.getmin = getmin
        self.a_inc = a_inc
        self.x_inc = x_inc
        self.ab = deque()
        self.inf = (1 << 63)-1
        self.px = self.inf*(-2*x_inc+1)
 
    def add_line(self, a, b):
        if not self.ab:
            self.ab.append((a, b))
            return
        if self.getmin ^ self.a_inc:
            if self.ab and a == self.ab[-1][0]:
                if (b-self.ab[-1][1])*(self.getmin*2-1) >= 0: return
                self.ab.pop()
            if self.x_inc:
                while len(self.ab) > 1:
                    c, d = self.ab.pop()
                    e, f = self.ab[-1]
                    if (b-d)*(e-c) > (d-f)*(c-a):
                        self.ab.append((c, d))
                        break
            else:
                (c, d), x = self.ab[-1], self.px
                if (a*x+b-c*x-d)*(self.getmin*2-1) >= 0: return
            self.ab.append((a, b))
        else:
            if self.ab and a == self.ab[0][0]:
                if (b-self.ab[0][1])*(self.getmin*2-1) >= 0: return
                self.ab.popleft()
            if not self.x_inc:
                while len(self.ab) > 1:
                    c, d = self.ab.popleft()
                    e, f = self.ab[0]
                    if (b-d)*(e-c) < (d-f)*(c-a):
                        self.ab.append((c, d))
                        break
            else:
                (c, d), x = self.ab[0], self.px
                if (a*x+b-c*x-d)*(self.getmin*2-1) >= 0: return
            self.ab.appendleft((a, b))
 
    def __getitem__(self, x):
        self.px = x
        if not self.ab: return self.inf*(self.getmin*2-1)
        if self.x_inc:
            while len(self.ab) > 1:
                a, b = self.ab.popleft()
                c, d = self.ab[0]
                if (a*x+b-c*x-d)*(self.getmin*2-1) < 0:
                    self.ab.appendleft((a, b))
                    return a*x+b
        else:
            while len(self.ab) > 1:
                a, b = self.ab.pop()
                c, d = self.ab[-1]
                if (a*x+b-c*x-d)*(self.getmin*2-1) < 0:
                    self.ab.append((a, b))
                    return a*x+b
        a, b = self.ab[0]
        return a*x+b

class CHT: #傾きの単調性と最小値クエリのxの単調性を仮定
    def __init__(self):
        self.q = deque([]) #直線群
        
    def f(self, f1, x): #直線f1のxの値
        return f1[0]*x+f1[1]
    
    def check(self, f1, f2, f3): #f2を削除しても良いかの判定
        return (f2[0]-f1[0])*(f3[1]-f2[1])>=(f2[1]-f1[1])*(f3[0]-f2[0])
    
    def add_line(self, a, b): #傾きa, 切片bの直線を追加
        while len(self.q)>=2 and self.check(self.q[-2], self.q[-1], (a, b)):
            self.q.pop()
        
        self.q.append((a, b))
    
    def get(self, x): #xでの最小値
        while len(self.q)>=2 and self.f(self.q[0], x)>=self.f(self.q[1], x):
            self.q.popleft()
        
        return self.f(self.q[0], x)

def dfs(l, r): # [l, r)
    global ans
    
    if r-l==1:
        ans[l] = min(ans[l], 1+A[l])
        return

    m = (l+r)//2
    # [l, m)
    cht = ConvexHullTrick()
    left = defaultdict(lambda: 10**18)
    
    for i in range(m+1, r+1):
        cht.add_line(-2*i, i*i+acc[i])
    
    for i in range(l, m):
        left[i] = cht[i]+i*i-acc[i]
        left[i] = min(left[i], left[i-1])
    
    # [m, r)
    cht = ConvexHullTrick()
    right = defaultdict(lambda: 10**18)
    
    for i in range(l, m+1):
        cht.add_line(-2*i, i*i-acc[i])
    
    for i in range(r, m, -1):
        right[i] = cht[i]+i*i+acc[i]
        right[i] = min(right[i], right[i+1])
    
    for i in range(l, r):
        ans[i] = min(ans[i], left[i], right[i+1])
    
    dfs(l, m)
    dfs(m, r)
    
N = int(input())
A = list(map(int, input().split()))
acc = [0]

for Ai in A:
    acc.append(acc[-1]+Ai)

ans = [10**18]*N
dfs(0, N)

for i in range(N):
    print(ans[i])
0