划定农田_纪中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 255
using 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;
}