답안 #60965

# 제출 시각 아이디 문제 언어 결과 실행 시간 메모리
60965 2018-07-25T04:32:38 Z tugushka Simurgh (IOI17_simurgh) C++14
30 / 100
3000 ms 20432 KB
#include "simurgh.h"
#include<bits/stdc++.h>
#define mp make_pair
#define pb push_back
using namespace std;

typedef pair < int, int > pii;

map < pii, int > id, gold;

vector < set < int > > v;
vector < int > cyc, exc;
int nn, mm, p[505], d[505], cyc_len;

void find_cycle(){
	cyc.clear();
	memset( p, -1, sizeof(p) );
	memset( d, -1, sizeof(d) );
	queue < int > q;
	q.push(0);
	d[0] = 0;
	int dep = 0, t1 = -1, t2, c;
	while( !q.empty() ){
		int sz = q.size();
		dep++;
		while(sz--){
			int u = q.front();
//			cout << u << ' ' << d[u] << ' ' << p[u] << endl;
			q.pop();
			for( set < int > :: iterator it = v[u].begin() ; it != v[u].end() ; it++ ){
				if( p[u] == *it ) continue;
//				cout << u << " --> " << *it << endl;
				if( d[*it] == -1 ){
					p[*it] = u;
					d[*it] = dep;
					q.push(*it);
				}
				else{
					t1 = *it;
					t2 = u;
					
//					cout << "Tails : " << t1 << ' ' << t2 << endl;
					
					break;
				}
			}
			if( t1 != -1 ) break;
		}
		if( t1 != -1 ) break;
	}
	
	vector < int > cw, ccw;
	while( t1 != t2 ){
		if( d[t1] < d[t2] ) ccw.pb(t2), t2 = p[t2];
		else if( d[t1] > d[t2] ) cw.pb(t1), t1 = p[t1];
		else ccw.pb(t2), cw.pb(t1), t2 = p[t2], t1 = p[t1];
	}
	
	cyc.pb(t1);
	for(int i = ccw.size()-1 ; i >= 0 ; i--) cyc.pb(ccw[i]);
	for(int i = 0 ; i < cw.size() ; i++) cyc.pb(cw[i]);
	cyc_len = cyc.size();
	cyc.pb(t1);
	
	/*
	cout << "Found\n";
	for(int i = 0 ; i < cyc.size() ; i++) cout << cyc[i] << ' ';
	cout << endl;
	*/
}

void eval(){
	exc.clear();
	memset(p, -1, sizeof(p));
	queue < int > q;
	for(int i = 0 ; i < cyc.size() ; i++) q.push(cyc[i]), p[cyc[i]] = cyc[i];
	while( !q.empty() ){
		int u = q.front();
		q.pop();
		for( set < int > :: iterator it = v[u].begin() ; it != v[u].end() ; it++ ){
			if( p[*it] == -1 ){
				p[*it] = u;
//				cout << "Excess : " << u << ' ' << *it << endl;
				exc.pb( id[ mp(u, *it) ] );
				q.push(*it);
			}
		}
	}
}

