FZU 2144 —— Shooting Game （贪心）

题意比较简单，就是给N个蚊子，主人公站在（0，0，0）这个位置上，然后N个蚊子有各自的起始点和移动方向，主人公可以在任意时间攻击，攻击的时候那些跟他距离不超过R的蚊子就算被击中，问最多能击中多少个以及最少要攻击几次。

做法就是对每个蚊子，求出它落在攻击范围内的时间段，列个关于时间t的距离不等式求解即可。

然后按照结束时间的顺序排序，贪心一下即可。

#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
struct Event{
    double st, ed;
    bool operator < (const Event& a)const{
        return ed<a.ed;
    }
}ev[100000];
int t, n, i, j, m, cnt;
double r, ax, ay, az, dx, dy, dz, a, b, c, d;
inline void getnum(int& x){
    char ch=getchar();
    bool mk=0;
    x = 0;
    while(ch<48 || ch>57){
        if(ch=='-')    mk=1;
        ch=getchar();
    }
    while(ch>=48 && ch<=57){
        x = x*10+ch-48;
        ch = getchar();
    }
    if(mk)    x=-x;
}
inline void getnum(double& x){
    int y;
    getnum(y);
    x = y*1.0;
}
int main(){
    getnum(t);
    for(int ct=1; ct<=t; ct++){
        m=0;
        getnum(n);
        getnum(r);
        while(n--){
            getnum(ax); getnum(ay); getnum(az);
            getnum(dx); getnum(dy); getnum(dz);
            a = dx*dx+dy*dy+dz*dz;
            b = 2*(ax*dx+ay*dy+az*dz);
            c = ax*ax+ay*ay+az*az-r*r;
            d = b*b-4*a*c;
            if(d<0)    continue;
            ev[m].st = (b*(-1)-sqrt(d))/(a*2);
            ev[m].ed = (b*(-1)+sqrt(d))/(a*2);
            if(ev[m].st<0 && ev[m].ed<0)    continue;//两个同时小于0表明不可能进入攻击范围
            m++;
        }
        printf("Case %d: %d ", ct, m);
        cnt=0;
        sort(ev, ev+m);
        for(i=0; i<m; i++){
            cnt++;
            a = ev[i].ed;
            j=i+1;
            while(j<m){
                if(ev[j].st<=a && ev[j].ed>=a){
                    j++;
                }
                else    break;
            }
            i = j-1;
        }
        printf("%d\n", cnt);
    }
    return 0;
}