poj3277 线段树 + 离散化 WA2次 long long freopen

 

 

 

首先说一下这题的想法，其实就是求总面积，我对x轴进行建线段树，需要离散化。

 

建完线段树以后从新扫描输入的数据，把对应的区域里的高度置成相应的高度，这里需要比较，如果 此区域的之前的高度，比将要加入的高度要低，那么就替换成新高度，相当于覆盖了旧的不能做的高度。

 

 

这里需要注意的一个问题，由于数据值很大，所以需要使用 long long 存储结果，还有对于long long 型第一次使用，哎，竟然要把它强制转化一下，哎，悲剧呀。还有杯具的是，我的freopen没有去掉。wa2次

 

 

#include<memory.h> #include<stdio.h> #include<algorithm> #include<iostream> using namespace std; const int ST_M = 100005 ; struct dnode { int f,t ; int ff,tt; int h; } data[ST_M]; struct node { int x; int i; } sdata[ST_M] ; struct stnode { int left,right; int f,t; int h ; } stree[ST_M*4] ; void create(int st,int ed,int r=1) { stree[r].left= st; stree[r].right= ed; stree[r].h= 0 ; if( st == ed ) { stree[r].f=sdata[st].x ; stree[r].t=sdata[st+1].x ; return ; } int mid=(st+ed)>>1 ; if( ed <=mid ) { create(st,ed,(r<<1)); } else if( st > mid ) { create(st,ed,(r<<1)+1); } else { create(st,mid,(r<<1)); create(mid+1,ed,(r<<1)+1); } } void change(int st,int ed,int di,int r=1) { if( st == stree[r].left && ed == stree[r].right ) { if( stree[r].h < data[di].h ) { stree[r].h = data[di].h; } return ; } int mid = (stree[r].left+stree[r].right)>>1; if( ed <=mid ) { change(st,ed,di,(r<<1)); } else if( st > mid ) { change(st,ed,di,(r<<1)+1); } else { change(st,mid,di,(r<<1)); change(mid+1,ed,di,(r<<1)+1); } } long long output(int loc,int r=1) { int mid; int mh =0; while( stree[r].left != stree[r].right) { if( mh < stree[r].h ) mh = stree[r].h; mid = (stree[r].left+stree[r].right)>>1; if( loc <= mid ) { r = r*2 ; } else { r=r*2+1; } } if( mh < stree[r].h ) mh = stree[r].h; //cout<<stree[r].f<<"_"<<stree[r].t<<"_"<<mh<<endl; long long tem = (long long)mh*(long long)(stree[r].t-stree[r].f) ; //printf("%lld#",tem); return tem; } int cmp(struct node a,struct node b) { return a.x < b.x ; } int main() { freopen("data.txt","r",stdin); int n; int i,j,k; while(cin>>n) { memset(stree,0,sizeof(stree)); for(i=1; i<=n; i++) { scanf("%d%d%d",&data[i].f , &data[i].t ,&data[i].h); sdata[i].x=data[i].f; sdata[2*n+1-i].x=data[i].t; sdata[i].i=-i; sdata[2*n+1-i].i=i; } sort(sdata+1,sdata+2*n+1,cmp) ; k=0; for(i=1; i<=2*n; i++) { j=i; ++k; while( sdata[i].x==sdata[j].x ) { if( sdata[j].i < 0 ) { data[-sdata[j].i ].ff = k ; } else { data[ sdata[j].i ].tt = k ; } j++; } sdata[k].x = sdata[i].x ; i=j-1; } /* for(i=1; i<=n; i++) { cout<<"#"<<data[i].ff<<" "<<data[i].tt<<endl; } */ //cout<<"@"<<k<<endl; create(1,k-1); // here is k-1 ,is the interval not the node //cout<<"create over"<<endl; int di; for(di=1;di<=n;di++) change(data[di].ff,(data[di].tt-1),di); //cout<<"change over"<<endl; long long sum=0; for(i=1;i<=k-1;i++) // k-1 interval { sum+=output(i); } printf("%lld/n",sum); } return 0; }