HDU 2469 Fire-control System解题报告

这题是一题计算几何题，第一次做计算几何题，出了很多错，WA了半天，最后搜到PROXIMA的解题代码，读懂后仿照着写了下，终于AC了，这题主要的难点在于要从点下手，而不是从角度下手

1、两条边扫一圈的时候，最好以点为前进单位，而不要以角度为单位。

    虽然这两种方法是等价的，但是后者需要先把角度相同的点都找出来，还要统计各个角度上分别有几个点，这要花费一些初始化的时间。然而，在最差情况下，每个点所在的角度都可以是不同的，那么以角度为前进单位就完全没有必要了，还白白浪费了初始化的时间。

2、最好把所有点复制一份，把第二份的所有点的角度都加上 2π。

    也就是说，下标为 0 ~ N-1 的点是第一圈，下标为 N ~ 2N-1 的点是第二圈。这样一来，当扫完了一圈回到开头时，下标不必回到 0，直接走到 N 继续下去就好了。这可以节省很多 MOD 运算的时间。#include<iostream> #include<algorithm> #include<cmath> const int MAXNUM=20000; const double PI=atan2(0.0,-1.0); using namespace std; int N,K; struct Point { int x,y; double angle,dist; }ps[MAXNUM]; double minangle; bool found[MAXNUM]; int cmp(Point a,Point b) { if(a.angle==b.angle) return a.dist<b.dist; else return a.angle<b.angle; } void calMin() { if(K==0||K==1) { minangle=0.00; return; } minangle=1e100; int end,scount,start; int size=N+N; double newm; for(int i=0;i<N;i++) { end=-1; scount=0;start=0; for(int j=0;j<size;j++) if(ps[j].dist<=ps[i].dist) { scount++; if(scount==K) { end=j; break; } } if(end!=-1) { scount--;end--; for(int q=end+1;q<size;q++) { if(ps[q].dist<=ps[i].dist) { scount++; while(ps[start].dist>ps[i].dist)start++; if(scount>K) { scount--; start++; while(ps[start].dist>ps[i].dist)start++; } } if((q-start+1)<=N) { newm=(ps[q].angle-ps[start].angle)*ps[i].dist; if(newm<minangle) minangle=newm; } } } } } int main() { int casenum=1; freopen("c://myin.txt","r",stdin); scanf("%d%d",&N,&K); while(N||K) { for(int i=0;i<N;i++) { scanf("%d%d",&ps[i].x,&ps[i].y); ps[i].angle=atan2((double)(ps[i].y),(double)(ps[i].x)); ps[i].dist=ps[i].x*ps[i].x+ps[i].y*ps[i].y; //printf("%d,%d的角度是%.2f/n",ps[i].x,ps[i].y,ps[i]); } sort(ps,ps+N,cmp); for(int j=0;j<N;j++)//将所有点复制一份副本以便节省扫描使用%的开销 { ps[N+j]=ps[j]; ps[N+j].angle+=PI+PI; } //for(int u=0;u<N;u++) //printf("%.2f ",ps[u].angle); //printf("/n"); calMin(); printf("Case #%d: %.2lf/n",casenum++,minangle/2.0); scanf("%d%d",&N,&K); } return 0; }