bzoj 1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居（切比雪夫距离+multiset贪心+并查集）

1604: [Usaco2008 Open]Cow Neighborhoods 奶牛的邻居

Time Limit: 5 Sec  Memory Limit: 64 MB
Submit: 1092  Solved: 441
[Submit][Status][Discuss]

Description

了解奶牛们的人都知道，奶牛喜欢成群结队．观察约翰的N(1≤N≤100000)只奶牛，你会发现她们已经结成了几个“群”．每只奶牛在吃草的时候有一个独一无二的位置坐标Xi，Yi(l≤Xi，Yi≤[1．.10^9]；Xi，Yi∈整数．当满足下列两个条件之一，两只奶牛i和j是属于同一个群的：

  1．两只奶牛的曼哈顿距离不超过C(1≤C≤10^9)，即lXi - xil+IYi - Yil≤C.

  2．两只奶牛有共同的邻居．即，存在一只奶牛k，使i与k，j与k均同属一个群．

    给出奶牛们的位置，请计算草原上有多少个牛群，以及最大的牛群里有多少奶牛

Input

   第1行输入N和C，之后N行每行输入一只奶牛的坐标．

Output

仅一行，先输出牛群数，再输出最大牛群里的牛数，用空格隔开．

Sample Input

4 2

1 1

3 3

2 2

10 10

Sample Output

2 3

和http://blog.csdn.net/jaihk662/article/details/77918172这题处理方法一样，贪心+multiset

①考虑将曼哈顿距离转化为切比雪夫距离

对于点(x, y)，令(x', y') = (x+y, x-y)，那么有|x1-x2|+|y1-y2|=max(|x1'-x2'|, |y1'-y2'|)<=C

②然后将所有点按x排序，依次加入multiset中，中间保证multiset中任意两点都满足|xi'-xj'|<=C（其中multiset的顺序为y从小到大），这样就只用考虑|y'-yi'|是否<=C了，因为multiset中y有序，所以可以直接二分

④对于第k个点(xk, yk)，将它加入multiset之后判断与它相邻的两个点是否满足|yk'-yi'|<=C，满足就并查集并一起

⑤所有点遍历完毕之后就可以统计答案了

#include<stdio.h>#include<set>#include<algorithm>using namespace std;typedef struct Res{int id;int x, y;bool operator < (const Res &b) const{if(y<b.y)return 1;return 0;}}Res;multiset<Res> st;Res s[100005];int ufs[100005], sum[100005];int Find(int x){if(ufs[x]==0)return x;return ufs[x] = Find(ufs[x]);}bool comp(Res a, Res b){if(a.x<b.x)return 1;return 0;}int main(void){int C, n, i, p, t1, t2, ans, cnt;multiset<Res>::iterator it;scanf("%d%d", &n, &C);for(i=1;i<=n;i++){s[i].id = i;scanf("%d%d", &s[i].x, &s[i].y);s[i].x += s[i].y;s[i].y = s[i].x-2*s[i].y;}sort(s+1, s+n+1, comp);p = 1;st.insert(s[1]);for(i=2;i<=n;i++){while(s[p].x+C<s[i].x && p<i){st.erase(st.find(s[p]));p++;}it = st.lower_bound(s[i]);if(it!=st.end() && (*it).y-C<=s[i].y){t1 = Find((*it).id);t2 = Find(s[i].id);if(t1!=t2)ufs[t1] = t2;}if(it!=st.begin() && (*(--it)).y+C>=s[i].y){t1 = Find((*it).id);t2 = Find(s[i].id);if(t1!=t2)ufs[t1] = t2;}st.insert(s[i]);}ans = cnt = 0;for(i=1;i<=n;i++)sum[Find(i)]++;for(i=1;i<=n;i++){if(sum[i])cnt++;ans = max(ans, sum[i]);}printf("%d %d\n", cnt, ans);return 0;}