BZOJ 1199: [HNOI2005]汤姆的游戏计算几何暴力

1199: [HNOI2005]汤姆的游戏

Time Limit: 1 Sec

Memory Limit: 256 MB

题目连接

http://www.lydsy.com/JudgeOnline/problem.php?id=1199

Description

汤姆是个好动的孩子,今天他突然对圆规和直尺来了兴趣。于是他开始在一张很大很大的白纸上画很多很多的矩形和圆。画着画着,一不小心将他的爆米花弄撒了,于是白纸上就多了好多好多的爆米花。汤姆发现爆米花在白纸上看起来就像一个个点,有些点落在矩形或圆内部,而有些则在外面。于是汤姆开始数每个点在多少个矩形或圆内部。毕竟汤姆还只是个孩子,而且点、矩形和圆又非常多。所以汤姆数了好一会都数不清,于是就向聪明的你求助了。你的任务是：在给定平面上N个图形(矩形或圆)以及M个点后，请你求出每个点在多少个矩形或圆内部(这里假设矩形的边都平行于坐标轴)。

Input

第一行为两个正整数N和M，其中N表示有多少个图形(矩形或圆)，M表示有多少个点。接下来的N行是对每个图形的描述，具体来说，第i+1行表示第i个图形。先是一个字母,若该字母为“r”，则表示该图形是一个矩形,这时后面将有4个实数x1,y1,x2,y2，表示该矩形的一对对角顶点的坐标分别为(x1,y1)和(x2,y2)；若该字母为“c”，则表示该图形是一个圆,这时后面将有3个实数x,y,r,表示该圆以(x,y)为圆心并以r为半径。最后M行是对每个点的描述，其中每行将有两个实数x,y,表示一个坐标为(x,y)的点。

Output

包含M行,每行是一个整数,其中第i行的整数表示第i个点在多少个图形内部(当某点在一个图形的边界上时，我们认为该点不在这个图形的内部)。

Sample Input

3 4
r 1.015 0.750 5.000 4.000
c 6.000 5.000 2.020
r 6.500 7.200 7.800 9.200
3.500 2.500
4.995 3.990
2.300 8.150
6.900 8.000

Sample Output

1
2
0
1

HINT

题意

题目没给取值范围，大概是25W的数据范围

给你n个点，和一些矩形和圆

问你每个点分别在多少个图形内呢

题解：

暴力做就好了，n^2可过……

代码：

#include<math.h>

#include<iostream>

#include<stdio.h>

#include<algorithm>

using namespace std;

struct node

{

    double x,y;

    int z;

    bool operator <(const node &a)const {

        return x < a.x;

    }

}Node[];

struct cir

{

    node Mid;

    double r;

}Cir[];

struct squ

{

    node a,b;

}Squ[];

int tot1=,tot2=;

int ans[];

double dis(node a,node b)

{

    return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);

}

int main()

{

    int n,m;

    char c[];

    scanf("%d%d",&n,&m);

    for(int i=;i<=n;i++)

    {

        scanf("%s",c);

        if(c[]=='c')

            scanf("%lf %lf %lf",&Cir[tot1].Mid.x,&Cir[tot1].Mid.y,&Cir[tot1].r),tot1++;

        if(c[]=='r')

        {

            double x1,y1,x2,y2;

            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);

            Squ[tot2].a.x = min(x1,x2);

            Squ[tot2].a.y = min(y1,y2);

            Squ[tot2].b.x = max(x1,x2);

            Squ[tot2].b.y = max(y1,y2);

            tot2++;

        }

    }

    for(int i=;i<=m;i++)

        scanf("%lf %lf",&Node[i].x,&Node[i].y),Node[i].z = i;

    sort(Node+,Node++m);

    for(int i=;i<tot1;i++)

    {

        int l = upper_bound(Node + ,Node + m + ,node{Cir[i].Mid.x - Cir[i].r,0.0,0.0}) - Node;

        int r = lower_bound(Node + ,Node + m + ,node{Cir[i].Mid.x + Cir[i].r,0.0,0.0}) - Node - ;

        for(int j = l; j <= r ; j++)

        {

            if(dis(Cir[i].Mid,Node[j])<Cir[i].r*Cir[i].r)

                ans[Node[j].z]++;

        }

    }

    for(int i=;i<tot2;i++)

    {

        int l = upper_bound(Node + ,Node + m + ,node{Squ[i].a.x ,0.0,0.0}) - Node;

        int r = lower_bound(Node + ,Node + m + ,node{Squ[i].b.x ,0.0,0.0}) - Node - ;

        for(int j = l; j <= r ; j++)

        {

            if(Node[j].y<Squ[i].b.y&&Squ[i].a.y<Node[j].y)

                ans[Node[j].z]++;

        }

    }

    for(int i=;i<=m;i++)

        printf("%d\n",ans[i]);

}

汤姆