TOJ4701 求阴影部分面积

能不能不爬我博客啊，MMP
设正方形ABCD边长为a，分别以B和D为圆心，绘制半径为a的1/4圆，与正方形的内切圆相交与M、N、O、P点，求这些交点围成的阴影部分面积。

Input

输入数据有多组，每组输入一个正整数a(a<=20)，表示正方形的边长，输入以EOF结束。

Output

每组输出一个小数，表示阴影部分的面积，保留为2位小数

Sample Input

Sample Output

Hint
圆周率PI取acos(-1.0)

求三角形面积
{
double a=1.0,b=2.0,c=sqrt(2.0),s,area;
s=(a+b+c)/2;
area=sqrt(s*(s-a)(s-b)(s-c));
printf(“area=%.12f”,area);
}
求一份的面积
{double PI=acos(-1.0);
double a1=acos(1.25/sqrt(2.0));
double a2=acos(-0.5/sqrt(2.0))-PI/2.0;
double s=a1*2+a2/2-0.661437827766;
printf(“%.12f”,s);
}

#include <stdio.h>
int main()
{