「个人训练」Radar Installation（POJ-1328）

这条题目A了十次。。。emmmmm 其实不难就是一个贪心。。。。

先说下算法（之前的和现在的）

之前考虑的其实很简单。用平面几何即可将雷达可以放置的区域转化为区间（顺便判断是否无解。问题就比较简单了：先排序。然后我们第一个区间我们肯定是放在区间的最右边，然后对后面的我们如何处理呢？若后一个区间的左边在目前区间的右边，我们肯定是不用动这个雷达的摆放位置的；如果不满足，那么我们就要看看后一个区间的右边和这个区间的摆放位置了。如果后一个区间的右边在目前区间的右边的左边，那么我们就应当放在后一个区间的右边。按照这个思路就可以了。

再谈下自己错误的地方

1. 浮点数，我前面几次全部都是用int保存的，省略了很多精度，导致了错误还死活调不出来。（这个应该是主要错误。）
2. 这题有大量i/o，应当采用c-style的输入输出或者关掉tie。
3. 各种代码细节问题，这个是编程基本功，以后慢慢提升吧。

代码

原题的代码套用了一个网上的模板，看起来有点恐怖233

#include <cctype>
#include <numeric>
#include <climits>
#include <cmath>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <string>
#include <sstream>
#include <algorithm>
#include <complex>
#include <deque>
#include <list>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <vector>
#include <utility>
#include <bitset>
using namespace std;

#define REP(i,n) for(int i=0;i!=static_cast<int>(n);++i)
#define REP_1(i,n) for(int i=1;i<=static_cast<int>(n);++i)
#define FOR(i,a,b) for(int i=static_cast<int>(a);i<static_cast<int>(b);++i)
#define FOR_1(i,a,b) for(int i=static_cast<int>(a);i<=static_cast<int>(b);++i)
#define DWN(i,b,a) for(int i=static_cast<int>(b);i>static_cast<int>(a);--i)
#define DWN_1(i,b,a) for(int i=static_cast<int>(b);i>=static_cast<int>(a);--i)
#define ECH(it,A) for(auto it=A.begin();it!=A.end();++it)

#define REP_2(i,j,n,m) REP(i,n) REP(j,m)
#define REP_2_1(i,j,n,m) REP_1(i,n) REP_1(j,m)
#define REP_3(i,j,k,n,m,l) REP(i,n) REP(j,m) REP(k,l)
#define REP_3_1(i,j,k,n,m,l) REP_1(i,n) REP_1(j,m) REP_1(k,l)
#define ALL(A) A.begin(), A.end()
#define LLA(A) A.rbegin(), A.rend()
#define CPY(A,B) memcpy(A,B,sizeof(A))
#define INS(A,P,B) A.insert(A.begin()+P,B)
#define ERS(A,P) A.erase(A.begin()+P)
#define BSC(A,X) find(ALL(A),X)// != A.end()
#define CTN(T,X) (T.find(X)!=T.end())

#define SZ(A) static_cast<int>(A.size())
#define PB push_back
#define MP(A,B) make_pair(A,B)
#define P(T) pair<T,T>
#define fi first
#define se second
#define D(v) (cout<< #v <<" is equal to "<<v<<endl)
#define DP(v) (cout<< #v << " is equal to "<<v.fi<<" "<<v.se<<endl)
#define ZERO(v) memset(v,0,sizeof(v))
#define MIO(v) memset(v,-1,sizeof(v))

#define IOS std::ios::sync_with_stdio(false); std::cin.tie(nullptr); std::cout.tie(nullptr);

const int inf=0x3f3f3f3f;
//const int MAXN;
const double eps=1e-6;
const double pi=acos(-1.0);
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
inline int sgn(double x)
{
    return (x>eps)-(x<-eps);
}
int d;
struct Interval
{
    double l,r;
    Interval(double _l=0,double _r=0):l(_l),r(_r) {}
    Interval(P(double) p):l(p.fi-sqrt(d*d-p.se*p.se)),r(p.fi+sqrt(d*d-p.se*p.se)) {}
    bool operator < (const Interval& rhs) const
    {
        if (r==rhs.r) return l<rhs.l;
        else return r<rhs.r;
    }
};

int main()
{
    int n;
    int kase=0;
    while(scanf("%d%d",&n,&d)==2)
    {
        if(!n&&!d) break;
        int ans=1;
        bool isVaild=true;
        if(d<0)//d能小于0
            isVaild=false;
        vector<Interval> vec;
        REP(i,n)
        {
            double x,y; scanf("%lf%lf",&x,&y);
            if(double(d)<fabs(y)) { isVaild=false;ans=-1; }
            else vec.PB(MP(x,y));
            //cout<<vec[i].l<<" "<<vec[i].r<<endl;
        }

        if(ans!=-1)
        {
            sort(ALL(vec));
            int tmp=0;
            int len=vec.size();
            for(int i=1;i!=len;++i)
            {
                if(vec[i].l>vec[tmp].r)
                {
                    ans++;
                    tmp=i;
                }
                else if(vec[i].r<vec[tmp].r)
                    tmp=i;
            }
        }
        printf("Case %d: %d\n",++kase,ans);
    }
    return 0;
}