鱼C论坛

 找回密码
 立即注册
查看: 1723|回复: 13

[技术交流] FCOI 13 官方题解

[复制链接]
发表于 2024-8-27 22:00:00 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
本帖最后由 zhangjinxuan 于 2024-8-30 12:07 编辑



更好的阅读体验:https://www.luogu.com.cn/paste/vmp3ru1q

FCOI 13 题解

问题 A - 平均

平均数即总和与个数的比,总和显然为 $\frac{N(N+1)}{2}$,比上 $N$ 即 $\frac{N+1}{2}$,乘十后输出答案就是 $5N+5$。


  1. print(int(input())*5+5)
复制代码


问题 B - 区间公因数

相邻数字的 $\gcd$ 一定为 $1$,因此当 $l\neq r$ 时,一定包含至少一个相邻数字,就是 $1$。当 $l=r$ 时,就直接时为 $l$,记得取模。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. signed main() {
  5.         int l, r;
  6.         cin >> l >> r;
  7.         if (l != r) cout << 1;
  8.         else cout << r% 998244353;
  9.         return 0;
  10. }
复制代码

问题 C - 二号线

可以结合 `if-else` 和逻辑或实现。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. signed main() {
  5.         string s;
  6.         cin >> s;
  7.         if (s == "1" || s == "4" || s == "5" || s == "6" || s == "9" || s == "10" || s == "18" || s == "LP" || s == "JT") puts("DT");
  8.         else if (s == "2" || s == "3") puts("DG");
  9.         else if (s == "YB") puts("QG");
  10.         else puts("ER");
  11.         return 0;
  12. }
复制代码


彩蛋解析:因为标题为二号线,所以彩蛋和二号线相关,因此考虑故意将二号线认作为其他类别的类型,尝试 $\text {DT, ER,QG}$ 后发现是 $\text{QG}$。

问题 D - 无限水

注意到自动灌水的格子可以无限取水的,所以若能实现无限灌水就代表答案为 $2$,否则就是空格子数量总和,也就是一个一个格子倒水。

而可以实现无限灌水的唯一条件就是存在至少两个空地相邻,即题目中给定的条件。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. int n;
  5. char s[2003];

  6. signed main() {
  7.         scanf("%d%s", &n, s + 1);
  8.         int cnt = 0;
  9.         for (int i = 1; i <= n; ++i) {
  10.                 if (s[i] == '.') {
  11.                         if (s[i - 1] == '.' && s[i + 1] == '.') return puts("2"), 0;
  12.                         ++cnt;
  13.                 }
  14.         }
  15.         cout << cnt;
  16.         return 0;
  17. }
复制代码


问题 E - AC 时间

找到前驱很简单的方法就是排序,然而这样会影响到整个序列的顺序。为了使得答案以正常的顺序输出,我们可以对于每一个数字再记录一个信息表示他在序列中原来的位置,排序完后直接根据位置信息复原答案即可。

对于 C/C++ 可以使用结构体实现,对于 C++ 可以直接使用 `std::pair`。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. int n, ans[200005];
  5. pair<int, int> a[200005];

  6. signed main() {
  7.         cin >> n;
  8.         for (int i = 1; i <= n; ++i) {
  9.                 cin >> a[i].first;
  10.                 a[i].second = i;
  11.         }
  12.         if (n <= 2000) {
  13.                 for (int i = 1; i <= n; ++i) {
  14.                         int ans = a[i].first;
  15.                         for (int j = 1; j <= n; ++j) {
  16.                                 if (a[i].first > a[j].first) ans = min(ans, a[i].first - a[j].first);
  17.                         }
  18.                         cout << ans << " \n"[i == n];
  19.                 }
  20.         } else {
  21.                 sort(a + 1, a + n + 1);
  22.                 for (int i = 1; i <= n; ++i) {
  23.                         ans[a[i].second] = a[i].first - a[i - 1].first;
  24.                 }
  25.                 for (int i = 1; i <= n; ++i) cout << ans[i] << " \n"[i == n];
  26.         }
  27.         return 0;
  28. }
复制代码



问题 F - 真算

Sub1
同时是 $4,6$ 的倍数一定要满足是 $12$ 的倍数,题目的代码却判断的是 $24$ 的倍数,数据 `12` 即可 Hack。

Sub2

注意到 $ab$ 的最大值比 $9.5\times10^{8}$ 大,因此 $ab \times cd$ 可以超过 `long long` 的存储范围,考虑令 $a=b=\lfloor \sqrt{998244353} \rfloor$,这样 $a,b$ 就能在值域较小的时候取到较大的模意义下的乘积。数据 `31595 31595 100000 100000` 即可。
Sub3
如果是随机生成每一次是 $O(\log n)$ 的,考虑这样构造,每一次使用一个大盒子取拼接小盒子,这样每一次遍历一次大盒子,单次会达到 $O(n)$,构造仅需要在第 $i$ 次操作用 $i$ 盒子去拼接 $i+1$ 即可。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. int n, ans[200005];
  5. pair<int, int> a[200005];

  6. signed main() {
  7.         ios::sync_with_stdio(0);
  8.         cin >> n;
  9.         if (n == 1) cout << 12;
  10.         else if (n == 2) cout << "31595 31595 100000 100000";
  11.         else if (n == 3) {
  12.                 cout << 50000 << " " << 49999 << "\n";
  13.                 for (int i = 1; i <= 50000; ++i) {
  14.                         cout << i << " \n"[i == 50000];
  15.                 }
  16.                 for (int i = 1; i <= 49999; ++i) {
  17.                         cout << i << " " << i + 1 << "\n";
  18.                 }
  19.         }
  20.         return 0;
  21. }
复制代码


问题 G - 连除

尝试在 $A_l$ 后面添加小括号,那么答案就是 $\displaystyle \frac{ A_l }{ \prod_{i=l+1}^{r} A_i}$,就转换成了两数模意义下除法和区间连乘的问题。

对于一个序列,我们可以预处理他们的前缀乘积,即对于任意 $S_i$ 都等于 $\displaystyle \prod_{j=1} ^{i} A_i$,这里是需要取模的。

然后查询 $l+1\sim r$ 的连乘就是 $\frac{S_r}{S_{l}}$,这一步同样也是一次模意义下除法。所以一段区间的连除就是 $\frac{A_l}{\frac{S_r}{S_{l}}}$。

模意义下的除法题目已经给出定义,即标准的乘法逆元,求解这个可以使用快速幂实现,即费马小定理求解逆元。


  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define int long long

  4. int n, q, a[200005], pres[200005];
  5. const int p = 998244353;

  6. long long qpow(long long a, long long b) {
  7.         long long res = 1;
  8.         while (b) {
  9.                 if (b & 1) res = res * a % p;
  10.                 a = a * a % p;
  11.                 b >>= 1;
  12.         }
  13.         return res;
  14. }

  15. signed main() {
  16.         ios::sync_with_stdio(0);
  17.         cin >> n >> q;
  18.         pres[0] = 1;
  19.         for (int i = 1; i <= n; ++i) {
  20.                 cin >> a[i];
  21.                 pres[i] = (pres[i - 1] * a[i]) % p;
  22.         }
  23.         while (q--) {
  24.                 int l, r;
  25.                 cin >> l >> r;
  26.                 cout << a[l] * qpow(pres[r] * qpow(pres[l], p - 2) % p, p - 2) % p << "\n";
  27.         }
  28.         return 0;
  29. }
复制代码


问题 H - 数字 LIS

考虑如何使用最少的量来维护一个最长不降子序列,显然需要三种东西给,一个是上一个元素,一个是以本身结尾的答案,一个是全局的答案,因此就能用 $O(N^3)$ 的复杂度表示出整个问题,而不是枚举数字的 $O(10^N)$,这样,只需要把这 $O(N^3)$ 的状态作为 DP 的维度进行 DP 即可。

但是这个算法只有在 $n\le100$ 的时候才能勉强跑一跑,对于 $n\le1000$,几乎不可能优化,所以考虑使用打表的技巧,由于空间过大,所以考虑使用滚动数组优化,即可在本地运行。

打表程序:

  1. #include <bits/stdc++.h>
  2. using namespace std;
  3. #define lim 1000
  4. long long n, f[2][11][lim+2][lim+2];
  5. int main() {
  6.         freopen("ans.cpp", "w", stdout);
  7.         cin >> n;
  8.         for (int i = 0; i < 11; ++i) for (int j = 0; j < lim; ++j) for (int k = 0; k < lim; ++k) f[0][i][j][k] = k;
  9.         cout << "#include <bits/stdc++.h>\nusing namespace std;\n";
  10.         cout << "const int ans[1004] = {-1, ";
  11.         for (int i = 1; i <= n; ++i) {
  12.                 for (int lst = 0; lst <= 10; ++lst)
  13.                 for (int now = 0; now <= n + 1; ++now) for (int ans = 0; ans <= n + 1; ++ans) {
  14.                         int sum = 0;
  15.                         for (int j = 0; j <= 9; ++j)
  16.                                 if (j >= lst) (sum += f[(i - 1) & 1][j][now + 1][max(ans, now + 1)]) %= 998244353;
  17.                                 else (sum += f[(i - 1) & 1][j][1][max(ans, 1)]) %= 998244353;               
  18.                         f[i & 1][lst][now][ans] = sum;
  19.                 }
  20.                 cout << f[i & 1][10][1][1] << ", " << flush;
  21.                 cerr << i << endl;
  22.         }
  23.         cout << "};\n";
  24.         cout << "\nint main() {\n\tint t;\n\tcin >> t;\n\twhile(t--) {int n; cin >> n; cout << ans[n] << "\\n"; }\n\treturn 0;\n}";
  25.         return 0;
  26. }
复制代码


至少需要 $5\text{min}$,但应该是已知的最优解,如果您发现了时间更低的做法,欢迎告诉我。

本帖被以下淘专辑推荐:

小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

发表于 2024-8-27 23:46:29 From FishC Mobile | 显示全部楼层
求py代码
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 08:33:36 | 显示全部楼层
打表???!!!
小甲鱼最新课程 -> https://ilovefishc.com
回复

使用道具 举报

发表于 2024-8-28 08:50:40 | 显示全部楼层
这非常规过分了吧。。。
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 09:22:01 | 显示全部楼层
我有一种最后一道题的想法。。。
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 09:25:24 | 显示全部楼层
如果序列足够大,将输入序列分段,将每一段分配给一个线程,每个线程计算出该段所有可能的最长不降子序列,后将所有返回,由主线程再次将其分配,将每个不降序列进行排序,即可
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2024-8-28 18:15:25 | 显示全部楼层
本帖最后由 zhangjinxuan 于 2024-8-28 18:31 编辑
zhangchenyvn 发表于 2024-8-28 09:25
如果序列足够大,将输入序列分段,将每一段分配给一个线程,每个线程计算出该段所有可能的最长不降子序列, ...


无意义,因为多线程相当于使用超级计算机。

而且总序列量级是 $10^{1000}$ 的。
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 18:41:24 | 显示全部楼层
懂了,非常规=逆天
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2024-8-28 18:42:16 | 显示全部楼层

G:


  1. n, q = map(int, input().split())
  2. p = 998244353
  3. a = [0] + list(map(int, input().split()))
  4. pres = [1]
  5. for i in range(1, n + 1):
  6.     pres.append(pres[-1] * a[i] % p)
  7. for i in range(q):
  8.     l, r = map(int, input().split())
  9.     print(a[l] * pow(pres[r] * pow(pres[l], p - 2, p) % p, p - 2, p) % p)
复制代码


记得用 PyPy3。

H 我觉得转换也不复杂(
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2024-8-28 18:44:57 | 显示全部楼层
sfqxx 发表于 2024-8-28 18:41
懂了,非常规=逆天

来我们重庆育才,虽然不保证一定金牌,但是绝对保证提升你的假算和打表能力(
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 19:10:14 | 显示全部楼层
zhangjinxuan 发表于 2024-8-28 18:15
无意义,因为多线程相当于使用超级计算机。

而且总序列量级是 $10^{1000}$ 的。

有道理
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 19:13:18 | 显示全部楼层
本帖最后由 zhangchenyvn 于 2024-8-28 19:15 编辑


H题py文件我也转了一下
  1. # 定义常量
  2. lim = 1000
  3. n = int(input())
  4. f = [[[[0] * (lim + 2) for _ in range(lim + 2)] for _ in range(11)] for _ in range(2)]

  5. for i in range(11):
  6.     for j in range(lim):
  7.         for k in range(lim):
  8.             f[0][i][j][k] = k

  9. with open("ans.py", "w") as output_file:
  10.     output_file.write("# Converted Python Code\n")
  11.     output_file.write("ans = [-1, ")

  12.     for i in range(1, n + 1):
  13.         for lst in range(11):
  14.             for now in range(n + 2):
  15.                 for ans in range(n + 2):
  16.                     sum_value = 0
  17.                     for j in range(10):
  18.                         if j >= lst:
  19.                             sum_value += f[(i - 1) & 1][j][now + 1][max(ans, now + 1)]
  20.                         else:
  21.                             sum_value += f[(i - 1) & 1][j][1][max(ans, 1)]
  22.                         sum_value %= 998244353
  23.                     f[i & 1][lst][now][ans] = sum_value
  24.         
  25.   
  26.         output_file.write(f"{f[i & 1][10][1][1]}, ")
  27.         print(i)  

  28.     output_file.write("]\n\n")
  29.     output_file.write("def main():\n")
  30.     output_file.write("\tt = int(input())\n")
  31.     output_file.write("\twhile t > 0:\n")
  32.     output_file.write("\t\tn = int(input())\n")
  33.     output_file.write("\t\tprint(ans[n])\n")
  34.     output_file.write("\t\tt -= 1\n")
  35.     output_file.write("\nif __name__ == '__main__':\n\tmain()\n")
复制代码
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-8-28 19:30:35 | 显示全部楼层
小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2024-9-6 13:55:07 | 显示全部楼层
zhangchenyvn 发表于 2024-8-28 19:13
H题py文件我也转了一下

小甲鱼最新课程 -> https://ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-5-10 07:58

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表