oj.uz

Submission #547553

# Submission time Handle Problem Language Result Execution time Memory
547553 2022-04-11T00:36:18 Z dolphingarlic Grad (COI14_grad) C++14
100 / 100
 328 ms 122700 KB 
#include <bits/stdc++.h>
using namespace std;
 
struct Point { int x, y; } p[101010];
 
struct SP;
 
struct Node { // Node in tree decomposition
	int idx[2], depth;
	double dist;
 
	Node(int _u, int _v) : dist(hypot(p[_u].x - p[_v].x, p[_u].y - p[_v].y)) {
		idx[0] = _u;
		idx[1] = _v;
	}
 
	Node *complement;
	vector<Node *> anc;
	vector<SP> anc_sp;
} *has_city[101010];
 
struct SP { // Shortest path between two nodes in tree decomposition
	int idx[2][2];
	double dist[2][2];
 
	SP(Node *u) {
		idx[0][0] = idx[1][0] = u->idx[0];
		idx[0][1] = idx[1][1] = u->idx[1];
		dist[0][0] = dist[1][1] = 0;
		dist[0][1] = dist[1][0] = u->dist;
	}
 
	SP(Node *u, Node *v) {
		assert(u->anc[0] == v);
		idx[0][0] = u->idx[0], idx[0][1] = u->idx[1];
		idx[1][0] = v->idx[0], idx[1][1] = v->idx[1];
		if (u->idx[0] == v->idx[0]) {
			dist[0][0] = 0;
			dist[0][1] = v->dist;
			dist[1][0] = u->dist;
			dist[1][1] = u->complement->dist;
		} else {
			dist[0][0] = v->dist;
			dist[0][1] = 0;
			dist[1][0] = u->complement->dist;
			dist[1][1] = u->dist;
		}
	}
 
	SP(SP a, SP b) {
		assert(a.idx[1][0] == b.idx[0][0] && a.idx[1][1] == b.idx[0][1]);
		idx[0][0] = a.idx[0][0], idx[0][1] = a.idx[0][1];
		idx[1][0] = b.idx[1][0], idx[1][1] = b.idx[1][1];
		for (int i : {0, 1}) for (int j : {0, 1})
			dist[i][j] = min(a.dist[i][0] + b.dist[0][j], a.dist[i][1] + b.dist[1][j]);
	}
};
 
map<pair<int, int>, Node*> mp;
 
int main() {
	cin.tie(0)->sync_with_stdio(0);
	int n;
	cin >> p[1].x >> p[1].y >> p[2].x >> p[2].y >> n;
	has_city[1] = has_city[2] = mp[{1, 2}] = new Node(1, 2);
 
	for (int cnt = 2; n; n--) {
		char c;
		cin >> c;
		if (c == 'd') {
			cnt++;
			int a, b;
			cin >> p[cnt].x >> p[cnt].y >> a >> b;
			if (a > b) swap(a, b);
			Node *par = mp[{a, b}];
 
			Node *u = new Node(a, cnt), *v = new Node(b, cnt);
			u->complement = v, v->complement = u;
			u->depth = v->depth = par->depth + 1;
			has_city[cnt] = u;
 
			u->anc.push_back(par);
			u->anc_sp.push_back(SP(u, par));
			for (int i = 0; i < u->anc.size() && i < u->anc[i]->anc.size(); i++) {
				u->anc.push_back(u->anc[i]->anc[i]);
				u->anc_sp.push_back(SP(u->anc_sp[i], u->anc[i]->anc_sp[i]));
			}
			mp[{a, cnt}] = u;
 
			v->anc.push_back(par);
			v->anc_sp.push_back(SP(v, par));
			for (int i = 0; i < v->anc.size() && i < v->anc[i]->anc.size(); i++) {
				v->anc.push_back(v->anc[i]->anc[i]);
				v->anc_sp.push_back(SP(v->anc_sp[i], v->anc[i]->anc_sp[i]));
			}
			mp[{b, cnt}] = v;
 
			has_city[cnt] = u;
		} else {
			int a, b;
			cin >> a >> b;
			if (has_city[a]->depth < has_city[b]->depth) swap(a, b);
			double ans = 1e18;
			Node *anode = has_city[a], *bnode = has_city[b];
			SP asp = SP(anode), bsp = SP(bnode);
 
			for (int k = 0, diff = anode->depth - bnode->depth; diff; k++, diff >>= 1) {
				if (diff & 1) {
					asp = SP(asp, anode->anc_sp[k]);
					anode = anode->anc[k];
				}
			}
			for (int k = anode->anc.size() - 1; ~k; k--) {
				if (k < anode->anc.size() && anode->anc[k] != bnode->anc[k]) {
					asp = SP(asp, anode->anc_sp[k]);
					anode = anode->anc[k];
					bsp = SP(bsp, bnode->anc_sp[k]);
					bnode = bnode->anc[k];
				}
			}
 
			for (int k : {0, 1}) for (int l : {0, 1}) if (asp.idx[0][k] == a && bsp.idx[0][l] == b)
				for (int m : {0, 1}) for (int o : {0, 1}) if (asp.idx[1][m] == bsp.idx[1][o])
					ans = min(ans, asp.dist[k][m] + bsp.dist[l][o]);
			if (anode != bnode) {
				asp = SP(asp, anode->anc_sp[0]), bsp = SP(bsp, bnode->anc_sp[0]);
				for (int k : {0, 1}) for (int l : {0, 1}) if (asp.idx[0][k] == a && bsp.idx[0][l] == b)
					for (int m : {0, 1}) for (int o : {0, 1}) if (asp.idx[1][m] == bsp.idx[1][o])
						ans = min(ans, asp.dist[k][m] + bsp.dist[l][o]);
			}
			cout << fixed << setprecision(6) << ans << '\n';
		}
	}
	return 0;
}

