結果

問題 No.2955 Pizza Delivery Plan
ユーザー twooimp2
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

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(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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0