题目大意:在平面直角坐标系的第一象限中,给出n个矩形(可能重叠)。有m次询问,每次询问点(t,t)的左下方的正方形区域中矩形的总面积(重叠部分重叠几次就得统计几次)。

题目分析:线段树的叶子节点x维护的是点(x,x)左下方的正方形区域中矩形的总面积。每添加一个矩形,便做一次区间更新。

代码如下:

# include<bits/stdc++.h>
using namespace std;
# define LL long long
# define mid (l+(r-l)/2) const int N=200000; LL lazy[N*4+1];
LL lazy_a0[N*4+1],lazy_d[N*4+1];
LL lazy_a[N*4+1],lazy_b[N*4+1],lazy_c[N*4+1]; void read(int &x)
{
x=0;
char c;
while((c=getchar())&&(c<'0'||c>'9'));
x=c-'0';
while(c=getchar()){
if(c<'0'||c>'9') break;
x=x*10+c-'0';
}
} void pushDown(int rt,int l,int r)
{
lazy[rt<<1]+=lazy[rt];
lazy[rt<<1|1]+=lazy[rt];
lazy[rt]=0; lazy_a0[rt<<1]+=lazy_a0[rt];
lazy_a0[rt<<1|1]+=lazy_a0[rt]+(LL)(mid-l+1)*lazy_d[rt];
lazy_d[rt<<1]+=lazy_d[rt];
lazy_d[rt<<1|1]+=lazy_d[rt];
lazy_a0[rt]=lazy_d[rt]=0; lazy_a[rt<<1]+=lazy_a[rt];
lazy_a[rt<<1|1]+=lazy_a[rt];
lazy_a[rt]=0; lazy_b[rt<<1]+=lazy_b[rt];
lazy_b[rt<<1|1]+=lazy_b[rt];
lazy_b[rt]=0; lazy_c[rt<<1]+=lazy_c[rt];
lazy_c[rt<<1|1]+=lazy_c[rt];
lazy_c[rt]=0;
} void build(int rt,int l,int r)
{
lazy[rt]=0;
lazy_a0[rt]=lazy_d[rt]=0;
lazy_a[rt]=lazy_b[rt]=lazy_c[rt]=0;
if(l==r) return ;
build(rt<<1,l,mid);
build(rt<<1|1,mid+1,r);
} void update1(int rt,int l,int r,int L,int R,LL x)
{
if(L<=l&&r<=R){
lazy[rt]+=x;
}else{
if(L<=mid) update1(rt<<1,l,mid,L,R,x);
if(R>mid) update1(rt<<1|1,mid+1,r,L,R,x);
}
} void update2(int rt,int l,int r,int L,int R,LL a0,int d)
{
if(L<=l&&r<=R){
LL a1=(LL)a0+(LL)d*(LL)(l-L);
lazy_a0[rt]+=(LL)a1;
lazy_d[rt]+=(LL)d;
}else{
if(L<=mid) update2(rt<<1,l,mid,L,R,a0,d);
if(R>mid) update2(rt<<1|1,mid+1,r,L,R,a0,d);
}
} void update3(int rt,int l,int r,int L,int R,int lowx,int lowy)
{
if(L<=l&&r<=R){
++lazy_a[rt];
lazy_b[rt]-=(LL)(lowx+lowy);
lazy_c[rt]+=(LL)lowx*(LL)lowy;
}else{
if(L<=mid) update3(rt<<1,l,mid,L,R,lowx,lowy);
if(R>mid) update3(rt<<1|1,mid+1,r,L,R,lowx,lowy);
}
} LL query(int rt,int l,int r,int x)
{
if(l==r){
LL ans=0;
ans+=lazy[rt];
LL last=lazy_a0[rt]+(r-l)*lazy_d[rt];
ans+=(r-l+1)*(lazy_a0[rt]+last)/2;
ans+=lazy_a[rt]*l*l+lazy_b[rt]*l+lazy_c[rt];
return ans;
}else{
pushDown(rt,l,r);
if(x<=mid) return query(rt<<1,l,mid,x);
return query(rt<<1|1,mid+1,r,x);
}
} int main()
{
int T,n;
read(T);
while(T--)
{
build(1,0,N);
scanf("%d",&n);
int x1,x2,y1,y2;
while(n--){
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
int dltx=x2-x1;
int dlty=y2-y1;
update1(1,0,N,max(x2,y2)+1,N,(LL)dltx*(LL)dlty);
if(y1>=x2)
update2(1,0,N,y1,y2,0,dltx);
else if(y2<=x1)
update2(1,0,N,x1,x2,0,dlty);
else if(y2>x1&&y1<x2){
if(y2<x2&&y1<=x1){
update3(1,0,N,x1,y2-1,x1,y1);
update2(1,0,N,y2,x2,(LL)dlty*(y2-x1),dlty);
}else if(y2>=x2&&y1<=x1){
update3(1,0,N,x1,x2-1,x1,y1);
update2(1,0,N,x2,y2,(LL)dltx*(LL)(x2-y1),dltx);
}else if(y2>=x2&&y1>x1){
update3(1,0,N,y1,x2-1,x1,y1);
update2(1,0,N,x2,y2,(LL)dltx*(x2-y1),dltx);
}else if(y2<x2&&y1>x1){
update3(1,0,N,y1,y2-1,x1,y1);
update2(1,0,N,y2,x2,(LL)dlty*(LL)(y2-x1),dlty);
}
}
}
scanf("%d",&n);
while(n--)
{
scanf("%d",&x1);
printf("%lld\n",query(1,0,N,x1));
}
}
return 0;
}

  

