poj 1263 反射问题

//题意：一条射线在n个圆间反射，求依次经过哪些圆。

//用了一些自己写的简单的函数，这道题要注意第11个圆也得判断。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const double EP=1e-8;
const double INF=1e9;
struct Point{double x, y;};
struct Circle{Point c; double r;}sp[26];
struct Line{double a, b, c;};
struct Ray{Point s, e;}ray1, ray2;
double dist(Point p1, Point p2){
    return ((p1.x-p2.x)*(p1.x-p2.x)+(p1.y-p2.y)*(p1.y-p2.y));
}
double dot_mult(Point sp, Point ep, Point op){
    return (sp.x-op.x)*(ep.x-op.x)+(sp.y-op.y)*(ep.y-op.y);
}
Line makeline(Point p1, Point p2){
    Line t1;
    int sign=1;
    t1.a=p2.y-p1.y;
    if(t1.a<0){
        t1.a*=-1;
        sign=-1;
    }
    t1.b=sign*(p1.x-p2.x);
    t1.c=sign*(p1.y*p2.x-p1.x*p2.y);
    return t1;
}
double dist(Point p1, Ray ray){
    Line l1=makeline(ray.s, ray.e);
    return fabs(l1.a*p1.x+l1.b*p1.y+l1.c)/sqrt(l1.a*l1.a+l1.b*l1.b);
}
double dist(Ray ray, Circle cir){
    double di=dist(cir.c, ray);
    if(di-cir.r>EP)
        return INF;
    return sqrt(dist(ray.s, cir.c)-di*di)-sqrt(cir.r*cir.r-di*di);
}
Point vert_p(Point p1, Line l1){
    Point p2;
    p2.x=(l1.b*(l1.b*p1.x-l1.a*p1.y)-l1.a*l1.c)/(l1.a*l1.a+l1.b*l1.b);
    p2.y=(l1.a*(l1.a*p1.y-l1.b*p1.x)-l1.b*l1.c)/(l1.a*l1.a+l1.b*l1.b);
    return p2;
}
Point sym_p(Point sp, Point op){
    Point p1;
    p1.x=2*op.x-sp.x;
    p1.y=2*op.y-sp.y;
    return p1;
}
Point sym_line(Point p1, Line l1){
    return sym_p(p1, vert_p(p1, l1));
}

int main(){
    //freopen("1.txt", "r", stdin);
    int n, i, t, tt, sum, test=0;
    double d, tmp;
    while(scanf("%d", &n)&&n){
        sum=0;tt=-1;test++;
        printf("Scene %d/n", test);
        for(i=0; i<n; i++)
            scanf("%lf %lf %lf", &sp[i].c.x, &sp[i].c.y, &sp[i].r);
        scanf("%lf %lf %lf %lf", &ray1.s.x, &ray1.s.y, &ray1.e.x, &ray1.e.y);
        ray1.e.x+=ray1.s.x;ray1.e.y+=ray1.s.y;
        while(sum<=10){
            t=-1;d=INF;
            for(i=0; i<n; i++){
                if(tt==i)continue;
                if(dot_mult(sp[i].c, ray1.e, ray1.s)<EP)    //判断射线方向是否射向此圆
                    continue;
                tmp=dist(ray1, sp[i]);
                if(tmp-d<-EP){
                    t=i;
                    d=tmp;
                }
            }
            if(t==-1){
                printf("inf/n/n");
                break;
            }
            if(sum==10)break;
            sum++;tt=t;
            ray2.s.x=ray1.s.x+d/sqrt(dist(ray1.s,ray1.e))*(ray1.e.x-ray1.s.x);
            ray2.s.y=ray1.s.y+d/sqrt(dist(ray1.s,ray1.e))*(ray1.e.y-ray1.s.y);
            ray2.e=sym_line(ray1.s, makeline(sp[t].c, ray2.s));
            printf("%d ", t+1);
            ray1=ray2;
        }
        if(t!=-1)
        printf(".../n/n");
    }
    return 0;
}