結果
問題 | No.2955 Pizza Delivery Plan |
ユーザー |
![]() |
提出日時 | 2024-11-09 03:10:27 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 330 ms / 2,000 ms |
コード長 | 1,027 bytes |
コンパイル時間 | 434 ms |
コンパイル使用メモリ | 82,560 KB |
実行使用メモリ | 81,280 KB |
最終ジャッジ日時 | 2024-11-09 03:10:36 |
合計ジャッジ時間 | 8,363 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
import sysinput = sys.stdin.readlineN,K=map(int,input().split())P=[list(map(int,input().split())) for i in range(N)]DP=[[1<<60]*(1<<N) for i in range(N)]DIS=[[0]*N for i in range(N)]for i in range(N):for j in range(N):x,y=P[i]z,w=P[j]DIS[i][j]=((x-z)**2 + (y-w)**2)**(1/2)D=[0]*Nfor i in range(N):x,y=P[i]D[i]=(x**2+y**2)**(1/2)for i in range(N):DP[i][1<<i]=D[i]for j in range(1<<N):for i in range(N):for k in range(N):if j & (1<<k) != 0:continueto=j | (1<<k)if to.bit_count()>K:continueDP[k][to]=min(DP[k][to],DP[i][j]+DIS[i][k])DIS_COME=[1<<60]*(1<<N)for i in range(N):for j in range(1<<N):DIS_COME[j]=min(DIS_COME[j],DP[i][j]+D[i])DP2=[1<<60]*(1<<N)DP2[0]=0# 3^NのDP(部分集合を列挙)for i in range(1<<N):p=iwhile p!=0:p=(p-1) & iDP2[i]=min(DP2[i],DP2[p]+DIS_COME[p^i])print(DP2[-1])