結果

問題 No.2376 障害物競プロ
ユーザー mkawa2mkawa2
提出日時 2023-07-07 23:08:10
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,941 bytes
コンパイル時間 1,824 ms
コンパイル使用メモリ 85,924 KB
実行使用メモリ 126,864 KB
最終ジャッジ日時 2023-09-29 00:51:53
合計ジャッジ時間 54,736 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 76 ms
71,876 KB
testcase_01 AC 76 ms
71,664 KB
testcase_02 AC 78 ms
71,920 KB
testcase_03 AC 78 ms
71,396 KB
testcase_04 AC 364 ms
79,520 KB
testcase_05 AC 658 ms
84,124 KB
testcase_06 AC 591 ms
84,052 KB
testcase_07 AC 1,168 ms
85,828 KB
testcase_08 AC 1,155 ms
85,992 KB
testcase_09 AC 1,517 ms
93,320 KB
testcase_10 AC 1,741 ms
95,532 KB
testcase_11 AC 3,242 ms
126,864 KB
testcase_12 AC 2,857 ms
122,264 KB
testcase_13 AC 1,133 ms
86,812 KB
testcase_14 AC 1,251 ms
85,924 KB
testcase_15 AC 1,591 ms
96,276 KB
testcase_16 AC 1,715 ms
98,052 KB
testcase_17 AC 3,125 ms
124,304 KB
testcase_18 AC 3,052 ms
123,352 KB
testcase_19 TLE -
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 -- -
testcase_42 -- -
testcase_43 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys

# sys.setrecursionlimit(200005)
# sys.set_int_max_str_digits(1000005)
int1 = lambda x: int(x)-1
pDB = lambda *x: print(*x, end="\n", file=sys.stderr)
p2D = lambda x: print(*x, sep="\n", end="\n\n", file=sys.stderr)
def II(): return int(sys.stdin.readline())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def LI1(): return list(map(int1, sys.stdin.readline().split()))
def LLI1(rows_number): return [LI1() for _ in range(rows_number)]
def SI(): return sys.stdin.readline().rstrip()

dij = [(0, 1), (-1, 0), (0, -1), (1, 0)]
# dij = [(0, 1), (-1, 0), (0, -1), (1, 0), (1, 1), (1, -1), (-1, 1), (-1, -1)]
inf = (1 << 63)-1
# md = 10**9+7
md = 998244353

def ccw(s, t, p):
    x1 = p[0]-s[0]
    y1 = p[1]-s[1]
    x2 = t[0]-s[0]
    y2 = t[1]-s[1]
    if x1*y2-x2*y1 > 0: return 1
    if x1*y2-x2*y1 < 0: return -1
    if x1*x2+y1*y2 < 0: return -2
    if x1**2+y1**2 > x2**2+y2**2: return 2
    return False

# 共有点の有無のみ
def common_point(a, b, s, t):
    cc = []
    cc.append(ccw(a, b, s))
    cc.append(ccw(a, b, t))
    cc.append(ccw(s, t, a))
    cc.append(ccw(s, t, b))
    if any(c == 0 for c in cc): return True
    if any(abs(c) == 2 for c in cc): return False
    if cc[0] == cc[1] or cc[2] == cc[3]: return False
    return True

co = 10**10

n, m = LI()
xy, st = [], []
for _ in range(n):
    x, y, s, t = LI()
    xy.append((x*co, y*co))
    st.append((s*co, t*co))

dist = [[inf]*2*n for _ in range(2*n)]

def can(a, b, i, j):
    for k in range(n):
        if k == i or k == j: continue
        if common_point(a, b, xy[k], st[k]): return False
    return True

for i in range(n):
    dist[i*2][i*2] = 0
    dist[i*2+1][i*2+1] = 0
    if can(xy[i], st[i], i, i):
        x, y = xy[i]
        s, t = st[i]
        dist[i*2][i*2+1] = dist[i*2+1][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
    for j in range(i):
        if can(xy[i], xy[j], i, j):
            x, y = xy[i]
            s, t = xy[j]
            dist[i*2][j*2] = dist[j*2][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(xy[i], st[j], i, j):
            x, y = xy[i]
            s, t = st[j]
            dist[i*2][j*2+1] = dist[j*2+1][i*2] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(st[i], xy[j], i, j):
            x, y = st[i]
            s, t = xy[j]
            dist[i*2+1][j*2] = dist[j*2][i*2+1] = round(((s-x)**2+(t-y)**2)**0.5)
        if can(st[i], st[j], i, j):
            x, y = st[i]
            s, t = st[j]
            dist[i*2+1][j*2+1] = dist[j*2+1][i*2+1] = round(((s-x)**2+(t-y)**2)**0.5)
# p2D(dist)

for j in range(n*2):
    for i in range(n*2):
        ij = dist[i][j]
        if ij == inf: continue
        for k in range(n*2):
            jk = dist[j][k]
            if ij == inf: continue
            dist[i][k] = min(dist[i][k], ij+jk)

for _ in range(m):
    i, b, j, d = LI1()
    ans = dist[i*2+b][j*2+d]
    print(ans/co)
0