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]);
}