#1025096 (C++23) No.2955 Pizza Delivery Plan

• 問題: No.2955 Pizza Delivery Plan
• ユーザー: nouka28
• 提出日時: 2024-11-08 21:59:32
• 言語: C++23 (gcc 12.3.0 + boost 1.83.0)
• 結果: TLE
• 実行時間: -
• コード長: 1,246 bytes
• コンパイル時間: 3,089 ms
• コンパイル使用メモリ: 257,804 KB
• 実行使用メモリ: 10,624 KB
• 最終ジャッジ日時: 2024-11-08 21:59:40
• 合計ジャッジ時間: 6,785 ms
• ジャッジサーバーID （参考情報）: judge5 / judge2

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 10,624 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	2 ms; 5,248 KB
testcase_03	AC	2 ms; 5,248 KB
testcase_04	AC	2 ms; 5,248 KB
testcase_05	AC	2 ms; 5,248 KB
testcase_06	AC	2 ms; 5,248 KB
testcase_07	AC	2 ms; 5,248 KB
testcase_08	TLE	-
testcase_09	--	-
testcase_10	--	-
testcase_11	--	-
testcase_12	--	-
testcase_13	--	-
testcase_14	--	-
testcase_15	--	-
testcase_16	--	-
testcase_17	--	-
testcase_18	--	-
testcase_19	--	-
testcase_20	--	-
testcase_21	--	-
testcase_22	--	-
testcase_23	--	-
testcase_24	--	-
testcase_25	--	-
testcase_26	--	-
testcase_27	--	-
testcase_28	--	-
testcase_29	--	-
testcase_30	--	-

ソースコード

#include<bits/stdc++.h>
using namespace std;

#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")

#define int long long

using ld=long double;

signed main(){
	int N,K;cin>>N>>K;
	vector<int> X(N+1),Y(N+1);
	for(int i=0;i<N;i++){
		cin>>X[i]>>Y[i];
	}

	vector<vector<ld>> d(N+1,vector<ld>(N+1));

	for(int i=0;i<=N;i++){
		for(int j=0;j<=N;j++){
			d[i][j]=hypotl(X[i]-X[j],Y[i]-Y[j]);
		}
	}

	X[N]=0,Y[N]=0;

	vector<ld> cost(1<<N);

	// cout<<"!"<<endl;

	for(int S=1;S<(1<<N);S++){
		vector<int> v={N};
		for(int i=0;i<N;i++){
			if(S>>i&1)v.push_back(i);
		}

		int n=v.size();

		vector<vector<ld>> dp(1<<n,vector<ld>(n,1e18));

		dp[0][0]=0;

		for(int s=0;s<(1<<n);s++){
			for(int j=0;j<n;j++){
				for(int k=0;k<n;k++){
					if(s>>k&1)continue;

					dp[s|(1<<k)][k]=min(dp[s|(1<<k)][k],dp[s][j]+d[v[j]][v[k]]);					
				}
			}
		}

		cost[S]=dp[(1<<n)-1][0];

		// cout<<S<<" : "<<dp[(1<<n)-1][0]<<endl;
	}

	
	// cout<<"!"<<endl;

	vector<ld> dp(1<<N,1e18);

	dp[0]=0;

	for(int S=0;S<(1<<N);S++){
		int al=(1<<N)-1;

		int st=al^S;
		int sub=st;

		while(sub){
			if(__popcount(sub)<=K){
				dp[S|sub]=min(dp[S|sub],dp[S]+cost[sub]);
			}
			sub=(sub-1)&st;
		}
	}

	cout<<fixed<<setprecision(20)<<dp[(1<<N)-1]<<endl;
}