『模拟赛』CSP-S模拟7(更新 T4

Rank

image

A. median

签。

你说得对,但是赛时嗯打 150 行 5k 代码超级分讨过了。

因为容斥做的不好,求总的然后减总会差点东西,所以直接分着加。发现如果中位数在这五个数中不止出现一次那么就会算重,所以分三种大情况考虑。

一,中位数只有一个。那么此时我们需要找到另外两个小于它的和另外两个大于它的,符合乘法原理。

二,中位数有两个。那么此时要分为一个数比它们小和两个数比它们小的情况。

三,中位数三个及以上。此时无论其它数选什么都不会改变贡献。

考虑实现以上操作需要维护什么。首先值域范围 \(10^9\),硬维护数量铁开不下,但 \(n\) 只有 \(10^5\),也就是说最多有 \(5\times 10^5\) 种不同的数,读入之后离散化一下即可,然后处理出每个数组的前缀和,再根据上述情况做就行了。外层枚举中位数在离散化后的编号,内层找含有该数的数组,然后进 dfs 找到方案即可。复杂度看起来是 \(\mathcal{O(5n\times 2^5)}\) 的,但实际上跑不满,效率还是很高的。

赛时由于是一点一点打的,可能常数有点大,很多赘余可以合并起来,但我不想重构了,所以就放这个删了注释的赛时代码了,理解之后还是很好理解的。

点击查看代码
#include<bits/stdc++.h>
#define fo(x,y,z) for(register int (x) = (y);(x)<=(z);(x)++)
#define fu(x,y,z) for(register int (x) = (y);(x)>=(z);(x)--)
using namespace std;
typedef long long ll;
#define lx ll
inline lx qr()
{char ch = getchar();lx x = 0 , f = 1;for(;ch<'0'||ch>'9';ch = getchar()) if(ch == '-') f = -1;for(;ch>= '0' && ch<= '9';ch = getchar()) x = (x<<3) + (x<<1) + (ch^48);return x*f;
}
#undef lx
#define qr qr()
#define pii pair<int , int>
#define fi first
#define se second
const int Ratio = 0;
const int N = 1e5 + 5;
const int mod = 998244353;
int n, m;
int a[6][N], al[N * 5], tot, cnt[6][N * 5];
bool yz1[6], yz2[6], yzz[6], yz3[6], yz4[6];
ll ans;
namespace Wisadel
{void Wdfs1(int wc, int targ, int now, int rsum){if(rsum == 2){ll res = cnt[wc][targ] - cnt[wc][targ - 1];ll lca = 1, rca = 1;fo(i, 1, 5){if(i == wc) continue;if(yz1[i]) lca = lca * cnt[i][targ - 1] % mod;if(yz2[i]) rca = rca * (cnt[i][tot] - cnt[i][targ]) % mod;}res = res * lca % mod * rca % mod;ans = (ans + al[targ] * res) % mod;return;}if(now > 5) return ;if(now == wc || yz1[now]){Wdfs1(wc, targ, now + 1, rsum); return;}if(cnt[now][tot] - cnt[now][targ] > 0) yz2[now] = 1, Wdfs1(wc, targ, now + 1, rsum + 1), yz2[now] = 0;Wdfs1(wc, targ, now + 1, rsum);}void Wdfs(int wc, int targ, int now, int lsum){if(lsum == 2){Wdfs1(wc, targ, 1, 0);return ;}if(now > 5) return ;if(now == wc){Wdfs(wc, targ, now + 1, lsum); return;}if(cnt[now][targ - 1] > 0) yz1[now] = 1, Wdfs(wc, targ, now + 1, lsum + 1), yz1[now] = 0;Wdfs(wc, targ, now + 1, lsum);}void Wdfss(int now, int targ, int ned){if(ned <= 0){ll res = 1, cc = 1;fo(i, 1, 5){if(yzz[i]) res = res * (cnt[i][targ] - cnt[i][targ - 1]) % mod;else cc = cc * (cnt[i][tot] - (cnt[i][targ] - cnt[i][targ - 1])) % mod;}ans = (ans + res * cc % mod * al[targ] % mod) % mod;return;}if(now > 5) return;if(cnt[now][targ] - cnt[now][targ - 1] > 0) yzz[now] = 1, Wdfss(now + 1, targ, ned - 1), yzz[now] = 0;Wdfss(now + 1, targ, ned);}void Wdfsss(int now, int targ, int ned){if(ned <= 0){ll res = 1, lsu = 1, rsu = 1;fo(i, 1, 5){if(yz3[i]) res = res * (cnt[i][targ] - cnt[i][targ - 1]) % mod;else if(yz4[i]) lsu = lsu * cnt[i][targ - 1] % mod;else rsu = rsu * (cnt[i][tot] - cnt[i][targ]) % mod;}ans = (ans + res * lsu % mod * rsu % mod * al[targ] % mod) % mod;return;}if(now > 5) return;if(yz3[now]){Wdfsss(now + 1, targ, ned); return ;}if(cnt[now][targ - 1] > 0) yz4[now] = 1, Wdfsss(now + 1, targ, ned - 1), yz4[now] = 0;Wdfsss(now + 1, targ, ned);}void Wdfs2(int now, int targ, int ned){if(ned <= 0){Wdfsss(1, targ, 1);Wdfsss(1, targ, 2);return;}if(now > 5) return;if(cnt[now][targ] - cnt[now][targ - 1] > 0) yz3[now] = 1, Wdfs2(now + 1, targ, ned - 1), yz3[now] = 0;Wdfs2(now + 1, targ, ned);}short main(){freopen("median.in", "r", stdin) , freopen("median.out", "w", stdout);n = qr;fo(i, 1, 5)fo(j, 1, n) a[i][j] = qr, al[(i - 1) * n + j] = a[i][j];tot = n * 5;sort(al + 1, al + tot + 1);tot = unique(al + 1, al + tot + 1) -al - 1;fo(i, 1, 5)fo(j, 1, n)cnt[i][lower_bound(al + 1, al + tot + 1, a[i][j]) -al] ++;fo(i, 1, 5) fo(j, 1, tot) cnt[i][j] += cnt[i][j - 1];fo(i, 1, tot)fo(j, 1, 5)if(cnt[j][i] - cnt[j][i - 1] > 0)Wdfs(j, i, 1, 0);fo(i, 1, tot){if(i > 1 && i < tot) Wdfs2(1, i, 2);fo(tim, 3, 5)Wdfss(1, i, tim);}printf("%lld\n", ans);return Ratio;}
}
int main(){return Wisadel::main();}

B. travel

也算是签吧,虽然改了数据绑包之后保龄了。

赛时完全没想到自环对答案也是有影响的,就直接删了没管,然后 hack 全输出 Yes 的人的时候顺手给我 hack 了。

Upd:题意解释,任选起点和终点后,若不存在一个边界值,使得之后的 \(f(x,y,k)\) 都为定值,即 \(f(x,y,k)\) 的值是在定义域上不断变化的,那么它就是有趣的。

借助样例解释理解一下,就是当选择一个起点和终点后能一直跑下去并且存在不一样的路径数。事实上就是找环,显然当取环上某一点作为起点和终点时,如果 \(k\) 是环大小的倍数,则答案为 1,否则为 0,无论 \(k\) 取多大都不会成为一个定值。再考虑我漏掉的自环的影响,发现如果一条路径上存在两个及以上的自环,它的每一个自环都可以走任意次,以相连的两个点各有一个自环为例,根据 \(k\) 值的不同,它所得到的答案的数量为 \(k\),显然不存在一个极限值使得其成为定值。

那么就考虑如何实现了。其实是非常简单的,一开始仍然不把自环连进去,重边用 set 处理一下,找环用 dfs 就行,找不到环再去枚举每一条路径求其上自环数量,复杂度 \(\mathcal{O(n+m)}\)

点击查看代码
#include<bits/stdc++.h>
#define fo(x,y,z) for(register int (x) = (y);(x)<=(z);(x)++)
#define fu(x,y,z) for(register int (x) = (y);(x)>=(z);(x)--)
using namespace std;
typedef long long ll;
#define lx ll
inline lx qr()
{char ch = getchar();lx x = 0 , f = 1;for(;ch<'0'||ch>'9';ch = getchar()) if(ch == '-') f = -1;for(;ch>= '0' && ch<= '9';ch = getchar()) x = (x<<3) + (x<<1) + (ch^48);return x*f;
}
#undef lx
#define qr qr()
#define pii pair<int , int>
#define fi first
#define se second
const int Ratio = 0;
const int N = 5e5 + 5;
const int mod = 998244353;
int n, m;
int hh[N], to[N << 1], ne[N << 1], cnt, tot[N], rd[N], zh[N];
bool vis[N], bb[N], huan, yz[N];
set<int> st[N];
namespace Wisadel
{void Wadd(int u, int v){to[++cnt] = v;ne[cnt] = hh[u];hh[u] = cnt;}void Wdfs(int u){vis[u] = 1,yz[u] = 1;for(int i = hh[u]; i != -1; i = ne[i]){int v = to[i];if(yz[v]){huan = 1; return;}Wdfs(v);if(huan) return;}yz[u] = 0;}bool Wbfs(int x){queue<int> q;int sum = 0;q.push(x);while(q.size()){int u = q.front(); q.pop();bb[u] = 1;sum += zh[u];if(sum > 1) return 1;for(int i = hh[u]; i != -1; i = ne[i])q.push(to[i]);}return 0;}short main(){freopen("travel.in", "r", stdin) , freopen("travel.out", "w", stdout);n = qr, m = qr; huan = 0;memset(hh, -1, sizeof hh);fo(i, 1, m){int a = qr, b = qr;if(a == b){zh[a] += 1; continue;}st[a].insert(b);}fo(i, 1, n) for(auto j : st[i]) Wadd(i, j);fo(i, 1, n){if(!vis[i]) Wdfs(i);if(huan) break;}if(huan) printf("Yes\n");else{bool can = 0;fo(i, 1, n)if(!bb[i])if(Wbfs(i)){can = 1; break;}if(can) printf("Yes\n");else printf("No\n");}return Ratio;}
}
int main(){return Wisadel::main();}

C. game

挺简单的博弈论,不过赛时唐了错解没注释掉,然后又顺手给我 hack 了。

首先最基本的,偶数个 1 同时存在的局面一定是必败局面;额外很显然的,两个相同的数的局面也一定是必败局面;那么给出推论:偶数个相同的数存在的局面一定是必败局面。

那么现在只需证明两件事。其一是这些局面无论如何操作结果均不在这个集合内。比较显然的是,后首跟着先手做相同操作,一定能使结果局面仍为必败,因为这样操作的终点是剩余两个 1。此时我们发现可以将必败局面集合扩展到偶数个并且可分为两两相同的若干组的局面。

其二是不在集合内的局面有可以一步操作达到这个局面的操作。我们将个数按升序排序。可以发现此时只要我们选取最大的无法配对为两两相同的一组的一个数,一定存在一组分配方案使得结果局面在集合内。这个证明比较容易,对于当前有奇数的数的情况,达到必败局面所需的个数为 \(\sum_{i=1}^n\ a_{i+1}-a_i\ [i \mod 2 = 0]\),显然其值是小于最大值的;对于偶数的情况,我们让最大的无法配对的数与最小的无法配对的数配对,所需结果仍然是小于拿去后可操作数的,所以结论得证。

然后就非常简单了,对于奇数直接 Yes,偶数就排个序,判一下能否两两配对即可。

点击查看代码
#include<bits/stdc++.h>
#define fo(x,y,z) for(register int (x) = (y);(x)<=(z);(x)++)
#define fu(x,y,z) for(register int (x) = (y);(x)>=(z);(x)--)
using namespace std;
typedef long long ll;
#define lx ll
inline lx qr()
{char ch = getchar();lx x = 0 , f = 1;for(;ch<'0'||ch>'9';ch = getchar()) if(ch == '-') f = -1;for(;ch>= '0' && ch<= '9';ch = getchar()) x = (x<<3) + (x<<1) + (ch^48);return x*f;
}
#undef lx
#define qr qr()
#define pii pair<int , int>
#define fi first
#define se second
const int Ratio = 0;
const int N = 2e5 + 5;
const int mod = 998244353;
int n;
int a[N];
namespace Wisadel
{short main(){freopen("game.in", "r", stdin) , freopen("game.out", "w", stdout);int T = qr;while(T--){n = qr;fo(i, 1, n) a[i] = qr;if(n & 1){printf("Yes\n"); continue;}bool can = 0;sort(a + 1, a + 1 + n);fo(i, 1, n){if(i & 1) continue;if(a[i] != a[i - 1]){can = 1; break;}}if(can) printf("Yes\n");else printf("No\n");}return Ratio;}
}
int main(){return Wisadel::main();}

D. counter

二分的题,看着像 ABC 常出的那种头脑风暴题,改出来再写。

改出来了。

简而言之,首先说明一个性质,一个数 \(x\) 扩展到另一个数的范围是 \([x + 1, x + 9]\),也就是说,每一个区间中的连续九个数一定至少有一个是必经的。

然后我们考虑猫树,其实就是一个不完整的线段树,把每段中间的九个数扒下来,剩下的继续向下递归建树,深度仍然是 \(\log n\) 的。

我们分别处理出每一段上每个点到中间九个点的最短路和中间九个点到段上每个点的最短路,直接 \(\mathcal{O(n)}\) bfs 就行。查询时我们要先递归找到一段使得 \(x\)\(y\) 在中间九个数的两边或者在其中,也就是点 \(x\) 和点 \(y\) 在猫树上的 lca 那一段,然后遍历中间 9 个点作为中途点时 \(x\)\(y\) 的距离,取 min 就做完了。

这道题时间消耗主要在建树上,我们 bfs 的常数要尽可能的小,反向建边(每个点到中间九个点)时最好提前把边连上。还有就是一种最短路的更新可能超出当前段的范围,分别向左右扩展 30 即可。

细节还挺多的,在不影响观感的前提下加了尽量多的注释,代码看不懂可以评论问我。

点击查看代码
#include<bits/stdc++.h>
#define fo(x,y,z) for(register int (x) = (y);(x)<=(z);(x)++)
#define fu(x,y,z) for(register int (x) = (y);(x)>=(z);(x)--)
using namespace std;
typedef long long ll;
#define lx ll
inline lx qr()
{char ch = getchar();lx x = 0 , f = 1;for(;ch<'0'||ch>'9';ch = getchar()) if(ch == '-') f = -1;for(;ch>= '0' && ch<= '9';ch = getchar()) x = (x<<3) + (x<<1) + (ch^48);return x*f;
}
#undef lx
#define qr qr()
#define pii pair<int , int>
#define fi first
#define se second
const int Ratio = 0;
const int N = 1e5 + 35, M = 1e5 + 30;
const int maxn = 1e8;
int n, m, res;
int f[18][10][N], g[18][10][N], dis[N], zc[N], vis[N];
bool yz[N];
vector<int>e[N];
namespace Wisadel
{void Wbfs(int x, int l, int r){// 正着的,x -> vqueue<int> q;int LL = max(l - 30, 0), RR = min(r + 30, M);fill(dis + LL, dis + RR + 1, maxn);fill(yz + LL, yz + RR + 1, 0);// 扩大边界 + 初始化dis[x] = 0, yz[x] = 1;q.push(x);while(q.size()){int u = q.front(); q.pop();int tt = u, zct = 0;while(tt){// 暴力出发点的每一位if(tt % 10) zc[++zct] = tt % 10;tt /= 10;}sort(zc + 1, zc + 1 + zct);zct = unique(zc + 1, zc + 1 + zct) - zc - 1;// 离散化一下(好像是负提升?)fo(i, 1, zct){// bfs 最短路,可以往大可以往小,注意边界if(u + zc[i] <= RR && !yz[u + zc[i]])yz[u + zc[i]] = 1, dis[u + zc[i]] = dis[u] + 1, q.push(u + zc[i]);if(u - zc[i] >= LL && !yz[u - zc[i]])yz[u - zc[i]] = 1, dis[u - zc[i]] = dis[u] + 1, q.push(u - zc[i]);}}}void Wbfs1(int x, int l, int r){// 反,x <- vqueue<int> q;int LL = max(l - 30, 0), RR = min(r + 30, M);fo(i, LL, RR) e[i].clear();fo(i, LL, RR){// 提前连好边int x = i;dis[i] = maxn, yz[i] = 0;while(x){if(x % 10 && vis[x % 10] != i){vis[x % 10] = i;if(i - x % 10 >= LL)e[i - x % 10].push_back(i);if(i + x % 10 <= RR)e[i + x % 10].push_back(i);}x /= 10;}}dis[x] = 0, yz[x] = 1;q.push(x);while(q.size()){int u = q.front(); q.pop();for(auto v : e[u])if(!yz[v]) yz[v] = 1, dis[v] = dis[u] + 1, q.push(v);}}#define mid ((l + r) >> 1)void Wbuild(int dep, int l, int r){fo(i, -4, 4){// 处理相关距离if(mid + i < l || mid + i > r) continue;Wbfs(mid + i, l, r);fo(j, l, r) f[dep][i + 5][j] = dis[j];Wbfs1(mid + i, l, r);fo(j, l, r) g[dep][i + 5][j] = dis[j];}if(r - l + 1 <= 9) return;// 不大于 9 就不用向下再分了Wbuild(dep + 1, l, mid - 5), Wbuild(dep + 1, mid + 5, r);}int Wq(int dep, int l, int r, int x, int y){if(r - l + 1 <= 9) return f[dep][x - mid + 5][y];// x,y 在一块内if(min(x, y) <= mid + 4 && mid - 4 <= max(x, y)){int ans = maxn;fo(i, 1, 9)// 枚举每个中途点ans = min(ans, f[dep][i][y] + g[dep][i][x]);return ans;}if(max(x, y) < mid - 4) return Wq(dep + 1, l, mid - 5, x, y);if(min(x, y) > mid + 4) return Wq(dep + 1, mid + 5, r, x, y);// 判断该往哪边递归}short main(){// freopen(".in", "r", stdin) , freopen(".out", "w", stdout);freopen("counter.in", "r", stdin) , freopen("counter.out", "w", stdout);Wbuild(1, 0, M);int T = qr;while(T--){n = qr ^ (res + 1), m = qr ^ (res + 1);printf("%d\n", (res = Wq(1, 0, M, n, m)) == maxn ? (res = -1) : res);}return Ratio;}
}
int main(){return Wisadel::main();}

模拟赛之决战数据之巅

赛后一看人都傻了,一车 300+,要是 noip 也能这样就好了。

然后就是各种群魔乱舞,T2 输出 Yes 拿 90pts 都算好的,T3 多测输出 Yes 也能拿 90pts 就很逆天了,T4 下午发现没人打到的点也锅了几个,牛牛牛!

收到了高三放六天假的孬消息,然而我们只有不到 30h 的 假,友好学校不够友好,下次记得把我们 OIer 带上。

想玩炉石了。


完结撒花~

image

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/806370.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

微积分快速入门5部分:基本算术、规律及花式算术

12 微积分的基本算术 12.1 加法12.2 乘法12.3 简单除法(倒数)你们原来的份额是 1/x(当 x=2 时,你有 1/2)。 有人进来 你的新份额变成1/(x+1)你的蛋糕数量是如何变化的?在求出总变化(及其恼人的代数)后,我们除以 dx,就得到了 “每 dx ”的变化:现在,我们去掉剩余的 d…

pbootcms常用的13个IF判断语句大全汇总

PBootCMS 提供了丰富的模板标签和条件判断功能,帮助开发者实现各种动态效果。以下是常用的 13 个 IF 判断语句及其具体应用示例。 1. 导航高亮 用途: 用于非首页的导航高亮。 语法:html{pboot:if([nav:scode]=={sort:tcode})}class="active"{/pboot:if}完整示例:…

残基和原子

从您提供的 aa_feature 类的截图信息来看,以下是对 aa_feature 类中各个属性的整理: 主要属性说明aa_embedding:residue_embedding: 一个嵌入层,形状为 (25, 64),用于表示氨基酸残基的嵌入。 res_pos_embedding: 一个嵌入层,形状为 (192, 64),用于表示氨基酸残基的位置嵌…

Windows下安装Nessus 10.8.3安装破解教程

1、下载: 下载地址:https://www.tenable.com/downloads/nessus 浏览器访问 https://127.0.0.1:8834 重点:Register offline,选择“Managed Scanner”, 再选择 “Tenable security center”,最后一步设置账号密码,账号密码没要求。 ​​ 2、获取插件包 2.1在命令行模式下(…

Web入门 ——生成验证码

<!DOCTYPE html> <html lang="zh-CN"><head><meta charset="UTF-8"><meta name="viewport" content="width=device-width, initial-scale=1.0"><title>验证码</title><script>let cu…

[OI] 猫树

喵~>_<概述 猫树是一种与线段树类似的数据结构 线段树在解决不带修区间问题时,对于单次查询是 \(O(\log)\) 的,此时运用猫树就可以将单次查询的复杂度降到 \(O(1)\) 合理运用分块的思想,在区间内分块,思考我们在块上查询的过程设黑色的线段为整个区间,红色的线段为查…

往来现金中五花八门的单据

导语: T+往来现金里有很多单据,收款单、收入单、付款单、费用单、其他应收单、其他应付单…… 每张单据还有不同的单据类型,比如收款单这一种里就有普通收款、预收款、直接收款。 新人难免懵,比如下面这位:很多小伙伴会认为收入单就是收款,其实不是这样的原理,看完下面你…

【Ehviewer绿色版】2.0.8.4最新版本下载2024安卓苹果

Ehviewer开发应用程序(App)是一项综合性的工作,涉及从构思到发布等多个环节。以下是开发一个基本应用程序的教程,Ehviewer包括从概念设计到发布的完整流程。Ehviewer本教程将分别介绍 iOS 和 Android 平台的开发过程。 ehviewer官方安装包下载: http://ez.oubaidu.com/ 一…

SqlSugar 基础知识

1、实体特性[SugarColumn(IsPrimaryKey = true)] 标识是否为主键[SugarColumn(IsIdentity = true)] 是否为自增长[SugarColumn(ColumnName = "id")] 对应数据库表里面的某列[SugarColumn(IsIgnore = true)] 忽略熟悉,在ORM会过滤掉[SugarColumn(ColumnDescripti…