結果

問題 No.2955 Pizza Delivery Plan
コンテスト
ユーザー 沙耶花
提出日時 2024-11-08 21:37:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 165 ms / 2,000 ms
コード長 938 bytes
コンパイル時間 4,384 ms
コンパイル使用メモリ 258,752 KB
最終ジャッジ日時 2025-02-25 02:37:38
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 28
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <atcoder/all>
#include <bits/stdc++.h>
using namespace std;
using namespace atcoder;
using mint = modint998244353;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf32 1000000001
#define Inf64 1000000000000000000LL


int main(){
	
	int N,K;
	cin>>N>>K;
	vector<double> x(N),y(N);
	rep(i,N) cin>>x[i]>>y[i];
	x.insert(x.begin(),0);
	y.insert(y.begin(),0);
	N++;
	vector dp(1<<N,vector(N,vector<double>(K+1,1e100)));
	dp[1][0][0] = 0.0;
	rep(i,1<<N){
		for(int k=K;k>=0;k--){
		rep(j,N){
			
				if(dp[i][j][k]==1e100)continue;
				rep(l,N){
					if(k==K&&l!=0)continue;
					if(l!=0&&((i>>l)&1))continue;
					double d = sqrt((x[j]-x[l])*(x[j]-x[l])+(y[j]-y[l])*(y[j]-y[l]));
					int ni = i|(1<<l);
					int nj = l;
					int nk = k+1;
					if(l==0)nk = 0;
					dp[ni][nj][nk] = min(dp[ni][nj][nk],dp[i][j][k]+d);
				}
			}
		}
	}
	cout<<fixed<<setprecision(20)<<dp[(1<<N)-1][0][0]<<endl;
	return 0;
}
0