算法-程序设计课程week16-CSP-M4

ZJM要抵御宇宙射线

描述

据传，2020年是宇宙射线集中爆发的一年，这和神秘的宇宙狗脱不了干系！但是瑞神和东东忙于正面对决宇宙狗，宇宙射线的抵御工作就落到了ZJM的身上。假设宇宙射线的发射点位于一个平面，ZJM已经通过特殊手段获取了所有宇宙射线的发射点，他们的坐标都是整数。而ZJM要构造一个保护罩，这个保护罩是一个圆形，中心位于一个宇宙射线的发射点上。同时，因为大部分经费都拨给了瑞神，所以ZJM要节省经费，做一个最小面积的保护罩。当ZJM决定好之后，东东来找ZJM一起对抗宇宙狗去了，所以ZJM把问题扔给了你~

输入
输入 第一行一个正整数N，表示宇宙射线发射点的个数

接下来N行，每行两个整数X,Y，表示宇宙射线发射点的位置

输出
输出包括两行

第一行输出保护罩的中心坐标x,y 用空格隔开

第二行输出保护罩半径的平方

（所有输出保留两位小数，如有多解，输出x较小的点，如扔有多解，输入y较小的点）

无行末空格

输入样例 1

5
0 0
0 1
1 0
0 -1
-1 0

输出样例 1

0.00 0.00
1.00

思路

1. 以每个点为圆心，计算最长距离，这个距离就是该点为圆心的圆的半径。
2. 看这些半径哪个小，最小半径可以形成最小圆。

心得

暴力解法，这个我确实想不出来，这题目是什么类型？考数学吗？
我一开始还真以为这是求最小覆盖圆来着…没想到只是一个暴力。

代码

#include <bits/stdc++.h>
#define LL long long

using namespace std;
struct point {
	int x, y;
	bool operator<(const point &p) const {
		if (x != p.x)
			return x < p.x;
		return y < p.y;
	}
}


int n;
LL ans_r = LONG_LONG_MAX;
int ans_i;
point a[1005];

LL getMaxR(int i) {
	LL res = 0;
	for (int j = 0; j < n; j++) {
		if (i == j) {
			continue;
		}
		LL tmp = (LL)(a[i].y - a[j].y) * (a[i].y - a[j].y) + (LL)(a[i].x - a[j].x) * (a[i].x - a[j].x);
		res = max(tmp, res);
	}
	return res;
}

int main() {
	cin >> n;
	for (int i = 0; i < n; i++){
		 cin >> a[i].x >> a[i].y;
	}
	sort(a, a + n);
	for (int i = 0; i < n; i++) {
		LL temp_r = getMaxR(i);
		if (temp_r < ans_r) {
			ans_r = temp_r;
			ans_i = i;
		}
	}
	printf("%.2f %.2f\n", (float)a[ans_i].x, (float)a[ans_i].y);
	printf("%.2lf", (double)ans_r);
	return 0;
}

TT数鸭子

描述

这一天，TT因为疫情在家憋得难受，在云吸猫一小时后，TT决定去附近自家的山头游玩。

TT来到一个小湖边，看到了许多在湖边嬉戏的鸭子，TT顿生羡慕。此时他发现每一只鸭子都不一样，或羽毛不同，或性格不同。TT在脑子里开了一个map<鸭子，整数> tong，把鸭子变成了一些数字。现在他好奇，有多少只鸭子映射成的数的数位中不同的数字个数小于k。

输入
输入第一行包含两个数n,k，表示鸭子的个数和题目要求的k。

接下来一行有n个数,a_ia
i
​
，每个数表示鸭子被TT映射之后的值。

输出
输出一行，一个数，表示满足题目描述的鸭子的个数。

无行末空格

输入样例 1

6 5
123456789 9876543210 233 666 1 114514

输出样例 1

4

思路

读入字符串，统计字符串中数字出现次数，再统计出现了多少类数字。

心得

这个题巨简单，但是我一开始不理解为什么一共10个阿拉伯数字，K确能大于10！
真的是 起 夜 级 李 姐
这里卡住了，我根本没做这道题…

代码

#include <bits/stdc++.h>
using namespace std;
int counti[10];
char s[20];

