什么,扫描线需要线段树?

那我第一个不干啊(其实是不会写)

这里介绍一种裸的扫描线:

我们根据x排序,对于相等的 \(x\) ,将 \(y\) 进入和退出分类讨论,然后全部放进set里面.每次 \(x\) 不相等的时候,答案就是 (现在y覆盖的乘以(现在的x-以前的x))

具体的判断方法:

1.y的判断:

将每个长方形的上方点记做出口,下方点记做入口.用一个set记录在某区间内所有的 \(y\) 值.每次从下往上扫,如果某个y是入口就将 \(sz\) +1,否则就将 \(sz\) -1.如果 \(sz\)==0 的时候就将你整个区间的值加到re里面



如图,每次x变的时候只需要保存中间y的值,然后用移动的x乘y就好了

注:我保存y进入和退出状态的原因就是为了记录中间是否有空位.可以发现,如果中间有空位,那么证明所有进入的点已经退出了,所以那一段不需要加上去(不懂可以画一下图)

long long query_up(){
long long re = 0,prev = -1,sz = 0;
for (multiset<pair<long long,bool> >::iterator i=se.begin();i!=se.end();++i){
pair<long long,bool> now = *i;
if (sz==0) {prev = now.f;sz++;}
else if (now.s) sz++;
else sz--;
if (sz==0) re+=(now.f-prev);//如果现在所以的y都已经出去了,那么答案就是最后一个y的出口-第一个y的入口
}
return re;
}

想到了这点以后这题的难点基本上就解决了.

因为有负数,我将每个数都加上了1e8,这样就可以完全不管负数了

不开long long见祖宗

完整代码:

#include <iostream>
#include <set>
#include <algorithm>
#include <iterator>
#include <utility>
using namespace std;
const long long MAXN = 1e5+5;
#define pp pair<long long,long long>
#define f first
#define s second
long long n,ans = 0;
multiset<pair<long long,bool> > se;
struct Edge{
long long x,y; bool in,in2;
}edge[MAXN*4];
Edge add_edge(long long a, long long b, bool bo,bool bo2){
Edge tmp;
tmp.x = a;
tmp.y = b;
tmp.in = bo;
tmp.in2 = bo2;
return tmp;
}
bool sorted(Edge a, Edge b){
return a.x<b.x;
}
long long query_up(){
long long re = 0,prev = -1,sz = 0;
for (multiset<pair<long long,bool> >::iterator i=se.begin();i!=se.end();++i){
pair<long long,bool> now = *i;
if (sz==0) {prev = now.f;sz++;}
else if (now.s) sz++;
else sz--;
if (sz==0) re+=(now.f-prev);
}
return re;
}
int main(){
cin >> n;
for (long long i=0;i<n;i++){
long long a,b,c,d; cin >> a >>b >> c >> d;
a+=1e8;b+=1e8;c+=1e8;d+=1e8;
edge[4*i] = add_edge(a,b,1,0);
edge[4*i+1] = add_edge(a,d,1,1);
edge[4*i+2] = add_edge(c,b,0,0);
edge[4*i+3] = add_edge(c,d,0,1);
//两种状态,第一种表示x变不变,第二种表示y变不变
}
sort(edge,edge+4*n,sorted);
long long prev = 0;
for (long long i=0;i<4*n;i++){
if (edge[i].x!=prev){//如果x变了
if (se.size()) ans += (edge[i].x-prev)*query_up();
prev = edge[i].x;
}
if (edge[i].in==1) se.insert(make_pair(edge[i].y,edge[i].in2));//如果这点是进入的点,就将y加入set
else se.erase(se.find(make_pair(edge[i].y,edge[i].in2)));
//否则将y扔出set
}
cout << ans;
}

留一组测试数据造福后人

3
3 7 7 3
1 5 5 1
2 2 7 -2

答案45

