程序设计思维与实践 CSP-M4

TT数鸭子

题意

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

解题思路

暴力对每个数的每一位进行判断，计算出不同数字个数即可。

代码

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
int n, k;
int cntall;
void diff(long long temp){
	if(k > 10){
		cntall++;
		return;
	}
	bool cnt[10];
	memset(cnt, 0, sizeof(cnt));
	int cnt1 = 0;
	while(temp){
		int a = temp % 10;
		if(!cnt[a]){
			cnt[a] = true;
			cnt1++;
			if(cnt1 >= k)
				return;
		}
		temp /= 10;
	}
	cntall++;
}
int main(){
	ios::sync_with_stdio(false);
	cin >> n >> k;
	for(int i = 0; i < n; ++i){
		long long temp;
		cin >> temp;
		diff(temp);
	}
	cout << cntall << endl;
	return 0;
}

ZJM要抵御宇宙射线

题意

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

解题思路

依然是道暴力，将每个点都枚举为圆心，计算到其余所有点的距离的最大值即为该圆的半径，取所有圆中半径的最小值（注意平方）

代码

#include <iostream>
#include <cstdio>
#include <cmath>
#include <algorithm>
using namespace std;
int n;
double ans = 1e11;
int mini = -1;
struct P{
	double x, y;
}point[1010];
bool cmp(P a, P b){
	if(a.x != b.x) return a.x < b.x;
	return a.y < b.y;
}
int main() {
	cin >> n;
	for(int i = 0; i < n; ++i)
		cin >> point[i].x >> point[i].y;
	sort(point, point + n, cmp);
	for(int i = 0; i < n; ++i){
		double max1 = 0;
		for(int j = 0; j < n; ++j){
			if(j == i) continue;
			double t = 1.0 * (point[j].x - point[i].x) * (point[j].x - point[i].x) + 1.0 * (point[j].y - point[i].y) * (point[j].y - point[i].y);
			max1 = max(t, max1);
		}
		if(max1 < ans){
			mini = i;
			ans = max1;
		}
	}
	printf("%.2lf %.2lf\n", point[mini].x, point[mini].y);
	printf("%.2lf\n", ans);
    return 0;
}

宇宙狗的危机

题意

在瑞神大战宇宙射线中我们了解到了宇宙狗的厉害之处，虽然宇宙狗凶神恶煞，但是宇宙狗有一个很可爱的女朋友。
最近，他的女朋友得到了一些数，同时，她还很喜欢树，所以她打算把得到的数拼成一颗树。
这一天，她快拼完了，同时她和好友相约假期出去玩。贪吃的宇宙狗不小心把树的树枝都吃掉了。所以恐惧包围了宇宙狗，他现在要恢复整棵树，但是它只知道这棵树是一颗二叉搜索树，同时任意树边相连的两个节点的gcd(greatest common divisor)都超过1。
但是宇宙狗只会发射宇宙射线，他来请求你的帮助，问你能否帮他解决这个问题。

解题思路

这是一道区间dp的问题，我们枚举区间$[l,r]$，对任意一个区间$[l,r]$来说,其上一层的父亲节点必然是$l-1$或$r+1$
对于每一个区间，枚举一个根节点，将区间分为左半部分和右半部分两个子区间，那么这两个子区间的父亲节点即为我们枚举出来的这个点。
我们用$dp[i][j][0]$表示区间$[l,r]$构成的子树是否有一种方案满足其根节点的值与$a_{l-1}$不互质
$dp[i][j][1]$表示区间$[l,r]$构成的子树是否有一种方案满足其根节点的值与$a_{r+1}$不互质
转移就是枚举$k(l\le k\le r)$
如果$dp[l][k-1][1]=1$且$dp[k+1][r][0]=1$且$a_k$与$a_{l-1}$互质，则$dp[l][r][0]=1$;
如果$dp[l][k-1][1]=1$且$dp[k+1][r][0]=1$且$a_k$与$a_{r+1}$互质，则$dp[l][r][1]=1$;
防止重复计算gcd可以预处理所有节点两两间的gcd存在一个二维数组中。

代码

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
int a[705], dp[705][705][2], GCD[705][705];
int gcd(int a,int b){return b == 0 ? a : gcd(b,a%b);}
int main() {
	int T;
	cin >> T;
	while(T--){
		int n;
		cin >> n;
		memset(dp, 0, sizeof(dp));
		for(int i = 1; i <= n; i++)
			cin >> a[i];
		for(int i = 1; i <= n; i++)
			for(int j = 1; j <= n; j++)
				GCD[i][j] = gcd(a[i], a[j]);
		for(int i = 1; i <= n; i++){
			if(GCD[i][i + 1] != 1){
				dp[i + 1][i + 1][0] = 1;
				dp[i][i][1] = 1;
			}
		}
		for(int i = 1; i < n; i++){
			for(int l = 1; l + i <= n; l++){
				int r = l + i;
				for(int k = l; k <= r; k++){
					if((k == l || dp[l][k-1][1]) && (dp[k+1][r][0] || k == r)){
						if(l != 1 && GCD[l - 1][k] != 1)
							dp[l][r][0] = 1;
						if(r != n && GCD[k][r + 1] != 1)
							dp[l][r][1] = 1;
					}
				}
			}
		}
		bool flag = false;
		if(dp[2][n][0] || dp[1][n-1][1]) flag = true;
		for(int i = 2; i < n - 1; i++){
			if(dp[1][i - 1][1] && dp[i + 1][n][0])
				flag = true;
		}
		if(!flag) cout << "No\n";
		else cout << "Yes\n";
	}
    return 0;
}