結果

問題 No.913 木の燃やし方
ユーザー vwxyzvwxyz
提出日時 2024-04-17 15:54:25
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 911 ms / 3,000 ms
コード長 2,047 bytes
コンパイル時間 293 ms
コンパイル使用メモリ 82,460 KB
実行使用メモリ 136,900 KB
最終ジャッジ日時 2024-10-09 09:58:18
合計ジャッジ時間 28,035 ms
ジャッジサーバーID
(参考情報)
judge4 / judge3
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
53,504 KB
testcase_01 AC 42 ms
53,760 KB
testcase_02 AC 41 ms
53,760 KB
testcase_03 AC 163 ms
77,780 KB
testcase_04 AC 168 ms
77,888 KB
testcase_05 AC 154 ms
77,908 KB
testcase_06 AC 149 ms
77,760 KB
testcase_07 AC 122 ms
76,728 KB
testcase_08 AC 824 ms
109,180 KB
testcase_09 AC 775 ms
108,096 KB
testcase_10 AC 787 ms
107,732 KB
testcase_11 AC 798 ms
108,784 KB
testcase_12 AC 843 ms
112,220 KB
testcase_13 AC 810 ms
109,468 KB
testcase_14 AC 777 ms
107,940 KB
testcase_15 AC 767 ms
107,704 KB
testcase_16 AC 831 ms
112,540 KB
testcase_17 AC 799 ms
111,128 KB
testcase_18 AC 806 ms
110,736 KB
testcase_19 AC 546 ms
133,200 KB
testcase_20 AC 793 ms
113,076 KB
testcase_21 AC 884 ms
110,464 KB
testcase_22 AC 898 ms
110,464 KB
testcase_23 AC 911 ms
109,408 KB
testcase_24 AC 830 ms
120,424 KB
testcase_25 AC 526 ms
123,380 KB
testcase_26 AC 813 ms
112,176 KB
testcase_27 AC 575 ms
135,764 KB
testcase_28 AC 857 ms
110,848 KB
testcase_29 AC 861 ms
110,208 KB
testcase_30 AC 891 ms
110,720 KB
testcase_31 AC 690 ms
133,628 KB
testcase_32 AC 702 ms
135,740 KB
testcase_33 AC 695 ms
136,900 KB
testcase_34 AC 858 ms
111,104 KB
testcase_35 AC 870 ms
110,976 KB
testcase_36 AC 845 ms
110,080 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

from collections import deque

class Convex_Hull_Trick_deque:
    """
    f_i = a_ix + b_i とする。f_i の追加および、min_i f(x) の取得ができるデータ構造。
    ただし、傾き a_i は降順に追加されなければならない。
    また、クエリ x も昇順に実行されなければならない。
    """

    def __init__(self):
        self.lines=deque()

    def add_line(self,a,b):
        lines=self.lines
        while len(lines) >= 2:
            a1,b1=lines[-2]
            a2,b2=lines[-1]
            if (a2-a1)*(b-b2)<(b2-b1)*(a-a2):
                break
            lines.pop()
        lines.append((a, b))

    def __call__(self, x):
        lines=self.lines
        a,b=lines[0]
        y=a*x+b
        while len(lines) >= 2:
            a2, b2 = lines[1]
            y2 = a2 * x + b2
            if y < y2:
                break
            y = y2
            lines.popleft()
        return y

N=int(input())
A=list(map(int,input().split()))
C=[0]+A
for i in range(1,N+1):
    C[i]+=C[i-1]
inf=1<<60
ans_lst=[inf]*N
def solve(l,r):
    if r-l==1:
        ans_lst[l]=min(ans_lst[l],1+A[l])
    elif r-l==2:
        ans_lst[l]=min(ans_lst[l],1+A[l],4+A[l]+A[l+1])
        ans_lst[l+1]=min(ans_lst[l+1],1+A[l+1],4+A[l]+A[l+1])
    else:
        m=(l+r)//2
        mi=[inf]*(r-l)
        CHT_left=Convex_Hull_Trick_deque()
        for i in range(l,m):
            CHT_left.add_line(-2*i,i*i-C[i])
        for j in range(m+1,r+1):
            mi[j-1-l]=min(mi[j-1-l],CHT_left(j)+j*j+C[j])
        CHT_right=Convex_Hull_Trick_deque()
        for j in range(m+1,r+1):
            CHT_right.add_line(-2*j,j*j+C[j])
        for i in range(l,m):
            mi[i-l]=min(mi[i-l],CHT_right(i)+i*i-C[i])

        for i in range(l+1,m):
            mi[i-l]=min(mi[i-l],mi[i-l-1])
        for j in range(r-2,m-1,-1):
            mi[j-l]=min(mi[j-l],mi[j-l+1])
        for i in range(l,r):
            ans_lst[i]=min(ans_lst[i],mi[i-l])
        solve(l,m)
        solve(m,r)
solve(0,N)
print(*ans_lst,sep="\n")
0