제출 #552705

#제출 시각아이디문제언어결과실행 시간메모리
552705Zanite가로등 (APIO19_street_lamps)C++17
20 / 100
5068 ms266624 KiB
// "I assure you that you guys won't make it to the top 4"
// - Tzaph, 21 December 2021

#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

#define ll long long
#define ld long double
#define si short int
#define i8 __int128
#define pii pair<int, int>
#define pll pair<ll, ll>
#define pld pair<ld, ld>
#define psi pair<si, si>
#define pi8 pair<i8, i8>
#define pq priority_queue
#define fi first
#define se second

#define sqr(x) ((x)*(x))
#define pow2(x) (1ll << (x))
#define debug(x) cout << #x << " = " << (x) << '\n'
#define debugV(x, a) cout << #x << "[" << (a) << "] = " << (x[a]) << '\n'

#define yume using
#define wo namespace
#define kanaeyo std
#define nazotte __gnu_pbds

yume wo kanaeyo;
yume wo nazotte;

template<typename T>
using ordered_set = tree<T, null_type, less<T>, rb_tree_tag,tree_order_statistics_node_update>;

template<typename T> void chmin(T &a, T b) {a = min(a, b);}
template<typename T> void chmax(T &a, T b) {a = max(a, b);}
 
template<typename T> void maddto(T &a, T b, T mod) {a += b; a %= mod;}
template<typename T> void msubto(T &a, T b, T mod) {a -= b; while (a < 0) a += mod; a %= mod;}
template<typename T> void mmulto(T &a, T b, T mod) {a *= b; a %= mod;}
 
template<typename T> T madd(T a, T b, T mod) {a += b; a %= mod; return a;}
template<typename T> T msub(T a, T b, T mod) {a -= b; while (a < 0) a += mod; return a;}
template<typename T> T mmul(T a, T b, T mod) {a *= b; a %= mod; return a;}

const ll ModA = 998244353;
const ll ModC = 1e9 + 7;
const ll INF = 1e18;
const ll iINF = 1e9;

const ld EPS = 1e-9;
const ld iEPS = 1e-6;

const int maxN = 3e5 + 5;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());

inline pll add(pll a, pll b) {
	return {a.fi + b.fi, a.se + b.se};
}

struct Treap {
	struct Node {
		pll val, sum;
		int prior, pos, mn, mx;
		Node *left, *right;

		Node(int _pos, pll _val) {
			pos = mn = mx = _pos;
			val = sum = _val;
			prior = rng();
			left = right = nullptr;
		}
		
		void update() {
			sum = val;
			if (left) {sum = add(sum, left->sum);}
			if (right) {sum = add(sum, right->sum);}

			mn = (left ? left->mn : pos);
			mx = (right ? right->mx : pos);
		}
	};

	Node *root;

	Treap() : root(nullptr) {}

	void split(Node *T, Node *&L, Node *&R, int sep) {
		// splits T into L = [mn..pos] and R = [pos+1..mx]
		if (!T) {
			L = R = nullptr;
			return;
		}

		if (T->pos <= sep) {
			split(T->right, T->right, R, sep);
			L = T;
		} else {
			split(T->left, L, T->left, sep);
			R = T;
		}

		T->update();
	}

	void merge(Node *&T, Node *L, Node *R) {
		// merges L and R into T, assuming that all elements in L < R
		if (!L || !R) {
			T = (L ? L : R);
			return;
		}

		if (L->prior >= R->prior) {
			merge(L->right, L->right, R);
			T = L;
		} else {
			merge(R->left, L, R->left);
			T = R;
		}

		T->update();
	}

	bool find(int pos) {
		// returns TRUE if the node with position pos is present in the treap
		Node *cur = root;
		while (cur) {
			if (cur->pos == pos) {return true;}
			if (cur->pos < pos) {cur = cur->right;}
			else {cur = cur->left;}
		}
		return false;
	}

	void updateNode(Node *cur, int pos, pll val) {
		// updates node at position pos with val, assuming the node exists
		if (cur->pos == pos) {
			cur->val.fi += val.fi;
			cur->val.se += val.se;
			cur->update();
			return;
		}
		
		if (cur->pos < pos) {
			updateNode(cur->right, pos, val);
		} else {
			updateNode(cur->left, pos, val);
		}
		cur->update();
	}

