poj City Horizon (线段树+二分离散)

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

City Horizon

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 15255		Accepted: 4111

Description

Farmer John has taken his cows on a trip to the city! As the sun sets, the cows gaze at the city horizon and observe the beautiful silhouettes formed by the rectangular buildings.

The entire horizon is represented by a number line with N (1 ≤ N ≤ 40,000) buildings. Building i's silhouette has a base that spans locations A_i through B_i along the horizon (1 ≤ A_i < B_i ≤ 1,000,000,000) and has height H_i(1 ≤ H_i ≤ 1,000,000,000). Determine the area, in square units, of the aggregate silhouette formed by all N buildings.

Input

Line 1: A single integer: N  Lines 2..N+1: Input line i+1 describes building i with three space-separated integers: A_i, B_i, and H_i

Output

Line 1: The total area, in square units, of the silhouettes formed by all N buildings

Sample Input

4
2 5 1
9 10 4
6 8 2
4 6 3

Sample Output

16

【题解】：
　　　　这题是二分离散线段树，因为1 ≤Ai<Bi≤ 1,000,000,000，1 ≤Hi≤ 1,000,000,000显然直接建树是不可能的，因为N<40000，数据量不大可以先离散化再建树
　　　　先将输入的Ai Bi Hi 保存在结构体node中，将Ai,Bi保存到一个index数组里面，排序，去重，这样每个Ai,Bi对应一个index数组的位置；
　　　　结构体node按照高度Hi排序，排序之后就不用考虑覆不覆盖的问题了，统一覆盖就行了；
　　　　
　　　　【注】：被坑了，还是自己不够严谨的问题，居然数据类型以为只要用来求和的sum是__int64就可以了，可是在进行乘法的时候height*（r-l）溢出了，果断height__int64过了
　　　　坑惨了，有图有真相：
　　　　

【code】：

 /**
 status:Accepted        memory:7424K
 time:407MS        language:C++
 code length:2904B   author:cj
 oj: poj   submittime:2013-08-06 12:36:08
 */

 #include<iostream>
 #include<stdio.h>
 #include<algorithm>
 #include<string.h>

 using namespace std;

 #define N 40010
 #define lson p<<1
 #define rson p<<1|1

 struct Nod
 {
     int from,to,height;
 }node[N];

 struct TNod
 {
     int l,r;
     __int64 height;  //我以为只要sum值设为__int64就行了，height值也应该是__int64的，因为后面要进行乘法运算
     __int64 sum;
 }tnode[N<<];

 bool cmp(Nod a,Nod b) //按高度排序
 {
     return a.height<b.height;
 }

 int temp[N<<],index[N<<],brush[N<<];

 int findPos(int a,int k)  //二分查找离散
 {
     int l,r,mid=-;
     l=;
     r=k;
     while(l<=r)
     {
         mid = (l+r)>>;
         if(a<index[mid])  r=mid-;
         else if(a>index[mid])   l=mid+;
         else return mid;
     }
     return mid;
 }

 void building(int l,int r,int p)  //建树
 {
     tnode[p].l = l;
     tnode[p].r = r;
     tnode[p].height = ;
     tnode[p].sum = ;
     if(l+==r)    return;
     int mid = (l+r)>>;
     building(l,mid,lson);
     building(mid,r,rson);
 }

 void update(int l,int r,int p,int h)
 {
     if(tnode[p].l==l&&r==tnode[p].r)
     {
         tnode[p].height = h;
         tnode[p].sum = (index[r]-index[l])*tnode[p].height;  //这里的height是64位int，相乘得64为int
         return;
     }
     if(tnode[p].height>)　　//向下更新
     {
         tnode[lson].sum = (index[tnode[lson].r] - index[tnode[lson].l])*tnode[p].height;
         tnode[rson].sum = (index[tnode[rson].r] - index[tnode[rson].l])*tnode[p].height;
         tnode[lson].height = tnode[rson].height = tnode[p].height;
         tnode[p].height = ;
     }
     int mid = (tnode[p].l+tnode[p].r)>>;
     if(r<=mid)  update(l,r,lson,h);
     else if(l>=mid)  update(l,r,rson,h);
     else
     {
         update(l,mid,lson,h);
         update(mid,r,rson,h);
     }
     if(tnode[lson].height&&tnode[rson].height&&tnode[lson].height==tnode[rson].height)  //向上更新
     {
         tnode[p].height = tnode[lson].height;
     }
     tnode[p].sum = tnode[lson].sum + tnode[rson].sum;   //向上更新
 }

 int main()
 {
     int n;
     scanf("%d",&n);
     int i,cnt=;
     for(i=;i<n;i++)
     {
         scanf("%d%d%d",&node[i].from,&node[i].to,&node[i].height);
         temp[cnt++]=node[i].from;
         temp[cnt++]=node[i].to;
     }
     sort(temp,temp+cnt);  //排序
     int j=;
     index[]=temp[];
     for(i=;i<cnt;i++)  //去掉重复值
     {
         if(index[j]!=temp[i])
         {
             index[++j]=temp[i];
         }
     }
     int k = j;
     building(,k,);
     sort(node,node+n,cmp);  //按照高度排序
     for(i=;i<n;i++)
     {
         int from = findPos(node[i].from,k);
         int to = findPos(node[i].to,k);
         update(from,to,,node[i].height);
     }
     printf("%I64d\n",tnode[].sum);  //直接输出，头结点的sum值
     return ;
 }