POJ 3004 && HDU 1922 Subway planning

题目大意:

平面直角坐标系中有一些点代表着一些城市。国家要在(0,0)点设置中心车站并向建设地铁线路,但要求线路是直线。某个城市可以使用地铁线路的前提是它与线路的直线距离不超过d。问最少建设多少条线路才能让所有的城市都能使用地铁。


解题思路:

扫描线+最小区间覆盖。

POJ 3004 && HDU 1922 Subway planning_第1张图片

如图所示:对于任何一个点,都有一个铁路允许的设立的角度区间,在d固定的情况下,点距离原点越远区间越小。

因为是一个圆圈,所以要枚举任一点为区间覆盖的起点。


下面是代码:

#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 
#include 

#define eps 1e-8
#define pi acos(-1.0)
#define inf 107374182
#define inf64 1152921504606846976
#define lc l,m,tr<<1
#define rc m + 1,r,tr<<1|1
#define iabs(x)  ((x) > 0 ? (x) : -(x))
#define clear1(A, X, SIZE) memset(A, X, sizeof(A[0]) * (SIZE))
#define clearall(A, X) memset(A, X, sizeof(A))
#define memcopy1(A , X, SIZE) memcpy(A , X ,sizeof(X[0])*(SIZE))
#define memcopyall(A, X) memcpy(A , X ,sizeof(X))
#define max( x, y )  ( ((x) > (y)) ? (x) : (y) )
#define min( x, y )  ( ((x) < (y)) ? (x) : (y) )


using namespace std;
struct node1
{
    double x,y,l,r;
    bool operator <(const node1 a)const
    {
        return l0)now=0.5*pi;
                    else now=1.5*pi;
                }
                else if(Point[i].y==0)
                {
                    if(Point[i].x>0)now=0;
                    else now=pi;
                }
                else
                {
                    now=asin(iabs(Point[i].y)/temp);
                    if(Point[i].x<0&&Point[i].y>=0)now=pi-now;
                    else if(Point[i].x<=0&&Point[i].y<0)now+=pi;
                    else if(Point[i].x>0&&Point[i].y<=0)now=2*pi-now;
                }
                drift=asin(d/temp);
                //printf("*%lf %lf\n",now,drift);
                Point[i].l=now-drift;
                Point[i].r=now+drift;
            }
        }
        if(n==0)
        {
            puts("0");
            continue;
        }
        sort(Point,Point+n);
        for(int i=n; i<2*n; i++)
        {
            Point[i]=Point[i-n];
            Point[i].l+=2*pi;
            Point[i].r+=2*pi;
        }
        ans=inf;
        for(int i=0; iPoint[j].r)r=Point[j].r;
                else if(r


转载于:https://www.cnblogs.com/lin375691011/p/3996618.html

你可能感兴趣的:(POJ 3004 && HDU 1922 Subway planning)