POJ 1375 圆的切线

题意：

求从光源透过圆形的阴影的线段，并排序输出。

题解：

求圆的切线，利用勾股定理和坐标旋转即可。

处理出所有与地平线的交点，标号正负1，排序，然后扫一遍就好了~

#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>

#define N 2222
#define BUG system("pause")

using namespace std;

struct PO
{
    double x,y,r;
    int fg;
    void prt() {printf("%lf    %lf\n",x,y);}
}p[N],s;

struct FP
{
    double x;
    int fg;
}fp[N];

struct LI
{
    PO a,b;
    int fg;
}li[N],flr;

int n,cnt,gs;
double lt;

inline PO operator +(PO a,PO b)
{
    PO as;
    as.x=a.x+b.x; as.y=a.y+b.y;
    return as;
}

inline PO operator -(PO a,PO b)
{
    PO as;
    as.x=a.x-b.x; as.y=a.y-b.y;
    return as;
}

inline PO operator /(PO a,double k)
{
    a.x/=k; a.y/=k;
    return a;
}

inline PO operator *(PO a,double k)
{
    a.x*=k; a.y*=k;
    return a;
}

inline void read()
{
    scanf("%lf%lf",&s.x,&s.y);
    for(int i=1;i<=n;i++) scanf("%lf%lf%lf",&p[i].x,&p[i].y,&p[i].r);
}

inline PO rotate(PO a,double sss,double ccc)
{
    PO ans;
    ans.x=a.x*ccc-a.y*sss;
    ans.y=a.x*sss+a.y*ccc;
    return ans;
}

inline double getdis(PO &a,PO &b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

inline void getli(PO &s,PO &c,PO &a,PO &b)
{
    double sc=getdis(s,c);
    double sp=sqrt(sc*sc-c.r*c.r);
    double sss=c.r/sc;
    double ccc=sp/sc;
    PO len=(c-s)*sp/sc;
    a=s+rotate(len,sss,ccc);
    b=s+rotate(len,-sss,ccc);
}

inline double cross(PO &a,PO &b,PO &c)
{
    return (b.x-a.x)*(c.y-a.y)-(c.x-a.x)*(b.y-a.y);
}

inline PO getpoint(PO &a,PO &b,PO &c,PO &d)
{
    double k1=cross(a,d,c),k2=cross(b,c,d);
    PO res=b-a;
    res=a+res*k1/(k1+k2);
    return res;
}

inline bool cmp(const FP &a,const FP &b)
{
    return a.x<b.x;
}

inline void go()
{
    cnt=gs=0;PO a,b;
    for(int i=1;i<=n;i++)
    {
        getli(s,p[i],a,b);
        if(a.x>b.x) swap(a,b);
        li[++cnt].a=s; li[cnt].b=a; li[cnt].fg=1;
        li[++cnt].a=s; li[cnt].b=b; li[cnt].fg=-1;
    }
    flr.a.x=flr.a.y=flr.b.y=0; flr.b.x=1;
    for(int i=1;i<=cnt;i++)
    {
        PO tmp=getpoint(flr.a,flr.b,li[i].a,li[i].b);
        fp[i].x=tmp.x;
        fp[i].fg=li[i].fg;
    }
    sort(fp+1,fp+1+cnt,cmp);
    for(int i=1;i<=cnt;i++)
    {
        gs+=fp[i].fg;
        if(gs==1&&fp[i].fg==1) lt=fp[i].x;//注意缺一不可！ 
        else if(gs==0) printf("%.2lf %.2lf\n",lt,fp[i].x);
    }
    puts("");
}

int main()
{
    while(scanf("%d",&n),n) read(),go();
    return 0;
}