C++の練習を兼ねて, AtCoder Regular Contest 042 の 問題A (掲示板) ~ 問題B (アリの高橋くん) を解いてみた.
■感想.
1. 問題B は, 図形の問題に見え, ややこしい感じがしたが, 何とかAC版となった.
2. 時間を見つけて, 引き続き, 過去問を振り返っていきたいと思う.
本家のサイトARC 042 解説をご覧下さい.
■C++版プログラム(問題A/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 |
#include <bits/stdc++.h> using namespace std; using P = 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 pb push_back #define a first #define b second vector<vector<int>> t(101010); int main(){ // 1. 入力情報. int N, M, a; scanf("%d %d", &N, &M); rep(i, M){ scanf("%d", &a); t[a].pb(i + 1); } // 2. 情報を移し替える. priority_queue<P> pq; // スレッド書き込み有り. vector<int> v; // スレッド書き込み無し. repx(i, 1, N + 1){ if(t[i].size()) pq.push({t[i].back(), i}); else v.pb(i); } // 3. 出力. while(!pq.empty()){ P p = pq.top(); pq.pop(); printf("%d\n", p.b); } for(auto &p : v) printf("%d\n", p); 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 |
[入力例] 3 3 2 3 1 [出力例] 1 3 2 ※AtCoderのテストケースより [入力例] 3 3 1 1 1 [出力例] 1 2 3 ※AtCoderのテストケースより [入力例] 10 10 3 1 4 1 5 9 2 6 5 3 [出力例] 3 5 6 2 9 1 4 7 8 10 ※AtCoderのテストケースより [入力例] 20 50 15 6 9 15 9 13 6 15 8 6 4 13 3 12 6 5 2 16 9 6 19 18 6 13 12 13 20 10 1 1 5 13 15 20 14 11 2 9 16 1 6 8 5 5 4 2 20 12 11 16 [出力例] 16 11 12 20 2 4 5 8 6 1 9 14 15 13 10 18 19 3 7 17 |
■C++版プログラム(問題B/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 |
#include <bits/stdc++.h> using namespace std; using P = pair<double, double>; #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 P p[22]; int main(){ // 1. 入力情報. double x, y, a, b; int N; scanf("%lf %lf %d", &x, &y, &N); rep(i, N){ scanf("%lf %lf", &a, &b); p[i].a = a; p[i].b = b; } // 2. (N + 1)番目は, 1番目の頂点を埋める. p[N].a = p[0].a; p[N].b = p[0].b; // 3. 凸多角形を一周して確認. // (0, 0) から A = (ax, ay), B = (bx, by) について, // 垂線を引いたとして, H = (hx, hy) となったとすると, // k = (|B| * |B| - A * B) / (|A| * |A| - A * B) // とすると, // hx = (k * ax + bx) / (k + 1), hy = (k * ay + by) / (k + 1) となるように見える. // 0.0 <= k <= 1.0 に注意する必要がありそう. double ans = 123456789.0; rep(i, N){ // 3-1. vector A, B を 設定. double ax = p[i].a - x, ay = p[i].b - y; double bx = p[i + 1].a - x, by = p[i + 1].b - y; // 3-2. k の 値は? double mol = bx * bx + by * by - ax * bx - ay * by; double den = ax * ax + ay * ay - ax * bx - ay * by; double k = mol / den; // printf("k=%.12lf\n", k); // 3-3. k の 値が, [0, 1] の 範囲外であれば, Skip. if(k < 0.0 && k > 1.0) continue; // 3-4. hx, hy の 値は? double hx = (k * ax + bx) / (k + 1.0); double hy = (k * ay + by) / (k + 1.0); // 3-5. 距離は? ans = min(ans, hypot(hx, hy)); } // 4. 出力. printf("%.12lf\n", ans); 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 |
[入力例] 0 0 4 100 100 -100 100 -100 -100 100 -100 [出力例] 100 ※AtCoderのテストケースより ※ 但し, 上記のプログラムでは, 以下の内容が出力される. 100.000000000000 [入力例] 10 10 3 0 100 -100 -100 100 -100 [出力例] 31.3049516850 ※AtCoderのテストケースより ※ 但し, 上記のプログラムでは, 以下の内容が出力される. 31.304951684997 [入力例] 34 6 7 -43 -65 -23 -99 54 -68 65 92 16 83 -18 43 -39 2 [出力例] 25.0284205314 ※AtCoderのテストケースより ※ 但し, 上記のプログラムでは, 以下の内容が出力される. 25.028420531444 [入力例] 13 23 10 12 97 35 88 53 67 85 10 98 -10 75 -35 45 -55 25 -72 -23 -48 -15 25 [出力例] 26.919463855110 |
■参照サイト
AtCoder Regular Contest 042