FZU 2187 扫描线

来源:互联网 发布:spss mac 23 序列号 编辑:程序博客网 时间:2024/06/06 18:28

和HDU 1255差不多

这次是求只被覆盖一次的矩形面积和

修改callen函数即可

data[k].len表示被覆盖的纵长度

data[k].key表示被只被覆盖一次的纵长度


#include "stdio.h"#include "string.h"#include "stdlib.h"#include "math.h"#include "iostream"#include "algorithm"using namespace std;struct comp{int flag; // 1代表矩阵的起始边,-1代表矩阵的终边int  x,y1,y2;} node[200010];  // 记录每条与Y轴平行的边struct comp1{int l,r,mid,s;int  len,ml,mr,key;} data[800010];bool cmp(comp a,comp b){return a.x<b.x;}int y[200010];void build(int l,int r,int k){data[k].l=l;data[k].r=r;data[k].mid=(l+r)/2;data[k].ml=y[l];data[k].mr=y[r];data[k].s=data[k].len=0;if (l+1==r) return ;build(l,data[k].mid,k*2);build(data[k].mid,r,k*2+1); // 注意应该从 mid开始}void callen(int k) // 计算纵长度{    if (data[k].s>1)    {        data[k].len=data[k].mr-data[k].ml;        data[k].key=0;    }    else        if (data[k].s==1)        {            data[k].len=data[k].mr-data[k].ml;            if (data[k].l+1==data[k].r) data[k].key=data[k].len;            else data[k].key=data[k].len-data[k*2].len-data[k*2+1].len;        }        else            if (data[k].l+1==data[k].r)            {                data[k].len=0;                data[k].key=0;            }            else            {                data[k].len=data[k*2].len+data[k*2+1].len;                data[k].key=data[k*2].key+data[k*2+1].key;            }}void updata(int k,comp b){if (b.y1==data[k].ml && b.y2==data[k].mr){data[k].s+=b.flag;callen(k);return ;}if (b.y2<=data[k*2].mr) updata(k*2,b);elseif (b.y1>=data[k*2+1].ml) updata(k*2+1,b);else{comp temp;temp=b;temp.y2=data[k*2].mr;updata(k*2,temp);temp=b;temp.y1=data[k*2+1].ml;updata(k*2+1,temp);}callen(k);}int main(){int n,i,t,p,Case;int x1,x2,y1,y2;long long  sum;scanf("%d",&Case);for (p=1;p<=Case;p++){    scanf("%d",&n);t=1;for (i=0;i<n;i++){scanf("%d%d%d%d",&x1,&y1,&x2,&y2);node[t].x=x1;node[t].y1=y1;node[t].y2=y2;node[t].flag=1;y[t++]=y1;node[t].x=x2;node[t].y1=y1;node[t].y2=y2;node[t].flag=-1;y[t++]=y2;}sort(node+1,node+t,cmp);sort(y+1,y+t);build(1,t-1,1);sum=0;updata(1,node[1]);for (i=2;i<t;i++){sum+=(long long )data[1].key*(node[i].x-node[i-1].x);updata(1,node[i]);}printf("Case %d: %lld\n",p,sum);}return 0;}


0 0
原创粉丝点击