結果
問題 | No.2955 Pizza Delivery Plan |
ユーザー |
👑 |
提出日時 | 2024-05-28 19:08:09 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 738 ms / 2,000 ms |
コード長 | 1,343 bytes |
コンパイル時間 | 360 ms |
コンパイル使用メモリ | 82,444 KB |
実行使用メモリ | 124,724 KB |
最終ジャッジ日時 | 2024-10-25 19:33:01 |
合計ジャッジ時間 | 13,361 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 28 |
ソースコード
import sysimport mathINF = 8000000000000000.0def main():input = sys.stdin.readdata = input().split()N = int(data[0])K = int(data[1])x = [0] * (N + 1)y = [0] * (N + 1)index = 2for i in range(N):x[i] = int(data[index])y[i] = int(data[index + 1])index += 2x[N] = 0y[N] = 0dist = [[0] * (N + 1) for _ in range(N + 1)]for i in range(N + 1):for j in range(N + 1):res = (x[i] - x[j]) ** 2 + (y[i] - y[j]) ** 2dist[i][j] = math.sqrt(res)dp = [[[INF] * (K + 1) for _ in range(N + 1)] for _ in range(1 << N)]dp[0][N][K] = 0.0for bit in range(1 << N):for from_idx in range(N + 1):if from_idx < N and not (bit >> from_idx) & 1:continuefor k in range(1, K + 1):for to in range(N):if (bit >> to) & 1:continuedp[bit | (1 << to)][to][k - 1] = min(dp[bit | (1 << to)][to][k - 1], dp[bit][from_idx][k] + dist[from_idx][to])for k in range(K + 1):dp[bit][N][K] = min(dp[bit][N][K], dp[bit][from_idx][k] + dist[from_idx][N])print("{:.15f}".format(dp[(1 << N) - 1][N][K]))if __name__ == "__main__":main()