ural(Timus) 1333. Genie Bomber 2

几何题

题意：给出n个圆的圆心坐标和半径，给出一个单位方格，在（0,0），（0,1），问这些圆覆盖的面积占方格的百分比，圆超出方格的部分不计算。另外算百分比，答案精确到1%即可，即整数部分正确即可，小数部分不要求

这题，可以想按照题意直接下手，未免太难了，另外注意到答案的输出，其实对精度的要求很低（对小数都没要求）。我们可以用一直近似的算法来解决

将方格分割为一个一个的小格子，当分割得很小的时候，格子可以看做是一个点，然后看这个点在不在圆上或圆内，在的话相当于圆覆盖了这个点，覆盖了这个格子。所以我们将1*1方格分割为1000*1000的格子，然后逐一去判断，很暴力的方法，但就是这样过了

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;

const double PI = acos(-1.0);
const double E = 1e-8;
const int N = 15;
const int M = 1000;
const double B = 0.001;
struct crc
{
    double x,y,r;
}a[N];
int n;

double dis(double x ,double y , struct crc p)
{
    return sqrt( (x-p.x)*(x-p.x) + (y-p.y)*(y-p.y));
}

int judge(double x , double y)
{
    for(int i=0; i<n; i++)
    {
        double d = dis(x,y,a[i]);
        if( d < a[i].r || fabs(d - a[i].r) < E)
            return 1;
    }
    return 0;
}

int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        int res = 0;
        for(int i=0; i<n; i++)
            scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].r);
        for(int i=0; i<=M; i++)
            for(int j=0; j<=M; j++)
            {
                double x = B * 1. * i;
                double y = B * 1. * j;
                int ok = judge(x,y);
                res += ok;
            }
        printf("%.6lf\n",1.*res/(1.*M*M)*100.0);
    }
    return 0;
}