题解 P1884 【[USACO12FEB]过度种植(银)Overplanting 】的更多相关文章

  1. 线段树+扫描线【p1884】[Usaco12FEB]过度种植(银)Overplanting …

    Description 在一个笛卡尔平面坐标系里(则X轴向右是正方向,Y轴向上是正方向),有\(N(1<=N<=1000)\)个矩形,第i个矩形的左上角坐标是\((x1, y1)\),右下 ...

  2. 「洛谷1884」「USACO12FEB」过度种植【离散化扫描线】

    题目链接 [洛谷传送门] 题解 矩阵面积的并模板.(请求洛谷加为模板题) 很明显是要离散化的. 我们将矩阵与\(x\)轴平行的两个线段取出来.并且将这两个端点的\(x1\)和\(x2\)进行离散化. ...

  3. 洛谷P1884 [USACO12FEB]Overplanting S (矩形切割)

    一种矩形切割的做法: 1 #include<bits/stdc++.h> 2 using namespace std; 3 typedef long long LL; 4 const in ...

  4. [luogu] P2519 [HAOI2011]problem a (贪心)

    P2519 [HAOI2011]problem a 题目描述 一次考试共有n个人参加,第i个人说:"有ai个人分数比我高,bi个人分数比我低."问最少有几个人没有说真话(可能有相同 ...

  5. usaco 2010年3月银组题解

    usaco银组解题报告 一.石子游戏如果把‘O’当作0,‘X’当做1,则N个洞的每一种状态都可以看做是一个N位二进制数.于是,这个问题就变成了求环绕的N位格雷码.幸运的是,这个结构很容易就能够用一个简 ...

  6. 洛谷 P1821 [USACO07FEB]银牛派对Silver Cow Party 题解

    P1821 [USACO07FEB]银牛派对Silver Cow Party 题目描述 One cow from each of N farms (1 ≤ N ≤ 1000) conveniently ...

  7. 【luogu P1821 [USACO07FEB]银牛派对Silver Cow Party】 题解

    题目链接:https://www.luogu.org/problemnew/show/P1821 反向多存一个图,暴力跑两遍 #include <cstdio> #include < ...

  8. USACO[19-20]Dec银组题解

    1,MooBuzz 这题其实是道数学题. 我们先找找符合要求的数:1,2,4,7,8,11,13,14…… 我们发现再往后找都是这8个数中的一个加15k如:16……19……29…… 找规律发现k=n/ ...

  9. 【题解】Luogu p3047 [USACO12FEB]附近的牛Nearby Cows 树型dp

    题目描述 Farmer John has noticed that his cows often move between nearby fields. Taking this into accoun ...

随机推荐

  1. vue学习(十)mixin 偷懒

    一 mixin混入偷懒技术 架子 <div id="app"> {{msg}} </div> <script> let app = new Vu ...

  2. 修改虚拟机ip

    在命令行里输入:more  /etc/sysconfig/network-scripts/ifcfg-eth0 (注意more后要空一格,还有eth0,最后是数字零).  然后再输入:ifconfig ...

  3. 【STM32H7教程】第48章 STM32H7的FMC总线应用之是32路高速IO扩展

    完整教程下载地址:http://www.armbbs.cn/forum.php?mod=viewthread&tid=86980 第48章       STM32H7的FMC总线应用之是32路 ...

  4. 常用ES6-ES10知识点总结

    在工作中我们会常用到的一些es6-es10的一些特性还记得多少,今天就让我们重新复习一遍 ES6语法 1.Let 1.let声明的变量具有块级作用域, { let a = 1 } console.lo ...

  5. 005、mysql查询表的结构

    EXPLAIN hcc_ip 如果有一个表,表明为hcc_ip,使用以上语句可以得到下图 不忘初心,如果您认为这篇文章有价值,认同作者的付出,可以微信二维码打赏任意金额给作者(微信号:38247724 ...

  6. slf4j NoSuchMethodError 错误 ---- 版本冲突

    java.lang.NoSuchMethodError: org.slf4j.spi.LocationAwareLogger.log(Lorg/slf4j/Marker;Ljava/lang/Stri ...

  7. POJ 1422 DAG最小路径覆盖

    求无向图中能覆盖每个点的最小覆盖数 单独的点也算一条路径 这个还是可以扯到最大匹配数来,原因跟上面的最大独立集一样,如果某个二分图(注意不是DAG上的)的边是最大匹配边,那说明只要取两个端点只要一条边 ...

  8. leetcode - 两数之和Ⅳ 输入BST(653)

    题目描述:给定一个二叉搜索树和一个目标结果,如果 BST 中存在两个元素且它们的和等于给定的目标结果,则返回 true. 解题思路:根据二叉搜索树的特点,对二叉搜索树进行中序遍历可以得到一个从小到达排 ...

  9. SAP_BASIS常用事务代码

    1.SM66:监控当前系统的所有进程: 2.SM50/SM51:监视当前客户端的所有进程: 3.AL08:查看系统当前用户登录情况: 4.ST22:查看系统发生的DUMP: 5.ST03:查看系统当前 ...

  10. Arduino LiquidCrystal库函数中文对照

    我之所以会试试看LCD+Shield,一是因为我本来就有这块板,但一直不会用,第二个原因则是Arduino+Cookbook这本书.O’Reilly的这本手册对你在Arduino实践中的各种问题(几乎 ...