結果

問題 No.2012 Largest Triangle
ユーザー googol_S0googol_S0
提出日時 2022-07-15 23:03:57
言語 Python3
(3.12.2 + numpy 1.26.4 + scipy 1.12.0)
結果
TLE  
実行時間 -
コード長 1,716 bytes
コンパイル時間 307 ms
コンパイル使用メモリ 11,192 KB
実行使用メモリ 97,172 KB
最終ジャッジ日時 2023-09-10 04:32:49
合計ジャッジ時間 6,002 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 22 ms
13,760 KB
testcase_01 AC 23 ms
9,356 KB
testcase_02 AC 22 ms
9,204 KB
testcase_03 AC 22 ms
9,316 KB
testcase_04 AC 22 ms
9,248 KB
testcase_05 AC 22 ms
9,384 KB
testcase_06 AC 22 ms
9,360 KB
testcase_07 AC 22 ms
9,400 KB
testcase_08 AC 23 ms
9,232 KB
testcase_09 AC 22 ms
9,252 KB
testcase_10 AC 23 ms
9,428 KB
testcase_11 AC 24 ms
9,412 KB
testcase_12 AC 24 ms
9,268 KB
testcase_13 AC 24 ms
9,316 KB
testcase_14 AC 24 ms
9,404 KB
testcase_15 AC 24 ms
9,424 KB
testcase_16 TLE -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

from decimal import *
def sgn(x):
  if x==0:
    return 0
  if x>0:
    return 1
  return -1

def check(a,b,c):
  #if a==b or b==c:
  #  return sgn(b[0]-a[0])*sgn(c[1]-b[1])>=sgn(c[0]-b[0])*sgn(b[1]-a[1])
  return (b[0]-a[0])*(c[1]-b[1])>=(c[0]-b[0])*(b[1]-a[1])

def get_y(a,x):
  return a[0]*x+a[1]

from collections import deque
class CHT:
  __slots__=['H','F']
  def __init__(self,f):
    self.H=deque()
    self.F=f
  
  def add(self,x):
    h=self.H
    if self.F:
      x=(-x[0],-x[1])
    if len(h)==0:
      h.append(x)
      return 0
    if h[0][0]<=x[0]:
      if h[0][0]==x[0]:
        if h[0][1]<=x[1]:
          return 0
        h.popleft()
      while(len(h)>=2 and check(x,h[0],h[1])):
        h.popleft()
      h.appendleft(x)
    else:
      if(h[-1][0]==x[0]):
        if(h[-1][1]<=x[1]):
          return 0
        h.pop()
      while(len(h)>=2 and check(h[-2],h[-1],x)):
        h.pop()
      h.append(x)
  
  def query_i(self,x):
    h=self.H
    a=h[0]
    if len(h)>=2:
      b=h[1]
      while(len(h)>=2 and a[0]*x+a[1]>=b[0]*x+b[1]):
        h.popleft()
        a=b
        if len(h)<2:
          break
        b=h[1]
    if self.F:
      return -a[0]*x-a[1]
    else:
      return a[0]*x+a[1]

getcontext().prec=50
N=int(input())
eps=Decimal('1e-18')
X=[tuple(map(lambda x:Decimal(int(x)+eps),input().split())) for i in range(N)]
A=CHT(True)
B=CHT(False)
X.sort(key=lambda x:x[1])
for i in range(N):
  a,b=X[i][0],X[i][1]
  A.add((b,-a))
  B.add((b,-a))
X.sort(key=lambda x:x[0]/x[1])
ANS=0
for i in range(N):
  a,b=X[i][0],X[i][1]
  x,y=B.query_i(a/b),A.query_i(a/b)
  ANS=max(ANS,-x*b,-y*b,x*b,y*b)
  if b>=0:
    ANS=max(ANS,x*b,-y*b)
  else:
    ANS=max(ANS,-x*b,y*b)
print(round(ANS))
0