LOJ2790&DTOJ2809 [CEOI2015 Day2]核能国度(Nuclearia)

最新推荐文章于 2022-08-11 18:38:53 发布
最新推荐文章于 2022-08-11 18:38:53 发布
阅读量219 收藏
点赞数
分类专栏： 题解 # DTOJ # 其他题解
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43849488/article/details/104443878
版权
题解 同时被 3 个专栏收录
66 篇文章 0 订阅
订阅专栏
DTOJ
58 篇文章 0 订阅
订阅专栏
其他题解
7 篇文章 0 订阅
订阅专栏

LOJ2790&DTOJ2809 [CEOI2015 Day2]核能国度

题目

题目描述

原题
核能国可以看作一个由 W × H W \times H W×H的方格组成的矩形。核能国有 N N N个核电站，每个核电站占用一个方格。不幸的是，核能国遭遇了百年一遇的特大地震，导致所有的核电站都发生了核泄漏。
每个核电站的核泄漏程度可以用两个整数 a , b a, b a,b来表示。如果位于 P = [ x P , y P ] P=[x_P,y_P] P=[xP​,yP​]的核电站爆炸，方格 C = [ x C , y C ] C=[x_C,y_C] C=[xC​,yC​]会增加 m a x ( 0 , a − b × d ( P , C ) ) \mathrm{max}(0,a-b\times d(P,C)) max(0,a−b×d(P,C))贝克的辐射（贝克是单位），其中 d ( P , C ) d(P,C) d(P,C)是两个方格的切比雪夫距离，即 d ( P , C ) = m a x ( ∣ x P − x C ∣ , ∣ y P − y C ∣ ) d(P,C) =\mathrm{max}(|x_P - x_C|,|y_P - y_C|) d(P,C)=max(∣xP​−xC​∣,∣yP​−yC​∣)。
一个方格可能会受到多处核泄漏的影响。
例如，如果一个 a = 7 , b = 3 a = 7, b = 3 a=7,b=3的核电站爆炸了，所在的方格 X X X会受到 7 7 7贝克辐射（贝克是单位），满足 d ( X , Y ) = 1 d(X,Y) = 1 d(X,Y)=1的 8 8 8个方格 Y Y Y会受到 4 4 4贝克辐射，满足 d ( X , Z ) = 2 d(X,Z) = 2 d(X,Z)=2的 16 16 16个方格 Z Z Z会受到 1 1 1贝克辐射。
环保部门给了你 Q Q Q组询问，每组询问会划定核能国领土中的一个矩形，请回答：矩形区域内（每个方格）所受的平均辐射量为多少。

输入格式

第一行，两个正整数 W W W和 H ( W × H ⩽ 2.5 × 1 0 6 ) H(W × H \leqslant 2.5×10^6) H(W×H⩽2.5×106)，分别表示核能国的宽度与高度。
第二行，一个正整数 N N N，表示核电站的个数 ( 1 ≤ N ⩽ 2 × 1 0 5 ) (1 \leq N \leqslant 2×10^5) (1≤N⩽2×105)。
在接下来的 N N N行中，每行四个正整数 x i , y i , a i , b i ( 1 ⩽ x i ⩽ W , 1 ⩽ y i ⩽ H , 1 ⩽ a i , b i ⩽ 1 0 9 ) x_i,y_i,a_i,b_i(1 \leqslant x_i \leqslant W,1 \leqslant y_i \leqslant H,1 \leqslant a_i,b_i \leqslant 10^9) xi​,yi​,ai​,bi​(1⩽xi​⩽W,1⩽yi​⩽H,1⩽ai​,bi​⩽109)，表示有一个核电站位于方格 [ x i , y i ] [x_i,y_i] [xi​,yi​]，它的参数为 a i a_i ai​与 b i ​ b_i​ bi​​。每个格子最多有一个核电站。
第 N + 3 N+3 N+3行，一个正整数 Q Q Q，表示询问的次数 ( 1 ≤ Q ≤ 2 × 1 0 5 ) (1 \leq Q \leq 2×10^5) (1≤Q≤2×105)。
在接下来的 Q Q Q行中，每行四个 正整数 x 1 j , y 1 j , x 2 j , y 2 j ( 1 ⩽ x 1 j ⩽ x 2 j ⩽ W , 1 ⩽ y 1 j ⩽ y 2 j ⩽ H ) x_{1j},y_{1j},x_{2j},y_{2j}(1 \leqslant x_{1j} \leqslant x_{2j} \leqslant W,1 \leqslant y_{1j} \leqslant y_{2j} \leqslant H) x1j​,y1j​,x2j​,y2j​(1⩽x1j​⩽x2j​⩽W,1⩽y1j​⩽y2j​⩽H)，表示该询问给出的矩形区域的左上角在 [ x 1 j , y 1 j ] [x_{1j},y_{1j}] [x1j​,y1j​]且它的右下角在 [ x 2 j , y 2 j ] [x_{2j},y_{2j}] [x2j​,y2j​]。
你可以假设核能国内的总辐射量少于 2 63 2^{63} 263。

