#423045 (C++14) No.961 Vibrant Fillumination

提出ソース

結果

問題	No.961 Vibrant Fillumination
ユーザー	square1001
提出日時	2020-01-31 18:05:50
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	2,280 ms / 5,000 ms
コード長	4,762 bytes
コンパイル時間	1,862 ms
コンパイル使用メモリ	106,684 KB
実行使用メモリ	98,564 KB
最終ジャッジ日時	2023-10-17 08:24:46
合計ジャッジ時間	83,087 ms
ジャッジサーバーID （参考情報）	judge12 / judge15

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 4,348 KB
testcase_01	AC	2 ms 4,348 KB
testcase_02	AC	1,170 ms 23,748 KB
testcase_03	AC	1,122 ms 24,008 KB
testcase_04	AC	1,169 ms 24,008 KB
testcase_05	AC	1,155 ms 24,008 KB
testcase_06	AC	1,154 ms 23,484 KB
testcase_07	AC	1,163 ms 24,008 KB
testcase_08	AC	1,170 ms 24,012 KB
testcase_09	AC	1,158 ms 24,008 KB
testcase_10	AC	1,845 ms 88,796 KB
testcase_11	AC	1,841 ms 88,796 KB
testcase_12	AC	1,837 ms 88,796 KB
testcase_13	AC	1,842 ms 88,796 KB
testcase_14	AC	1,864 ms 88,796 KB
testcase_15	AC	1,864 ms 88,796 KB
testcase_16	AC	1,837 ms 88,796 KB
testcase_17	AC	1,897 ms 88,796 KB
testcase_18	AC	1,859 ms 88,796 KB
testcase_19	AC	1,829 ms 88,796 KB
testcase_20	AC	1,846 ms 88,796 KB
testcase_21	AC	1,857 ms 88,796 KB
testcase_22	AC	1,848 ms 88,796 KB
testcase_23	AC	1,837 ms 88,796 KB
testcase_24	AC	1,971 ms 88,796 KB
testcase_25	AC	1,842 ms 88,796 KB
testcase_26	AC	1,724 ms 11,972 KB
testcase_27	AC	1,562 ms 98,564 KB
testcase_28	AC	1,355 ms 12,008 KB
testcase_29	AC	1,626 ms 98,564 KB
testcase_30	AC	1,380 ms 98,192 KB
testcase_31	AC	1,269 ms 97,840 KB
testcase_32	AC	1,216 ms 97,772 KB
testcase_33	AC	2 ms 4,348 KB
testcase_34	AC	1,666 ms 83,084 KB
testcase_35	AC	1,834 ms 69,260 KB
testcase_36	AC	2,116 ms 78,808 KB
testcase_37	AC	2,280 ms 65,720 KB
testcase_38	AC	1,806 ms 65,036 KB
testcase_39	AC	1,903 ms 61,024 KB
testcase_40	AC	1,618 ms 51,308 KB
testcase_41	AC	1,665 ms 47,876 KB
testcase_42	AC	1,521 ms 40,720 KB
testcase_43	AC	1,295 ms 35,972 KB
testcase_44	AC	1,333 ms 35,724 KB
testcase_45	AC	1,378 ms 34,764 KB
testcase_46	AC	1,220 ms 32,656 KB
testcase_47	AC	1,238 ms 29,924 KB
testcase_48	AC	1,303 ms 32,036 KB
testcase_49	AC	1,187 ms 29,660 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

