CSP——M4

A——TT数鸭子

题意：
题目描述 这一天，TT因为疫情在家憋得难受，在云吸猫一小时后，TT决定去附近自家的山头游玩。
TT来到一个小湖边，看到了许多在湖边嬉戏的鸭子，TT顿生羡慕。此时他发现每一只鸭子都不 一样，或羽毛不同，或性格不同。TT在脑子里开了一个map<鸭子，整数> tong，把鸭子变成了 一些数字。现在他好奇，有多少只鸭子映射成的数的数位中不同的数字个数小于k。
输入描述
输入第一行包含两个数n,k，表示鸭子的个数和题目要求的k。
接下来一行有n个数,aiai a_iai​，每个数表示鸭子被TT映射之后的值。
输出描述
输出一行，一个数，表示满足题目描述的鸭子的个数。
无行末空格
样例输入
6 5
123456789 9876543210 233 666 1 114514
样例输出
4

思路
这题比较简单，开一个字符数组存数，对于每一个数（字符串),记录0-9数字出现的个数即可。

代码

#include<iostream>
#include<cstring>
#include<string> 
using namespace std;
string s[1000010];
int num[10];
int n,k;
int ans=0,ans1;
int main(){
	cin>>n>>k;
	for(int i=1;i<=n;i++){
		cin>>s[i];
		memset(num,0,sizeof(num));
		ans1=0;
		for(int j=0;j<s[i].size();j++){
			int idx=s[i][j]-'0';
			num[idx]++;
		}
		for(int l=0;l<=9;l++){
			if(num[l]!=0)
			ans1++;
		} 
		if(ans1<k)
		ans++;
	}
	cout<<ans<<endl;
} 

B——ZJM要抵御宇宙射线

题目描述:
据传，2020年是宇宙射线集中爆发的一年，这和神秘的宇宙狗脱不了干系！但是瑞神和东东忙 于正面对决宇宙狗，宇宙射线的抵御工作就落到了ZJM的身上。假设宇宙射线的发射点位于一个 平面，ZJM已经通过特殊手段获取了所有宇宙射线的发射点，他们的坐标都是整数。而ZJM要构 造一个保护罩，这个保护罩是一个圆形，中心位于一个宇宙射线的发射点上。同时，因为大部分 经费都拨给了瑞神，所以ZJM要节省经费，做一个最小面积的保护罩。当ZJM决定好之后，东东 来找ZJM一起对抗宇宙狗去了，所以ZJM把问题扔给了你~
输入描述: 输入 第一行一个正整数N，表示宇宙射线发射点的个数
接下来N行，每行两个整数X,Y，表示宇宙射线发射点的位置
输出描述: 输出包括两行
第一行输出保护罩的中心坐标x,y 用空格隔开
第二行输出保护罩半径的平方
（所有输出保留两位小数，如有多解，输出x较小的点，如扔有多解，输入y较小的点）
无行末空格
样例输入
5
0 0
0 1
1 0
0 -1
-1 0
样例输出:
0.00 0.00
1.00

思路
这道题由于n很小，直接暴力求出每个点到其他点的距离
即可，复杂度O（n^2),求出每个点到其他点的最远距离中的最小的即可，其中排序的目的是因为题意有多解的情况返回x，y字典序最小的。

代码

#include <bits/stdc++.h>
using namespace std;
const double inf=DBL_MAX;
struct point{
 double x;
 double y;
}p[1010];
double len(point a,point b){
 return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
}
bool cmp(point a,point b){
	return a.x!=b.x? a.x<b.x:a.y<b.y;
}
int n; 
int main(){
 cin>>n;
 for(int i=0;i<n;i++){
  cin>>p[i].x>>p[i].y;
 }
 sort(p,p+n,cmp);
 double mind=inf;
 int ans;
 for(int i=0;i<n;i++){
  double maxd=0;
  for(int j=0;j<n;j++){
    maxd=max(maxd,len(p[i],p[j]));
  }
  if(maxd<mind){
   mind=maxd;
   ans=i;
  }
 }
 printf("%.2lf %.2lf\n%.2lf",p[ans].x,p[ans].y,mind);
 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);}
给了辗转相除的函数真的贴心了。
输入描述
输入第一行一个t，表示数据组数。
对于每组数据，第一行输入一个n，表示数的个数
接下来一行有n个数aiai a_iai​，输入保证是升序的。
输出描述
每组数据输出一行，如果能够造出来满足题目描述的树，输出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

思路
本题运用动态规划的方法。两个状态Left[i][j]和Right[i][j]，分别表示"[i-1,j]能成为i的左子树"、"[i+1,j]能成为i的右子树"。二维数组GCD[i][j]表示ai和aj的最小公约数大于一。
状态转移过程如下：对每一个区间[i,j]，检查区间[i,j]中是否存在节点k，如果满足Left[i,k]&&Right[k,j]&&GCD[i+1][k] 那么Left[i+1][j]=1，如果Left[i,k]&&Right[k,j]&&GCD[j-1][k]，那么Right[j-1][i]=1。最后检查Left[k][1]&&Right[k][n]。
核心代码为

for(int i=1;i<=n;i++){
for(int j=i;j>=1;j--){
for(int k=j;k<=i;k++){
if(Left[k][j]&&Right[k][i]){
if(GCD[i+1][k]）
Left[i+1][j]=1;
if(GCD[j-1][k])
Right[j-1][i]=1;
					}
				}
			}
		}

代码为

#include<iostream>
#include<string.h>
using namespace std;

int t,n,a[710];
bool flag;
bool GCD[710][710];
bool Left[710][710],Right[710][710];

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

int main(){
	cin>>t;
	while(t--){
		memset(Left,0,sizeof(Left));
		memset(Right,0,sizeof(Right));
		flag=false;
		cin>>n;
		for(int i=1;i<=n;i++){
			cin>>a[i];
			GCD[i][i]=false;
			Left[i][i]=Right[i][i]=true;
		}
		for(int i=1;i<=n;i++){
			for(int j=1;j<=n;j++){			
					if(gcd(a[i],a[j])>1)
					GCD[i][j]=1;
					else 
					GCD[i][j]=0;
			}
		}
		for(int i=1;i<=n;i++){
			for(int j=i;j>=1;j--){
				for(int k=j;k<=i;k++){
					if(Left[k][j]&&Right[k][i]){
						if(GCD[i+1][k])
						Left[i+1][j]=1;
						if(GCD[j-1][k])
						Right[j-1][i]=1;
					}
				}
			}
		}
		for(int i=1;i<=n;i++){
			if(Left[i][1]&&Right[i][n]){
				flag=1;
				break;
			}
		}
		if(flag)
		cout<<"Yes"<<endl;
		else 
		cout<<"No"<<endl;
	}
}