結果

問題 No.5007 Steiner Space Travel
ユーザー hirayuu_yc
提出日時 2023-07-20 17:54:58
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 297 ms / 1,000 ms
コード長 1,752 bytes
コンパイル時間 846 ms
コンパイル使用メモリ 87,088 KB
実行使用メモリ 86,316 KB
スコア 7,636,089
最終ジャッジ日時 2023-07-20 17:55:10
合計ジャッジ時間 10,302 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge14
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

from random import randint
inf=998244353

def calc(a,b):
    return ((star[a][0]-star[b][0])**2+(star[a][1]-star[b][1])**2)*star[a][2]*star[b][2]

def sell(strict):
    pos=0
    go={0}
    ans=[0]
    ansdis=0
    while len(go)<N:
        next=-1
        dis=inf
        for i in range(N):
            if i not in go:
                if dist[pos][i][0]<dis:
                    next=i
                    dis=dist[pos][i][0]
        ans.append(next)
        pos=next
        ansdis+=dis
        go.add(pos)
    ans.append(0)
    while strict:
        imp=0
        l=-1
        r=-1
        for i in range(1,N):
            for j in range(i+1,N):
                pre=(dist[ans[i-1]][ans[i]][0]+dist[ans[j]][ans[j+1]][0])-(dist[ans[i-1]][ans[j]][0]+dist[ans[i]][ans[j+1]][0])
                if pre>imp:
                    l=i
                    r=j
                    imp=pre
        if l==-1:
            break
        for i in range((r-l+1)//2):
            ans[l+i],ans[r-i]=ans[r-i],ans[l+i]
    return ansdis,ans

N,M=map(int,input().split())
star=[tuple(map(int,input().split()+["5"])) for i in range(N)]+[[randint(0,1000),randint(0,1000),1] for i in range(M)]
dist=[[[calc(i,j),[j]] for i in range(N+M)]for j in range(N+M)]
for i in range(N+M):
    for j in range(N+M):
        for k in range(N+M):
            if dist[j][i][0]+dist[i][k][0]<dist[j][k][0]:
                dist[j][k]=[dist[j][i][0]+dist[i][k][0],dist[j][i][1]+dist[i][k][1]]
for i in range(M):
    print(*star[N+i][:2])
ansdis,ans=sell(True)
finalans=[]
for i in range(len(ans)-1):
    for j in dist[ans[i]][ans[i+1]][1]:
        finalans.append(j)
finalans.append(0)
print(len(finalans))
for i in finalans:
    if i>=N:
        print(2,i-N+1)
    else:
        print(1,i+1)
0