输出格式

对于每个询问，输出一行表示给定矩形区域内所有方格的平均辐射量，四舍五入至整数。

样例

样例输入

4 3
2
1 1 7 3
3 2 4 2
4
1 2 2 3
1 1 4 3
4 2 4 2
1 3 4 3

样例输出

4
4
2
2

样例解释

以下为两次爆炸后对每个方格产生的辐射量：

7 6 3 2
4 6 5 2
1 3 3 2
  1. 2 2 2^2 22方形区域内的总辐射为 14 14 14，所以平均值为 14 ÷ 4 = 3.5 14\div 4=3.5 14÷4=3.5，四舍五入至 4 4 4。
  2. 整个核能国的总辐射为 44 44 44，所以平均值为 44 ÷ 12 ≈ 3.67 44\div 12 \approx 3.67 44÷12≈3.67，四舍五入至 4 4 4。
  3. 单个格子的平均辐射量就是它所受到的辐射量。
  4. 最后一行的平均辐射量为 9 ÷ 4 = 2.25 9\div 4=2.25 9÷4=2.25，四舍五入至 2 2 2。

数据范围与提示

有 14 14 14组测试数据。奇数的测试组只包含 a a a是 b b b的倍数的核电站。对每个子任务的进一步限制如下：

测试组进一步限制分数
1 H = 1 , N ⋅ W ⩽ 1 0 8 , Q ⋅ W ⩽ 1 0 8 H=1,N\cdot W \leqslant 10^8,Q \cdot W \leqslant 10^8 H=1,N⋅W⩽108,Q⋅W⩽1083
2 H = 1 , N ⋅ W ⩽ 1 0 8 , Q ⋅ W ⩽ 1 0 8 H=1,N\cdot W \leqslant 10^8,Q \cdot W \leqslant 10^8 H=1,N⋅W⩽108,Q⋅W⩽1082
3 N ⋅ W ⋅ H ≤ 1 0 8 , Q ⋅ W ⋅ H ≤ 1 0 8 N\cdot W \cdot H \leq 10^8,Q \cdot W \cdot H \leq 10^8 N⋅W⋅H≤108,Q⋅W⋅H≤1083
4 N ⋅ W ⋅ H ≤ 1 0 8 , Q ⋅ W ⋅ H ≤ 1 0 8 N\cdot W \cdot H \leq 10^8,Q \cdot W \cdot H \leq 10^8 N⋅W⋅H≤108,Q⋅W⋅H≤1082
5 H = 1 , N ⋅ W ≤ 1 0 8 H=1,N\cdot W \leq 10^8 H=1,N⋅W≤1086
6 H = 1 , N ⋅ W ≤ 1 0 8 H=1,N\cdot W \leq 10^8 H=1,N⋅W≤1084
7 N ⋅ W ⋅ H ≤ 1 0 8 N\cdot W \cdot H \leq 10^8 N⋅W⋅H≤1086
8 N ⋅ W ⋅ H ≤ 1 0 8 N\cdot W \cdot H \leq 10^8 N⋅W⋅H≤1084
9 H = 1 H=1 H=115
10 H = 1 H=1 H=110
11没有符合界限定义的爆炸事件15
12没有符合界限定义的爆炸事件10
13无12
14无8

如果核电站位于核能国的边境或是在离边境稍近的位置，那么爆炸可能也会影响到核能国之外的方格。影响到核能国外方格的爆炸被称作界限。

题解

暴力算法1

最暴力的算法你能想到什么？枚举！
对于每个核电站爆炸事件，枚举它周围的方格受到的影响
直接枚举整个国度显得太暴力了，我们能不能稍稍优化一下呢？
显然是可以的
我们发现，每个核电站爆炸事件有一个“势力范围”
只有在这个方格和核电站的切比雪夫距离 ⩽ a b \leqslant\frac{a}{b} ⩽ba​时，它才会受影响
所以，我们只需要枚举这个核电站的势力范围就可以了

//官方程序
#include <stdio.h>
#include <stdlib.h>

#define MAXWH 10000000
#define MAXN 1000000
#define MAXQ 1000000

typedef long long int huge;

int w, h, n, q;

typedef struct NUCLEARIA
{
	huge Info[MAXWH];
	huge& operator () (int x, int y) {return Info[(y * w) + x];}//其实就是个二维数组
}
NUCLEARIA;

typedef struct PLANT
{
	int x;
	int y;
	int a;
	int b;
}
PLANT;

