2020牛客暑期多校训练营（第二场）K

题目描述

Given three concentric circles whose radiuses are $r_1, r_2, r_3$ respectively, and $A,B,C$ are the moving points on the given three circles respectively. Determine the expected area of $\triangle ABC$ .

输入描述

The first line contains one integer $T$ $(1 \leqslant T \leqslant 1000)$ , denoting the number of test cases.
For each test case, one line containing three integers $r_1, r_2, r_3$ $(1\leqslant r_1,r_2,r_3 \leqslant 100)$ , denoting the radiuses of three given concentric circles.

输出描述

Print $T$ lines each containing one real number with one decimal places after the decimal point, denoting the answer to corresponding test case.
It's guaranteed that the second decimal place after the decimal point is neither $4$ nor $5$ .

示例1

输入

2
1 1 1
2 3 5

输出

0.5
5.5

说明

For the first text case, the accurate answer is $\frac{3}{2\pi}=0.47746482927568600730665129011754\cdots$ .

分析

由于 $S_{\triangle ABC}$ 只与 $A,B,C$ 三点的相对位置有关，可以将 $A$ 视作顶点， $B,C$ 为动点。不妨以同心圆的圆心作为原点建立平面直角坐标系，令 $A(r_1,0)$ ， $B(r_2\cos\alpha,r_2\sin\alpha)$ ， $C(r_3\cos\beta,r_3\sin\beta)$ ； $\alpha,\beta\in\mathbb R$ ，且 $0\leqslant\alpha,\beta\leqslant 2\pi$ 。
利用向量的叉积计算面积

$S(\alpha,\beta)=\frac{1}{2}|\overrightarrow{AB}\times\overrightarrow{AC}|=\frac{1}{2}|(r_2\cos\alpha-r_1)r_3\sin\beta-(r_3\cos\beta-r_1)r_2\sin\alpha|$

设出现极角 $\alpha,\beta$ 的概率为 $f(\alpha,\beta)$ ，则 $S_{\triangle ABC}$ 的期望为

$E= \int_0^{2\pi}\int_0^{2\pi}S(\alpha,\beta)f(\alpha,\beta)d\alpha d\beta$

此题精度要求较低，不妨将 $2\pi$ 分成 $t$ 份，每份的角度为 $\frac{2\pi}{t}$ ，那么 $f(\alpha,\beta)=\frac{1}{t^2}$ 。可直接用矩形法求积分的近似数值解，设置步长为 $\frac{2\pi}{t}$ ，枚举 $[0,2\pi]$ 内的所有角度即可。

代码

/******************************************************************
Copyright: 11D_Beyonder All Rights Reserved
Author: 11D_Beyonder
Problem ID: 2020牛客暑期多校训练营（第二场） Problem K
Date: 8/4/2020
Description: Expectation and Integral
*******************************************************************/
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int t=400;
const double pi=acos(-1);
double _sin[t+1],_cos[t+1];
int main()
{
    int i,j,_;
    const double step=2*pi/t;
    double theta=0;
    for(i=1;i<=t;i++)
    {
        _sin[i]=sin(theta);
        _cos[i]=cos(theta);
        theta+=step;
    }
    for(cin>>_;_;_--)
    {
        double r1,r2,r3;
        scanf("%lf%lf%lf",&r1,&r2,&r3);
        //=====================
        //排序
        if(r1>r2) swap(r1,r2);
        if(r2>r3) swap(r2,r3);
        if(r1>r3) swap(r1,r3);
        //=====================
        //矩形法 
        //枚举 t*t 个角度组合
        double ans=0;
        for(i=1;i<=t;i++)
        {
            for(j=1;j<=t;j++)
            {
                ans+=fabs((r2*_cos[i]-r1)*r3*_sin[j]-(r3*_cos[j]-r1)*r2*_sin[i]);
            }
        }
        printf("%.1lf\n",ans/2/t/t);
    }
    return 0;
}