Problem Description
因为马拉松初赛中吃鸡腿的题目让不少人抱憾而归,威威猫一直觉得愧对大家,这几天他悄悄搬到直角坐标系里去住了。
生活还要继续,太阳也照常升起,今天,威威猫在第一象限晒了N条矩形的被子,被子的每条边都和坐标轴平行,不同被子的某些部分可能会叠在一起。这时候,在原点处突然发了场洪水,时间t的时候,洪水会蔓延到( t, t ),即左下角为( 0, 0 ) ,右上角为( t, t )的矩形内都有水。
悲剧的威威猫想知道,在时间t1, t2, t3 ... tx 的时候,他有多少面积的被子是湿的?
Input
输入数据首先包含一个正整数T,表示有T组测试数据;
每组数据的第一行首先是一个整数N,表示有N条被子;
接下来N行,每行包含四个整数x1, y1, x2, y2,代表一条被子的左下角和右上角的坐标;
然后接下来一行输入一个整数x,表示有x次询问;
再接下来x行,输入x个严格单调递增的整数,每行一个,表示威威猫想知道的时间ti。
[Technical Specification]
T <= 5
0 < N <= 20000
1 <= x1 < x2 <= 200000
1 <= y1 < y2 <= 200000
1 <= x <= 20000
1 <= ti <= 200000 (1 <= i <= x )
Output
对于每次询问,请计算并输出ti时有多少面积的被子是湿的,每个输出占一行。
Sample Input
1 2 1 1 3 3 2 2 4 4 5 1 2 3 4 5
Sample Output
0 1 5 8 8 这题里求的面积是各个矩形的面积和,并不是矩形面积的并,即相互不影响,因为对于各个时段,矩形的面积都可以用A*t^2+B*t+C表示,所以我们可以用树状数组或者线段树分别维护A,B,C,每次读入一个矩形,就记录它对不同时间段t的影响,然后把相应的系数加到线段树中。 当一个矩形读入时(左下角坐标(x1,y1),右上角坐标(x2,y2)),有四种影响情况: 1.0~max(x1,y1) 这个时期因为t时间形成的矩形在该矩形下面,对面积没有影响,所以不用考虑。 2.如果存在t时间形成的矩形部分覆盖该矩形,但没有碰到或者超过该矩形的上边或者右边,那么要同时更新A,B,C,表达式为(t-x1)*(t-y1)。这里的判断表达式为如果(max(x1,y1)<min(x2,y2)),则更新max(x1,y1)~min(x2,y2)。 3.t时间内形成的矩形恰好碰到或者越过右边或上边的一条,即这段时间内的矩形面积变化中(x2-x1)或(y2-y1)为常数,那么要更新,注意,这里更新的时候要分类讨论,具体看代码。 4.max(x2,y2)~maxn,这个时期因为t时间形成的矩形已经完全覆盖矩形,所以要更新C,加上面积常数。 写线段树的时候注意,三个变量要一起存,不能开三个线段树,会超内存的。 线段树代码:#include<iostream> #include<stdio.h> #include<stdlib.h> #include<string.h> #include<math.h> #include<vector> #include<map> #include<set> #include<queue> #include<stack> #include<string> #include<algorithm> using namespace std; #define ll __int64 #define maxn 200050 ll A,B,C; ll a[10]; struct node{ ll l,r,cnt1,cnt2,cnt3,flag; }b[4*maxn]; void build(ll l,ll r,ll i) { ll mid; b[i].l=l;b[i].r=r;b[i].cnt1=b[i].cnt2=b[i].cnt3=0;b[i].flag=1; if(l==r)return; mid=(l+r)/2; build(l,mid,i*2); build(mid+1,r,i*2+1); } void update(ll l,ll r,ll a[],ll i) { ll mid; if(b[i].l==l && b[i].r==r){ b[i].cnt1+=a[1];b[i].cnt2+=a[2];b[i].cnt3+=a[3];return; } mid=(b[i].l+b[i].r)/2; if(r<=mid)update(l,r,a,i*2); else if(l>mid)update(l,r,a,i*2+1); else{ update(l,mid,a,i*2); update(mid+1,r,a,i*2+1); } } void question(ll id,ll i) { ll mid; if(b[i].l==id && b[i].r==id){ A=b[i].cnt1;B=b[i].cnt2;C=b[i].cnt3;return; } if(b[i].cnt1){ b[i*2].cnt1+=b[i].cnt1; b[i*2+1].cnt1+=b[i].cnt1; b[i].cnt1=0; } if(b[i].cnt2){ b[i*2].cnt2+=b[i].cnt2; b[i*2+1].cnt2+=b[i].cnt2; b[i].cnt2=0; } if(b[i].cnt3){ b[i*2].cnt3+=b[i].cnt3; b[i*2+1].cnt3+=b[i].cnt3; b[i].cnt3=0; } mid=(b[i].l+b[i].r)/2; if(id<=mid)question(id,i*2); else question(id,i*2+1); } int main() { ll i,j,x1,x2,y1,y2,t,t1,t2; int T,n,m; scanf("%d",&T); while(T--) { scanf("%d",&n); build(1,maxn,1); for(i=1;i<=n;i++){ scanf("%I64d%I64d%I64d%I64d",&x1,&y1,&x2,&y2); if(max(x1,y1)<min(x2,y2)){ a[1]=1;a[2]=-(x1+y1);a[3]=x1*y1; update(max(x1,y1),min(x2,y2)-1,a,1); } if(x2<y2){ a[1]=0;a[2]=-x1+x2;a[3]=y1*(x1-x2); update(max(x2,y1),y2-1,a,1); } if(y2<x2){ a[1]=0;a[2]=-y1+y2;a[3]=x1*(y1-y2); update(max(y2,x1),x2-1,a,1); } a[1]=0;a[2]=0;a[3]=(x2-x1)*(y2-y1); update(max(x2,y2),maxn,a,1); } scanf("%d",&m); for(i=1;i<=m;i++){ scanf("%I64d",&t); A=B=C=0; question(t,1); printf("%I64d\n",A*t*t+B*t+C); } } return 0; }树状数组代码://速度快了一倍= =
#include <iostream> #include <cstdio> #include <cstring> #define maxn 200010 #define ll __int64 using namespace std; ll max(ll a,ll b){ return a>b?a:b; } ll min(ll a,ll b){ return a<b?a:b; } struct segment{ ll b[maxn],ans; void clear(){ memset(b,0,sizeof(b)); } int lowbit(ll x){ return x&(-x); } void update(ll pos,ll num){ while(pos<=maxn){ b[pos]+=num;pos+=lowbit(pos); } } ll getsum(ll pos){ ans=0; while(pos>0){ ans+=b[pos];pos-=lowbit(pos); } return ans; } }A,B,C; void fun(ll l,ll r,ll num1,ll num2,ll num3){ A.update(l,num1);A.update(r+1,-num1); B.update(l,num2);B.update(r+1,-num2); C.update(l,num3);C.update(r+1,-num3); } int main() { int t; scanf("%d",&t); while(t--) { A.clear(); B.clear(); C.clear(); int n,m; scanf("%d",&n); for(int i=0;i<n;i++) { ll x1,y1,x2,y2; scanf("%I64d%I64d%I64d%I64d",&x1,&y1,&x2,&y2); if(max(x1,y1)<min(x2,y2)) fun(max(x1,y1),min(x2,y2),1,-(x1+y1),x1*y1); if(x2<y2) fun(max(x2,y1)+1,y2,0,-x1+x2,y1*(x1-x2)); if(y2<x2) fun(max(y2,x1)+1,x2,0,-y1+y2,x1*(y1-y2)); fun(max(x2,y2)+1,maxn,0,0,(x2-x1)*(y2-y1)); } scanf("%d",&m); while(m--) { ll t; scanf("%I64d",&t); ll ans=0; ans+=A.getsum(t)*t*t; ans+=B.getsum(t)*t; ans+=C.getsum(t); printf("%I64d\n",ans); } } return 0; }
版权声明:本文为博主原创文章,未经博主允许不得转载。
时间: 2024-08-24 15:36:43