hdu3124——Moonmist

最近圆对。

思路：二分+扫描。

用两根线扫描，保证l<r。

参考：http://hi.baidu.com/bobobry/blog/item/daef6438fd21e03096ddd81d.html

附set用法：http://www.cnblogs.com/motadou/archive/2009/12/17/1626127.html

#include<iostream> #include<cstring> #include<cstdio> #include<cmath> #include<algorithm> #include<set> using namespace std; #define N 51000 #define esp 1e-8 int n; double add; class point { public: double x,y,r; }; point p[N]; class yline { public: double x,y; int id; bool operator <(const yline & a)const { return x <a.x ||(x==a.x &&y>a.y ); return false; } }; yline ll[N],rl[N]; class setnode { public: int id; bool operator<(const setnode &a )const { return (p[id].y >p[a.id ].y ||(p[id].y ==p[a.id ].y&&p[id].x <p[a.id ].x )) ; } }; set<setnode> line; set<setnode>::iterator it,pre,later;//迭代器，能够遍历容器 double sqr(double a) { return a*a; } bool crosspoint(int a,int b)//判断两圆是否相交 { if((sqr(p[a].x -p[b].x )+sqr(p[a].y-p[b].y ))<=sqr(add+add+p[a].r +p[b].r) ) { return true; } return false; } bool insert(int id)//加入圆 { setnode temp; temp.id =id; it=line.insert (temp).first; pre=later=it; later++; if(it==line.begin ()) { if(later==line.end ()) return false; else { if(crosspoint(it->id ,later->id )) return true; return false; } } else { pre--; if(later==line.end()) { if(crosspoint(it->id ,pre->id )) return true; return false; } else { if(crosspoint(it->id,pre->id)||crosspoint(it->id,later->id )) return true; return false; } } } void dele(int id)//删除圆 { setnode temp; temp.id =id; line.erase (temp); } bool judge(double mid)//判断是否相交 { line.clear (); int i=0,j=0; while(i<n||j<n) { if(i==n) dele(rl[j++].id ); else if(j==n) { if(insert(ll[i++].id )) return true; } else if(i<n&&ll[i].x -mid<rl[j].x +mid) { if(insert(ll[i++].id )==true) return true; } else { dele(rl[j++].id ); } } return false; } int main() { int i,ca,j; cin>>ca; while(ca--) { cin>>n; for(i=0;i<n;i++) { scanf("%lf%lf%lf",&p[i].x ,&p[i].y ,&p[i].r ); ll[i].x=p[i].x -p[i].r ; ll[i].y=p[i].y; ll[i].id =i; rl[i].x =p[i].x +p[i].r ; rl[i].y=p[i].y; rl[i].id =i; } sort(ll,ll+n); sort(rl,rl+n); double l=0; double r=sqrt(sqr(p[0].x-p[1].x )+sqr(p[0].y-p[1].y ))-p[0].r -p[1].r ; double mid=0; while(r-l>esp) { mid=(r+l)*0.5; add=mid; if(judge(mid))//出现相交 r=mid; else l=mid; } printf("%.6f\n",mid+mid);//即l+r } return 0; }