vector<int> find_roads(int _n, vector<int> _u, vector<int> _v) {
	nn = _n;
	v.resize(nn);
	mm = _u.size();
	for(int i = 0 ; i < mm ; i++){
		v[_u[i]].insert(_v[i]);
		v[_v[i]].insert(_u[i]);
		id[ mp(_u[i], _v[i]) ] = i;
		id[ mp(_v[i], _u[i]) ] = i;
	}
	
	while( mm >= nn ){
		find_cycle();
		eval();
		
		vector < int > now;
		for(int i = 0 ; i < cyc_len ; i++){
			vector < int > query;
			for(int j = 0 ; j < exc.size() ; j++) query.pb(exc[j]);
			for(int j = 0 ; j < cyc_len ; j++){
				if( i == j ) continue;
				query.pb( id[mp( cyc[j], cyc[j+1] )] );
			}
			
			now.pb( count_common_roads(query) );
			
//			cout << "Query : "; for(int i = 0 ; i < query.size() ; i++) cout << query[i];
//			cout << endl << now.back() << endl;;
		}
		
		int f = -1;
		for(int i = 1 ; i < cyc_len ; i++){
			if( now[i] != now[i-1] ){
				if( now[i] < now[i-1] ) f = i;
				else f = i-1;
				break;
			}
		}
		
//		cout << "! " << f << endl;
		
		if( f == -1 ){
			for(int i = 0 ; i < cyc_len ; i++){
				v[ cyc[i] ].erase( cyc[i+1] );
				v[ cyc[i+1] ].erase( cyc[i] );
				mm--;
			}
			continue;
		}
		
		bool g = 1;
		for(int i = f-1 ; i >= 0 ; i--){
			if( now[i] != now[i+1] ) g = !g;
			if( !g ){
				v[ cyc[i] ].erase( cyc[i+1] );
				v[ cyc[i+1] ].erase( cyc[i] );
				mm--;
			}
			else{
				gold[ mp( cyc[i], cyc[i+1] ) ] = 1;
				gold[ mp( cyc[i+1], cyc[i] ) ] = 1;
			}
		}
		g = 1;
		for(int i = f+1 ; i < cyc_len ; i++){
			if( now[i] != now[i-1] ) g = !g;
			if( !g ){
				v[ cyc[i] ].erase( cyc[i+1] );
				v[ cyc[i+1] ].erase( cyc[i] );
				mm--;
			}
			else{
				gold[ mp( cyc[i], cyc[i+1] ) ] = 1;
				gold[ mp( cyc[i+1], cyc[i] ) ] = 1;
			}
		}		
	}
	
	
	vector < int > res;
	for(int i = 0 ; i < nn ; i++){
		for( set < int > :: iterator it = v[i].begin() ; it != v[i].end() ; it++ ){
			if( *it > i ) res.pb( id[ mp(*it, i) ] );
		}
	}
	
//	for(int i = 0 ; i < res.size() ; i++) cout << res[i] << ' ';
	
	return res;
}
/*
6 6
0 1
1 2
2 3
3 4
0 4
0 5
0 1 2 3 5

4 6
0 1
0 2
0 3
1 2
1 3
2 3
0 1 5

*/

Compilation message

simurgh.cpp: In function 'void find_cycle()':
simurgh.cpp:61:20: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for(int i = 0 ; i < cw.size() ; i++) cyc.pb(cw[i]);
                  ~~^~~~~~~~~~~
simurgh.cpp:22:28: warning: unused variable 'c' [-Wunused-variable]
  int dep = 0, t1 = -1, t2, c;
                            ^
simurgh.cpp: In function 'void eval()':
simurgh.cpp:76:20: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
  for(int i = 0 ; i < cyc.size() ; i++) q.push(cyc[i]), p[cyc[i]] = cyc[i];
                  ~~^~~~~~~~~~~~
simurgh.cpp: In function 'std::vector<int> find_roads(int, std::vector<int>, std::vector<int>)':
simurgh.cpp:109:22: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
    for(int j = 0 ; j < exc.size() ; j++) query.pb(exc[j]);
                    ~~^~~~~~~~~~~~
