hdu 3642 Get The Treasury(扫描线、立方体交)
来源:互联网 发布:c语言 函数库啥意思- 编辑:程序博客网 时间:2024/05/17 20:28
hdu 3642 Get The Treasury
和此题做对比
这个题相比之有不同的就是需要求覆盖两次以上部分的体积,而且不一定底面都是在xy平面上,但处理思路依然是一样的
主要问题是如果沿用上次求面积交的方法没办法处理处覆盖多次的情况,于是又回到了老办法上,没有多想先写下试试,结果没有超时,就这样了
不过还是要研究下更快的求面积交的算法,感觉挺重要
#include<cstdio>#include<cstring>#include<algorithm>#include<iostream>#include<map>using namespace std;#define MAXN 2005#define ll long long#define lch p<<1#define rch p<<1|1#define mid (t[p].l+t[p].r)>>1struct node{int l,r;int cover;ll dy;}t[MAXN<<2];struct segment{ll x,yd,yu;ll zd,zu;int lr;segment(){}segment(ll x,ll yd,ll yu,ll zd,ll zu,int lr):x(x),yd(yd),yu(yu),zd(zd),zu(zu),lr(lr){}bool operator < (const segment &a) const{return x<a.x;}}seg[MAXN],tseg[MAXN];ll yy[MAXN],zz[MAXN];void construct(int l,int r,int p){t[p].l=l,t[p].r=r;t[p].cover=t[p].dy=0;if(l+1>=r) return ;int m=mid;construct(l,m,lch);construct(m,r,rch);}void modify(int l,int r,int op,int p){if(t[p].l==l&&t[p].r==r){t[p].cover+=op;return ;}int m=mid;if(r<=m) modify(l,r,op,lch);else if(l>=m) modify(l,r,op,rch);else modify(l,m,op,lch),modify(m,r,op,rch);}void query(ll &ans,int p){if(t[p].cover>2){ans+=yy[t[p].r]-yy[t[p].l];}else if(t[p].l+1<t[p].r){t[lch].cover+=t[p].cover;t[rch].cover+=t[p].cover;query(ans,lch);query(ans,rch);t[lch].cover-=t[p].cover;t[rch].cover-=t[p].cover;}}map<ll,int> rec;int main(){int cas,n;ll x1,y1,z1,x2,y2,z2;scanf("%d",&cas);for(int kcas=1;kcas<=cas;kcas++){scanf("%d",&n);int i;for(i=0;i<(n<<1);i+=2){scanf("%I64d%I64d%I64d%I64d%I64d%I64d",&x1,&y1,&z1,&x2,&y2,&z2);seg[i]=segment(x1,y1,y2,z1,z2,1);seg[i+1]=segment(x2,y1,y2,z1,z2,-1);yy[i]=y1,yy[i+1]=y2;zz[i]=z1,zz[i+1]=z2;}sort(seg,seg+i);sort(yy,yy+i);sort(zz,zz+i);int up=unique(yy,yy+i)-yy;int upz=unique(zz,zz+i)-zz;rec.clear();for(int i=0;i<up;i++)rec[yy[i]]=i;ll ans=0;for(int j=0;j<upz-1;j++){ll area=0;int tn=0;for(int i=0;i<(n<<1);i++)if(seg[i].zd<=zz[j]&&seg[i].zu>=zz[j+1])tseg[tn++]=seg[i];if(tn==0) continue; construct(0,up-1,1);for(int i=0;i<tn-1;i++){modify(rec[tseg[i].yd],rec[tseg[i].yu],tseg[i].lr,1);ll tans=0;query(tans,1);area+=(tseg[i+1].x-tseg[i].x)*tans;}ans+=area*(zz[j+1]-zz[j]);}printf("Case %d: %I64d\n",kcas,ans);}return 0;}
0 0
- hdu 3642 Get The Treasury(扫描线、立方体交)
- HDU 3642 Get The Treasury (线段树、扫描线求立方体体积交)
- hdu 3642 Get The Treasury(扫描线)
- 【hdu3642】Get The Treasury (立方体体积交+线段树+扫描线)
- hdu 3642 Get The Treasury 扫描线求长方体体积交
- HDU 3642——Get The Treasury(线段树+扫描线+离散化+体积交多次)
- hdu 3642 Get The Treasury (三维的扫描线)
- HDU 3642 Get The Treasury 【线段树】【扫描线】
- hdoj 3642 Get The Treasury 【线段树 扫描线 求立方体积交】
- Hdu 3642 Get The Treasury(长方体的体积交)
- hdu 3642 Get The Treasury
- HDU 3642 Get The Treasury
- HDU-3642-Get The Treasury
- HDU 3642 Get The Treasury
- HDU 3642 Get The Treasury (线段树 体积交[转成面积交])
- HDU 3642 Get The Treasury(离散化+线段树:扫描线)
- Get The Treasury - HDU 3642 扫描线 重复三次的体积
- HDU 3642 Get The Treasury[离散化 + 扫描线 + 线段树]
- Travelling - HDU 3001 状压dp
- GetMemory错误讲解(指针练习)
- linux下运行netca不能启动监听
- C++/mfc错误总结
- 51单片机的复位电路原理
- hdu 3642 Get The Treasury(扫描线、立方体交)
- 字符串的一些处理
- struts2 session过期或无操作时自动返回登录
- 一些iphone 开源项目
- NeHe学习之基础二
- Sip基础四
- SpringMVC实战
- POJ 3268 Silver Cow Party(Dijkstra)
- linux tr命令详解