	void update(int pos, pll val) {
		// if the node exists, update it
		if (find(pos)) {updateNode(root, pos, val);}
		else {
			// insert a new node
			Node *tmp = new Node(pos, val);
			Node *t1, *t2;
			split(root, t1, t2, pos);
			merge(root, t1, tmp);
			merge(root, root, t2);
		}
	}

	pll query(Node *T, int L, int R) {
		// returns sum of all nodes in the range [L..R] under node T
		if (!T) return {0, 0};
		if (T->mx < L || R < T->mn) return {0, 0};
		if (L <= T->mn && T->mx <= R) return T->sum;

		pll ret = (L <= T->pos && T->pos <= R) ? (T->val) : make_pair(0ll, 0ll);
		if (T->left) {ret = add(ret, query(T->left, L, R));}
		if (T->right) {ret = add(ret, query(T->right, L, R));}
		return ret;
	}

	pll query(int L, int R) {
		// returns sum of all nodes in the range [L..R]
		return query(root, L, R);
	}

	void print(Node *T) {
		if (!T) return;
		printf("%d | {%lld, %lld}\n", T->pos, T->val.fi, T->val.se);
		print(T->left);
		print(T->right);
	}

	void print() {
		print(root);
	}
};

struct DynamicSegtreeBBST {
	struct Node {
		Treap inner;
		Node* left;
		Node* right;

		Node() {
			inner = Treap();
			left = right = nullptr;
		}
	};

	Node* root;
	int dim1, dim2;

	DynamicSegtreeBBST(int _dim1, int _dim2) {
		root = new Node();
		dim1 = _dim1, dim2 = _dim2;
	}

	void update(Node* cur, int l, int r, int d1, int d2, pll upd) {
		cur->inner.update(d2, upd);
		if (l == r) {return;}

		int m = (l + r) >> 1;
		if (d1 <= m) {
			if (!cur->left) {cur->left = new Node();}
			update(cur->left, l, m, d1, d2, upd);
		} else {
			if (!cur->right) {cur->right = new Node();}
			update(cur->right, m+1, r, d1, d2, upd);
		}
	}

	pll query(Node* cur, int l, int r, int d1l, int d1r, int d2l, int d2r) {
		if (l > d1r || d1l > r) {return {0, 0};}
		if (d1l <= l && r <= d1r) {
			return cur->inner.query(d2l, d2r);
		}

		int m = (l + r) >> 1;
		pll ret = {0, 0};
		if (cur->left) {ret = add(ret, query(cur->left, l, m, d1l, d1r, d2l, d2r));}
		if (cur->right) {ret = add(ret, query(cur->right, m+1, r, d1l, d1r, d2l, d2r));}
		return ret;
	}

	void print(Node* cur, ll l, ll r) {
		if (!cur) return;
		if (l == r) {
			debug(l);
			cur->inner.print();
			return;
		}

		ll m = (l + r) >> 1;
		print(cur->left, l, m);
		print(cur->right, m+1, r);
	}

	void print() {
		print(root, 1, dim1);
	}
};

struct PairTreap {
	struct Node {
		pll val;
		int prior;
		Node *left, *right;

		Node(pll _val) {
			val = _val;
			prior = rng();
			left = right = nullptr;
		}
	};

	Node *root;

	PairTreap() : root(nullptr) {}

	void split(Node *T, Node *&L, Node *&R, ll sep) {
		// splits T into L = [mn..sep] and R = [sep+1..mx]
		if (!T) {
			L = R = nullptr;
			return;
		}

		if (T->val.fi <= sep) {
			split(T->right, T->right, R, sep);
			L = T;
		} else {
			split(T->left, L, T->left, sep);
			R = T;
		}
	}

	void merge(Node *&T, Node *L, Node *R) {
		// merges L and R into T, assuming that all elements in L < R
		if (!L || !R) {
			T = (L ? L : R);
			return;
		}

		if (L->prior >= R->prior) {
			merge(L->right, L->right, R);
			T = L;
		} else {
			merge(R->left, L, R->left);
			T = R;
		}
	}

	pll find(Node *cur, int x) {
		// returns the largest node in the treap whose first element is <= x
		if (cur->val.fi <= x) {
			pll ret = cur->val;
			if (cur->right) {ret = max(ret, find(cur->right, x));}
			return ret;
		} else {
			if (cur->left) return find(cur->left, x);
			else return {-1, -1};
		}
	}

	pll find(int x) {
		return find(root, x);
	}

	void insert(pll val) {
		// insert a new node
		Node *tmp = new Node(val);
		Node *t1, *t2;
		split(root, t1, t2, val.fi);
		merge(root, t1, tmp);
		merge(root, root, t2);
	}

