ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

[USACO12FEB] Overplanting S - 二维差分,离散化

2020-03-14 10:03:26  阅读:347  来源: 互联网

标签:USACO12FEB ++ px py 差分 int my mx Overplanting


给定 \(N\leq 1000\) 个矩形,坐标绝对值 \(\leq 10^8\),求矩形面积并。

Solution

难度:L3

离散化后用二维前缀和处理出每个小方格被覆盖的次数,最后扫一遍所有小方格,如果覆盖次数 \(\geq 1\),就算上这个小方格的面积

#include <bits/stdc++.h>
using namespace std;

#define int long long
const int N = 4005;
struct rect {
    int x1,y1,x2,y2;
} r[N];
int n,px[N],py[N],s[N][N],idx,idy;
map <int,int> mx,my;

signed main() {
    ios::sync_with_stdio(false);
    cin>>n;
    for(int i=1;i<=n;i++) {
        cin>>r[i].x1>>r[i].y2>>r[i].x2>>r[i].y1;
        mx[r[i].x1]++;
        mx[r[i].x2]++;
        my[r[i].y1]++;
        my[r[i].y2]++;
    }
    for(auto i=mx.begin();i!=mx.end();i++) i->second=++idx, px[idx]=i->first;
    for(auto i=my.begin();i!=my.end();i++) i->second=++idy, py[idy]=i->first;
    ++idx;
    ++idy;
    for(int i=1;i<=n;i++) {
        r[i].x1=mx[r[i].x1]+1;
        r[i].x2=mx[r[i].x2]+1;
        r[i].y1=my[r[i].y1]+1;
        r[i].y2=my[r[i].y2]+1;
        s[r[i].x1][r[i].y1]++;
        s[r[i].x2][r[i].y2]++;
        s[r[i].x1][r[i].y2]--;
        s[r[i].x2][r[i].y1]--;
    }
    for(int i=1;i<=idx;i++) {
        for(int j=1;j<=idy;j++) {
            s[i][j]+=s[i-1][j]+s[i][j-1]-s[i-1][j-1];
        }
    }
    int ans=0;
    for(int i=1;i<=idx;i++) {
        for(int j=1;j<=idy;j++) {
            if(s[i][j]>0) ans+=(px[i]-px[i-1])*(py[j]-py[j-1]);
        }
    }
    cout<<ans;

}

标签:USACO12FEB,++,px,py,差分,int,my,mx,Overplanting
来源: https://www.cnblogs.com/mollnn/p/12490665.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有