typedef struct QUERY
{
	int x1;
	int y1;
	int x2;
	int y2;
}
QUERY;

NUCLEARIA Nuclearia;
PLANT Plant[MAXN];
QUERY Query[MAXQ];

int min(int a, int b)
{
	return (a < b) ? a : b;
}

int max(int a, int b)
{
	return (a > b) ? a : b;
}

int abs(int a)
{
	return (a < 0) ? -a : a;
}

int dist(int x1, int y1, int x2, int y2)
{
	return max(abs(x1 - x2), abs(y1 - y2));
}

void Print(huge sum, int area)//四舍五入
{
	huge rsl = sum / area;
	if((sum % area) * 2 >= area)
	{
		rsl++;
	}
	
	printf("%lld\n", rsl);
}

int main()
{
	scanf("%d%d", &w, &h);
	
	scanf("%d", &n);
	for(int i = 0; i < n; i++)
	{
		scanf("%d%d%d%d", &(Plant[i].x), &(Plant[i].y), &(Plant[i].a), &(Plant[i].b));
		Plant[i].x--;
		Plant[i].y--;
	}
	
	scanf("%d", &q);
	for(int i = 0; i < q; i++)
	{
		scanf("%d%d%d%d", &(Query[i].x1), &(Query[i].y1), &(Query[i].x2), &(Query[i].y2));
		Query[i].x1--;
		Query[i].y1--;
	}
	
	for(int i = 0; i < n; i++)
	{
		PLANT& P = Plant[i];
		int d = (P.a - 1) / P.b;
		int x1 = max(0, P.x - d);
		int y1 = max(0, P.y - d);
		int x2 = min(w, P.x + d + 1);
		int y2 = min(h, P.y + d + 1);
		//势力范围
		for(int x = x1; x < x2; x++)
		{
			for(int y = y1; y < y2; y++)
			{
				int d = dist(x, y, P.x, P.y);
				Nuclearia(x, y) += P.a - (d * P.b);
			}
		}
	}
	
	for(int i = 0; i < q; i++)
	{
		QUERY& Q = Query[i];
		huge rsl = 0;
		
		for(int x = Q.x1; x < Q.x2; x++)
		{
			for(int y = Q.y1; y < Q.y2; y++)
			{
				rsl += Nuclearia(x, y);
			}
		}
		
		Print(rsl, (Q.x2 - Q.x1) * (Q.y2 - Q.y1));
	}
	
	return 0;
}

暴力算法2

我们发现上面的算法在询问次数很多时是很耗时间的，因为每次我们都需要把矩形的辐射量加一遍
所以，我们可以在询问前先预处理一下
大家应该知道激光炸弹吧？我们可以借用一下这种矩阵的预处理方法，就是计算出左上角在 [ 1 , 1 ] [1,1] [1,1]且它的右下角在 [ x , y ] [x,y] [x,y]的矩形的总辐射，假设是 N u c l e a r i a ( x , y ) Nuclearia(x,y) Nuclearia(x,y)
最终的答案就是 N u c l e a r i a ( x 2 j , y 2 j ) − N u c l e a r i a ( x 1 j − 1 , y 2 j ) − N u c l e a r i a ( x 2 j , y 1 j − 1 ) + N u c l e a r i a ( x 1 j − 1 , y 1 j − 1 ) Nuclearia(x_{2j},y_{2j})-Nuclearia(x_{1j}-1,y_{2j})-Nuclearia(x_{2j},y_{1j}-1)+Nuclearia(x_{1j}-1,y_{1j}-1) Nuclearia(x2j​,y2j​)−Nuclearia(x1j​−1,y2j​)−Nuclearia(x2j​,y1j​−1)+Nuclearia(x1j​−1,y1j​−1)

//官方程序
#include <stdio.h>
#include <stdlib.h>

#define MAXWH 10000000
#define MAXN 1000000
#define MAXQ 1000000

typedef long long int huge;

int w, h, n, q;

typedef struct NUCLEARIA
{
	huge Info[MAXWH];
	huge& operator () (int x, int y) {return Info[(y * w) + x];}//其实就是个二维数组
}
NUCLEARIA;

typedef struct PLANT
{
	int x;
	int y;
	int a;
	int b;
}
PLANT;

typedef struct QUERY
{
	int x1;
	int y1;
	int x2;
	int y2;
}
QUERY;

NUCLEARIA Nuclearia;
PLANT Plant[MAXN];
QUERY Query[MAXQ];

int min(int a, int b)
{
	return (a < b) ? a : b;
}

int max(int a, int b)
{
	return (a > b) ? a : b;
}

int abs(int a)
{
	return (a < 0) ? -a : a;
}

int dist(int x1, int y1, int x2, int y2)
{
	return max(abs(x1 - x2), abs(y1 - y2));
}

