結果
問題 | No.2955 Pizza Delivery Plan |
ユーザー |
|
提出日時 | 2024-12-29 19:53:04 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,336 bytes |
コンパイル時間 | 5,516 ms |
コンパイル使用メモリ | 273,896 KB |
実行使用メモリ | 67,456 KB |
最終ジャッジ日時 | 2024-12-29 19:53:16 |
合計ジャッジ時間 | 11,308 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 11 WA * 17 |
ソースコード
#pragma GCC target("avx2")#pragma GCC optimize("O3")#pragma GCC optimize("unroll-loops")#include<bits/stdc++.h>using namespace std;using ll=long long;using ld=long double;void IO(){ios::sync_with_stdio(false);std::cin.tie(nullptr);cout<<fixed<<setprecision(10);}ll n,K;vector<ld> x,y;ld d1(ll i,ll j){ld n2=(x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j]);return sqrt(n2);}ld d2(ll i){ld n2=x[i]*x[i]+y[i]*y[i];return sqrt(n2);}int main(){IO();cin>>n>>K;x.resize(n);y.resize(n);for(ll i=0;i<n;i++){cin>>x[i]>>y[i];}vector<vector<vector<ld>>> dp(1<<n,vector<vector<ld>>(K+1,vector<ld>(n,1e18)));for(ll i=0;i<n;i++){dp[1<<i][K-1][i]=d2(i);}for(ll b=1;b<(1<<n);b++){for(ll i=0;i<n;i++){for(ll j=0;j<n;j++){if(b!=0&&!(b&(1<<j))){continue;}if((b&(1<<i))==0){if(i!=j){for(ll k=0;k<=K;k++){if(k>0){dp[b|(1<<i)][k-1][i]=min(dp[b|(1<<i)][k-1][i],dp[b][k][j]+d1(j,i));}dp[b|(1<<i)][K][i]=min(dp[b|(1<<i)][K][i],dp[b][k][j]+d2(j)+d2(i));}}}}}}ld ans=1e18;for(ll i=0;i<n;i++){for(ll j=0;j<=K;j++){ans=min(ans,dp[(1<<n)-1][j][i]+d2(i));}}cout<<ans<<endl;}