首页 > 代码库 > 【bzoj1176】[Balkan2007]Mokia/【bzoj2683】简单题
【bzoj1176】[Balkan2007]Mokia/【bzoj2683】简单题
bzoj1176
题目描述
维护一个W*W的矩阵,初始值均为S(题目描述有误,这里的S没有任何作用!).每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.
输入
第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小
接下来每行为一下三种输入之一(不包含引号):
1 x y a
2 x1 y1 x2 y2
3
输入1:你需要把(x,y)(第x行第y列)的格子权值增加a
输入2:你需要求出以左下角为(x1,y1),右上角为(x2,y2)的矩阵内所有格子的权值和,并输出
输入3:表示输入结束
输出
对于每个输入2,输出一行,即输入2的答案
样例输入
0 4
1 2 3 3
2 1 1 3 3
1 2 2 2
2 2 2 3 4
3
样例输出
3
5
bzoj2683
题目描述
同上,只是少了一个烦人的S
题解
CDQ分治
先读入所有修改和询问操作,将每个修改拆成4个。
然后按照CDQ分治的套路,按x排序,按id比较,按y查询。
(好像这道题我写的和别的CDQ分治题不一样。。。)
代码是1176的,2683直接去掉%*d即可。
(就算S真的有用,可以直接改ans的初始值,也是一样的)
#include <cstdio> #include <algorithm> using namespace std; struct data { int x , y , v , p , id , opt; }a[800000] , tmp[800000]; int n , tot , cnt , f[2000010] , ans[200000]; bool cmp(data a , data b) { return a.x == b.x ? (a.y == b.y ? a.opt < b.opt : a.y < b.y) : a.x < b.x; } void add(int b , int c , int d , int t) { a[++tot].x = b , a[tot].y = c , a[tot].v = d , a[tot].id = tot , a[tot].opt = t; if(t) a[tot].p = cnt; } void update(int x , int a) { int i; for(i = x ; i <= n ; i += i & -i) f[i] += a; } int query(int x) { int i , ans = 0; for(i = x ; i ; i -= i & -i) ans += f[i]; return ans; } void solve(int l , int r) { if(l == r) return; int mid = (l + r) >> 1 , i , p1 = l , p2 = mid + 1; for(i = l ; i <= r ; i ++ ) { if(a[i].id <= mid && !a[i].opt) update(a[i].y , a[i].v); if(a[i].id > mid && a[i].opt) ans[a[i].p] += a[i].v * query(a[i].y); } for(i = l ; i <= r ; i ++ ) if(a[i].id <= mid && !a[i].opt) update(a[i].y , -a[i].v); for(i = l ; i <= r ; i ++ ) { if(a[i].id <= mid) tmp[p1 ++ ] = a[i]; else tmp[p2 ++ ] = a[i]; } for(i = l ; i <= r ; i ++ ) a[i] = tmp[i]; solve(l , mid) , solve(mid + 1 , r); } int main() { int i , k , b , c , d , e; scanf("%*d%d" , &n); while(scanf("%d" , &k) && k != 3) { if(k == 1) scanf("%d%d%d" , &b , &c , &d) , add(b , c , d , 0); else scanf("%d%d%d%d" , &b , &c , &d , &e) , cnt ++ , add(d , e , 1 , 1) , add(b - 1 , e , -1 , 1) , add(d , c - 1 , -1 , 1) , add(b - 1 , c - 1 , 1 , 1); } sort(a + 1 , a + tot + 1 , cmp); solve(1 , tot); for(i = 1 ; i <= cnt ; i ++ ) printf("%d\n" , ans[i]); return 0; }
【bzoj1176】[Balkan2007]Mokia/【bzoj2683】简单题
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。