zoukankan      html  css  js  c++  java
  • 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[250060];
    struct cir
    {
        node Mid;
        double r;
    }Cir[250060];
    struct squ
    {
        node a,b;
    }Squ[250060];
    int tot1=0,tot2=0;
    int ans[250060];
    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[10];
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            scanf("%s",c);
            if(c[0]=='c')
                scanf("%lf %lf %lf",&Cir[tot1].Mid.x,&Cir[tot1].Mid.y,&Cir[tot1].r),tot1++;
            if(c[0]=='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=1;i<=m;i++)
            scanf("%lf %lf",&Node[i].x,&Node[i].y),Node[i].z = i;
        sort(Node+1,Node+1+m);
    
        for(int i=0;i<tot1;i++)
        {
            int l = upper_bound(Node + 1,Node + m + 1,node{Cir[i].Mid.x - Cir[i].r,0.0,0.0}) - Node;
            int r = lower_bound(Node + 1,Node + m + 1,node{Cir[i].Mid.x + Cir[i].r,0.0,0.0}) - Node - 1;
            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=0;i<tot2;i++)
        {
            int l = upper_bound(Node + 1,Node + m + 1,node{Squ[i].a.x ,0.0,0.0}) - Node;
            int r = lower_bound(Node + 1,Node + m + 1,node{Squ[i].b.x ,0.0,0.0}) - Node - 1;
            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=1;i<=m;i++)
            printf("%d
    ",ans[i]);
    }
  • 相关阅读:
    VMware虚拟机下网络配置模式
    2021考研规划(持续更新)
    20210326 名词解释及常用下载地址(持续更新)
    20210326学习笔记1---java及hadoop组件最新版本汇总
    20210326日记
    20210325一天总结--进步最大的一天
    20210325学习感悟--学习是开始几小时烦躁,越学习越上瘾;暴食看剧打游戏打飞机,最开始爽,时间越长越厌倦甚至感觉痛苦。
    20210326继续解决----20210325学习笔记2--运行MapReduce Jar(我为什么这样起标题,因为结构化数据才好搜索)
    20210325学习笔记1--解决了打包不生成jar文件的问题
    20210325日记--加油,相信只要基础够扎实熟练,就能找到转职成功。
  • 原文地址:https://www.cnblogs.com/qscqesze/p/4814168.html
Copyright © 2011-2022 走看看