void Summarize()//预处理
{
	for(int x = 0; x < w; x++)
	{
		for(int y = 0; y < h; y++)
		{
			if(x) Nuclearia(x, y) += Nuclearia(x - 1, y);
			if(y) Nuclearia(x, y) += Nuclearia(x, y - 1);
			if((x) && (y)) Nuclearia(x, y) -= Nuclearia(x - 1, y - 1);
		}
	}
}

void Print(huge sum, int area)//四舍五入
{
	huge rsl = sum / area;
	if((sum % area) * 2 >= area)
	{
		rsl++;
	}
	
	printf("%lld\n", rsl);
}

int main()
{
	scanf("%d%d", &w, &h);
	
	scanf("%d", &n);
	for(int i = 0; i < n; i++)
	{
		scanf("%d%d%d%d", &(Plant[i].x), &(Plant[i].y), &(Plant[i].a), &(Plant[i].b));
		Plant[i].x--;
		Plant[i].y--;
	}
	
	scanf("%d", &q);
	for(int i = 0; i < q; i++)
	{
		scanf("%d%d%d%d", &(Query[i].x1), &(Query[i].y1), &(Query[i].x2), &(Query[i].y2));
		Query[i].x1 -= 2;
		Query[i].y1 -= 2;
		Query[i].x2--;
		Query[i].y2--;
	}
	
	for(int i = 0; i < n; i++)
	{
		PLANT& P = Plant[i];
		int d = (P.a - 1) / P.b;
		int x1 = max(0, P.x - d);
		int y1 = max(0, P.y - d);
		int x2 = min(w, P.x + d + 1);
		int y2 = min(h, P.y + d + 1);
		//势力范围
		for(int x = x1; x < x2; x++)
		{
			for(int y = y1; y < y2; y++)
			{
				int d = dist(x, y, P.x, P.y);
				Nuclearia(x, y) += P.a - (d * P.b);
			}
		}
	}
	
	Summarize();
	
	for(int i = 0; i < q; i++)
	{
		QUERY& Q = Query[i];
		huge rsl = Nuclearia(Q.x2, Q.y2);
		if(Q.x1 >= 0) rsl -= Nuclearia(Q.x1, Q.y2);
		if(Q.y1 >= 0) rsl -= Nuclearia(Q.x2, Q.y1);
		if((Q.x1 >= 0) && (Q.y1 >= 0)) rsl += Nuclearia(Q.x1, Q.y1);
		
		Print(rsl, (Q.x2 - Q.x1) * (Q.y2 - Q.y1));
	}
	
	return 0;
}

正解

