专题十六：模拟 & implementation problem

cf

1.Problem L: Logic Puzzle

这个扫雷的题太漂亮了。每次只看左上角是不是1。是的话就全部减1。最后检查一遍，如果出现非0的情况，那么这个扫雷题图本身就是错误的。

#include<cstdio>
int a[110][110];
char maze[110][110];
void sub(int x, int y) {
	for (int i = -1; i <= 1; i++) {
		for (int j = -1; j <= 1; j++) {
			a[x + i][y + j] -= 1;
		}
	}
}
int main() {
	int N, M;
	scanf("%d%d", &N, &M);
	for (int i = 0; i <= N + 1; i++) {
		for (int j = 0; j <= M + 1; j++) {
			scanf("%d", &a[i][j]);
		}
	}
	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= M; j++) {
			maze[i][j] = '.';
		}
		maze[i][M + 1] = '\0';
	}

	for (int i = 1; i <= N; i++) {
		for (int j = 1; j <= M; j++) {
			if (a[i - 1][j - 1]) {
				maze[i][j] = 'X';
				sub(i, j);
			}
		}
	}
	int flag = 1;
	for (int i = 0; i <= N + 1; i++) {
		for (int j = 0; j <= M + 1; j++) {
			if (a[i][j]) {
				flag = 0;
				break;
			}
		}
	}
	if (!flag) printf("impossible\n");
	else {
		for (int i = 1; i <= N; i++) {
			printf("%s\n", maze[i] + 1);
		}
	}
	return 0;
}

2.Problem G. Invited Speakers

这道题叫构造统一的解题方法。其实，这道题当时应该时想到可以绕圈啊，每次都逆时针绕圈，每次绕的圈都比之前大1，这样一定不会相交的。

• 李鼎懋证明了一个结论，平面上的2n个点，总有一种连线段的方式，是他们两两都不相交。

这里还是用那个绕圈圈的办法，献上代码。

#include<cstdio>
#include<algorithm>
using namespace std;
typedef pair<int, int> P;
P a[10], b[10];
int main() {
	int T;
	scanf("%d", &T);
	while (T--) {
		int N;
		scanf("%d", &N);
		for (int i = 0; i < N; i++) {
			scanf("%d%d", &a[i].first, &a[i].second);
		}
		for (int i = 0; i < N; i++) {
			scanf("%d%d", &b[i].first, &b[i].second);
		}
		sort(a, a + N);
		sort(b, b + N);
		int C = 101;
		for (int i = 0; i < N; i++) {
			printf("6\n");
			printf("%d %d\n", a[i].first, a[i].second);
			printf("%d %d\n", a[i].first, C + i);
			printf("%d %d\n", -C - i, C + i);
			printf("%d %d\n", -C - i, -C - i);
			printf("%d %d\n", b[i].first, -C - i);
			printf("%d %d\n", b[i].first, b[i].second);
		}
	}
	return 0;
}

3.C. Three displays

这道题难度标号1400，但是方法还不是那么的好想。看到 n = 3000，感觉到平方复杂度是可以满足要求的。实际上是这样写的：

• 让 j 从 0 循环到 N - 1，这是第一层循环，然后第二层循环，搜索 i 和 k。
• 只要s[ i ] < s[ j ]，就用min函数更新c1的值。同理，只要s[ j ] < s[ k ]，就更新c2的值。这样可以找到c1 与 c2的最小值。

4.F Finding Lines

• 题意：Given n points in  , is there a straight line passing through at least a p% of them?
• 解法：If line exists it is uniquely determined by any two of its points.
• Repeat 250 times: 1.Pick two distinct points uniformly at random；2.Check if line dened by the points has enough points.
• Pr[false negative] ≈  < 5e-5
• 由于RAND_MAX很有可能只有32767这么小，可以采用相乘的形式进行放大。
• 可以用cstdlib中的srand函数初始化随机数种子。如果需要程序每次执行时使用一个不同的种子，可以用ctime中的time（NULL）为参数调用srand。一般来说，只在程序执行的开头调用一次srand。

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<ctime>
using namespace std;
typedef long long ll;
const int maxn = 100010;
ll x[maxn], y[maxn];
int N, P;
int check(int a, int b) {
	int cnt = 0;
	for (int i = 0; i < N; i++) {
		if ((y[i] - y[a]) * (x[a] - x[b]) == (x[i] - x[a]) * (y[a] - y[b])) {
			cnt++;
		}
	}
	return cnt;
}
int main() {
	scanf("%d%d", &N, &P);
	for (int i = 0; i < N; i++) {
		scanf("%d%d", &x[i], &y[i]);
	}
	bool flag = false;
	int ans = 0;
	if (N <= 2) flag = true;
	else {
		srand(time(NULL));
		
		for (int i = 0; i < 250; i++) {
			int a = (ll)rand() * rand() % N;
			int b = (ll)rand() * rand() % N;
			if (a == b) {
				continue;
			}
			ans = max(ans, check(a, b));
		}
	}
	if ((N * P + 99) / 100 <= ans) flag = true;
	if (flag) printf("possible\n");
	else printf("impossible\n");
	return 0;
}