int main() {
	int n, k;
	int sum = 0;
	scanf("%d %d", &n, &k);
	for (int i = 1; i <= n; i++) {
		memset(counti, 0, sizeof(counti));	//每次清零
		scanf("%s", s);
		int len = strlen(s);
		for (int j = 0; j < len; j++) {
			counti[s[j] - '0']++;  //标记 出现的数
		}
		int total = 0;
		for (int j = 0; j < 10; j++) {
			if (counti[j]) {
				total++;
			}
		}
		if (total < k) {
			sum++;
		}
	}

	printf("%d\n", sum);
	return 0;
}

宇宙狗的危机

描述

在瑞神大战宇宙射线中我们了解到了宇宙狗的厉害之处，虽然宇宙狗凶神恶煞，但是宇宙狗有一个很可爱的女朋友。

最近，他的女朋友得到了一些数，同时，她还很喜欢树，所以她打算把得到的数拼成一颗树。

这一天，她快拼完了，同时她和好友相约假期出去玩。贪吃的宇宙狗不小心把树的树枝都吃掉了。所以恐惧包围了宇宙狗，他现在要恢复整棵树，但是它只知道这棵树是一颗二叉搜索树，同时任意树边相连的两个节点的gcd(greatest common divisor)都超过1。

但是宇宙狗只会发射宇宙射线，他来请求你的帮助，问你能否帮他解决这个问题。

补充知识：

GCD：最大公约数，两个或多个整数共有约数中最大的一个 ，例如8和6的最大公约数是2。

一个简短的用辗转相除法求gcd的例子：

int gcd(int a,int b){return b == 0 ? a : gcd(b,a%b);}

输入
输入第一行一个t，表示数据组数。

对于每组数据，第一行输入一个n，表示数的个数

接下来一行有n个数​，输入保证是升序的。

输出
每组数据输出一行，如果能够造出来满足题目描述的树，输出Yes，否则输出No。

无行末空格。

输入样例 1

1
6
3 6 9 18 36 108

输出样例 1

Yes

输入样例 2

2
2
7 17
9
4 8 10 12 15 18 33 44 81

输出样例 2

No
Yes

思路

这道题目的题面很简单，需要数组符合以下要求：

1. 可构成一颗二叉搜索树；
2. 任意树边相连的两个节点的gcd(greatest common divisor)都超过1。

要用区间dp来做，用动态规划的方式构建这棵树，dp定义如下：

m[i][j]：i到j能不能构造成二叉搜索树。
GCD[i][j]：两个数gcd是不是大于1
l[i][j]：i可以是j的左子树
r[i][j]：j可以是i的右子树

依次枚举区间长度len，m中长度为len的区间，根节点，转移函数为

if (lt[l][root] && rt[root][r]) {
	ok[l][r] = 1;
	if (GCD[l - 1][root]) {
		rt[l - 1][r] = 1;
	}
	if (GCD[root][r + 1]) {
		lt[l][r + 1] = 1;
	}
}

心得

1. 小写的gcd编译会报错，看来是STL保留字

代码

#include <bits/stdc++.h>
#define LL long long
using namespace std;

int t, n;
int a[700 + 5], ok[700 + 5][700 + 5], GCD[700 + 5][700 + 5];
int lt[700 + 5][700 + 5];
int rt[700 + 5][700 + 5];

int main() {
	cin >> t;
	while (t--) {
		memset(ok, 0, sizeof(ok));
		memset(lt, 0, sizeof(lt));
		memset(rt, 0, sizeof(rt));
		scanf("%d", &n);
		for (int i = 1; i <= n; i++) {
			scanf("%d", &a[i]);
			lt[i][i] = rt[i][i] = ok[i][i] = 1;
			for (int j = 1; j <= n; j++) {
				GCD[i][j] = __gcd(a[i], a[j]) > 1;
			}
		}

		for (int i = 1; i <= n; i++) {							//枚举长度
			for (int l = 1, r = l + i - 1; r <= n; l++, r++) {	//枚举区间
				for (int root = l; root <= r; root++) {			//枚举根节点
					if (lt[l][root] && rt[root][r]) {
						ok[l][r] = 1;
						if (GCD[l - 1][root]) {
							rt[l - 1][r] = 1;
						}
						if (GCD[root][r + 1]) {
							lt[l][r + 1] = 1;
						}
					}
				}
			}
		}
		if (ok[1][n]) {	 //代表从1到n整个序列都可以构成二叉树
			printf("Yes\n");
		} else {
			printf("No\n");
		}
	}
	return 0;
}