通过 cf Round745-C Train Maintenance 浅谈分块

最新推荐文章于 2021-12-14 13:14:04 发布

South-twilight

最新推荐文章于 2021-12-14 13:14:04 发布

阅读量207

点赞数

分类专栏：总结文章标签：算法

本文链接：https://blog.csdn.net/qq_35776579/article/details/120602035

版权

总结专栏收录该内容

20 篇文章 0 订阅

订阅专栏

题目传送门

题意：有 $n$ 类火车，每类火车在工作 $x_i$ 天后，休息 $y_i$ 天。有 $m$ 天，给定每天的操作 $op_i$ ， $k_i$ ： $o p = 1$ ，则在第 $i$ 天将第 $k_i$ 类车投入使用； $o p = 2$ 则回收它。保证一类车只有一辆在运行。求每天有多少车辆在维修。

思考：这题感觉很熟悉，但是不会

题解：
听队友说，这道题是个标准的分块。原谅我分块基本没写过 -_-
但是根据我看过的分块题，大致总结一下分块的类型：

大于 $\sqrt n$ 的操作次数有限（小于 $\sqrt n$ ）；小于 $\sqrt n$ 种类有限（ $\sqrt n$ 类）
大于 $\sqrt n$ 的拆到各个块整体做（类似打标记的操作），对于边界部分暴力处理；小于 $\sqrt n$ 的部分直接暴力做。

这两类的 $\sqrt n$ 指代的东西都不同，前者是这道题，后者可以参考线段树区间修改，区间查询的题目。

对这道题，首先，对答案统计可以差分。然后我们考虑对 $x_i + y_i$ 分块，对于大于 $s q r t (m)$ 的部分，直接用一个数组 $a n s 2 [i]$ 暴力差分维护。对于小于的部分，我们考虑一下怎么把它进行转换为同一类。对所有的相同的 $x_i + y_i$ ，由于 $x_i + y_i$ 是一个循环节，也就是对开始的那一天 $t$ ，如果满足 $t \mod (x_i+y_i)$ 相同，那么他们在撤走之前，相同时间点的运行状态是一致的。

我们考虑用一个数组 $c [u] [v]$ 维护当前局面下 $x_i + y_i = u$ ，初始时间满足 $\mod u = v$ 的数的个数，这个可以进行暴力维护，然后每次把当前的状态取出来 $a n s 1 [i]$ 。 $O(\sqrt m)$

最后对 $a n s 2 [i]$ 统计前缀和后， $a n s 1 [i] + a n s 2 [i]$ 就是答案。

#include <bits/stdc++.h>
#define For(aa, bb, cc) for(int aa = (bb); aa <= (int)(cc); ++aa)
using namespace std;
const int maxn = 2e5 + 10;
const int bc = 450;
int n, m, bl;
int x[maxn], y[maxn];
int d[maxn];
int d2[bc][bc];
int ans[maxn], be[maxn];

void block_update(int p, int id, int v){
	int t = x[id] + y[id];
	int l = (p + x[id]) % t, r = (p - 1) % t;
	if(l <= r) {
		for(int i = l; i <= r; ++i) d2[t][i] += v;
	}
	else {
		for(int i = 0; i <= r; ++i) d2[t][i] += v;
		for(int i = l; i < t; ++i) d2[t][i] += v;
	}
}

int block_query(int p) { 
	int sum = 0;
	for(int i = 2; i <= bl; ++i) sum += d2[i][p % i];
	return sum;
}

int main() {
#ifndef ONLINE_JUDGE
	freopen("in.txt", "r", stdin);
	freopen("out.txt", "w", stdout);
#endif
	ios::sync_with_stdio(false);
	cin.tie(NULL);
	cin >> n >> m;
	bl = min((int)sqrt(m) + 1, bc);
	For(i, 1, n) cin >> x[i] >> y[i];
	For(i, 1, m) {
		int op, k;
		cin >> op >> k;
		if(op == 1) {
			if(x[k] + y[k] >= bl) {
				for(int j = i + x[k]; j <= m; j += x[k] + y[k]){
					++d[j];
					if(j + y[k] <= m) --d[j + y[k]];
				}
			}
			else block_update(i, k, 1);
			be[k] = i;
		}
		else {
			if(x[k] + y[k] >= bl) {
				for(int j = be[k] + x[k]; j <= m; j += x[k] + y[k]){
					if(j <= i && j + y[k] > i) { 
						--d[i];
						if(j + y[k] <= m) ++d[j + y[k]];
					}
					else if(j > i) {
						--d[j];
						if(j + y[k] <= m) ++d[j + y[k]];
					}
				}
			}
			else block_update(be[k], k, -1);
		}
		ans[i] = block_query(i);
	}
	For(i, 1, m) { 
		d[i] += d[i - 1];
		printf("%d\n", d[i] + ans[i]);
	}
	return 0;
}