2016HUAS_ACM暑假集训2L - Points on Cycle（圆上的点）

一个简单的几何题，自己在纸上列出方程解出结果的表达式，再用程序表达出来就行了。

不过老司机（老司机的woodcoding）说用旋转向量法比较简单，有时间要去看一看。

大致题意：一个圆心在原点的圆，半径未知，现在给你圆上的一点，让你在这个圆上找到另外两点，使得这三点构成的三角形的周长最长。

样例输入：（第一行为一个整数N，表示后面有N组案例，每个案例给出一组圆上点的坐标）

2

1.500        2.000

563.585    1.251

样例输出：（其他两个点的坐标）

0.982 -2.299 -2.482 0.299

-280.709  -488.704  -282.876  487.453

主要思路：首先，一个圆的内接三角形周长最长的是正三角形。设给出的为点A(x,y)，圆的半径为R，连接AO，反向延长AO与圆交于点O'，以O'为圆心，R为半径画圆，两个圆的交点坐标就是所求的答案。

记O'(a,b)，（a=-x,b=-y）目的只是为了防止混淆。

联立x²+y²=R² ① 和 (x-a)²+(y-b)²=R² ②，然后自己打草稿即可。

 #include<stdio.h>
 #include<math.h>
 int T;
 int main()
 {
     double x,y,a,b,x1,y1,x2,y2,R,A,B,C;//R是圆半径的平方
     scanf("%d",&T);
     while(T--)
     {
         scanf("%lf %lf",&x,&y);
         a=-x;b=-y;//(x,y)的对称点(a,b)
         R=x*x+y*y;
         A=R;
         B=-R*b;
         C=R*R/-a*a*R;
         y1=(-B-sqrt(B*B-*A*C))/(*A);
         y2=(-B+sqrt(B*B-*A*C))/(*A);//求出y1,y2后还不能马上算x1,x2，防止x作除数
         if(x==)
         {
             x1=-sqrt(R-y1*y1);
             x2=sqrt(R-y2*y2);
         }
         else
         {
             x1=(R/-b*y1)/a;
             x2=(R/-b*y2)/a;
         }
         printf("%.3lf %.3lf %.3lf %.3lf\n",x1,y1,x2,y2);
     }
     return ;
 }