結果

問題 No.5007 Steiner Space Travel
ユーザー hirayuu_yc
提出日時 2023-07-20 16:20:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 938 ms / 1,000 ms
コード長 1,278 bytes
コンパイル時間 294 ms
コンパイル使用メモリ 87,004 KB
実行使用メモリ 144,800 KB
スコア 7,280,096
最終ジャッジ日時 2023-07-20 16:20:35
合計ジャッジ時間 26,782 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge13
純コード判定しない問題か言語
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
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():
    pos=0
    go={0}
    ans=[]
    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+=dist[pos][next][1]
        pos=next
        ansdis+=dis
        go.add(pos)
    ans+=dist[pos][0][1]
    return ansdis,ans+[0]

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)]
prest=[]
preans=[]
predis=inf
t=7
for _ in range(t):
	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]]
	ansdis,ans=sell()
	if predis>ansdis:
		predis=ansdis
		preans=ans[:]
		prest.clear()
		for i in range(M):
			prest.append(star[N+i][:2])
for i in range(M):
	print(*prest[i])
print(len(preans))
for i in preans:
	if i>=N:
		print(2,i-N+1)
	else:
		print(1,i+1)
0