HDU 5858 Hard problem 2016多校10 1002 [相交圆面积]【计算几何】

最新推荐文章于 2018-04-20 10:07:17 发布

Tabris_

最新推荐文章于 2018-04-20 10:07:17 发布

阅读量688

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.youkuaiyun.com/qq_33184171/article/details/52251839

===== 模板 ===== 同时被 2 个专栏收录

15 篇文章

订阅专栏

==== 计算几何 ====

15 篇文章

订阅专栏

该博客介绍了如何解决一道涉及计算几何的数学问题，具体是求解正方形内切圆与四分之一圆相交形成的阴影部分面积。博主分享了解题思路，从最初尝试割补法失败，到后来利用相交圆面积公式解决问题。题目要求给出正方形边长并保留两位小数。博主提醒，等比例放大图形时，面积会按比例平方增长，并提供了相应代码。

题目连接：就当前HDU的情况就先不给了

——————————————————————–.
Hard problem

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 195 Accepted Submission(s): 148

Problem Description
cjj is fun with math problem. One day he found a Olympic Mathematics problem for primary school students. It is too difficult for cjj. Can you solve it?
这里写图片描述

Give you the side length of the square L, you need to calculate the shaded area in the picture.

The full circle is the inscribed circle of the square, and the center of two quarter circle is the vertex of square, and its radius is the length of the square.

Input
The first line contains a integer T(1<=T<=10000), means the number of the test case. Each case contains one line with integer l(1<=l<=10000).

Output
For each test case, print one line, the shade area in the picture. The answer is round to two digit.

Sample Input
1
1

Sample Output
0.29

———————————————-.
题目大意：就是正方形边长为n时求阴影面积保留小数后两位

题目分析：
当天不想做多校然后队友给看了这题说是小学奥数题目所以就一直割补法然后 GG
今天偶然看了看过去整理的模板猛然看到一个叫做相交圆面积的东西然后就想到了这道题目。。

然后大致是这样的。
这里写图片描述

小的圆是蓝的打的圆是红的绿色部分是相交的面积（红色+蓝色==绿色？就当是对的了：不可能错）

然后根据相交圆面积求绿色部分最后拿小圆的面积减去它最后在乘个二就行了。。。

之后做题的时候只求一次边长为1的阴影部分面价就好了因为平面图形等比放大N倍面积就是原来的n*n倍

就是这么水然而被第一句话误导了一直在寻找巧妙的办法…
什么小学奥数题目。。。小学能学这个？？！！！

附本题代码
———————————————–.

#include <bits/stdc++.h>
using namespace std;
#define lalal puts("*******");
const double PI  =  acos(-1);

struct Round {
    double x, y;
    double r;
}c1,c2;

double dis(Round a, Round b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

double solve(Round a, Round b)
{
    double d = dis(a, b);
    if (d >= a.r + b.r)  return 0;

    if (d <= fabs(a.r - b.r))
    {
        double r = a.r < b.r ? a.r : b.r;
        return PI * r * r;
    }

    double ang1 = acos((a.r * a.r + d * d - b.r * b.r) / 2. / a.r / d);
    double ang2 = acos((b.r * b.r + d * d - a.r * a.r) / 2. / b.r / d);
    double ret = ang1 * a.r * a.r + ang2 * b.r * b.r - d * a.r * sin(ang1);
    return ret;
}

int main()
{
    c1.x=0.0,c1.y=0.0,c1.r=1.0;
    c2.x=0.5,c2.y=0.5,c2.r=0.5;

    double s = solve(c1,c2);
    s = PI*0.5-s*2;
  //  printf("%.15lf\n",s);

    int _,n;
    scanf("%d",&_);
    while(_--)
    {
        scanf("%d",&n);
        printf("%.2lf\n",s*n*n);
    }
    return 0;
}