【CSP - M4】A - TT数鸭子、B - ZJM要抵御宇宙射线、C - 宇宙狗的危机

A - TT数鸭子

题意：

这一天，TT因为疫情在家憋得难受，在云吸猫一小时后，TT决定去附近自家的山头游玩。
TT来到一个小湖边，看到了许多在湖边嬉戏的鸭子，TT顿生羡慕。此时他发现每一只鸭子都不一样，或羽毛不同，或性格不同。TT在脑子里开了一个map<鸭子，整数> tong，把鸭子变成了一些数字。现在他好奇，有多少只鸭子映射成的数的数位中不同的数字个数小于k。

Input

输入第一行包含两个数n,k，表示鸭子的个数和题目要求的k。
接下来一行有n个数,$a_i$，每个数表示鸭子被TT映射之后的值。

Output

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

Sample Input

6 5
123456789 9876543210 233 666 1 114514

Sample Output

4

思路做法：

将数字按字符串形式读入，遍历字符串计算数位上不同数字个数，小于k就将答案加一

总结:

题目的范围里包含k很大的情况，当k大于10时，输入的数字个数即为答案，因为cnt < k才算答案数，做题时多用了k的判断，但是将>10写成了>=10就错了，因为k==10时可允许9个不同的数，但是总数为0~9共10个数。

代码:

#include <stdio.h>
#include <string.h>
bool occur[15];
const int N = 20;
char s[N];
int n, k;
bool isOk(){
	for(int j = 0; j < 10; ++j) occur[j] = false;
	int len = strlen(s), cnt = 0;
	for(int j = 0; j < len; ++j){
		if(!occur[s[j]-'0']){
			cnt++; occur[s[j]-'0'] = true;
		}
		if(cnt >= k) return false;
	}
	return true;
}
int main(){
	int ans = 0; scanf("%d%d", &n, &k);
	for(int i = 0; i < n; ++i){
		scanf("%s", s);
		if(k > 10) ans++;
		if(isOk()) ans++;
	}
	printf("%d\n", ans);
	return 0;
}

B - ZJM要抵御宇宙射线

题意：

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

Input

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

Output

输出包括两行
第一行输出保护罩的中心坐标x,y 用空格隔开
第二行输出保护罩半径的平方
（所有输出保留两位小数，如有多解，输出x较小的点，如扔有多解，输入y较小的点）
无行末空格

Sample Input

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

Sample Output

0.00 0.00
1.00

思路做法：

数据量不大，可以遍历每个点，求最大的半径可以包括所有点，再不停更新最小的半径即可

总结:

T2没那么复杂，而且要是中心不在发射点上的情况我也不会

代码:

#include <stdio.h>
#include <algorithm>
using namespace std;
#define ll long long
const int N = 1e3+20;
struct Point{
	int x, y;
	Point(){}
	Point(int _x, int _y) : x(_x), y(_y) {}
	bool operator < (const Point& p) const {
		if(x != p.x) return x < p.x;
		return y < p.y;
	}
}pos[N];
int n; 
ll calcu(int j){
	ll R = 0;
	for(int i = 0; i < n; ++i){
		if(i == j) continue;
		ll r = (1ll)*(pos[j].x-pos[i].x)*(pos[j].x-pos[i].x)+(1ll)*(pos[j].y-pos[i].y)*(pos[j].y-pos[i].y);
		if(r > R) R = r;
	}
	return R;
}
int main(){
	scanf("%d", &n);
	for(int i = 0; i < n; ++i){
		scanf("%d%d", &pos[i].x, &pos[i].y);
	}
	sort(pos, pos+n);
	int x, y;
	ll R = 1e14;
	for(int i = 0; i < n; ++i){
		ll r = calcu(i);
		if(r < R){
			R = r; x = pos[i].x; y = pos[i].y;
		}
	}
	printf("%.2f %.2f\n%.2f\n", 1.0*x, 1.0*y, 1.0*R);
	return 0;
}

C - 宇宙狗的危机

题意：

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

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

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

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

补充知识：

GCD：最大公约数，两个或多个整数共有约数中最大的一个 ，例如8和6的最大公约数是2。
一个简短的用辗转相除法求gcd的例子：

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

Input

输入第一行一个t，表示数据组数。
对于每组数据，第一行输入一个n，表示数的个数
接下来一行有n个数$a_i$，输入保证是升序的。

Output

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

Sample Input

1
6
3 6 9 18 36 108

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

Sample Output

Yes

No
Yes

思路做法：

区间DP，从小区间到大区间不停更新，由于二叉搜索树对于左子树和右子树的大小关系严格，因此用2个二维数组le和ri分别记录可不可以成树的情况，$le[i][j]$表示一个bool量，代表j为根时的左子树能不能到i，类似的，$ri[i][j]$代表i为根时的右子树能不能到j。对于一个较大的区间而言，$dp[i][j]$代表题意中的yes或no，它可以由更小的子区间更新而来，即遍历k，$i<=k<=j$，若$le[i][k]\&\&ri[k][j]$为真，那么2个小区间自然可以合并为1个大区间，同时le和ri也要更新，对于$dp[i][j]$为真的区间，遍历k，其中$i<=k<=j$，若$gcd(a[i-1],a[k])>1$即可以a[i-1]为根，右子树为[i, j]构建二叉搜索树，另一个更新方式同理。最后观察$dp[0][n-1]$的取值情况

总结:

做题时没想到区间DP这一点，只是用了gcd>1的条件构建了图，但是没有考虑二叉搜索树左右子树的限制

代码:

#include <stdio.h>
const int N = 7e2+20;
int a[N], gcd_value[N][N];
bool le[N][N], ri[N][N], dp[N][N];
/*
区间dp  dp[i][j]  [i, j]是否满足
左右子树必须划分出来 le[i][j] ri[i][j]
小区间到大区间 dp[i][j] 左 gcd[i-1][k] 右 gcd[k][j+1] i<=k<=j
*/
int gcd(int a,int b) { return b == 0 ? a : gcd(b, a % b); }
int main(){
	int t; scanf("%d", &t);
	while(t--){
		int n; scanf("%d", &n);
		for(int i = 0; i < n; ++i){
			scanf("%d", &a[i]);
		}
		for(int i = 0; i < n; ++i){
			for(int j = 0; j < n; ++j){
				gcd_value[i][j] = gcd(a[i], a[j]);
				gcd_value[j][i] = gcd_value[i][j];
				le[i][j] = 0; ri[i][j] = 0; dp[i][j] = 0;
			}
		}
		for(int i = 0; i < n; ++i){
			le[i][i] = 1; ri[i][i] = 1;
		}
		for(int len = 1; len <= n; ++len){  // 区间大小 
			for(int i = 0; i <= n-len+1; ++i){
				int j = i+len-1;
				for(int k = i; k <= j; ++k){
					if(le[i][k] && ri[k][j]){
						dp[i][j] = 1;
						if(i > 0 && gcd_value[i-1][k] > 1)
							ri[i-1][j] = 1;
						if(j < n-1 && gcd_value[k][j+1] > 1)
							le[i][j+1] = 1; 
					}
				}
			}
		}
		if(dp[0][n-1]) printf("Yes\n");
		else printf("No\n");
	}
	return 0;
}