接着，我们尝试在进一步的优化
回想一下我们在利用树状数组做区间修改时的方法：把区间的前面标记为 + a +a +a，最后面标记为 − a -a −a，这样，就可以进行区间修改了
我们可以借用一下这个思路，来标记一次核电站爆炸事件
首先，我们先找出这个核电站爆炸事件的“势力范围”，接着，把“势力范围”的左上角和右下角标记为 + a % b +a\% b +a%b，把左下角和右上角标记为 − a % b -a\% b −a%b
然后，我们开两个数组，在这个“势力范围”的对角线上（两条对角线，所以要两个数组）存储 b b b的值，为了节省时间，我们可以只在4个角上标记，最后再把对角线上的数全部加出来
因为每个核电站爆炸事件的影响相等方格的，都会围成一个正方形（因为取的是切比雪夫距离），也就是说，我们在统计时，只需要先把对角线上的数复制到 N u c l e a r i a Nuclearia Nuclearia数组中，再把左边的和上面的数相加再减去左上的数就可以了
看起来有点抽象，我们可以再总结一下：

  1. 对于每个核电站爆炸事件，计算出它的“势力范围”，左上角为 ( x 1 , y 1 ) (x1,y1) (x1,y1)，右下角为 ( x 2 , y 2 ) (x2,y2) (x2,y2)
  2. 标记： N u c l e a r i a ( x 1 , y 1 ) = N u c l e a r i a ( x 2 , y 2 ) = a % b Nuclearia(x1,y1) = Nuclearia(x2,y2) = a\% b Nuclearia(x1,y1)=Nuclearia(x2,y2)=a%b， N u c l e a r i a ( x 1 , y 2 ) = N u c l e a r i a ( x 2 , y 1 ) = a % b Nuclearia(x1,y2) = Nuclearia(x2,y1) = a\% b Nuclearia(x1,y2)=Nuclearia(x2,y1)=a%b
  3. 修改两条对角线上的值（两条对角线的数组分别为 P o s D i a g PosDiag PosDiag（主对角线）和 N e g D i a g NegDiag NegDiag（次对角线））： P o s D i a g ( x 1 + 1 , y 1 + 1 ) = b , P o s D i a g ( x 2 , y 2 ) = − b , N e g D i a g ( x 1 + 1 , y 2 − 1 ) = − b , N e g D i a g ( x 2 , y 1 ) = b PosDiag(x1+1,y1+1)=b,PosDiag(x2,y2)=-b,NegDiag(x1+1,y2-1)=-b,NegDiag(x2,y1)=b PosDiag(x1+1,y1+1)=b,PosDiag(x2,y2)=−b,NegDiag(x1+1,y2−1)=−b,NegDiag(x2,y1)=b
  4. 枚举整个国度的 P o s D i a g PosDiag PosDiag和 N e g D i a g NegDiag NegDiag， P o s D i a g ( x , y ) + = P o s D i a g ( x − 1 , y − 1 ) , N e g D i a g ( x , y ) + = N e g D i a g ( x − 1 , y + 1 ) PosDiag(x,y)+=PosDiag(x-1,y-1),NegDiag(x,y)+=NegDiag(x-1,y+1) PosDiag(x,y)+=PosDiag(x−1,y−1),NegDiag(x,y)+=NegDiag(x−1,y+1)
  5. 枚举整个国度的 N u c l e a r i a Nuclearia Nuclearia， N u c l e a r i a ( x , y ) + = P o s D i a g ( x , y ) + N e g D i a g ( x , y ) Nuclearia(x,y)+=PosDiag(x,y)+NegDiag(x,y) Nuclearia(x,y)+=PosDiag(x,y)+NegDiag(x,y)
  6. 重复两次：枚举整个国度的 N u c l e a r i a Nuclearia Nuclearia， N u c l e a r i a ( x , y ) + = N u c l e a r i a ( x − 1 , y ) + N u c l e a r i a ( x , y − 1 ) − N u c l e a r i a ( x − 1 , y − 1 ) Nuclearia(x,y)+=Nuclearia(x-1,y)+Nuclearia(x,y-1)-Nuclearia(x-1,y-1) Nuclearia(x,y)+=Nuclearia(x−1,y)+Nuclearia(x,y−1)−Nuclearia(x−1,y−1)
    这就是整个过程了，最后的答案计算方法和上面一样
    你以为这就完了吗？
    不！没完！
    注意题目的最后一句话：如果核电站位于核能国的边境或是在离边境稍近的位置，那么爆炸可能也会影响到核能国之外的方格。影响到核能国外方格的爆炸被称作界限
    也就是说 x 1 x1 x1和 y 1 y1 y1有可能 ⩽ 0 \leqslant 0 ⩽0， x 2 x2 x2也有可能 > W >W >W， y 2 y2 y2也有可能 > H >H >H！
    对于这个的处理，我们需要再开两个数组 C o l Col Col和 R o w Row Row，存储超出边界的部分，具体的实现就看程序吧，这些细节的处理，这里就不再细讲了
//官方程序
#include <stdio.h>
#include <stdlib.h>

#define MAXWH 10000000
#define MAXN 1000000
#define MAXQ 1000000

typedef long long int huge;

int w, h, n, q;

typedef struct NUCLEARIA
{
	huge Info[MAXWH];
	huge& operator () (int x, int y) {return Info[(y * w) + x];}
}
NUCLEARIA;

typedef struct PLANT
{
	int x;
	int y;
	int a;
	int b;
}
PLANT;

typedef struct QUERY
{
	int x1;
	int y1;
	int x2;
	int y2;
}
QUERY;

NUCLEARIA Nuclearia, PosDiag, NegDiag;
PLANT Plant[MAXN];
QUERY Query[MAXQ];
huge Row[MAXWH];
huge Col[MAXWH];

int max(int a, int b)
{
	return (a > b) ? a : b;
}

void UpdatePosDiag(int x1, int y1, int x2, int y2, int b)
{
	if((x1 < 0) && (y1 < 0))
	{
		int m = -max(x1, y1);
		Nuclearia(0, 0) += m * b;
		x1 += m;
		y1 += m;
	}
	
	if(x1 < 0)
	{
		int m = -x1;
		Col[y1] += b;
		Col[y1 + m] -= b;
		x1 += m;
		y1 += m;
	}
	
	if(y1 < 0)
	{
		int m = -y1;
		Row[x1] += b;
		Row[x1 + m] -= b;
		x1 += m;
		y1 += m;
	}
	
	PosDiag(x1, y1) += b;
	
	if((x2 + 1 < w) && (y2 + 1 < h)) PosDiag(x2 + 1, y2 + 1) -= b;
}

