結果

問題 No.5007 Steiner Space Travel
コンテスト
ユーザー Edomonndo365
提出日時 2023-04-25 12:27:46
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 307 ms / 1,000 ms
コード長 2,453 bytes
コンパイル時間 460 ms
コンパイル使用メモリ 87,180 KB
実行使用メモリ 81,568 KB
スコア 6,629,940
最終ジャッジ日時 2023-04-25 12:28:02
合計ジャッジ時間 11,689 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import sys
from time import time
from random import randrange

START = time()
INF=10**9

class Node:
    def __init__(self, isTerminal: bool, idx: int, pos: list):
        self.isTerminal = isTerminal
        self.idx = idx
        self.pos = pos
        self.dist = INF
    def __str__(self):
        return "{} {}".format(2-self.isTerminal, self.idx)

def get_time(START):
    return time() - START

def dist(p1,p2):
    x1,y1=p1
    x2,y2=p2
    return (x1-x2)**2 + (y1-y2)**2

def calc(idx: int):
    d = dist(ans[idx].pos, ans[idx+1].pos)
    if ans[idx].isTerminal and ans[idx+1].isTerminal:
        a = 25
    elif ans[idx].isTerminal or ans[idx+1].isTerminal:
        a = 5
    else:
        a = 1
    return a * d

N,M=map(int,input().split())
Terminals = [list(map(int,input().split())) for _ in range(N)]

ans = [Node(True,1,Terminals[0])]
candidates = set(range(1,N))
for _ in range(N-1):
    best_index = -1
    best_dist = INF
    for cand in candidates:
        d = dist(ans[-1].pos, Terminals[cand])
        if d < best_dist:
            best_dist = d
            best_index = cand
    assert best_index!=-1
    ans.append(Node(True, best_index+1, Terminals[best_index]))
    candidates.discard(best_index)

ans.append(Node(True,1,Terminals[0]))

n = len(ans)
loop_cnt1 = 0
while get_time(START) < 0.2:
    loop_cnt1+=1
    v1 = randrange(1,n-1)
    v2 = randrange(1,n-1)
    while v1==v2:
        v2 = randrange(1,n-1)
    current_score = calc(v1-1) + calc(v1) + calc(v2-1) + calc(v2)
    ans[v1],ans[v2] = ans[v2],ans[v1]
    changed_score = calc(v1-1) + calc(v1) + calc(v2-1) + calc(v2)
    if changed_score < current_score:
        continue
    else:
        ans[v1],ans[v2] = ans[v2],ans[v1]

candidates = []
for i in range(N):
    d = dist(ans[i].pos, ans[i+1].pos)
    candidates.append([i,d])
candidates.sort(key=lambda x:x[1], reverse=1)

Stations = []
for i,cand in enumerate(candidates[:M], start=1):
    idx = cand[0]
    Stations.append(Node(False, i, ans[idx].pos))


for station in Stations:
    for i in range(len(ans)):
        if ans[i].pos == station.pos:
            ans.insert(i+1, station)
            break

bestScore = 0
loop_cnt2 = 0


for node in Stations:
    print(*node.pos)
print(len(ans))
for node in ans:
    print(node)

print("Score", bestScore, file=sys.stderr)
print("Loop1", loop_cnt1, file=sys.stderr)
print("Loop2", loop_cnt2, file=sys.stderr)
print("time:",get_time(START) * 1000, file=sys.stderr)
0