Week16-CSP-M4

TT数鸭子

描述

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

输入

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

输出

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

思路：

• 需要判断不同数字的个数，因此直接开10个整数的数组cnt[10]，cnt[i]表示各数位上为i的数的个数。通过遍历string更新cnt数组，最后挑出数组中值>0的个数，与K进行比较即可。

实现：

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;

int cnt[10]={0};
int n,k,res,ans;
string a;

void solve()
{
	memset(cnt,0,sizeof(cnt));
	res = 0;
	int tmp = 0;
	for(int i=0;i<a.size();++i)
	{
		tmp = a[i]-'0';
		cnt[tmp]++;
	}
	for(int i=0;i<=9;++i)
		if(cnt[i]>0)res++;

	if(res < k)ans++;
}


int main()
{
	ios::sync_with_stdio(false);
	cin>>n>>k;
	for(int i=0;i<n;++i)
	{
		cin >> a;
		solve();
	}
	cout<<ans;
	
	return 0;
}

---

ZJM要抵御宇宙射线

描述

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

输入

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

输出

输出包括两行：第一行输出保护罩的中心坐标x,y 用空格隔开；第二行输出保护罩半径的平方

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

思路：

• 因为n <=1000因此n^2复杂度是允许的，所以直接遍历点，对每个点求出到其他点距离平方的最大值，最后进行排序即可。

实现：

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#include<utility>
using namespace std;

int n,x,y;
long long dis,dx,dy;
struct point
{
	double x,y;
	double rad;
	
	bool operator<(point & t)
	{
		if(rad != t.rad)return rad < t.rad;
		if(x != t.x)return x < t.x;
		if(y != t.y)return y < t.y;
	}
};
vector<point> res;
vector<pair<int,int> > p;

point st;
int main()
{
	ios::sync_with_stdio(false);
	cin>>n;
	st.x=st.y=st.rad = -1e9;
	for(int i=0;i<n;++i)
	{
		cin>>x>>y;
		p.push_back({x,y});		
	}
	for(int i=0;i<n;++i)
	{
		for(int j=0;j<n;++j)
		{
			if(i!=j)
			{
				dx = abs(p[j].first-p[i].first);
				dy = abs(p[j].second-p[i].second);
				dis = dx*dx + dy*dy;
				if(dis > st.rad)
					st.rad = dis; 
			}	
		}
		st.x=p[i].first,st.y=p[i].second;
		res.push_back(st);
		st.x=st.y=st.rad = -1e9;
	}	
	sort(res.begin(),res.end());
	printf("%.2f %.2f\n",res[0].x,res[0].y);
	printf("%.2f",res[0].rad);
	
	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。

无行末空格。

思路：

• 注意保证输入是升序的，所以考虑对一段区间[l,r]来说，如果其是子树，那么可以得到其父节点或者是l-1，或者是r+1；如果将root作为根，并且[l,root-1]可以作为root的左子树、[root+1,r]可以作为root的右子树，那么显然[l,r]可以构成符合要求的子树。
• 确定状态：fl[l][r]表示以r为根、[l,r-1]作为左子树是否符合要求；fr[l][r]表示以l为根、[l+1,r]为右子树是否符合要求；
• 对于[l,r]这样的子树，如果尚未覆盖整个区间，那么进行转移：
  • 如果root与l-1相连符合要求，那么l-1作为根、root作为右子树的根符合要求；
  • 如果root与r+1相连符合要求，那么r+1作为根、root作为左子树的根符合要求。
• 初始化过程：对于长度为1的区间，必然可以；对于长度为2的区间，判断gcd大小；

实现：

#include<iostream>
#include<string>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;

int t,n;
int a[710];
int fl[710][710],fr[710][710],dp[710][710];
int gcd(int a,int b){return b == 0 ? a : gcd(b,a%b);} 
bool flag = 0;
void init()
{
	memset(fl,0,sizeof(fl));
	memset(fr,0,sizeof(fr));
	memset(dp,0,sizeof(dp));
	
	for(int i=1;i<=n;++i)
		for(int j=i+1;j<=n;++j)
			if(gcd(a[i],a[j]) > 1)dp[i][j]=dp[j][i]=1;
			
	for(int i=1;i<=n;++i)fl[i][i]=fr[i][i]=1;
	for(int i=2;i<=n;++i)
		if(dp[i-1][i])fl[i-1][i]=fr[i-1][i]=1;
			
	flag = 0;
}
int main()
{
	ios::sync_with_stdio(false);
	cin >> t;
	while(t--)
	{
		cin >> n;
		for(int i=1;i<=n;++i)
		cin>>a[i];
		init();
		
		for(int len=2;len<=n && !flag;++len)
		{
			for(int l=1;l<=n-len+1 && !flag;++l)
			{
				int r = l+len-1;
				for(int k=l;k<=r && !flag;++k)
				{
					if(fl[l][k] && fr[k][r])
					{
						if(l==1 && r==n)
						{
							flag = 1;
							break;
						}
						if(dp[l-1][k])fr[l-1][r]=1;
						if(dp[r+1][k])fl[l][r+1]=1;
					}
				}	
			}
		} 
		if(t!=0)
		{
			if(flag)cout<<"Yes"<<endl;
			else cout<<"No"<<endl;
		}
		else
		{
			if(flag)cout<<"Yes";
			else cout<<"No";	
		}	
	}
	return 0;
}

错误记录：

• 每次都是实时地计算gcd导致TLE，这样可能会重复计算很多；提前计算结果使用dp数组储存就过了。