592
POJ 1375 过一点求圆切线极角
题意：给出一点，几个圆的圆心半径，问两条切线与x轴相交组成的区间，然后合并一下。
这题我用解析做，精度损失的太高wa，然后又用极角做，还是精度问题，不过好把数的值扩大一些加上了eps之后过了。
具体看代码。
#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
struct seg
{
double l,r;
};
double dis(double lx,double ly,double cx,double cy)
{
return sqrt((lx-cx)*(lx-cx)+(ly-cy)*(ly-cy));
}
bool cmp(seg a,seg b)
{
return a.l<b.l;
}
double getx(double angle,double lx,double ly)
{
return lx-((ly)/tan(angle));
}
int main()
{
int n,i;
seg g[505];
double eps=1e-8;
double lx,ly,cx,cy,r;
while(~scanf("%d",&n),n)
{
scanf("%lf%lf",&lx,&ly);
for(i=0; i<n; i++)
{
scanf("%lf%lf%lf",&cx,&cy,&r);
double a0=atan2(ly-cy,lx-cx),a1=asin(r/dis(lx,ly,cx,cy));
g[i].l=getx(a0-a1,lx,ly),g[i].r=getx(a0+a1,lx,ly);
if(g[i].l>0)
g[i].l+=eps;
else
g[i].l-=eps;
if(g[i].r>0)
g[i].r+=eps;
else
g[i].r-=eps;
if(g[i].l>g[i].r)swap(g[i].l,g[i].r);
}
sort(g,g+n,cmp);
i=0;
while(i<n)
{
cx=g[i].l;
cy=g[i++].r;
while(i<n&&g[i].l<=cy)cy=max(cy,g[i++].r);
printf("%.2f %.2f\n",cx,cy);
}
puts("");
}
return 0;
}
最后更新：2017-04-03 16:59:42