#815588 (C++17) No.1069 電柱 / Pole (Hard)

提出ソース

結果

問題	No.1069 電柱 / Pole (Hard)
ユーザー	kmjp
提出日時	2022-11-10 23:58:30
言語	C++17 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	698 ms / 2,000 ms
コード長	2,860 bytes
コンパイル時間	3,058 ms
コンパイル使用メモリ	234,028 KB
実行使用メモリ	38,080 KB
最終ジャッジ日時	2024-07-26 19:15:25
合計ジャッジ時間	7,853 ms
ジャッジサーバーID （参考情報）	judge1 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	4 ms 10,052 KB
testcase_01	AC	4 ms 9,924 KB
testcase_02	AC	4 ms 9,924 KB
testcase_03	AC	4 ms 9,800 KB
testcase_04	AC	441 ms 30,372 KB
testcase_05	AC	698 ms 38,080 KB
testcase_06	AC	28 ms 27,096 KB
testcase_07	AC	18 ms 26,524 KB
testcase_08	AC	18 ms 18,564 KB
testcase_09	AC	35 ms 26,276 KB
testcase_10	AC	41 ms 28,552 KB
testcase_11	AC	6 ms 12,752 KB
testcase_12	AC	6 ms 11,896 KB
testcase_13	AC	8 ms 14,384 KB
testcase_14	AC	8 ms 12,512 KB
testcase_15	AC	9 ms 14,372 KB
testcase_16	AC	17 ms 20,960 KB
testcase_17	AC	7 ms 12,044 KB
testcase_18	AC	15 ms 18,544 KB
testcase_19	AC	24 ms 20,976 KB
testcase_20	AC	6 ms 10,304 KB
testcase_21	AC	12 ms 12,356 KB
testcase_22	AC	5 ms 12,152 KB
testcase_23	AC	27 ms 10,308 KB
testcase_24	AC	55 ms 33,124 KB
testcase_25	AC	54 ms 33,152 KB
testcase_26	AC	26 ms 33,672 KB
testcase_27	AC	27 ms 33,112 KB
testcase_28	AC	26 ms 31,580 KB
testcase_29	AC	15 ms 18,860 KB
testcase_30	AC	5 ms 10,560 KB
testcase_31	AC	4 ms 9,800 KB
testcase_32	AC	4 ms 10,048 KB
testcase_33	AC	4 ms 9,924 KB
testcase_34	AC	4 ms 10,052 KB
testcase_35	AC	4 ms 10,056 KB
testcase_36	AC	4 ms 10,056 KB
testcase_37	AC	5 ms 12,184 KB
testcase_38	AC	6 ms 12,196 KB
testcase_39	AC	4 ms 10,052 KB
testcase_40	AC	5 ms 12,048 KB
testcase_41	AC	6 ms 12,164 KB
testcase_42	AC	6 ms 12,312 KB
testcase_43	AC	5 ms 12,276 KB
testcase_44	AC	20 ms 20,436 KB
testcase_45	AC	11 ms 16,548 KB
testcase_46	AC	15 ms 18,392 KB
testcase_47	AC	23 ms 23,076 KB
testcase_48	AC	41 ms 29,444 KB
testcase_49	AC	24 ms 22,752 KB
testcase_50	AC	6 ms 12,192 KB
testcase_51	AC	32 ms 26,580 KB
testcase_52	AC	7 ms 14,196 KB
testcase_53	AC	38 ms 28,512 KB
testcase_54	AC	4 ms 11,972 KB
testcase_55	AC	4 ms 12,104 KB
testcase_56	AC	3 ms 10,052 KB
testcase_57	AC	4 ms 10,056 KB
testcase_58	AC	4 ms 12,100 KB
testcase_59	AC	5 ms 11,972 KB
testcase_60	AC	5 ms 10,048 KB
testcase_61	AC	4 ms 9,928 KB
testcase_62	AC	4 ms 11,848 KB
testcase_63	AC	4 ms 10,052 KB
testcase_64	AC	5 ms 12,228 KB
testcase_65	AC	4 ms 9,928 KB
testcase_66	AC	4 ms 9,800 KB
testcase_67	AC	4 ms 10,056 KB
testcase_68	AC	4 ms 9,928 KB
testcase_69	AC	60 ms 35,096 KB
testcase_70	AC	70 ms 38,016 KB
testcase_71	AC	68 ms 37,352 KB
testcase_72	AC	21 ms 25,424 KB
testcase_73	AC	30 ms 26,812 KB
testcase_74	AC	57 ms 37,712 KB
testcase_75	AC	31 ms 37,960 KB
testcase_76	AC	31 ms 31,892 KB
testcase_77	AC	28 ms 29,320 KB
testcase_78	AC	5 ms 10,304 KB
testcase_79	AC	14 ms 14,592 KB
testcase_80	AC	18 ms 22,668 KB
testcase_81	AC	12 ms 14,460 KB
testcase_82	AC	9 ms 14,404 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <bits/stdc++.h>
using namespace std;
typedef signed long long ll;