Compilation message

grad.cpp: In function 'int main()':
grad.cpp:84:22: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Node*>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   84 |    for (int i = 0; i < u->anc.size() && i < u->anc[i]->anc.size(); i++) {
      |                    ~~^~~~~~~~~~~~~~~
grad.cpp:84:43: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Node*>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   84 |    for (int i = 0; i < u->anc.size() && i < u->anc[i]->anc.size(); i++) {
      |                                         ~~^~~~~~~~~~~~~~~~~~~~~~~
grad.cpp:92:22: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Node*>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   92 |    for (int i = 0; i < v->anc.size() && i < v->anc[i]->anc.size(); i++) {
      |                    ~~^~~~~~~~~~~~~~~
grad.cpp:92:43: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Node*>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
   92 |    for (int i = 0; i < v->anc.size() && i < v->anc[i]->anc.size(); i++) {
      |                                         ~~^~~~~~~~~~~~~~~~~~~~~~~
grad.cpp:114:11: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<Node*>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  114 |     if (k < anode->anc.size() && anode->anc[k] != bnode->anc[k]) {
      |         ~~^~~~~~~~~~~~~~~~~~~

Subtask #1
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 1 ms 340 KB 41 numbers

Subtask #2
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 468 KB 300 numbers

Subtask #3
7.0 / 7.0

# Verdict Execution time Memory Grader output
1 Correct 2 ms 1108 KB 500 numbers

Subtask #4
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 148 ms 76360 KB 15000 numbers
2 Correct 141 ms 44948 KB 15000 numbers
3 Correct 156 ms 65612 KB 30000 numbers

Subtask #5
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 276 ms 112944 KB 28333 numbers
2 Correct 206 ms 97736 KB 28333 numbers
3 Correct 249 ms 87520 KB 40000 numbers

Subtask #6
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 280 ms 109880 KB 50000 numbers
2 Correct 244 ms 57460 KB 50000 numbers
3 Correct 290 ms 109704 KB 50000 numbers

Subtask #7
6.0 / 6.0

# Verdict Execution time Memory Grader output
1 Correct 281 ms 99396 KB 55000 numbers
2 Correct 244 ms 86872 KB 55000 numbers
3 Correct 267 ms 109616 KB 50000 numbers

Subtask #8
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 295 ms 109800 KB 50000 numbers
2 Correct 274 ms 108888 KB 50000 numbers
3 Correct 254 ms 45476 KB 50000 numbers
4 Correct 302 ms 109712 KB 50000 numbers

Subtask #9
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 296 ms 122700 KB 44000 numbers
2 Correct 268 ms 119856 KB 44000 numbers
3 Correct 236 ms 69188 KB 44000 numbers
4 Correct 301 ms 122568 KB 44000 numbers

Subtask #10
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 301 ms 109956 KB 50000 numbers
2 Correct 276 ms 109680 KB 50000 numbers
3 Correct 251 ms 50556 KB 50000 numbers
4 Correct 295 ms 109732 KB 50000 numbers

Subtask #11
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 312 ms 119224 KB 45713 numbers
2 Correct 328 ms 100512 KB 54285 numbers
3 Correct 239 ms 46176 KB 58571 numbers
4 Correct 289 ms 94860 KB 57000 numbers

Subtask #12
11.0 / 11.0

# Verdict Execution time Memory Grader output
1 Correct 293 ms 111084 KB 49285 numbers
2 Correct 290 ms 111072 KB 49285 numbers
3 Correct 248 ms 50156 KB 49285 numbers
4 Correct 296 ms 120108 KB 45000 numbers