划定农田_纪中1764_暴力

题目描述

---

在寒冷的冬天过去以后，农夫约翰决定重新划定他的农田了。农田是由N个用栅栏围起来的封闭区域组成的，每个封闭区域都是一个矩形，并且平行于x轴和y轴。这些封闭区域可能会存在着包含关系，但是一定不存在相交关系（边相交）。

问题描述：

---

请帮助约翰计算没有被包含在其他的区域内的区域共有多少块。

输入

---

第一行一个整数N，表示封闭区域的个数。

下来N行，每行四个非负整数x1,y1,x2,y2，（x1，y1）表示矩形的左下角的点，（x2，y2）表示矩形右上角的点。

输出

---

输出没有被其他区域包含的区域的数量。

数据范围

---

1<=N<=50000，坐标的范围不超过1000000。

Analysis

---

赤裸裸的暴力了这就是，按左下角x坐标排序然后用vis记录是否被包含，如果被包含就跳过退
当然这样做仍不够我们切掉这题，于是需要寻求一些更强力的优化
如果当前矩形左下角x坐标已经大于枚举的右下角，那么由于是按照上升排序的，后面的矩形都不可能有包含关系，这么做的效率是接近O(n)的

Code

---

#include <cstdio>#include <cstdlib>#include <cstring>#include <cmath>#include <ctime>#include <iostream>#include <algorithm>#include <string>#include <vector>#include <deque>#include <list>#include <set>#include <map>#include <stack>#include <queue>#include <numeric>#include <iomanip>#include <bitset>#include <sstream>#include <fstream>#define debug puts("-----")#define rep(i, st, ed) for (int i = st; i <= ed; i += 1)#define drp(i, st, ed) for (int i = st; i >= ed; i -= 1)#define fill(x, t) memset(x, t, sizeof(x))#define min(x, y) x<y?x:y#define max(x, y) x>y?x:y#define PI (acos(-1.0))#define EPS (1e-8)#define INF (1<<30)#define ll long long#define db double#define ld long double#define N 50001#define E N * 8 + 1#define L 255using namespace std;struct rect{int l, d, r, u;}r[N];bool vis[N];__attribute__((optimize("O2")))inline int read(){    int x = 0, v = 1;    char ch = getchar();    while (ch < '0' || ch > '9'){        if (ch == '-'){            v = -1;        }        ch = getchar();    }    while (ch <= '9' && ch >= '0'){        x = (x << 1) + (x << 3) + ch - '0';        ch = getchar();    }    return x * v;}__attribute__((optimize("O2")))inline int cmp(rect a, rect b){    return a.l < b.l;}__attribute__((optimize("O2")))int main(void){    freopen("painting.in", "r", stdin);    freopen("painting.out", "w", stdout);    int n = read();    rep(i, 1, n){        r[i] = (rect){read(), read(), read(), read()};    }    sort(r + 1, r + n + 1, cmp);    int ans = 0;    rep(i, 1, n){        if (!vis[i]){            ans += 1;            rep(j, i + 1, n){                if (r[i].u > r[j].u && r[i].l < r[j].l && r[i].r > r[j].r && r[i].d < r[j].d){                    vis[j] = true;                }                if (r[j].l > r[i].r){                    break;                }            }        }    }    printf("%d\n", ans);    return 0;}