HDU-4533 威威猫系列故事——晒被子(区间更新)的更多相关文章

  1. hdu4533 威威猫系列故事——晒被子

    Problem Description 因为马拉松初赛中吃鸡腿的题目让不少人抱憾而归,威威猫一直觉得愧对大家,这几天他悄悄搬到直角坐标系里去住了. 生活还要继续,太阳也照常升起,今天,威威猫在第一象限 ...

  2. HDU 4533 威威猫系列故事——晒被子

    题目链接 扫描线可做,然后当时比赛后问虎哥,他说可以标记,然后拖了很久,今天从早上折腾到晚上,终于把两种情况写出来,分析太弱.改天扫描线,再来一次. 被子如果被y = x 穿过,可以分成两部分,上和下 ...

  3. HDU 4540 威威猫系列故事——打地鼠

    威威猫系列故事--打地鼠 Time Limit: 300/100 MS (Java/Others) Memory Limit: 65535/32768 K (Java/Others) Total Su ...

  4. hdu 4523 威威猫系列故事——过生日 小模拟

    威威猫系列故事——过生日 Time Limit: 500/200 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total ...

  5. hdu 4540 威威猫系列故事——打地鼠 dp小水题

    威威猫系列故事——打地鼠 Time Limit: 300/100 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total ...

  6. HDU 4540 威威猫系列故事——打地鼠 (状态压缩DP)

    威威猫系列故事——打地鼠 Time Limit: 300/100 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total ...

  7. HDU--杭电--4504--威威猫系列故事——篮球梦--DP

    威威猫系列故事——篮球梦 Time Limit: 300/100 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others) Total ...

  8. HDUOJ----4504 威威猫系列故事——篮球梦

    威威猫系列故事——篮球梦 Time Limit: 300/100 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total ...

  9. hdoj 4526 威威猫系列故事——拼车记

    威威猫系列故事——拼车记 Time Limit: 500/200 MS (Java/Others)    Memory Limit: 65535/32768 K (Java/Others)Total ...

随机推荐

  1. [转]source inslght使用指导

    作为一个开放源代码的操作系统,Linux附带的源代码库使得广大爱好者有了一个广泛学习.深入钻研的机会,特别是Linux内核的组织极为复杂,同时,又不能像windows平台的程序一样,可以使用集成开发环 ...

  2. asp.net页面使用doPostBack的后台取值

    前台页面(aspx文件): --伪装按钮 <span onclick='__doPostBack("lkSend","key")'>发送</s ...

  3. 有哪些 PHP 调试技巧?

    我目前遇到的最让我称赞的debug方式是:xdebug的 xdebug_start_trace(); /* 业务代码 */ xdebug_stop_trace(); 他解决了我长久以来一个代码调试问题 ...

  4. NSString的几种常用方法

    NSString的几种常用方法   要把 “2011-11-29” 改写成 “2011/11/29”一开始想用ios的时间格式,后来用NSString的方法搞定. [string stringByRe ...

  5. Python清理内存中的密码

    基本不太好搞.可以参考如下讨论: http://stackoverflow.com/questions/728164/securely-erasing-password-in-memory-pytho ...

  6. 从决策树学习谈到贝叶斯分类算法、EM、HMM --别人的,拷来看看

    从决策树学习谈到贝叶斯分类算法.EM.HMM     引言 最近在面试中,除了基础 &  算法 & 项目之外,经常被问到或被要求介绍和描述下自己所知道的几种分类或聚类算法(当然,这完全 ...

  7. 北大poj-1001

    Description Problems involving the computation of exact values of very large magnitude and precision ...

  8. hdu 2097

    ps:WA了两次好像....Sky数是三个进制下的各位数之和相等...而不是都等于22...我傻逼了... 代码: #include "stdio.h" int inp(int a ...

  9. php大力力 [032节] php设计时候遇见麻烦:XQB50-H8268 进水电磁阀

    海信洗衣机 无法进水,刚才写程序,洗衣机不进水,在叫唤,去看了看,上网查了查,估计是进水电磁阀坏了. 打算自己拆了查出型号,淘宝买,自己修. 想起以前洗衣机坏了,找人修,对方报价好几百,淘宝看洗衣机主 ...

  10. 在android的spinner中,实现取VALUE值和TEXT值。 ZT

    在android的spinner中,实现取VALUE值和TEXT值.   为了实现在android的 spinner实现取VALUE值和TEXT值,我尝试过好些办法,在网上查的资料,都是说修改适配器, ...