結果
問題 | No.705 ゴミ拾い Hard |
ユーザー | vwxyz |
提出日時 | 2023-05-28 04:17:32 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,335 ms / 1,500 ms |
コード長 | 2,170 bytes |
コンパイル時間 | 387 ms |
コンパイル使用メモリ | 87,072 KB |
実行使用メモリ | 169,052 KB |
最終ジャッジ日時 | 2023-08-27 03:28:57 |
合計ジャッジ時間 | 25,849 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge15 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 72 ms
71,344 KB |
testcase_01 | AC | 70 ms
71,524 KB |
testcase_02 | AC | 71 ms
71,216 KB |
testcase_03 | AC | 70 ms
71,156 KB |
testcase_04 | AC | 71 ms
71,308 KB |
testcase_05 | AC | 71 ms
71,308 KB |
testcase_06 | AC | 71 ms
71,196 KB |
testcase_07 | AC | 71 ms
71,340 KB |
testcase_08 | AC | 72 ms
71,540 KB |
testcase_09 | AC | 69 ms
71,296 KB |
testcase_10 | AC | 71 ms
71,492 KB |
testcase_11 | AC | 70 ms
71,096 KB |
testcase_12 | AC | 71 ms
71,272 KB |
testcase_13 | AC | 70 ms
71,388 KB |
testcase_14 | AC | 127 ms
78,168 KB |
testcase_15 | AC | 123 ms
77,640 KB |
testcase_16 | AC | 127 ms
78,024 KB |
testcase_17 | AC | 127 ms
77,760 KB |
testcase_18 | AC | 123 ms
78,080 KB |
testcase_19 | AC | 134 ms
78,512 KB |
testcase_20 | AC | 131 ms
78,224 KB |
testcase_21 | AC | 126 ms
77,776 KB |
testcase_22 | AC | 136 ms
78,516 KB |
testcase_23 | AC | 143 ms
79,600 KB |
testcase_24 | AC | 1,294 ms
162,452 KB |
testcase_25 | AC | 1,289 ms
162,888 KB |
testcase_26 | AC | 1,275 ms
162,736 KB |
testcase_27 | AC | 1,283 ms
162,620 KB |
testcase_28 | AC | 1,267 ms
162,856 KB |
testcase_29 | AC | 1,274 ms
162,624 KB |
testcase_30 | AC | 1,335 ms
162,668 KB |
testcase_31 | AC | 1,215 ms
162,544 KB |
testcase_32 | AC | 1,201 ms
162,340 KB |
testcase_33 | AC | 1,085 ms
147,396 KB |
testcase_34 | AC | 1,059 ms
147,432 KB |
testcase_35 | AC | 1,078 ms
139,068 KB |
testcase_36 | AC | 1,128 ms
143,580 KB |
testcase_37 | AC | 1,091 ms
169,052 KB |
testcase_38 | AC | 1,140 ms
163,160 KB |
testcase_39 | AC | 1,218 ms
163,344 KB |
testcase_40 | AC | 67 ms
71,276 KB |
testcase_41 | AC | 67 ms
71,044 KB |
testcase_42 | AC | 70 ms
71,044 KB |
testcase_43 | AC | 1,142 ms
165,004 KB |
ソースコード
import sys readline=sys.stdin.readline def monotone_minima(N,cost,cnt=0,inf=1<<60): dp=[inf]*(N+1) dp[0]=0 if cnt: for _ in range(cnt): prev=dp dp=[inf]*(N+1) idx=[None]*(N+1) dp[0]=0 idx[0]=0 for n in range(N+1): if prev[n]+cost(n,N)<dp[N]: dp[N]=prev[n]+cost(n,N) idx[N]=n while inf in dp: lst=[i for i in range(N+1) if dp[i]!=inf] for l,r in zip(lst,lst[1:]): mid=(l+r)//2 for n in range(idx[l],min(idx[r],mid)+1): if prev[n]+cost(n,mid)<dp[mid]: dp[mid]=prev[n]+cost(n,mid) idx[mid]=n else: def transition(L,R): if R-L<=1: return mid=(L+R)//2 transition(L,mid) update=[inf]*(R-mid) idx=[inf]*(R-mid) for n in range(L,mid): if dp[n]+cost(n,mid)<update[0]: update[0]=dp[n]+cost(n,mid) idx[0]=n if dp[n]+cost(n,R-1)<update[R-mid-1]: update[R-mid-1]=dp[n]+cost(n,R-1) idx[R-mid-1]=n while inf in update: lst=[i+mid for i in range(R-mid) if update[i]!=inf] for l,r in zip(lst,lst[1:]): m=(l+r)//2 for n in range(idx[l-mid],min(idx[r-mid],m)+1): if dp[n]+cost(n,m)<update[m-mid]: update[m-mid]=dp[n]+cost(n,m) idx[m-mid]=n for n in range(mid,R): dp[n]=min(dp[n],update[n-mid]) transition(mid,R) transition(0,N+1) return dp[N] N=int(readline()) A=list(map(int,readline().split())) X=list(map(int,readline().split())) Y=list(map(int,readline().split())) inf=1<<60 def cost(l,r): if r-l<0: return inf if l==r: return 0 return abs(X[l]-A[r-1])**3+abs(Y[l])**3 ans=monotone_minima(N,cost,inf=inf) print(ans)