poj3109 Inner Vertices 扫描线+树状数组

http://poj.org/problem?id=3109

无限大的平面上有n个黑整点，剩下的整点都是白的，如果一个白整点四周都有黑点，那么它就会变成黑点，求最后有多少个黑点。

首先白点的变化不会发生连锁反应，意味着这个变化过程一定能够终止。

变化的白点肯定是横向线段与纵向线段的交点，所以只需要找出交点的个数即可。

先离散化再找出所有的横向与纵向线段，然后扫描线+树状数组统计即可。

#include <iostream>
#include <cstring>
#include <cctype>
#include <cstdlib>
#include <cmath>
#include <vector>
#include <cstdio>
#include <set>
#include <map>
#include <algorithm>
#define N 100005
typedef long long LL;
#define fi first
#define se second
using namespace std;
int n , m;
pair<int , int> a[N];
int top[N] , bot[N] , d[N] , c[N];
struct opt
{
  int p , x;
  bool f;
  bool operator < (const opt& r) const {
    if (x == r.x)
      return f > r.f;
    return x < r.x;
  }
}t[N << 1];

void add(int x , int w)
{
  for (int i = x ; i <= m ; i += i & -i)
    c[i] += w;
}
int sum(int x)
{
  int _s = 0;
  for (int i = x ; i ; i &= i - 1)
    _s += c[i];
  return _s;
}


void work()
{
  int i , j , x , y;
  scanf("%d",&n);
  for (i = 1 ; i <= n ; ++ i)
    scanf("%d%d",&a[i].fi , &a[i].se) , d[i] = a[i].se;
  sort(a + 1 , a + 1 + n) , sort(d + 1 , d + 1 + n);
  m = unique(d + 1 , d + 1 + n) - d - 1;
  for (i = 1 ; i <= n ; ++ i)
    a[i].se = lower_bound(d + 1 , d + 1 + m , a[i].se) - d;
  for (i = 1 ; i <= m ; ++ i) top[i] = 1 << 30;
  for (i = 1 ; i <= n ; ++ i)
  {
    if (top[a[i].se] == 1 << 30)
      top[a[i].se] = a[i].fi;
    bot[a[i].se] = a[i].fi;
  }
  for (i = 1 ; i <= m ; ++ i)
    t[2 * i - 1] = (opt) {i , top[i] , 1} , t[2 * i] = (opt) {i , bot[i] , 0};
  sort(t + 1 , t + m + m + 1);
  LL ans = 0; j = 1;
  for (i = 1 ; i <= n ; )
  {
    x = a[i].fi;
    int l = 1 << 30 , r = -1 << 30;
    while (i <= n && a[i].fi == x)
      l = min(a[i].se , l) , r = max(a[i].se , r) , ++ i;
    while (j <= m + m && t[j].x < x) ++ j;
    while (j <= m + m && t[j].x == x && t[j].f)
      add(t[j].p , 1) , ++ j;
    ans += sum(r) - sum(l -1);
    while (j <= m + m && t[j].x == x && !t[j].f)
      add(t[j].p , -1) , ++ j;
  }
  cout << ans << endl;
}

int main()
{
  work();
  return 0;
}