zoukankan      html  css  js  c++  java
  • poj 1106 Transmitters (叉乘的应用)

    http://poj.org/problem?id=1106

    Time Limit: 1000MS   Memory Limit: 10000K
    Total Submissions: 4488   Accepted: 2379

    Description

    In a wireless network with multiple transmitters sending on the same frequencies, it is often a requirement that signals don't overlap, or at least that they don't conflict. One way of accomplishing this is to restrict a transmitter's coverage area. This problem uses a shielded transmitter that only broadcasts in a semicircle. 

    A transmitter T is located somewhere on a 1,000 square meter grid. It broadcasts in a semicircular area of radius r. The transmitter may be rotated any amount, but not moved. Given N points anywhere on the grid, compute the maximum number of points that can be simultaneously reached by the transmitter's signal. Figure 1 shows the same data points with two different transmitter rotations. 

    All input coordinates are integers (0-1000). The radius is a positive real number greater than 0. Points on the boundary of a semicircle are considered within that semicircle. There are 1-150 unique points to examine per transmitter. No points are at the same location as the transmitter. 

    Input

    Input consists of information for one or more independent transmitter problems. Each problem begins with one line containing the (x,y) coordinates of the transmitter followed by the broadcast radius, r. The next line contains the number of points N on the grid, followed by N sets of (x,y) coordinates, one set per line. The end of the input is signalled by a line with a negative radius; the (x,y) values will be present but indeterminate. Figures 1 and 2 represent the data in the first two example data sets below, though they are on different scales. Figures 1a and 2 show transmitter rotations that result in maximal coverage.

    Output

    For each transmitter, the output contains a single line with the maximum number of points that can be contained in some semicircle.

    Sample Input

    25 25 3.5
    7
    25 28
    23 27
    27 27
    24 23
    26 23
    24 29
    26 29
    350 200 2.0
    5
    350 202
    350 199
    350 198
    348 200
    352 200
    995 995 10.0
    4
    1000 1000
    999 998
    990 992
    1000 999
    100 100 -2.5

    Sample Output

    3
    4
    4

    Source

     
     
    -----------------------------------------------------------------------------
    思维不敏捷啊,还是后悔看了题解,不解释。自己先想想吧,就是叉乘的应用
     1 #include <stdio.h>
     2 #include <string.h>
     3 #include <stdlib.h>
     4 #include <iostream>
     5 #include <algorithm>
     6 #include <math.h>
     7 #define eps 1e-6
     8 typedef struct point
     9 {
    10     double x,y;
    11 }point;
    12 
    13 bool  dy(double x,double y){ return x>y+eps; }
    14 bool  xy(double x,double y){ return x<y-eps; }
    15 bool dyd(double x,double y){ return x>y-eps; }
    16 bool xyd(double x,double y){ return x<y+eps; }
    17 bool  dd(double x,double y){ return fabs(x-y)<eps; }
    18 
    19 double crossProduct(point a,point b,point c)
    20 {
    21     return (c.x-a.x)*(b.y-a.y)-(c.y-a.y)*(b.x-a.x);
    22 }
    23 
    24 double dist(point a,point b)
    25 {
    26     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
    27 }
    28 
    29 point c[1005];
    30 double st,en,ri;
    31 point tmp;
    32 int solve(int n)
    33 {
    34         int ans;
    35         int maxx=0;
    36         for(int i=0;i<n;i++)
    37         {
    38             ans=1;
    39             for(int j=0;j<n;j++)
    40             {
    41                 if(i!=j&&dyd(crossProduct(tmp,c[i],c[j]),0.0))
    42                 {
    43                        ans++;
    44                 }
    45             }
    46             if(ans>maxx)
    47             {
    48                 maxx=ans;
    49                 //ans=0;
    50             }
    51         }
    52         return maxx;
    53 }
    54 
    55 int main()
    56 {
    57     int n;
    58     double a,b;
    59     while(scanf("%lf%lf%lf",&st,&en,&ri)!=EOF&&ri>=0)
    60     {
    61         point p;
    62         tmp.x=st;
    63         tmp.y=en;
    64         scanf("%d",&n);
    65         int cas=0;
    66         for(int i=0;i<n;i++)
    67         {
    68             scanf("%lf%lf",&p.x,&p.y);
    69             if(xyd(dist(tmp,p),ri))
    70             {
    71                 c[cas++]=p;
    72             }
    73         }
    74         printf("%d
    ",solve(cas));
    75     }
    76 }
  • 相关阅读:
    零售行业解决方案一
    N-Tier Entity Framework开源项目介绍
    软件代码生成之Codesmith模板.netTiers
    ActiveMQ消息队列介绍
    EntityFramework动态多条件查询与Lambda表达式树
    GitLab版本管理
    REST服务介绍
    2014年物联网Internet of Things应用简介
    wordpress如何屏蔽wp-json(禁用REST API)
    bootstrap tab切换如何让鼠标移动自动切换内容
  • 原文地址:https://www.cnblogs.com/ccccnzb/p/3871706.html
Copyright © 2011-2022 走看看