# 결과 실행 시간 메모리 Grader output
1 Correct 3 ms 376 KB correct
2 Correct 3 ms 512 KB correct
3 Correct 2 ms 512 KB correct
4 Correct 2 ms 512 KB correct
5 Correct 3 ms 512 KB correct
6 Correct 3 ms 516 KB correct
7 Correct 4 ms 520 KB correct
8 Correct 3 ms 520 KB correct
9 Correct 3 ms 528 KB correct
10 Correct 3 ms 640 KB correct
11 Correct 2 ms 688 KB correct
12 Correct 3 ms 688 KB correct
13 Correct 3 ms 944 KB correct
# 결과 실행 시간 메모리 Grader output
1 Correct 3 ms 376 KB correct
2 Correct 3 ms 512 KB correct
3 Correct 2 ms 512 KB correct
4 Correct 2 ms 512 KB correct
5 Correct 3 ms 512 KB correct
6 Correct 3 ms 516 KB correct
7 Correct 4 ms 520 KB correct
8 Correct 3 ms 520 KB correct
9 Correct 3 ms 528 KB correct
10 Correct 3 ms 640 KB correct
11 Correct 2 ms 688 KB correct
12 Correct 3 ms 688 KB correct
13 Correct 3 ms 944 KB correct
14 Correct 27 ms 956 KB correct
15 Correct 34 ms 964 KB correct
16 Correct 34 ms 972 KB correct
17 Correct 21 ms 972 KB correct
18 Correct 9 ms 972 KB correct
19 Correct 27 ms 992 KB correct
20 Correct 26 ms 992 KB correct
21 Correct 20 ms 992 KB correct
22 Correct 14 ms 992 KB correct
23 Correct 16 ms 992 KB correct
24 Correct 9 ms 992 KB correct
25 Correct 3 ms 992 KB correct
26 Correct 12 ms 992 KB correct
27 Correct 11 ms 992 KB correct
28 Correct 7 ms 992 KB correct
29 Correct 4 ms 992 KB correct
30 Correct 12 ms 992 KB correct
31 Correct 15 ms 992 KB correct
32 Correct 11 ms 992 KB correct
33 Correct 13 ms 992 KB correct
# 결과 실행 시간 메모리 Grader output
1 Correct 3 ms 376 KB correct
2 Correct 3 ms 512 KB correct
3 Correct 2 ms 512 KB correct
4 Correct 2 ms 512 KB correct
5 Correct 3 ms 512 KB correct
6 Correct 3 ms 516 KB correct
7 Correct 4 ms 520 KB correct
8 Correct 3 ms 520 KB correct
9 Correct 3 ms 528 KB correct
10 Correct 3 ms 640 KB correct
11 Correct 2 ms 688 KB correct
12 Correct 3 ms 688 KB correct
13 Correct 3 ms 944 KB correct
14 Correct 27 ms 956 KB correct
15 Correct 34 ms 964 KB correct
16 Correct 34 ms 972 KB correct
17 Correct 21 ms 972 KB correct
18 Correct 9 ms 972 KB correct
19 Correct 27 ms 992 KB correct
20 Correct 26 ms 992 KB correct
21 Correct 20 ms 992 KB correct
22 Correct 14 ms 992 KB correct
23 Correct 16 ms 992 KB correct
24 Correct 9 ms 992 KB correct
25 Correct 3 ms 992 KB correct
26 Correct 12 ms 992 KB correct
27 Correct 11 ms 992 KB correct
28 Correct 7 ms 992 KB correct
29 Correct 4 ms 992 KB correct
30 Correct 12 ms 992 KB correct
31 Correct 15 ms 992 KB correct
32 Correct 11 ms 992 KB correct
33 Correct 13 ms 992 KB correct
34 Execution timed out 3045 ms 7716 KB Time limit exceeded
35 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 2 ms 7716 KB correct
2 Correct 3 ms 7716 KB correct
3 Execution timed out 3050 ms 20432 KB Time limit exceeded
4 Halted 0 ms 0 KB -
# 결과 실행 시간 메모리 Grader output
1 Correct 3 ms 376 KB correct
2 Correct 3 ms 512 KB correct
3 Correct 2 ms 512 KB correct
4 Correct 2 ms 512 KB correct
5 Correct 3 ms 512 KB correct
6 Correct 3 ms 516 KB correct
7 Correct 4 ms 520 KB correct
8 Correct 3 ms 520 KB correct
9 Correct 3 ms 528 KB correct
10 Correct 3 ms 640 KB correct
11 Correct 2 ms 688 KB correct
12 Correct 3 ms 688 KB correct
13 Correct 3 ms 944 KB correct
14 Correct 27 ms 956 KB correct
15 Correct 34 ms 964 KB correct
16 Correct 34 ms 972 KB correct
17 Correct 21 ms 972 KB correct
18 Correct 9 ms 972 KB correct
19 Correct 27 ms 992 KB correct
20 Correct 26 ms 992 KB correct
21 Correct 20 ms 992 KB correct
22 Correct 14 ms 992 KB correct
23 Correct 16 ms 992 KB correct
24 Correct 9 ms 992 KB correct
25 Correct 3 ms 992 KB correct
26 Correct 12 ms 992 KB correct
27 Correct 11 ms 992 KB correct
28 Correct 7 ms 992 KB correct
29 Correct 4 ms 992 KB correct
30 Correct 12 ms 992 KB correct
31 Correct 15 ms 992 KB correct
32 Correct 11 ms 992 KB correct
33 Correct 13 ms 992 KB correct
34 Execution timed out 3045 ms 7716 KB Time limit exceeded
35 Halted 0 ms 0 KB -