結果
| 問題 |
No.2955 Pizza Delivery Plan
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2024-11-08 23:14:21 |
| 言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 276 ms / 2,000 ms |
| コード長 | 1,024 bytes |
| コンパイル時間 | 5,778 ms |
| コンパイル使用メモリ | 314,240 KB |
| 実行使用メモリ | 35,328 KB |
| 最終ジャッジ日時 | 2024-11-08 23:14:32 |
| 合計ジャッジ時間 | 10,564 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 28 |
ソースコード
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using namespace atcoder;
using ll=long long;
int main()
{
int n,k;
cin>>n>>k;
vector<double> x(n+1),y(n+1);
for(int i=0;i<n;i++){
cin>>x[i+1]>>y[i+1];
}
auto dist=[&](int i,int j)->double{
return sqrtf64((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]));
};
vector<vector<vector<double>>> dp(n+1,vector<vector<double>>(k+1,vector<double>(1<<n,-1)));
auto rec=[&](auto rec,int v,int m,int s)->double{
if(dp[v][m][s]!=-1)return dp[v][m][s];
else if((s+1)==1<<n)return dp[v][m][s]=dist(v,0);
double res=1e40;
if(v!=0)res=rec(rec,0,k,s)+dist(v,0);
if(m>0){
for(int i=1;i<=n;i++){
int u=i-1;
if(s>>u&1)continue;
res=min(res,rec(rec,i,m-1,s|(1<<u))+dist(v,i));
}
}
return dp[v][m][s]=res;
};
cout<<fixed<<setprecision(10);
cout<<rec(rec,0,k,0)<<endl;
return 0;
}