几何HDU5858

这里只我的
思路说一下。

思路：给出正方形边长L。BD为圆的半径为L/2，BA为正方形对角线的1/2，AD为1/4圆AEF的半径为L，则能根据余弦定理求出∠a，进而求出∠b，进而求出扇形ACD的面积以及三角形ABD的面积。用扇形ACD的面积减去两倍的三角形ABD的面积为饼BCD的面积。然后根据三角形ABD的三边求得∠c，进而求得∠d，进而求得以B为圆心的扇形BCD的面积，然后用扇形BCD的面积减去饼BCD的面积即为阴影部分的面积，然后乘以2即可。

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<algorithm>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<map>
#include<string>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f;
const ll INFF=0x3f3f3f3f3f3f3f3f;
const double pi=acos(-1.0);
const double eps=1e-9;
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        double l;
        scanf("%lf",&l);
        double a=l/2;
        double b=l*sqrt(2)/2;
        double c=l;
        double angle=(b*b+c*c-a*a)/(2*b*c);
        double p=(a*a+b*b-c*c)/(2*a*b);
        double angle_d=2*pi-2*(acos(p));
        double s1=l*l*angle_d/8;
        double s3=b*c*sin(acos(angle))/2;
        double s2=acos(angle)*l*l;
        double ans=s1-(s2-2*s3);
        printf("%.2lf\n",ans*2);
    }
    return 0;
}