C++の練習を兼ねて, AtCoder Beginner Contest 253 の 問題F (Operations on a Matrix) を解いてみた.
■感想.
1. 問題Fは, 方針が見えなかったので, 解説を参考にして, AC版に到達できたと思う.
2. 実装に苦労したものの, Binary Indexed Tree の 復習が出来たので, 非常に良かったと思う.
※ 公式のライブラリを拝借させて頂いてます.
3. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト AtCoder Beginner Contest 253 解説 の 各リンク を ご覧下さい.
■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 101 102 103 104 105 106 107 108 109 110 |
// 解き直し. // https://atcoder.jp/contests/abc253/editorial/4025 // C++(GCC 9.2.1) #include <bits/stdc++.h> using namespace std; using LL = long long; using vi = vector<int>; using vvi = vector<vi>; #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 const int MAX = 202020; LL rMatrix[MAX], qSub[MAX]; int qLatest[MAX]; int qf[MAX], ql[MAX], qr[MAX], qx[MAX], qi[MAX], qj[MAX]; // Binary Indexed Tree (Fenwick Tree) // https://youtu.be/lyHk98daDJo?t=7960 template<typename T> struct BIT{ int n; vector<T> d; BIT(int n = 0) : n(n), d(n + 1) {} void add(int i, T x = 1){ for(i++; i <= n; i += i & -i) d[i] += x; } T sum(int i){ T x = 0; for(i++; i; i -= i & -i) x += d[i]; return x; } T sum(int l, int r){ return sum(r - 1) - sum(l - 1); } }; int main(){ // 1. 入力情報. int N, M, Q; scanf("%d %d %d", &N, &M, &Q); vvi q3(Q); rep(q, Q){ // 1-1. クエリ入力情報. scanf("%d", &qf[q]); // 1-2. クエリ 1. if(qf[q] == 1) scanf("%d %d %d", &ql[q], &qr[q], &qx[q]); // 1-3. クエリ 2. if(qf[q] == 2){ scanf("%d %d", &qi[q], &qx[q]); // 直近 の クエリ 2 番号 を 保存. qLatest[qi[q]] = q; } // 1-4. クエリ 3. if(qf[q] == 3){ scanf("%d %d", &qi[q], &qj[q]); // 直近 の クエリ 2 番号 に対して, クエリ 3 番号 を 紐づける. q3[qLatest[qi[q]]].pb(q); } } // 2. 事前準備. BIT<LL> pre(M + 2); rep(q, Q){ // 2-1. クエリ 1. if(qf[q] == 1){ // 位置 l で, x を 加算. pre.add(ql[q], qx[q]); // 位置 r + 1 で, x を 減算. pre.add(qr[q] + 1, -qx[q]); } // 2-2. クエリ 2. if(qf[q] == 2){ // クエリ 3 で 減算分 を 保存. for(auto &x : q3[q]) qSub[x] = pre.sum(1, qj[x] + 1); } } // 3. クエリ回答. BIT<LL> cur(M + 2); rep(q, Q){ // 3-1. クエリ 1. if(qf[q] == 1){ // 位置 l で, x を 加算. cur.add(ql[q], qx[q]); // 位置 r + 1 で, x を 減算. cur.add(qr[q] + 1, -qx[q]); } // 3-2. クエリ 2. if(qf[q] == 2) rMatrix[qi[q]] = qx[q]; // 3-3. クエリ 3. if(qf[q] == 3) printf("%lld\n", cur.sum(1, qj[q] + 1) + rMatrix[qi[q]] - qSub[q]); } 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 |
[入力例] 3 3 9 1 1 2 1 3 2 2 2 3 2 3 3 3 3 3 1 1 2 3 3 3 3 2 3 2 3 3 1 2 [出力例] 1 2 2 5 3 4 ※AtCoderテストケースより [入力例] 1 1 10 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 1 1 1 1000000000 3 1 1 [出力例] 9000000000 ※AtCoderテストケースより [入力例] 10 10 10 1 1 8 5 2 2 6 3 2 1 3 4 7 1 5 9 7 3 3 2 3 2 8 2 8 10 3 8 8 3 1 10 [出力例] 6 5 5 13 10 0 ※AtCoderテストケースより [入力例] 30 30 30 1 3 22 77 2 3 111 3 1 5 3 3 1 3 3 3 1 8 27 33 1 3 15 30 2 5 222 1 9 30 55 3 10 15 3 12 17 1 2 12 1 3 8 8 3 15 15 2 12 333 2 17 444 3 17 8 3 21 18 1 7 27 345 1 1 11 123 3 25 15 3 10 10 2 26 666 2 17 777 3 17 15 1 3 23 334 3 17 15 1 5 25 889 3 17 1 3 17 15 [出力例] 77 111 111 195 165 141 195 444 165 540 664 777 1111 777 2000 |
■参照サイト
NOMURA プログラミングコンテスト2022(AtCoder Beginner Contest 253)