void UpdateNegDiag(int x1, int y1, int x2, int y2, int b)
{
	if(y2 > h - 1)
	{
		int m = y2 - (h - 1);
		x1 += m;
		y2 -= m;
	}
	
	if(x1 < 0)
	{
		int m = -x1;
		Col[(y2 - m) + 1] -= b;
		if(y2 + 1 < h) Col[y2 + 1] += b;
		x1 += m;
		y2 -= m;
	}
	
	if((x1 < w) && (y2 >= 0)) NegDiag(x1, y2) -= b;
	
	if(x2 > w - 1)
	{
		int m = x2 - (w - 1);
		x2 -= m;
		y1 += m;
	}
	
	if(y1 < 0)
	{
		int m = -y1;
		Row[(x2 - m) + 1] -= b;
		if(x2 + 1 < w) Row[x2 + 1] += b;
		x2 -= m;
		y1 += m;
	}
	
	if((x2 + 1 >= 0) && (x2 + 1 < w) && (y1 - 1 >= 0) && (y1 - 1 < h)) NegDiag(x2 + 1, y1 - 1) += b;
}

void SummarizeDiags()
{
	for(int x = 0; x < w; x++)
	{
		for(int y = 0; y < h; y++)
		{
			if((x) && (y)) PosDiag(x, y) += PosDiag(x - 1, y - 1);
			if((x) && (y != h - 1)) NegDiag(x, y) += NegDiag(x - 1, y + 1);
		}
	}
}

void AddDiags()
{
	for(int x = 0; x < w; x++)
	{
		for(int y = 0; y < h; y++)
		{
			Nuclearia(x, y) += PosDiag(x, y);
			Nuclearia(x, y) += NegDiag(x, y);
		}
	}
}

void SummarizeLines()
{
	for(int x = 1; x < w; x++)
	{
		Row[x] += Row[x - 1];
	}
	
	for(int y = 1; y < h; y++)
	{
		Col[y] += Col[y - 1];
	}
}

void AddLines()
{
	for(int x = 0; x < w; x++)
	{
		Nuclearia(x, 0) += Row[x];
	}
	
	for(int y = 0; y < h; y++)
	{
		Nuclearia(0, y) += Col[y];
	}
}

void Summarize()
{
	for(int x = 0; x < w; x++)
	{
		for(int y = 0; y < h; y++)
		{
			if(x) Nuclearia(x, y) += Nuclearia(x - 1, y);
			if(y) Nuclearia(x, y) += Nuclearia(x, y - 1);
			if((x) && (y)) Nuclearia(x, y) -= Nuclearia(x - 1, y - 1);
		}
	}
}

void Print(huge sum, int area)
{
	huge rsl = sum / area;
	if((sum % area) * 2 >= area)
	{
		rsl++;
	}
	
	printf("%lld\n", rsl);
}

