【POJ】2318 TOYS（计算几何+暴力）

http://poj.org/problem?id=2318

第一次完全是$O(n^2)$的暴力为什么被卡了～QAQ（一定是常数太大了...）

后来排序了下点然后单调搞了搞。。（然而还是可以随便造出让我的code变成$O(n^2)$的23333）

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define pii pair<int, int>
#define mkpii make_pair<int, int>
#define pdi pair<double, int>
#define mkpdi make_pair<double, int>
#define pli pair<ll, int>
#define mkpli make_pair<ll, int>
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endl
inline const ll getint() { ll r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }

const int N=5015;
struct dat { ll x, y; }E[2], a[N][2], b[N];
ll Cross(dat &a, dat &b, dat &c) {
	static ll x1, x2, y1, y2;
	x1=a.x-c.x; y1=a.y-c.y;
	x2=b.x-c.x; y2=b.y-c.y;
	return x1*y2-x2*y1;
}
int ans[N];

bool check(dat &l1, dat &l2, dat &r1, dat &r2, dat &p) {
	if(Cross(p, l1, l2)*Cross(p, r1, r2)<0) return 1;
	return 0;
}

int n, m;
void work(int last, dat p) {
	int xxx=n+1;
	for1(now, last, xxx) {
		if(check(a[now-1][0], a[now-1][1], a[now][0], a[now][1], p)) { ++ans[now-1]; return; }
	}
}

bool cmp(const dat &a, const dat &b) { return a.x<b.x; }

int main() {
	while(read(n), n) {
		read(m);
		rep(i, 2) read(E[i].x), read(E[i].y);
		a[0][0].x=E[0].x; a[0][0].y=E[0].y;
		a[0][1].x=E[0].x; a[0][1].y=E[1].y;
		a[n+1][0].x=E[1].x; a[n+1][0].y=E[0].y;
		a[n+1][1].x=E[1].x; a[n+1][1].y=E[1].y;
		for1(i, 1, n) rep(k, 2) read(a[i][k].x), a[i][k].y=E[k].y;
		for1(i, 1, m) read(b[i].x), read(b[i].y);
		sort(b+1, b+1+m, cmp);
		int last=1;
		for1(i, 1, m) {
			while(last<=n && max(a[last][0].x, a[last][1].x)<b[i].x) ++last;
			work(last, b[i]);
		}
		rep(i, n+1) printf("%d: %d\n", i, ans[i]), ans[i]=0;
		puts("");
	}
	return 0;
}

5 6 0 10 60 0
3 1
4 3
6 8
10 10
15 30
1 5
2 1
2 8
5 5
40 10
7 9
4 10 0 10 100 0
20 20
40 40
60 60
80 80
 5 10
15 10
25 10
35 10
45 10
55 10
65 10
75 10
85 10
95 10
0

0: 2
1: 1
2: 1
3: 1
4: 0
5: 1

0: 2
1: 2
2: 2
3: 2
4: 2