題意：
有一些圓形的鏡子，一束光線射進(jìn)來(lái)，試模擬光線的反射過(guò)程
測(cè)試數(shù)據(jù)：
/Files/yzhw/reflect.rar解法：
這個(gè)問(wèn)題首先從整體上看應(yīng)該用迭代法，每次確定v向量和初始位置，然后循環(huán)10次即可
下面解決2個(gè)問(wèn)題：
1、求與某個(gè)圓的的第一個(gè)交點(diǎn)
2、求新的速度矢量
對(duì)以第一個(gè)問(wèn)題，我們?cè)O(shè)參數(shù)k，然后列方程:
設(shè)x',y'為圓心坐標(biāo)，r為半徑。x0,y0為初始點(diǎn)，(dx,dy)為速度矢量
那么有
(x0+kdx-x')
2+(y0+kdy-y')
2=r
2
如果delta<0，那么說(shuō)明直線與圓沒(méi)有交點(diǎn)
否則取k=(-b-sqrt(delta))/2a
如果k<=0，則說(shuō)明交點(diǎn)在射線的反方向或者交點(diǎn)就在原來(lái)的位置，不滿足題意
第一個(gè)問(wèn)題算比較好的解決了。
第二個(gè)問(wèn)題我覺(jué)得比較簡(jiǎn)單的方法是利用向量的旋轉(zhuǎn)矩陣
藍(lán)色的向量和紅色的向量在知道交點(diǎn)后是很容易算出的。cosθ利用向量夾角公式也不是問(wèn)題，由于θ肯定是銳角，sinθ也不是問(wèn)題了。
然后就可以利用旋轉(zhuǎn)矩陣將紅色的向量右旋得到目標(biāo)速度向量（不一定是右旋，如果光線從紅色向量右側(cè)射入，那么就要左旋了。判斷左右可以用叉積或者干脆直接將兩個(gè)向量都求出，選擇與原來(lái)那個(gè)不同的向量即可）。
旋轉(zhuǎn)矩陣(左旋)：
[cosθ -sinθ]
[sinθ cosθ]
有一點(diǎn)要注意：精度卡到1e-8就夠了，求sin的時(shí)候sqrt(1-cos*cos+
eps)，不然可能會(huì)出現(xiàn)-0.000000001這種情況。。。
代碼：
1
# include <cstdio>
2# include <cmath>
3# include <cstdlib>
4using namespace std;
5# define eps 1e-8
6struct circle
7
{
8
double x,y,r;
9
}cir[30];
10int n,id[30];
11double tmp[30][4];
12double dx,dy,x0,y00;
13inline double dis(double x1,double y1,double x2,double y2)
14{
15 return (((x1)-(x2))*((x1)-(x2))+((y1)-(y2))*((y1)-(y2)));
16}
17inline bool less(double n1,double n2)
18{
19 return (fabs((n1)-(n2))>eps&&(n1)<(n2));
20}
21inline bool gt(double n1,double n2)
22{
23 return (fabs((n1)-(n2))>eps&&(n1)>(n2));
24}
25void cal(circle &pos,int &count,int _id)
26{
27 double b=2.0*((x0-pos.x)*dx+(y00-pos.y)*dy),a=dx*dx+dy*dy,c=(x0-pos.x)*(x0-pos.x)+(y00-pos.y)*(y00-pos.y)-pos.r*pos.r;
28 double delta=b*b-4.0*a*c;
29 if(less(delta,0)) return;
30 else
31
{
32 double k=(-b-sqrt(delta))/2.0/a;
33 if(!gt(k,0)) return;
34 id[count]=_id;
35 tmp[count][0]=x0+dx*k;
36 tmp[count][1]=y00+dy*k;
37 double x2=tmp[count][0]-pos.x,y2=tmp[count][1]-pos.y;//連心線向量
38 double x1=x0-tmp[count][0],y1=y00-tmp[count][1];//原向量
39 double cs=(x1*x2+y2*y1)/sqrt(x1*x1+y1*y1)/sqrt(x2*x2+y2*y2);//計(jì)算cos
40 double sn=sqrt(1-cs*cs);//計(jì)算sin
41 if(less(x1*y2-x2*y1,0))//點(diǎn)積判斷左右情況
42 x1=cs*x2+sn*y2,y1=-sn*x2+cs*y2;
43 else if(gt(x1*y2-x2*y1,0))
44 x1=cs*x2-sn*y2,y1=sn*x2+cs*y2;
45 else x1=x2,y1=y2;
46 k=sqrt(x1*x1+y1*y1);
47 x1/=k;
48 y1/=k;
49 tmp[count][2]=x1;
50 tmp[count][3]=y1;
51 count++;
52
}
53}
54int main()
55{
56 //freopen("reflect.in","r",stdin);
57 //freopen("ans.txt","w",stdout);
58 int count=1;
59 while(true)
60 {
61 scanf("%d",&n);
62 if(!n) break;
63 for(int i=0;i<n;i++)
64 scanf("%lf%lf%lf",&cir[i].x,&cir[i].y,&cir[i].r);
65 scanf("%lf%lf%lf%lf",&x0,&y00,&dx,&dy);
66 printf("Scene %d\n",count++);
67 for(int t=1;t<=11;t++)
68 {
69 int c=0;
70 for(int i=0;i<n;i++)
71 cal(cir[i],c,i);
72 if(c==0)
73 {
74 printf("inf\n\n");
75 break;
76 }
77 else if(t!=11)
78 {
79 double len=1e16;
80 int target;
81 for(int i=0;i<c;i++)
82 {
83 if(less(dis(x0,y00,tmp[i][0],tmp[i][1]),len))
84 len=dis(x0,y00,tmp[i][0],tmp[i][1]),target=i;
85 }
86 printf("%d ",id[target]+1);
87 x0=tmp[target][0];
88 y00=tmp[target][1];
89 dx=tmp[target][2];
90 dy=tmp[target][3];
91 // printf("%.4f %.4f %.4f %.4f\n",x0,y00,dx,dy);
92 }
93 else printf("\n\n");
94 }
95 }
96 return 0;
97}
98
99
100