C++の練習を兼ねて, AtCoder Beginner Contest 181 の 問題F (Silver Woods) を解いてみた.
■感想.
1. 問題Fは, 方針が見えなかったので, 解説を参照して実装して, ようやく, AC版となった.
2. Union-Find木 の 復習が出来たので, 非常に良かったと思う.
※ 公式のライブラリを拝借させて頂いてます.
3. 個人的には, 非常に面白い問題に感じた.
4. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイト AtCoder Beginner Contest 181 解説 の 各リンク を ご覧下さい.
■C++版プログラム(問題F／AC版).
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 |
// 解き直し. // https://atcoder.jp/contests/abc181/editorial/267 // C++(GCC 9.2.1) #include <bits/stdc++.h> using namespace std; using P = pair<double, pair<int, int>>; #define repex(i, a, b, c) for(int i = a; i < b; i += c) #define repx(i, a, b) repex(i, a, b, 1) #define rep(i, n) repx(i, 0, n) #define repr(i, a, b) for(int i = a; i >= b; i--) #define a first #define b second double x[111], y[111]; // https://github.com/atcoder/live_library/blob/master/uf.cpp // UnionFind // coding: https://youtu.be/TdR816rqc3s?t=726 // comment: https://youtu.be/TdR816rqc3s?t=6822 // -> 一部改変. struct UnionFind{ vector<int> d; UnionFind(int n = 0): d(n, -1) {} int find(int x){ return (d[x] < 0) ? x : (d[x] = find(d[x])); } bool unite(int x, int y){ x = find(x); y = find(y); if(x == y) return false; if(d[x] > d[y]) swap(x, y); d[x] += d[y]; d[y] = x; return true; } bool same(int x, int y){ return find(x) == find(y); } int size(int x){ return -d[find(x)]; } }; int main(){ // 1. 入力情報. int N; scanf("%d", &N); rep(i, N) scanf("%lf %lf", &x[i], &y[i]); // 2. 2点間の距離, 頂点のペア を 保存. // -> 頂点 1 … 直線 y = -100, 頂点 N + 2 … 直線 y = 100 とする. priority_queue<P, vector<P>, greater<P>> pq; rep(i, N){ P p0, p1; // 頂点 1 と 他の頂点との辺 を 追加. p0.a = abs(-100.0 - y[i]); p0.b.a = 0; p0.b.b = i + 1; pq.push(p0); // 頂点 N + 2 と 他の頂点との辺 を 追加. p1.a = abs(100.0 - y[i]); p1.b.a = N + 1; p1.b.b = i + 1; pq.push(p1); } rep(i, N){ repx(j, i + 1, N){ double d = hypot(x[i] - x[j], y[i] - y[j]); P p; p.a = d; p.b.a = i + 1; p.b.b = j + 1; pq.push(p); } } // 3. 距離の短い方から順に, 辺を追加. double ans = 0.0; UnionFind uf(N + 2); while(!pq.empty()){ // 3-1. 辺を取り出す. P p = pq.top(); pq.pop(); // 3-2. 円の直径を更新. ans = p.a; // 3-3. グラフに辺を追加. uf.unite(p.b.a, p.b.b); // 3-4. 辺を追加した結果, 二直線が, 同じ連結成分に含まれてしまう場合は終了. if(uf.same(0, N + 1)) break; } // 4. 出力. printf("%.12lf\n", ans / 2.0); return 0; } |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 |
[入力例] 2 0 -40 0 40 [出力例] 40 ※AtCoderテストケースより ※但し, 上記のプログラムでは, 以下の内容が出力される. 40.000000000000 [入力例] 4 0 -10 99 10 0 91 99 -91 [出力例] 50.5 ※AtCoderテストケースより ※但し, 上記のプログラムでは, 以下の内容が出力される. 50.500000000000 [入力例] 10 -90 40 20 -30 0 -90 10 -70 80 70 -90 30 -20 -80 10 90 50 30 60 -70 [出力例] 33.541019662496845446 ※AtCoderテストケースより ※但し, 上記のプログラムでは, 以下の内容が出力される. 33.541019662497 [入力例] 10 65 -90 -34 -2 62 99 42 -13 47 -84 84 87 16 -78 56 35 90 8 90 19 [出力例] 35.003571246374276203 ※AtCoderテストケースより ※但し, 上記のプログラムでは, 以下の内容が出力される. 35.003571246374 [入力例] 30 -85 42 12 69 17 -3 22 -51 -86 21 74 -58 78 -41 43 64 21 3 30 -58 4 32 14 22 28 42 71 -69 78 57 60 -4 87 17 -75 13 -86 32 28 95 -34 -32 70 -67 27 94 38 -80 21 -4 -54 87 89 81 -12 -13 87 -55 -42 81 [出力例] 20.573040611441 [入力例] 50 7 72 65 -52 46 33 59 30 30 -53 31 42 5 52 -7 56 40 41 34 99 50 -35 98 44 52 -69 7 97 28 -85 -96 41 46 -31 8 95 19 37 9 32 17 32 61 -94 10 33 89 19 38 -64 31 88 -7 25 86 61 3 56 10 79 -55 32 44 27 -55 74 18 44 50 -75 19 58 17 43 -36 60 68 43 24 79 45 41 19 74 46 34 30 48 -48 27 79 28 20 34 59 -44 72 68 23 72 [出力例] 29.017236257094 |