	void erase(pll val) {
		Node *t1, *t2, *t3, *t4;
		split(root, t1, t2, val.fi);
		split(t1, t3, t4, val.fi - 1);
		merge(root, t3, t2);
	}

	void print(Node* cur) {
		if (!cur) return;
		print(cur->left);
		cout << cur->val.fi << ' ' << cur->val.se << '\n';
		print(cur->right);
	}

	void print() {
		print(root);
	}
};

int n, q;
bool lamps[maxN];

int main() {
	scanf("%d %d", &n, &q);

	char c;
	for (int i = 1; i <= n; i++) {
		scanf(" %c", &c);
		lamps[i] = (c == '1');
	}

	DynamicSegtreeBBST D = DynamicSegtreeBBST(n, n);
	PairTreap segments = PairTreap();

	for (int start = -1, i = 1; i <= n+1; i++) {
		if (lamps[i]) {
			if (start == -1) {start = i;}
		} else {
			if (start != -1) {
				segments.insert({start, (i-1)});
				D.update(D.root, 1, n, start, i-1, {q, 1});
				start = -1;
			}
		}
	}
	//segments.print();

	char typ[10];
	for (int a, b, i = 1; i <= q; i++) {
		scanf(" %s", typ);
		if (typ[0] == 't') {
			scanf("%d", &a);
			if (lamps[a]) {
				// toggle the lamp off

				// remove segment of current lamp
				pll cseg = segments.find(a);
				//cout << "cseg " << cseg.fi << ' ' << cseg.se << '\n';
				D.update(D.root, 1, n, cseg.fi, cseg.se, {-(q-i), -1});
				segments.erase(cseg);

				// insert new segments (except if it's of length 0)

				// left segment
				if (a != cseg.fi) {
					segments.insert({cseg.fi, (a-1)});
					D.update(D.root, 1, n, cseg.fi, a-1, {q-i, 1});
				}
				// right segment
				if (a != cseg.se) {
					segments.insert({(a+1), cseg.se});
					D.update(D.root, 1, n, a+1, cseg.se, {q-i, 1});
				}

			} else {
				// toggle the lamp on
				pll newseg = {a, a};

				// check if lamps next to it are on

				// left
				if (lamps[a-1]) {
					pll lseg = segments.find(a-1);
					D.update(D.root, 1, n, lseg.fi, lseg.se, {-(q-i), -1});
					newseg.fi = lseg.fi;
					segments.erase(lseg);
				}
				// right
				if (lamps[a+1]) {
					pll rseg = segments.find(a+1);
					D.update(D.root, 1, n, rseg.fi, rseg.se, {-(q-i), -1});
					newseg.se = rseg.se;
					segments.erase(rseg);
				}

				// insert new segment
				segments.insert(newseg);
				D.update(D.root, 1, n, newseg.fi, newseg.se, {q-i, 1});
			}

			lamps[a] = !lamps[a];

		} else {
			scanf("%d %d", &a, &b);

			pll ans = D.query(D.root, 1, n, 1, a, b-1, n);
			//cout << ans.fi << ' ' << ans.se << '\n';
			ll fin = ans.fi - (1ll * ans.se * (q-i));
			printf("%lld\n", fin);
		}

		/*
		debug(i);
		segments.print();

		//for (auto i : segments) {cout << "{" << i.fi << ", " << i.se << "} ";}
		cout << '\n';
		D.print();
		cout << '\n';
		*/
	}
}

컴파일 시 표준 에러 (stderr) 메시지

street_lamps.cpp: In function 'int main()':
street_lamps.cpp:363:7: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  363 |  scanf("%d %d", &n, &q);
      |  ~~~~~^~~~~~~~~~~~~~~~~
street_lamps.cpp:367:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  367 |   scanf(" %c", &c);
      |   ~~~~~^~~~~~~~~~~
street_lamps.cpp:389:8: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  389 |   scanf(" %s", typ);
      |   ~~~~~^~~~~~~~~~~~
street_lamps.cpp:391:9: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  391 |    scanf("%d", &a);
      |    ~~~~~^~~~~~~~~~
street_lamps.cpp:443:9: warning: ignoring return value of 'int scanf(const char*, ...)' declared with attribute 'warn_unused_result' [-Wunused-result]
  443 |    scanf("%d %d", &a, &b);
      |    ~~~~~^~~~~~~~~~~~~~~~~
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...