计算几何 求相交圆面积(余弦定理)

来源:互联网 发布:我知我师我爱我师ppt 编辑:程序博客网 时间:2024/05/01 02:48

求任意两圆相交的面积
注意PI : 最好用公式表示,否则精度可能不够会一直wa

#include<bits/stdc++.h>  using namespace std;  #define PI acos(-1)  int main()  {      int T;      double x1,x2,y1,y2,r1,r2;      double a,b,s1,s2,d,sum;      cin >> T;      while(T--)      {          sum = 0;          cin >> x1 >> y1 >> r1 >> x2 >> y2 >> r2;          d = sqrt( (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2) );          if(r1 + r2 <= d)              puts("0.00");          else if(r2 >= r1 && d + r1 <= r2)          {              printf("%.2lf\n", PI * r1 * r1);          }           else if(r1 > r2 && d + r2 <= r1)          {              printf("%.2lf\n", PI * r2 * r2);          }          else          {              a = (r1 * r1 + d * d - r2 * r2) / (2.0 * r1 * d);              b = acos(a) * 2.0;              s1 = r1 * r1 * sin(b) / 2;              s2 = b / (2 * PI) * PI * r1 * r1;              sum = s2 - s1;              a = (r2 * r2 + d * d - r1 * r1) / (2.0 * r2 * d);              b = acos(a) * 2.0;              s1 = r2 * r2 * sin(b) / 2;              s2 = b / (2 * PI) * PI * r2 * r2;              sum += s2 - s1;              printf("%.2lf\n", sum);          }      }      return 0;  }  
0 0