#undef _P
#define _P(...) (void)printf(__VA_ARGS__)
#define FOR(x,to) for(x=0;x<(to);x++)
#define FORR(x,arr) for(auto& x:arr)
#define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++)
#define ALL(a) (a.begin()),(a.end())
#define ZERO(a) memset(a,0,sizeof(a))
#define MINUS(a) memset(a,0xff,sizeof(a))
//-------------------------------------------------------

int N,M,K;
int S,T;
int X[202020],Y[202020];
bitset<2020> B[2020];
vector<int> E[202020];
double D[2020];
vector<int> P[2020];
double dist[2020][2020];

pair<double,vector<int>> hoge(int start,set<pair<int,int>> NG,vector<int> pref,double prel) {
	int i;
	
	FOR(i,N) D[i]=1e9, P[i].clear(), B[i].reset();
	FORR(p,pref) D[p]=-1;
	D[start]=prel;
	
	P[start]=pref;
	P[start].push_back(start);
	FORR(p,pref) B[start][p]=1;
	priority_queue<pair<double,int>> Q;
	Q.push({-prel,start});
	while(Q.size()) {
		double co=-Q.top().first;
		int cur=Q.top().second;
		Q.pop();
		if(D[cur]!=co) continue;
		FORR(e,E[cur]) {
			if(NG.count({e,cur})) continue;
			if(NG.count({cur,e})) continue;
			if(B[cur][e]) continue;
			double d=co+dist[cur][e];
			if(D[e]>d) {
				D[e]=d;
				P[e]=P[cur];
				P[e].push_back(e);
				B[e]=B[cur];
				B[e][e]=1;
				Q.push({-d,e});
			}
		}
	}
	
	return {D[T],P[T]};
	
}


void solve() {
	int i,j,k,l,r,x,y; string s;
	
	cin>>N>>M>>K>>S>>T;
	S--,T--;
	FOR(i,N) cin>>X[i]>>Y[i];
	FOR(x,N) FOR(y,N) dist[x][y]=hypot(X[x]-X[y],Y[x]-Y[y]);
	FOR(i,M) {
		cin>>x>>y;
		E[x-1].push_back(y-1);
		E[y-1].push_back(x-1);
	}
	
	vector<pair<double,vector<int>>> ret;
	vector<pair<double,vector<int>>> cand;
	set<pair<int,int>> NG;
	ret.push_back(hoge(S,NG,vector<int>(),0));
	while(ret.size()<K) {
		double sum=0;
		vector<int> pre;
		FOR(x,ret.back().second.size()-1) {
			if(x) {
				pre.push_back(ret.back().second[x-1]);
				sum+=dist[ret.back().second[x-1]][ret.back().second[x]];
			}
			NG.clear();
			FOR(i,ret.size()) {
				if(ret[i].second.size()<x+2) continue;
				FOR(y,x+1) if(ret[i].second[y]!=ret.back().second[y]) break;
				if(y!=x+1) continue;
				NG.insert({ret[i].second[x],ret[i].second[x+1]});
			}
			
			cand.push_back(hoge(ret.back().second[x],NG,pre,sum));
		}
		
		
		
		sort(ALL(cand));
		while(cand.size()) {
			int ng=0;
			FORR(r,ret) if(r.second==cand[0].second) ng=1;
			if(ng==0) break;
			cand.erase(cand.begin());
		}
		if(cand.empty()) break;
		if(cand[0].first>1e8) break;
		ret.push_back(cand[0]);
		cand.erase(cand.begin());
	}
	
	FOR(i,K) {
		if(i<ret.size()) {
			_P("%.12lf\n",ret[i].first);
		}
		else {
			_P("-1\n");
		}
	}
	
	
}


int main(int argc,char** argv){
	string s;int i;
	if(argc==1) ios::sync_with_stdio(false), cin.tie(0);
	FOR(i,argc-1) s+=argv[i+1],s+='\n'; FOR(i,s.size()) ungetc(s[s.size()-1-i],stdin);
	cout.tie(0); solve(); return 0;
}

yukicoder

結果

テストケース

ソースコード