SCUACM22暑假集训前劝退赛部分题解

最新推荐文章于 2024-07-22 09:38:37 发布
最新推荐文章于 2024-07-22 09:38:37 发布
阅读量126 收藏
点赞数 1
分类专栏： 基础dp&&矩阵快速幂 数论&组合数学 文章标签： 动态规划 c++ 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hans774882968/article/details/125139227
版权

文章目录


传送门

作者：hans774882968以及hans774882968

本文juejin传送门

A

简单几何题，用三角函数知识很容易推（对于考上SCU的同学）。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 1e6 + 5;

int r, a, b, h;

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

int main() {
    read (r, a, b, h);
    if (2 * r <= b) {
        puts ("Drop");
        return 0;
    }
    double H = 1.0 * h * a / (a - b);
    double theta = atan (1.0 * a / 2 / H);
    printf ("Stuck\n%.10lf\n", r / sin (theta) - b / (2 * tan (theta) ) );
    return 0;
}

B

一看就是博弈dp。先打个表

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 5e3 + 5;

int dp[N][N], a[N];
int n, m;

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

bool dfs (int x, int y) {
    if (~dp[x][y]) return dp[x][y];
    dp[x][y] = false;
    if (! (x + y) ) return dp[x][y] = false;
    if (!x || !y) return dp[x][y] = true;
    if (x == y) return dp[x][y] = true;
    if (x < y) return dp[x][y] = dfs (y, x);
    rep (k1, 1, x) {
        for (int k2 = 0; k2 <= y; k2 += k1) {
            dp[x][y] |= !dfs (x - k1, y - k2);
            if (dp[x][y]) break;
        }
    }
    rep (k1, 1, y) {
        for (int k2 = 0; k2 <= x; k2 += k1) {
            dp[x][y] |= !dfs (x - k2, y - k1);
            if (dp[x][y]) break;
        }
    }
    return dp[x][y];
}

void bf (int n) {
    memset (dp, -1, sizeof dp);
    rep (i, 0, n) rep (j, 0, n) dfs (i, j);
    // rep (i, 0, n) rep (j, 0, n) cout << dp[i][j] << " \n"[j == n];
    rep (i, 0, n) rep (j, 0, n) if (!dp[i][j]) dbg (i, j);
    puts ("");
    rep (i, 0, n) rep (j, 0, n) if (!dp[i][j] && i < j) dbg (i, j);
}

int main() {
    bf (300);
    // memset (dp, -1, sizeof dp);
    // int lim = 500;
    // rep (i, 1, lim) {
        // if (a[i]) continue;
        // rep (j, i + 1, lim) {
            // if (!dfs (i, j) ) {
                // a[i] = j;
                // break;
            // }
        // }
        // if (a[i]) a[a[i]] = i;
        // if (i % 100 == 0) cout << "i=" << i << endl; //dbg
    // }
    // rep (i, 1, lim) cout << a[i] << " \n"[i == lim];
    return 0;
}

发现每一行至多一个false，但没有发现其他规律。因此有一种做法是跑上面那段被注释的代码，把每一行的false的点（可能不存在）给找出来。但是要等很久，算了。

正解是：既然为false的点不超过5000，那么就直接用刷表法。根据博弈dp常识，值为true的点推不出任何点，值为false的点可以推出一系列值为true的点。根据调和级数，复杂度O(n^2*logn)。说实话这是我第一次看到这种套路。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 5e3 + 5;

bool dp[N][N];
int n, m;

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

void init (int n) {
    rep (i, 0, n) rep (j, 0, n) {
        if (dp[i][j]) continue;
        rep (k1, 1, n - i) {
            for (int k2 = 0; j + k2 <= n; k2 += k1) {
                dp[i + k1][j + k2] = true;
            }
        }
        rep (k1, 1, n - j) {
            for (int k2 = 0; i + k2 <= n; k2 += k1) {
                dp[i + k2][j + k1] = true;
            }
        }
    }
}

int main() {
    init (5000);
    int T;
    read (T);
    while (T--) {
        read (n, m);
        puts (dp[n][m] ? "Alice" : "Bob");
    }
    return 0;
}

E

设dp[i]为当前得到的数集的最大值为i，能走的步数的期望。再设ans[i]表示当前得到的数集的最大值为i，能获得的分数的期望。因为题目中，分数只获得1次，所以我们假设当前存在一个数集的最大值，但并没有拿过分数。显然这一构造性假设是合法的。

没有数，就用“有一个0”来表示，于是ans[0]即所求。

值得注意的是，我们也可以用数学形式E(x[i])和E(x[i] ^ 2)来表示dp[i]和ans[i]。

推dp

全期望公式：E(X) = sum(E(X|Y = y[i]) * P(Y = y[i]))

设下一个数为j。当下一个数小于i，能走的步数为1；大于等于i，能走的步数为E(1 + x[j]) = 1 + dp[j]。故由全期望公式：

dp[i] = sum(p[j])(j < i) + sum(p[j] * (1 + dp[j]))(j >= i)

移项得：

dp[i] = (1 + sum(p[j] * dp[j])) / (1 - p[i]) (j > i)

推ans

当下一个数小于i，能获得的分数为1。大于等于i，由上述数学形式得E((1 + x[j]) ^ 2) = E(x[j] ^ 2) + 2 * E(x[j]) + 1 = ans[j] + 2 * dp[j] + 1。故由全期望公式：

ans[i] = sum(p[j])(j < i) + sum(p[j] * (ans[j] + 2 * dp[j] + 1))(j >= i)

移项得：

ans[i] = (1 + 2 * dp[i] * p[i] + sum(p[j] * (ans[j] + 2 * dp[j]))) / (1 - p[i]) (j > i)

总结：关注期望dp的数学记号！

代码
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 5e3 + 5;
const int mod = 998244353;

LL dp[N], ans[N], p[N];
int a[N];
int n;

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

LL q_pow (LL a, LL b) {
    LL r = 1;
    for (; b; b >>= 1) {
        if (b & 1) r = r * a % mod;
        a = a * a % mod;
    }
    return r;
}

int main() {
    read (n);
    rep (i, 1, n) read (a[i]);
    LL isuma = q_pow (accumulate (a + 1, a + n + 1, 0LL), mod - 2);
    rep (i, 1, n) p[i] = a[i] * isuma % mod;
    dwn (i, n, 0) {
        LL s1 = 1, inv1_p = q_pow ( (1 - p[i] + mod) % mod, mod - 2);
        rep (j, i + 1, n)
            s1 = (s1 + p[j] * dp[j] % mod) % mod;
        dp[i] = s1 * inv1_p % mod;
        LL s2 = (1 + 2 * p[i] * dp[i] % mod) % mod;
        rep (j, i + 1, n)
            s2 = (s2 + p[j] * (ans[j] + 2 * dp[j] % mod) % mod) % mod;
        ans[i] = s2 * inv1_p % mod;
    }
    printf ("%lld\n", ans[0]);
    return 0;
}

H

a[i]*a[j] = a[k]，以及数据范围1e6都提示我们需要枚举倍数。开个桶计数，然后枚举值i和i的倍数j，则贡献为c[i] * c[j] * c[j/i]。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 1e6 + 5;

int n, c[N];

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

int main() {
    read (n);
    int mx = 0;
    rep (_, 1, n) {
        int x;
        read (x);
        c[x]++;
        mx = max (mx, x);
    }
    LL ans = 0;
    rep (i, 1, mx) {
        if (!c[i]) continue;
        for (int j = i; j <= mx; j += i) {
            ans += 1LL * c[i] * c[j] * c[j / i];
        }
    }
    printf ("%lld\n", ans);
    return 0;
}

I

签到。

  • x是s的子集，x为1的部分随便取，x为0但s为1的部分固定为1。
  • x不是s的子集，个数显然为0。

易错点：x == s时0也是答案，但题目问的是正整数，要减掉。

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
#define rep(i,a,b) for(int i = (a);i <= (b);++i)
#define re_(i,a,b) for(int i = (a);i < (b);++i)
#define dwn(i,a,b) for(int i = (a);i >= (b);--i)

const int N = 2e5 + 5;

int x, s;

void dbg() {
    puts ("");
}
template<typename T, typename... R>void dbg (const T &f, const R &... r) {
    cout << f << " ";
    dbg (r...);
}
template<typename Type>inline void read (Type &xx) {
    Type f = 1;
    char ch;
    xx = 0;
    for (ch = getchar(); ch < '0' || ch > '9'; ch = getchar() ) if (ch == '-') f = -1;
    for (; ch >= '0' && ch <= '9'; ch = getchar() ) xx = xx * 10 + ch - '0';
    xx *= f;
}
void read() {}
template<typename T, typename ...R>void read (T &x, R &...r) {
    read (x);
    read (r...);
}

int main() {
    read (x, s);
    if (! ( (x & s) == x && (x | s) == s) ) {
        puts ("0");
        return 0;
    }
    printf ("%lld\n", (1LL << __builtin_popcount (x) ) - (x == s) );
    return 0;
}
hans774882968
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
cug2015ACM集训数学进阶部分题解
07-17
ACM数学进阶题解 ，思路很好。保存起来，方便以后使用
2020牛客多校暑期集训营第一场题解.pdf
01-20
2020牛客多校暑期集训营第一场题解.pdf
暑假集训前部分题目题解（持续更新）
The_OIer的博客
08-05 296
8.5 毒瘤 [HNOI2016] 树 毒瘤 [CQOI2015] 任务查询系统
暑假集训模拟赛部分题题解
ezoixx118的博客
08-14 348
暑假集训=（狂补数学+狂打模拟赛），暑假作业也翘了。。。 由于我太菜，模拟赛几乎场场爆0，就写了写部分好题的题解如下： （有些题我至今没理解，所以没法题解） 模拟赛1~19：还未有写题解的计划（以后可能会补） 模拟赛20：B(Tree Game) 模拟赛21：C(Flags) 模拟赛22~24：暂无。。。 模拟赛25：A(RGB Sequence) C(Prime Flip)...
2021NEUACM暑假集训补题部分题解_20210724_Training
Passenger12138的博客
07-05 350
2021NEUACM暑假集训水赛补题部分题解
IOI2021 国家集训队作业部分题解
zx2003的博客
01-16 1446
读者本地应该都有由各个集训队员撰写的题解。本文旨在对其进行一些补充。 代码基本都能在cf gym的提交中找到，用户名zx2003（应该没人不知道吧。。。）WF2017-2019需要在loj上找，因为当时觉得loj交起来比较舒适。 WF2014 B 本题复杂度是可以优于yht给出 O(w2log⁡d)O(w^2\log d)O(w2logd) 的，虽然也没什么意思。 对于每种离散物品，将背包容量按模 wiw_iwi​ 的值分类后，每一类都是可以决策单调性优化转移的，写个分治就没了，对zx还写了这个，每一类都是
牛客国庆集训派对Day1 —— 题解
Scar_Halo
10-01 804
牛客国庆集训派对Day2 —— 题解
Scar_Halo
10-02 573
出题人期望的难度： A 矩阵乘法 B 字符串的幂 C 生命游戏 D 数格点 E 数据排序 F 平衡二叉树 G 数组合并 H 卡牌游戏 I 游戏 J 魔法阵 K 排队 ...
牛客国庆集训派对Day3 —— 题解
Scar_Halo
10-03 559
牛客国庆集训派对Day4 —— 题解
Scar_Halo
10-04 481
2022广西师范大学暑期训练赛部分题 题解
qq_64554898的博客
10-07 668
【代码】2022广西师范大学暑期训练赛部分题 题解。
CSDN::START 暑期集训总结，题解补完计划开启
fuathear的专栏
08-19 130
2018的暑假集训只剩下最后一周了，经过一个月的训练，除去hdu的多校自闭赛，开了至少150道题了。 自己粗略估计了下，150道题里面现场能做出来的也就只有50道左右，远远低于预期。虽然每次训练赛后自己也在补题，但还是仅限于提交代码并且AC。对于一些有价值的题目，并没有整理下思路、写出题解，以后遇到近似题目，仍不确定能不能写出来。 所以从今天开始，我将陆续把部分题解及思路整理下来，发到csdn...
2020牛客暑期多校集训营第二场题解.pdf
01-20
2020牛客暑期多校集训营第二场题解.pdf
C语言 之 理解指针（1）
fallzzzzz的博客
07-20 1131
就是先获取该数组的首元素的地址，从前面我们知道，指针的类型的作用是指针+1 所跳过的字节数，int为整形，大小为4个字节，所以p+1就跳过4个字节，p+2就跳过8个字节，计算机的CPU（中央处理器）在处理数据的时候，需要的数据是在内存中读取的，处理后数据也会放回内存中，如果想对这些空间进行高效的管理，就可以。并且我们知道，数组的元素存放是随着下标增长，内存地址变大的。上面的代码创建了一个整型变量n，因为一个整形的大小是4个字节，所以内存中申请了4个字节，⽤于存放整数10，其中每个字节都有地址。
C++二叉搜索树
Jk_Mr的博客
07-20 1127
树结构是一种重要的数据结构。树是一种非线性的数据结构，它是由n（n>=0）个有限结点组成一个具有层次关系的集合。把它叫做树是因为它看起来像一棵倒挂的树，也就是说它是根朝上，而叶朝下的。
C++笔试强训7
m0_74189279的博客
07-21 972
这意味着编译器会尝试将函数的代码直接插入到每个调用该函数的地方，而不是像通常那样进行函数调用（即，生成调用指令，跳转到函数代码，执行函数体，然后返回）。友元函数是在类外部定义的普通函数，它只能通过传入的参数（如果有的话）来访问类的成员，或者通过类的对象显式地访问（如果该函数设计为接受类的对象作为参数）。函数参数默认值又叫缺省参数，缺省参数必须从右向左开始缺省，必须连续给定缺省值，也就是说，从左到右，一旦遇到一个参数是缺省的，那么之后的参数也必须都是缺省的。指针是类成员函数特有的，用于指向调用该函数的对象。
商品信息管理系统(C语言)
最新发布
zeyeqianli的博客
07-22 324
该案例使用了C语言中最具特色的结构体，将每个商品的所有信息存在结构体中，并且定义一个结构体类型的数组保存所有商品的信息，并且按照模块化的编程思想，将要实现的每个功能编写成独立的函数，这样即方便阅读同时也方便差错修改。该系统的主要功能如下：(1) 登录系统。(2) 创建商品信息。(3) 打印商品信息。(4) 查询商品信息。(5) 修改商品信息。(6) 删除商品信息。(7) 商品价格排名显示。(8) 退出系统。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值