int main()
{
	scanf("%d%d", &w, &h);
	
	scanf("%d", &n);
	for(int i = 0; i < n; i++)
	{
		scanf("%d%d%d%d", &(Plant[i].x), &(Plant[i].y), &(Plant[i].a), &(Plant[i].b));
		Plant[i].x--;
		Plant[i].y--;
	}
	
	scanf("%d", &q);
	for(int i = 0; i < q; i++)
	{
		scanf("%d%d%d%d", &(Query[i].x1), &(Query[i].y1), &(Query[i].x2), &(Query[i].y2));
		Query[i].x1 -= 2;
		Query[i].y1 -= 2;
		Query[i].x2--;
		Query[i].y2--;
	}
	
	for(int i = 0; i < n; i++)
	{
		PLANT& P = Plant[i];
		int d = (P.a - 1) / P.b;
		int x1 = P.x - d;
		int y1 = P.y - d;
		int x2 = P.x + d + 1;
		int y2 = P.y + d + 1;
		int r = P.a % P.b;
		
		if(r)
		{
			Nuclearia(max(0, x1), max(0, y1)) += r;
			if(x2 < w) Nuclearia(x2, max(0, y1)) -= r;
			if(y2 < h) Nuclearia(max(0, x1), y2) -= r;
			if((x2 < w) && (y2 < h)) Nuclearia(x2, y2) += r;
			
			x1++;
			y1++;
			x2--;
			y2--;
		}
		
		if(P.a >= P.b)
		{
			UpdatePosDiag(x1, y1, x2, y2, P.b);
			UpdateNegDiag(x1, y1, x2, y2, P.b);
		}
	}
	
	SummarizeDiags();
	AddDiags();
	SummarizeLines();
	AddLines();
	Summarize();
	Summarize();
	
	for(int i = 0; i < q; i++)
	{
		QUERY& Q = Query[i];
		huge rsl = Nuclearia(Q.x2, Q.y2);
		if(Q.x1 >= 0) rsl -= Nuclearia(Q.x1, Q.y2);
		if(Q.y1 >= 0) rsl -= Nuclearia(Q.x2, Q.y1);
		if((Q.x1 >= 0) && (Q.y1 >= 0)) rsl += Nuclearia(Q.x1, Q.y1);
		
		Print(rsl, (Q.x2 - Q.x1) * (Q.y2 - Q.y1));
	}
	
	return 0;
}
jacky0705
关注
专栏目录
【CEOI2015 Day2】世界冰球锦标赛（ 提高+/省选-） 题解
qq_35684989的博客
06-01 1467
给出 Bobek 的预算和每场比赛的票价，试求：如果总票价不超过预算，他有多少种观赛方案。如果存在以其中一种方案观看某场比赛而另一种方案不观看，则认为这两种方案不同。
BZOJ 2790「POI2012」Distance
Yucohny
10-19 214
#include &lt;cmath&gt; #include &lt;vector&gt; #include &lt;cstdio&gt; #include &lt;cstring&gt; #include &lt;iostream&gt; #include &lt;algorithm&gt; #define ll long long #define rel(i,x,y) for(ll i=(x...
【CF1025B】 Weakened Common Divisor
Brute♂force
12-19 322
题目 题目描述 与 GCDGCD （最大公约数）类似，我们引进 WCDWCD （弱公约数）的概念， WCDWCD 的定义如下： 给出几对数 \left( a_1,b_1 \right) ,\left( a_2,b_2 \right) ,\cdots ,\left( a_n,b_n \right)(a 1 ​ ,b 1 ​ ),(a 2 ​ ,b 2 ​ ),⋯,(a n ​ ,b n ​...
Codeforces Ozon Tech Challenge 2020 (Div1 + Div2) A ~ D
OlaMins
03-09 2208
A. Kuroni and the Gifts time limit per test 1 second memory limit per test 256 megabytes input standard input output ...
[CEOI2015]Calvinball championship
OZY的博客
10-31 377
CEOI
loj#2788-「CEOI2015 Day1」管道【树上差分】
QuantAsk
08-11 264
但是我们不能离线建生成树，因为我们存不下所有的边，考虑一下别的方向的优化。我们会发现对于非树边来说，如果这一条非树边能被其他非树边完全覆盖，那么说明这条边就没有用，所以我们对于非树边来说也只需要保留一棵最小生成树即可。一个朴素的想法是我们搞出一棵生成树来，然后对于非树边。路径上的边都标记成非割边，然后剩下的就是割边了。然后至于标记方面用树上差分来处理就好了。条边的一张图，求它的所有割边。我们存不下所有的边，但是。...
Loj#3005-「JOISC 2015 Day 4」Limited Memory【交互题】
QuantAsk
06-15 262
题目链接:https://loj.ac/p/3005有一个长度为nnn的括号串SSS，其中包括[]和两种括号类型，一个合法的括号串要求同类型的括号一一对应。你每次可以询问SSS中的一个字符并且传递一个[0,222)[0,2^{22})[0,222)的数字到下一次。你的程序每次只知道字符串长度nnn和上一次传递过来的数字（最开始时是000），你需要在150001500015000次内得到这个字符串是否为一个合法括号串。1≤n≤1001\leq n\leq 1001≤n≤100很有趣的题目。我们需要知道的信..
[LOJ6507]「雅礼集训 2018 Day7」A
时光真疯狂, 我一路执迷于匆忙.
11-03 1234
Description 给出一个长度为n的序列a，要求资瓷 区间或，区间与和求区间最小值 n&lt;=5e5,ai&lt;2^31 Solution 让我们先来考虑一个暴力： 如果一次操作对某个区间的影响是一样的（即最小值还是那个数），那么我们直接打标记退出；若一次操作对某个区间无影响就直接退出 只需要维护区间or和区间and 然后你就会发现你过了_ (:з」∠) _ 这是为什么呢？ 让我们来分析...
LOJ6435 & 洛谷5465 「PKUSC2018」星际穿越 倍增
The_OIer的博客
09-18 344
题目链接： loj6435 洛谷5465 蒟蒻zyd：这不是大水题吗？看我写个O(nlog2n)O(nlog^2n)O(nlog2n)的诡异ST表卡卡常数跑过去 题目：输出区间距离和 蒟蒻zyd：（笑容逐渐消失） 没想到一道倍增题能这么巧（毒）妙（瘤）…… 一些奇奇怪怪的性质 这里是需要用到的性质…… 为了方便，这里把“花费1单位时间进行传送”称为“走了1步”qwq Ps.一些类似的情况就不画图了...
LOJ6354 & 洛谷4366：[Code+#4]最短路——题解
最新发布
05-17
2. 创建一个集合 visited，用于存储已经确定最短路径的节点。 3. 将源点加入 visited 集合中。 4. 对于源点周围的每个节点，更新它们到源点的距离。如果距离更短，则更新 dist 数组。 5. 从 dist 数组中选择距离...
LOJ #3163. 「CEOI2019」动态直径（欧拉序+线段树）
Freopen的博客
10-24 621
重点在于实现。 #3163. 「CEOI2019」动态直径 #include<bits/stdc++.h> #define maxn 200005 #define pii pair<int,int> #define mp make_pair #define lc u<<1 #define rc u<<1|1 #define LL long long ...
【LOJ3276】「JOISC 2020 Day2」遗迹
cz_xuyixuan的博客
03-24 1060
题目链接 点击打开链接 题目解法 首先考虑对于确定的 hih_ihi​ ，判断其是否满足条件。显然，选出 AAA 集合的过程如下： 取出 hih_ihi​ 的两个最大值的下标，加入集合 SSS ，弹出 SSS 的最大值，加入集合 AAA ，重复 NNN 次。 也有这样一个等价的过程：从后往前考虑各个 hih_ihi​ ，若在最终状态中，仍然存在 1≤x≤hi1\leq x\leq h_i1≤x≤h...
【LOJ3275】「JOISC 2020 Day2」有趣的 Joitter 交友
cz_xuyixuan的博客
03-24 656
题目链接 点击打开链接 题目解法 问题可以转述为如下形式：在一张会自行补边的有向图上不断加边，若 xxx 连向了 yyy ，且 y,zy,zy,z 在一个二元环内， xxx 也会连向 zzz ，每次加入一条边，求出当前边数。 显然由二元环连接的点集中每一条可能的边都存在。 考虑将由二元环连接的点集缩点，则为了计算答案，需要维护指向该点集的入点集合 ineineine 。 由此，一个点集 SSS 对...
洛谷P1172 安全逃离题解
菜鸡Jacky0705的CSDN博客
12-05 1553
洛谷【P1172】安全逃离题解题目 题目 农夫john最近在研究如果发生重大事故，如何让农场里的奶牛逃离问题。他想要确信在紧急情况下，所有的奶牛都有一个安全逃离方案。因为在紧急情况下，奶牛们都会失去观察和判断能力，所以最近john一直在教奶牛们逃离的方法，他的方法很简单，就是任何时候都只向北方或东方逃离，北方是行坐标减1的方向，东方是列坐标加1的方向。奶牛们虽笨，不过这一点事关自己的生命，所以他们...
洛谷P7916 [CSP-S 2021]交通规划
菜鸡Jacky0705的CSDN博客
11-15 1337
洛谷P7916 [CSP-S 2021]交通规划Part1 k⩽2k\leqslant 2k⩽2的情况Step1 分析Step2 计算分界线上边权之和 如果你没有学习过对偶图等知识，那么对你来说其他题解可能有点难理解，但这篇题解不会！在阅读这篇题解之前，你不需要掌握任何高深的知识，就可以学会这道题的做法 Part1 k⩽2k\leqslant 2k⩽2的情况 Step1 分析 首先，k=1k=1k=1的情况非常简单，你只需要把所有点染成相同颜色即可，这时候答案为000，显然是最小的 而对于k=2k=2k=2
DTOJ2229 抢车位
菜鸡Jacky0705的CSDN博客
03-26 910
DTOJ2229 抢车位题目题目描述输入格式输出格式样例样例输入1样例输出1样例输入2样例输出2样例输入3样例输出3数据范围与提示样例1解释数据范围题解 题目 题目描述 很久以前 ，cxm做了一个题，叫“抢车位”，大意是让你调度的汽车 使得每个汽车都有位。AC以后，cxm去实地考察了这个游戏 ，发现最有意思的地方是“以旧换新”：你最多拥有101010辆汽车， 便宜的汽车换贵只用补差价 但是贵的汽车...
DTOJ3489 可怜与超市（supermarket）
菜鸡Jacky0705的CSDN博客
04-30 830
DTOJ3489 可怜与超市（supermarket）题解 题解 这道题主要的难点在于处理“使用第iii张优惠券时必须先使用第xix_ixi​张优惠券”这个问题 这非常容易让我们联想到树，我们把iii向xix_ixi​连一条边，就可以构成一棵树 然后就很容易想到可以使用树形DP了 我们用fi,j,kf_{i,j,k}fi,j,k​表示第iii个节点，买了jjj个，用(k=1k=1k=1)还是没用(...
DTOJ1984 花园
菜鸡Jacky0705的CSDN博客
03-23 596
DTOJ1984 花园题目题目描述输入格式输出格式样例样例输入样例输出数据范围与提示样例解释数据规模与约定题解 题目 题目描述 奇怪的大学有一座奇怪的花园，花园由N座温室组成，温室依次标号为1,2,⋯⋯ ,N1,2,\cdots \cdots ,N1,2,⋯⋯,N，温室之间由N−1N-1N−1条双向道路连接 每一座温室都种植这一种花，随着季节的变换，温室里的花的种类也在不断发生着变化 ShenX平...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值