#include <cmath>
#include <vector>
#include <cassert>
#include <iostream>
#include <algorithm>
using namespace std;
int main() {
	cin.tie(0);
	ios_base::sync_with_stdio(false);
	int N;
	cin >> N;
	vector<long long> H(N);
	for (int i = 0; i < N; ++i) cin >> H[i];
	vector<int> A(N), B(N), C(N), D(N), E(N);
	for (int i = 0; i < N; ++i) {
		cin >> A[i] >> B[i] >> C[i] >> D[i] >> E[i];
		--E[i];
		if (E[i] < 0 || N <= E[i]) {
			assert(false);
		}
	}
	vector<vector<int> > G(N);
	for (int i = 0; i < N; ++i) {
		G[E[i]].push_back(i);
	}
	vector<int> perm(N);
	int active = 0;
	for (int i = 0; i < N; ++i) {
		perm[i] = i;
		if (!G[i].empty()) ++active;
	}
	sort(perm.begin(), perm.end(), [&](int i, int j) { return G[i].size() > G[j].size(); });
	int threshold = int(sqrt(N)) + 1;
	int lp = 0, sum = 0, cols = 0;
	vector<vector<int> > group;
	vector<int> col;
	for (int i = 0; i < active; ++i) {
		sum += G[perm[i]].size();
		if (sum >= threshold || i == active - 1) {
			if (lp == i) {
				int divs = max(int(G[perm[i]].size()) / threshold, 1);
				for (int j = 0; j < divs; ++j) {
					int dl = j * G[perm[i]].size() / divs, dr = (j + 1) * G[perm[i]].size() / divs;
					group.push_back(vector<int>(G[perm[i]].begin() + dl, G[perm[i]].begin() + dr));
					col.push_back(cols);
				}
			}
			else {
				group.push_back(vector<int>());
				for (int j = lp; j <= i; ++j) {
					group.back().insert(group.back().end(), G[perm[j]].begin(), G[perm[j]].end());
				}
				col.push_back(cols);
			}
			sum = 0, lp = i + 1;
			++cols;
		}
	}
	int Q;
	cin >> Q;
	vector<int> X(Q), Y(Q);
	for (int i = 0; i < Q; ++i) {
		cin >> X[i] >> Y[i];
	}
	vector<vector<long long> > sub(Q, vector<long long>(cols));
	for (int i = 0; i < group.size(); ++i) {
		vector<int> cps, cx = { 0, 2 * N }, cy = { 0, 2 * N };
		for (int j : group[i]) {
			cps.push_back(E[j]);
			cx.push_back(A[j]);
			cx.push_back(C[j]);
			cy.push_back(B[j]);
			cy.push_back(D[j]);
		}
		sort(cx.begin(), cx.end());
		cx.erase(unique(cx.begin(), cx.end()), cx.end());
		sort(cy.begin(), cy.end());
		cy.erase(unique(cy.begin(), cy.end()), cy.end());
		sort(cps.begin(), cps.end());
		cps.erase(unique(cps.begin(), cps.end()), cps.end());
		vector<vector<int> > F(cps.size());
		for (int j : group[i]) {
			F[lower_bound(cps.begin(), cps.end(), E[j]) - cps.begin()].push_back(j);
		}
		int h = cx.size(), w = cy.size();
		vector<vector<long long> > subdat(h + 1, vector<long long>(w + 1));
		for (int j = 0; j < F.size(); ++j) {
			vector<int> csx, csy;
			for (int k : F[j]) {
				csx.push_back(A[k]);
				csx.push_back(C[k]);
				csy.push_back(B[k]);
				csy.push_back(D[k]);
			}
			sort(csx.begin(), csx.end());
			csx.erase(unique(csx.begin(), csx.end()), csx.end());
			sort(csy.begin(), csy.end());
			csy.erase(unique(csy.begin(), csy.end()), csy.end());
			vector<vector<int> > imos(csx.size(), vector<int>(csy.size()));
			for (int k : F[j]) {
				int xa = lower_bound(csx.begin(), csx.end(), A[k]) - csx.begin();
				int ya = lower_bound(csy.begin(), csy.end(), B[k]) - csy.begin();
				int xb = lower_bound(csx.begin(), csx.end(), C[k]) - csx.begin();
				int yb = lower_bound(csy.begin(), csy.end(), D[k]) - csy.begin();
				++imos[xa][ya];
				++imos[xb][yb];
				--imos[xa][yb];
				--imos[xb][ya];
			}
			for (int k = 1; k < csx.size(); ++k) {
				for (int l = 0; l < csy.size(); ++l) {
					imos[k][l] += imos[k - 1][l];
				}
			}
			for (int k = 0; k < csx.size(); ++k) {
				for (int l = 1; l < csy.size(); ++l) {
					imos[k][l] += imos[k][l - 1];
				}
			}
			for (int k = 0; k < csx.size(); ++k) {
				csx[k] = lower_bound(cx.begin(), cx.end(), csx[k]) - cx.begin();
			}
			for (int k = 0; k < csy.size(); ++k) {
				csy[k] = lower_bound(cy.begin(), cy.end(), csy[k]) - cy.begin();
			}
			for (int k = 0; k < csx.size() - 1; ++k) {
				for (int l = 0; l < csy.size() - 1; ++l) {
					if (imos[k][l] == 0) continue;
					subdat[csx[k]][csy[l]] ^= H[cps[j]];
					subdat[csx[k + 1]][csy[l + 1]] ^= H[cps[j]];
					subdat[csx[k]][csy[l + 1]] ^= H[cps[j]];
					subdat[csx[k + 1]][csy[l]] ^= H[cps[j]];
				}
			}
		}
		for (int j = 1; j < cx.size(); ++j) {
			for (int k = 0; k < cy.size(); ++k) {
				subdat[j][k] ^= subdat[j - 1][k];
			}
		}
		for (int j = 0; j < cx.size(); ++j) {
			for (int k = 1; k < cy.size(); ++k) {
				subdat[j][k] ^= subdat[j][k - 1];
			}
		}
		for (int j = 0; j < Q; ++j) {
			if (sub[j][col[i]] != 0) continue;
			int px = lower_bound(cx.begin(), cx.end(), X[j] + 1) - cx.begin() - 1;
			int py = lower_bound(cy.begin(), cy.end(), Y[j] + 1) - cy.begin() - 1;
			sub[j][col[i]] = subdat[px][py];
		}
	}
	for (int i = 0; i < Q; ++i) {
		long long ans = 0;
		for (int j = 0; j < cols; ++j) {
			ans ^= sub[i][j];
		}
		cout << ans << '\n';
	}
	return 0;
}

yukicoder

結果

テストケース

ソースコード