結果

問題 No.2438 Double Least Square
ユーザー chineristACchineristAC
提出日時 2023-08-19 00:51:04
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,026 bytes
コンパイル時間 196 ms
コンパイル使用メモリ 82,428 KB
実行使用メモリ 77,752 KB
最終ジャッジ日時 2024-05-06 07:53:56
合計ジャッジ時間 4,525 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 40 ms
52,992 KB
testcase_01 AC 41 ms
53,248 KB
testcase_02 AC 39 ms
53,248 KB
testcase_03 AC 41 ms
53,120 KB
testcase_04 AC 40 ms
52,864 KB
testcase_05 AC 39 ms
53,504 KB
testcase_06 AC 40 ms
52,992 KB
testcase_07 AC 39 ms
53,120 KB
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
testcase_20 AC 136 ms
77,280 KB
testcase_21 AC 138 ms
77,288 KB
testcase_22 AC 131 ms
77,256 KB
testcase_23 AC 138 ms
77,412 KB
testcase_24 AC 140 ms
77,420 KB
testcase_25 AC 134 ms
77,344 KB
testcase_26 AC 136 ms
77,408 KB
testcase_27 AC 137 ms
77,272 KB
testcase_28 AC 141 ms
77,408 KB
testcase_29 AC 141 ms
77,336 KB
testcase_30 AC 39 ms
53,248 KB
testcase_31 AC 39 ms
53,120 KB
testcase_32 AC 39 ms
53,504 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from itertools import permutations
from heapq import *
from math import acos

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

N = int(input())
H = int(input())
point = [tuple(mi()) for i in range(N)]

def calc(f,g):
    res = 0

    a,b,c = 0,0,0
    for i in f:
        x,y = point[i]
        a += x*x
        b += 2*x*(H-y)
        c += (H-y)*(H-y)
    if a!=0:
        t = (-b/(2*a))
        res += a * t * t + b * t + c

    a,b,c = 0,0,0
    for i in g:
        x,y = point[i]
        a += x*x
        b += 2 * x * y
        c += y*y
    if a!=0:
        t = (-b/(2*a))
        res += a * t * t + b * t + c

    return res

arg = []
for i in range(N):
    x,y = point[i]
    r = ((y-H/2)**2 + x**2)**.5
    arg.append(acos((H/2-y)/r))

idx = [i for i in range(N)]
idx.sort(key=lambda i:arg[i])

ans = calc([i for i in range(N)],[])
for i in range(N):
    ans = min(ans,calc([idx[j] for j in range(i+1,N)],[idx[j] for j in range(i+1)]))

print(ans)
0