結果

問題 No.2955 Pizza Delivery Plan
コンテスト
ユーザー twooimp2
提出日時 2024-12-29 19:57:38
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 340 ms / 2,000 ms
コード長 1,334 bytes
コンパイル時間 5,144 ms
コンパイル使用メモリ 274,444 KB
実行使用メモリ 67,328 KB
最終ジャッジ日時 2024-12-29 19:57:49
合計ジャッジ時間 10,869 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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(20);
}

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&(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-1][i]=min(dp[b|(1<<i)][K-1][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;
}
0