結果

問題 No.913 木の燃やし方
ユーザー roarisroaris
提出日時 2023-03-16 05:40:29
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 4,084 bytes
コンパイル時間 177 ms
コンパイル使用メモリ 81,772 KB
実行使用メモリ 240,408 KB
最終ジャッジ日時 2024-09-18 09:10:31
合計ジャッジ時間 47,809 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 42 ms
54,524 KB
testcase_01 AC 43 ms
54,656 KB
testcase_02 AC 42 ms
53,888 KB
testcase_03 AC 193 ms
79,028 KB
testcase_04 AC 200 ms
79,584 KB
testcase_05 WA -
testcase_06 AC 189 ms
79,024 KB
testcase_07 AC 158 ms
77,736 KB
testcase_08 AC 1,502 ms
191,420 KB
testcase_09 AC 1,433 ms
187,844 KB
testcase_10 AC 1,408 ms
186,860 KB
testcase_11 AC 1,455 ms
190,212 KB
testcase_12 AC 1,531 ms
198,592 KB
testcase_13 AC 1,425 ms
191,852 KB
testcase_14 AC 1,386 ms
185,928 KB
testcase_15 AC 1,385 ms
186,720 KB
testcase_16 AC 1,472 ms
213,876 KB
testcase_17 AC 1,457 ms
212,384 KB
testcase_18 AC 1,450 ms
210,112 KB
testcase_19 AC 1,268 ms
240,408 KB
testcase_20 AC 1,474 ms
220,000 KB
testcase_21 WA -
testcase_22 WA -
testcase_23 WA -
testcase_24 WA -
testcase_25 AC 1,241 ms
226,100 KB
testcase_26 AC 1,415 ms
198,360 KB
testcase_27 AC 1,188 ms
217,072 KB
testcase_28 AC 1,489 ms
182,900 KB
testcase_29 AC 1,547 ms
188,528 KB
testcase_30 WA -
testcase_31 AC 1,280 ms
210,812 KB
testcase_32 AC 1,288 ms
214,988 KB
testcase_33 AC 1,296 ms
218,316 KB
testcase_34 AC 1,598 ms
201,692 KB
testcase_35 AC 1,550 ms
202,540 KB
testcase_36 AC 1,582 ms
209,696 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