#889111 (PyPy3) No.2376 障害物競プロ

• 問題: No.2376 障害物競プロ
• ユーザー: mkawa2
• 提出日時: 2023-07-08 11:19:16
• 言語: PyPy3 (7.3.15)
• 結果: AC
• 実行時間: 3,968 ms / 4,000 ms
• コード長: 3,389 bytes
• コンパイル時間: 271 ms
• コンパイル使用メモリ: 82,240 KB
• 実行使用メモリ: 90,828 KB
• 最終ジャッジ日時: 2024-07-26 15:12:02
• 合計ジャッジ時間: 84,310 ms
• ジャッジサーバーID （参考情報）: judge3 / judge2

ファイルパターン	結果
sample	AC * 4
other	AC * 40

ソースコード

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
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))
    xy.append((x, y))
    st.append((s, t))
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)
        dist[i*2][i*2+1] = dist[i*2+1][i*2] = ((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)
            dist[i*2][j*2] = dist[j*2][i*2] = ((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)
            dist[i*2][j*2+1] = dist[j*2+1][i*2] = ((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)
            dist[i*2+1][j*2] = dist[j*2][i*2+1] = ((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)
            dist[i*2+1][j*2+1] = dist[j*2+1][i*2+1] = ((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 jk == 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)
    print(ans)
 
 
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX