fzu 2144 Shooting Game(贪心区间覆盖)

题目链接：fzu 2144 Shooting Game

题目大意：在立体的空间上，有n只蚊子，给出蚊子的坐标，以及蚊子的移动方向（向量的方式给出）， 然后人的攻击范围为半径为r的球体。问说人最多可以打几只蚊子，以及需要花多少时间。

解题思路：对于每只蚊子计算进入攻击范围的时间和离开攻击范围的时间，注意有可能不进入攻击空间，计算区间可以设一个k然后解方程求出区间，然后对求出的区间做区间覆盖问题。

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>

using namespace std;

const int N = 100005;
const double eps = 0;

inline bool scan_d(double &ret) {  
	char c; int sgn;  
	if(c = getchar(),c == EOF) return 0; //EOF  
	while(c != '-' && (c < '0' || c > '9')) c = getchar();  
	sgn = (c == '-') ? -1 : 1;  
	ret = (c == '-') ? 0 : (c - '0');  
	while(c=getchar(),c>='0'&&c<='9') ret=ret*10+(c-'0');  
	ret *= sgn;  
	return 1;  
} 

struct point {
	double x, y, z;
	void get() { scan_d(x); scan_d(y); scan_d(z); }
};

struct state {
	double x, y;
}s[N];
int n, tmp;
double r;

void add(point p, point q) {
	double a = q.x * q.x + q.y * q.y + q.z * q.z;
	double b = 2 * (p.x * q.x + p.y * q.y + p.z * q.z);
	double c = p.x * p.x + p.y * p.y + p.z * p.z - r * r;
	double t = b * b - 4 * a * c;

	if (t < 0) return;
	else if (t == 0) {
		double i = -b / (2.0 * a);
		if (i < 0) return;
		s[tmp].x = s[tmp].y = i;
	}
	else {
		double i = (-b + sqrt(t) ) / (2 * a);
		double j = (-b - sqrt(t) ) / (2 * a);
		if (i < 0 && j < 0) return;
		else if (i < 0) i = 0;
		else if (j < 0) j = 0;

		s[tmp].x = min(i, j);
		s[tmp].y = max(i, j);
	}
	tmp++;
}

void init() {
	tmp = 0;
	point p, q;

	scanf("%d%lf", &n, &r);
	for (int i = 0; i < n; i++) {
		p.get(); q.get();
		add(p, q);
	}
}

bool cmp(const state& a, const state& b) {
	if (fabs(a.y - b.y) > eps) return a.y - b.y < eps;
	return a.x - b.x < eps;
}

int solve() {
	int ans = 0;
	double c = -1;
	for (int i = 0; i < tmp; i++) {
		if (s[i].x - c > eps) {
			c = s[i].y;
			ans++;
		}
	}
	return ans;
}

int main () {
	int cas;
	scanf("%d", &cas);
	for (int i = 1; i <= cas; i++) {
		init();
		sort(s, s + tmp, cmp);
		printf("Case %d: %d %d\n", i, tmp, solve());
	}
	return 0;
}