程序设计思维与实践 Week16 CSP模测

A TT数鸭子

1. 题目大意

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

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

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

样例：

65
123456789 9876543210 233 666 1 114514

4

2. 思路历程

• 要求出多少个数的数位中不同的数字个数小于k，先要对各个数字的不同位数进行计数。
• 数字只有0-9，因此用一个长度为10的数组 a[i] 表示某个数中数字 i 有a[i] 个，计数后遍历数组a即可。
• 为了方便从数中得到单个的数字，把数字当作字符串读入，scanf遇到空格停止读入因此输入无需另外处理。

3. 反思

• 真不知道自己怎么想的，竟然用相邻两个数字是否相同来对不同的位数进行计数（唉有时候错的真的让人没话说
• 还是读题要仔细，思考不能太慌张吧，要谨记欧

4. 代码

#include <iostream>
#include <string.h>
using namespace std;
const int maxn = 1000011;

int n, k;
int cnt, ans = 0, a[11];
char duck[maxn][20];

int main()
{
    cin >> n >> k;
    for (int i = 0; i < n; i++) scanf("%s", duck[i]);
    
    for (int i = 0; i < n; i++)
    {
        memset(a, 0, sizeof(a));
        for (int j = 0; j < strlen(duck[i]); j++)
        {
            a[duck[i][j] - '0']++;
        }
        cnt = 0;
        for (int i = 0; i < 10; i++)
        {
            if (a[i]) cnt++;
        }
        if (cnt < k) ans++;
    }
    
    cout << ans;
    return 0;
}

B ZJM抵御宇宙射线

1. 题目大意

据传，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

2. 思路历程

• 刚开始的时候考虑根据平均值找到最中心的点，再用bfs的方法找到最远的点，最后求出中心点到其的距离，想的可太麻烦了。
• 这其实是个非常暴力的题，用枚举法解决，尝试以任何一个点为中心求得可以覆盖所有点的半径，再在所有半径中找到最短的一条即可。
• 由于多解时的优先顺序为先输出 x 较小的点再输出 y 较小的点，因此注意需要在点的结构体里重构比较函数。

3. 反思

• 我刚开始想到的先找中心再bfs方法用了各种各样的数据结构（集合、队列等等），越来越复杂，还是要学会根据题目的位置预判难度，然后尽早换掉繁琐的思路
• 发现自己思考问题还是按照原来的方式，比如这道题要求中心和半径就非常“规矩”地在找中心和半径，计算机思维里比较简单的枚举法也没有想到。（可能可以归结为没有见过很多题？
• 不管怎么说，还是勤加练习吧

4. 代码

#include <iostream>
#include <math.h>
#include <algorithm>
using namespace std;
const int maxn = 1011;

struct Pos
{
    int x, y;
    
} pos[maxn];

bool cmp(Pos a, Pos b)
{
    if (a.x != b.x) return a.x < b.x;
    return a.y < b.y;
}

long long get_dis(Pos a, Pos b)
{
    long long x = abs(a.x - b.x);
    long long y = abs(a.y - b.y);
    return x * x + y * y;
}

int main()
{
    int n;
    cin >> n;
    for (int i = 0; i < n; i++) cin >> pos[i].x >> pos[i].y;
    
    sort(pos, pos + n, cmp);
    Pos temp;
    long long min_dis = 1e15;
    for (int i = 0; i < n; i++)
    {
        long long max_dis = -1;
        for (int j = 0; j < n; j++)
        {
            if (i == j) continue;
            long long d = get_dis(pos[i], pos[j]);
            if (d > max_dis) max_dis = d;
        }
        
        if (max_dis < min_dis)
        {
            min_dis = max_dis;
            temp.x = pos[i].x;
            temp.y = pos[i].y;
        }
    }
    
    printf("%.2f %.2f\n%.2f", (double)temp.x, (double)temp.y, (double)min_dis);
    return 0;
}

C 宇宙狗的危机

1. 题目大意

在瑞神大战宇宙射线中我们了解到了宇宙狗的厉害之处，虽然宇宙狗凶神恶煞，但是宇宙狗有一 个很可爱的女朋友。
最近，他的女朋友得到了一些数，同时，她还很喜欢树，所以她打算把得到的数拼成一颗树。
这一天，她快拼完了，同时她和好友相约假期出去玩。贪吃的宇宙狗不小心把树的树枝都吃掉 了。所以恐惧包围了宇宙狗，他现在要恢复整棵树，但是它只知道这棵树是一颗二叉搜索树，同 时任意树边相连的两个节点的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个数$a_i$，输入保证是升序的。

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

样例：

1
6
3 6 9 18 36 108

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

Yes

No
Yes

2. 思路历程

• 题目要求二叉树，也给出了区间，是一个区间dp问题，也要按照dp的套路慢慢分析
• f[i][j]表示区间[i,j]合法
  L[i][k]表示[i, k - 1]能构成k的左子树
  R[k][j]表示[k + 1, j]能构成k的右子树
• 如果k可以与l - 1相连，说明以l - 1为根，[l, r]为右儿子的状态合法
  如果k可以与r + 1相连，说明以r + 1为根，[l, r]为左儿子的状态合法
• 当L[l][k] = 1且R[k][r] = 1，即[l, k]和[k, r]都是一个合法的二叉搜索树的时候，时，两子树合并后仍是二叉树，因此状态合法，可以转移
  最后情况考虑到l == 1且r == n时，说明找到了一棵符合要求的二叉树

3. 反思

• 做模测题的速度还是不行，这一题甚至都没能写完时间就到了（唉
• 我知道自己的dp学的不怎么扎实，所以拿到题目的时候第一反应是选择用之前朴素的树的方法做，构造了二叉搜索树也写了层次遍历的函数。
• 可惜调了很久也不可行，dp真的很重要欧，dp的转移条件，特殊情况判断都需要好好掌握，复习的时候一定要注意把所有的dp整理归类！

4. 代码

#include <iostream>
#include <algorithm>
#include <string.h>
using namespace std;
const int maxn = 707;

int t, n, a[maxn], L[maxn][maxn], R[maxn][maxn], f[maxn][maxn];
bool flag;

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

void solve()
{
    cin >> n;
    flag = false;
    memset(f, 0, sizeof(f));
    memset(L, 0, sizeof(L));
    memset(R, 0, sizeof(R));
    
    for (int i = 1; i <= n; i++)
    {
        cin >> a[i];
        L[i][i] = R[i][i] = 1;
    }
    
    for (int i = 1; i <= n; i++)
    {
        for (int j = i; j <= n; j++)
        {
            if (gcd(a[i], a[j]) > 1)
                f[i][j] = f[j][i] = 1;
        }
    }
    
    for (int l = n; l >= 1; l--)
    {
        for (int r = l; r <= n; r++)
        {
            for (int k = l; k <= r; k++)
            {
                if (L[l][k] && R[k][r])
                {
                    if (l == 1 && r == n)
                    {
                        
                        cout << "Yes" << endl;
                        return;
                    }
                    if (f[l - 1][k]) R[l - 1][r] = 1;
                    if (f[k][r + 1]) L[l][r + 1] = 1;
                }
            }
        }
    }
    cout << "No" << endl;
}

int main()
{
    cin>> t;
    while (t--)
    {
        solve();
    }
    
    return 0;
}