ZOJ3720 Magnet Darts(点在多边形内)

第一道点在多边形内的判断题，一开始以为是凸的。其实题意很简单的啦，最后转化为判断一个点是否在一个多边形内。

如果只是简单的凸多边形的话，我们可以枚举每条边算下叉积就可以知道某个点是不是在范围内了。但对于更一般的多边形，要采用的方法可能就要复杂一些，其中一种是叫做射线法的东西，在该点上画一条射线，如果射线与多边形偶数次相交，则说明点在多边形外，否则在多边形内。实践中可以取水平的一条，至于怎么判断就要看下代码了，判断的过程中如何发现点已经在多边形上了则直接返回。

计算几何做的少，则是算是补充了一个盲区吧。

#pragma warning(disable:4996)

#include<iostream>

#include<cstring>

#include<cstdio>

#include<algorithm>

#include<cmath>

#define ll long long

#define mod 1000000007

#define maxn 1000000

#define eps 1e-6

using namespace std;



struct Point

{

	int x, y;

	Point(int xi, int yi) :x(xi), y(yi){}

	Point(){}

	int dot(const Point &b){

		return x*b.x + y*b.y;

	}

	int det(const Point &b){

		return x*b.y - y*b.x;

	}

	Point operator - (const Point &b){

		return Point(x - b.x, y - b.y);

	}

}p[25];



int dcmp(double x){

	return (x > eps) - (x < -eps);

}



int n;

double A, B;

double lx, rx, ly, ry;



bool onLine(Point x, Point aa, Point bb)

{

	return (bb - aa).det(x-aa) == 0 && (bb - x).dot(aa - x) <= 0;

}



bool inside(Point x)

{

	int ret = 0;

	for (int i = 0; i < n; i++){

		if (onLine(x, p[i], p[(i + 1) % n])) return true;

		int k = (p[(i + 1) % n] - p[i]).det(x - p[i]);

		int d1 = x.y - p[i].y;

		int d2 = x.y - p[(i + 1) % n].y;

		if (k > 0 && d1 >= 0 && d2 < 0) ret++;

		if (k < 0 && d2 >= 0 && d1 < 0) ret--;

	}

	return ret != 0;

}



int main()

{

	while (cin>>lx>>ly>>rx>>ry){

		scanf("%d%lf%lf", &n, &A, &B);

		for (int i = 0; i < n; i++){

			scanf("%d%d", &p[i].x, &p[i].y);

		}

		p[n].x = p[0].x; p[n].y = p[0].y;

		double ans = 0;

		int llx = ceil(lx), rrx = floor(rx);

		int lly = ceil(ly), rry = floor(ry);



		for (int i = llx; i <= rrx; i++){

			for (int j = lly; j <= rry; j++){

				if (inside(Point(i, j))){

					double h = min(i + 0.5, rx) - max(i - 0.5, lx);

					double w = min(j + 0.5, ry) - max(j - 0.5, ly);

					ans += (A*i + B*j)*h*w;

				}

			}

		}

		ans /= (double)(rx - lx)*(ry - ly);

		printf("%.3lf\n", ans);

	}

	return 0;

}