中超8 hdu 7063 Square Card（计算几何，两圆交集面积，圆内正方形）

hdu 7063 Square Card

题意：

• 平面上有两个圆，一个圆得分，一个圆奖励

  现将一个边长为a的正方形扔到平面上，会绕得分圆的中心旋转

  若某一时刻，正方形严格处在圆里，将得分 or 奖励

  求同时得分和得奖励 与 得分的概率比

分析：

• 首先是要仔细读题，给我是读傻了

  然后，就是转换了

  要把概率比转换成面积比

  考虑一个正方形完全落入圆内，最极端的情况即两点在圆上，可求出正方形中心到圆心的距离 $r_a=\sqrt{r^2-(\frac{a}{2}{)}^2}-\frac{a}{2}$​ ，正方形的中心落在半径为r_a的圆形区域即得分，这样分别求出在两个圆内的有效区域，同时得分和得奖励的有效区域便是两个求出圆的交集部分（上板子）

• 注意：题目只保证了得分区域够大，并没有保证得奖励的区域，所以得奖励的区域还要特判一下

#include<bits/stdc++.h>
#define P pair<double, double>
using namespace std;
const double pi = acos(-1);
//求两圆相交面积
double area(P a, double r1, P b, double r2)
{
	double d = sqrt((a.first-b.first)*(a.first-b.first) + (a.second-b.second)*(a.second-b.second));
	if(d >= r1+r2) return 0;
	if(r1 > r2) swap(r1,r2);
	if(r2-r1>=d) return pi*r1*r1;
	double ang1 = acos((r1*r1+d*d-r2*r2)/(2*r1*d));
	double ang2 = acos((r2*r2+d*d-r1*r1)/(2*r2*d));
	return ang1*r1*r1 + ang2*r2*r2 - r1*d*sin(ang1);
}

int main()
{
    int T;  
    cin>>T;
    while(T--)
    {
        double ra,xa,ya,rb,xb,yb,a;
		cin>>ra>>xa>>ya>>rb>>xb>>yb>>a;
		if(rb*rb-a*a/4 < 0){ cout<<"0.000000"<<endl; continue; }
        P ap=P(xa,ya), bp=P(xb,yb);
		double r1 = sqrt(ra*ra-a*a/4)-a/2, r2 = sqrt(rb*rb-a*a/4)-a/2;
		double res = area(ap, r1, bp, r2);
		printf("%.6f\n", res/(pi*r1*r1));
    }
    return 0;
}