『模拟赛』多校A层冲刺NOIP2024模拟赛03

Rank

炸了,触底反弹

image

A. 五彩斑斓(colorful)

签,又没签上。

考虑如何一步步优化暴力。最暴力的思想 \(\mathcal{O(n^4)}\) 枚举每个矩形,判断四个顶点颜色。稍微优化些,两次 \(\mathcal{O(n^2)}\) 跑出对于行/列每个点下一个与之颜色相同的坐标,利用容斥全部减去不合法的方案数,然后再枚举每个点,随机数据下跑得很快,颜色数量越少效率越低。

考虑正解,依旧容斥,我们可以直接枚举两行,再枚举列,记录每一个相同颜色组成的点对出现的数量,简单组合计数可知对于 \(n\) 组点对能组成的方案数为 \(\sum_{i=1}^n\ i\),增加出现数量时直接在总 \(ans\) 里减一下就行了,每次统计后清空计数数组,这样复杂度为 \(\mathcal{O(n^3)}\),可以轻松过。

据说有 \(\mathcal{O(n^4)}\) 算法各种大力优化碾过去了。

点击查看代码
#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 = 400 + 5;
const int mod = 998244353;
int n, m;
int c[N][N], tim[N * N];
ll ans;
namespace Wisadel
{short main(){freopen("colorful.in", "r", stdin) , freopen("colorful.out", "w", stdout);n = qr, m = qr;fo(i, 1, n) fo(j, 1, m) c[i][j] = qr;ans = 1ll * n * (n + 1) * m * (m + 1) / 4;fo(i, 1, n) fo(j, i, n){fo(k, 1, m) if(c[i][k] == c[j][k]) tim[c[i][k]]++, ans -= tim[c[i][k]];fo(k, 1, m) if(c[i][k] == c[j][k]) tim[c[i][k]] = 0;}printf("%lld\n", ans);return Ratio;}
}
signed main(){return Wisadel::main();}

B. 错峰旅行(travel)

签,双没签上。

看到第一感觉是动态开点线段树,不过赛时怎么压都过不了最后一个大样例,卡到极限 70pts,自己画蛇添足唐错挂 10pts。

其实是可以用线段树做的。首先很显然的是,对于一段时间 \([s,t]\),如果不拥挤城市数都是 \(n\),那么总方案数为 \(\mathcal{n^{t-s+1}}\)。初始每一天合法城市数都是 \(n\),由于题目说明了同一城市没有交集,所以每个限制条件相当于做了一次区间减的操作。我们用标记永久化的思想,每次正常区间减,不下放标记,最后做一次最坏 \(\mathcal{O(n)}\) 的查询就做完了。

点击查看代码
#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 = 1e6 + 5;
const int mod = 1e9 + 7;
int n, m, s, t, tot, ans = 1;
int v[N * 20], son[N * 20][2];
namespace Wisadel
{int Wqp(ll x, int y){int res = 1;while(y){if(y & 1) res = 1ll * res * x % mod; x = x * x % mod; y >>= 1;}return res;}#define ls (son[rt][0])#define rs (son[rt][1])#define mid ((l + r) >> 1)void Wupd(int rt, int l, int r, int x, int y){if(x <= l && r <= y){v[rt]--;return;}if(x <= mid){if(!ls) ls = ++tot;Wupd(ls, l, mid, x, y);}if(y > mid){if(!rs) rs = ++tot;Wupd(rs, mid + 1, r, x, y);}}void Wq(int rt, int l, int r, int now){if(!rt){ans = 1ll * ans * Wqp(now, r - l + 1) % mod;return;}if(ls || rs) Wq(ls, l, mid, now + v[rt]), Wq(rs, mid + 1, r, now + v[rt]);else ans = 1ll * ans * Wqp(now + v[rt], r - l + 1) % mod;}short main(){freopen("travel.in", "r", stdin) , freopen("travel.out", "w", stdout);n = qr, m = qr, s = qr, t = qr;t -= s - 1;tot = 1, v[1] = n;fo(i, 1, m){int x = qr, l = qr - s + 1, r = qr - s + 1;Wupd(1, 1, t, l, r);}Wq(1, 1, t, 0);printf("%d\n", ans);return Ratio;}
}
signed main(){return Wisadel::main();}

正解给出的是差分的思想,也很简单,有用的信息最多 \(2\times m\) 个,我们记录每个限制信息的左右端点,到左端点处合法城市数减一,到右端点 +1 处合法城市数加一。实现也很简单,pair 记录后 sort 一遍,加入首末时间 \(\mathcal{O(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 = 1e6 + 5;
const int mod = 1e9 + 7;
int n, m, s, t, tot, ans = 1;
pii p[N << 1];
namespace Wisadel
{int Wqp(ll x, int y){int res = 1;while(y){if(y & 1) res = 1ll * res * x % mod; x = x * x % mod; y >>= 1;}return res;}short main(){freopen("travel.in", "r", stdin) , freopen("travel.out", "w", stdout);n = qr, m = qr, s = qr, t = qr;fo(i, 1, m){int x = qr, l = qr, r = qr;p[++tot] = {l, -1}, p[++tot] = {r + 1, 1};}sort(p + 1, p + 1 + tot);p[0] = {s, 0};p[++tot] = {t + 1, 0};fo(i, 0, tot - 1){n += p[i].se;if(p[i].fi != p[i + 1].fi)ans = 1ll * ans * Wqp(n, p[i + 1].fi - p[i].fi) % mod;}printf("%d\n", ans);return Ratio;}
}
signed main(){return Wisadel::main();}

C. 线段树(segment)

签(?),叒没签上。

还以为真是线段树,赛时脑子被禁锢了转不了,暴力都没写出来。

据 wang54321 说这是他们暑假第二场模拟赛原题。

正解区间 dp。考虑维护包含 \([i,j]\) 的询问数量 \(s_{i,j}\),设 \(f_{i,j}\) 表示所有询问在区间 \([i,j]\) 子树内的答案,我们枚举所有区间,然后枚举其中可行的分割点,则有状态转移方程:

\[f_{i,j} = min\left\{f_{i,k}+f_{k +1,r}-s_{i.j} \right\}\ k\in\left(l,r\right) \]

边界处理上,初始化 \(f_{i,i} = s_{i,i}\),其余赋为极大值即可。复杂度 \(\mathcal{O(n^3)}\)

点击查看代码
#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 = 500 + 5;
const int mod = 1e9 + 7;
int n, m;
int sum[N][N], f[N][N];
namespace Wisadel
{short main(){// freopen(".in", "r", stdin) , freopen(".out", "w", stdout);freopen("segment.in", "r", stdin) , freopen("segment.out", "w", stdout);n = qr, m = qr;memset(f, 0x3f, sizeof f);fo(i, 1, m){int a = qr, b = qr;sum[a][b]++;}fo(i, 1, n) fu(j, n, i) sum[i][j] += sum[i][j + 1];fo(j, 1, n) fo(i, 1, j) sum[i][j] += sum[i - 1][j];fo(i, 1, n) f[i][i] = sum[i][i];fo(len, 2, n) fo(l, 1, n - len + 1){int r = l + len - 1;fo(k, l, r - 1)f[l][r] = min(f[l][r], f[l][k] + f[k + 1][r] - sum[l][r]);}printf("%d\n", f[1][n]);return Ratio;}
}
signed main(){return Wisadel::main();}

D. 量子隧穿问题(experiment)

叕没做出来(?)

一个 trick 是计数转概率,答案 \(=\) 总方案数 \(\times\) 合法概率。

这样一来,盒子的状态可以得出有猫的概率,状态转移方程如下:

\[f_i'=f_i\times f_{to_i} \]

\[f_{to_i}'=f_{to_i}+f_i\times (1-f_{to_i}) \]

然后就可以做 dp 了。吗?

显然这么简单就不至于放 T4 了。发现每个隧道有可能连成环,整个图是一个基环树。当 \(n\) 存在与一个环内时转移时就会因为二者概率不独立而出错。因此我们需要断一条边,枚举两点的概率 1/0,然后再根据上面的转移方程就能无后效性地 dp 了。

注意由于 巴甫洛夫的狗 是从 1 到 \(n\) 运动的,我们断的必须是第一个与 \(n\) 在一个环内的盒子的边,可以用并查集辅助解决。

点击查看代码
#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 = 1e9 + 7;
int n, m;
int to[N], p[N], fx[N];
bool in[N];
string s;
namespace Wisadel
{int Wfind(int x){if(x == fx[x]) return x;return fx[x] = Wfind(fx[x]);}short main(){freopen("experiment.in", "r", stdin) , freopen("experiment.out", "w", stdout);int T = qr;while(T--){n = qr; cin >> s;memset(in, 0, sizeof in);fo(i, 1, n) to[i] = qr, fx[i] = i;int st = 0, ans = 0;fo(i, 1, n) fx[Wfind(i)] = Wfind(to[i]);fo(i, 1, n) in[i] = (Wfind(i) == Wfind(n));fo(i, 1, n) fx[i] = i;fu(i, n, 1){int _ = Wfind(i), __ = Wfind(to[i]);if(in[i] && _ != __) fx[_] = __;else if(in[i]) st = i;}fo(x, 0, 1) fo(y, 0, 1){fo(i, 0, n - 1)if(s[i] == 'C') p[i + 1] = 1;else if(s[i] == '.') p[i + 1] = 0;else p[i + 1] = (1 + mod) / 2;int res = 1;fo(i, 1, n){if(i == st){res = 1ll * (x ? p[i] : (1 - p[i])) * (y ? p[to[i]] : (1 - p[to[i]])) % mod;p[i] = x;p[to[i]] = y;}int xx = p[i], yy = p[to[i]];p[i] = 1ll * xx * yy % mod;p[to[i]] = (1ll * yy + 1ll * xx * (1 - yy) % mod) % mod;}ans = (ans + 1ll * res * p[n] % mod) % mod;}fo(i, 0, n - 1) if(s[i] == '?') ans = ans * 2 % mod;printf("%d\n", (ans + mod) % mod);}return Ratio;}
}
signed main(){return Wisadel::main();}

怀疑这个出题人看过 某人 和 某人 的唐博,开幕巴甫洛夫的狗雷击整个机房。

挺烂的,作为 17 岁 的开局有点糟糕,希望低开高走吧。

一整天状态不太好,好在下午及时调整过来了。

一些消极的话,有反转

今天的超链接好像放的过于多了


完结撒花~

image

最后再祝自己生日快乐!

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

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

相关文章

加装spark-3.5.3

集群版本 hadoop-3.4.0 hive-3.1.3 zookeeper-3.9.2 hbase-2.6.0(1.0.0以上需要zookeeper-3.4.0以上) spark-3.5.3(只能选2.13.0) scala-2.13.0(jdk8仅支持x.x.0系)总结一下:JDK8和scala-2.13.0必选。1.安装scala 1.1 下载解压 tar zxvf scala-2.13.0.tgz 1.2 配置环境变…

高级程序语言第二次个人作业

高级程序语言第二次作业这个作业属于哪个课程 https://edu.cnblogs.com/campus/fzu这个作业要求在哪里 https://edu.cnblogs.com/campus/fzu/2024C/homework/13282学号 222200424姓名 赵伟豪编程练习3.13.23.33.43.53.63.73.8示例程序3.13.23.33.43.53.63.73.83.93.10总结与收获…

浏览器的渲染原理

浏览器渲染原理 五个渲染流程Parse 阶段:解析 HTMLStyle 阶段:样式计算三个阶段:收集,划分和索引所有样式表中存在的样式规则 访问每个元素并找到适用于该元素的所有规则,CSS 引擎遍历 DOM 节点,进行选择器匹配,并且匹配的节点执行样式设置 结合层叠规则和其他信息为节点…

CSP2024 前集训:多校A层冲刺NOIP2024模拟赛03

前言T1 没想到正难则反,脑瘫了没敢用 bitset(复杂度擦边但卡常能过),T2 空间开大了挂了 \(100pts\),\(T3\) 是原。 T1 五彩斑斓部分分 \(20pts\):\(O(n^4)\) 暴力。部分分 \(20+?pts\):进行一些优化,极限数据下仍是 \(O(n^4)\)。部分分 \(60\sim 100pts\):bitset 优化…

在C#中使用适配器Adapter模式和扩展方法解决面向的对象设计问题

之前有阵子在业余时间拓展自己的一个游戏框架,结果在实现的过程中发现一个设计问题。这个游戏框架基于MonoGame实现,在MonoGame中,所有的材质渲染(Texture Rendering)都是通过SpriteBatch类来完成的。举个例子,假如希望在屏幕的某个地方显示一个图片材质(imageTexture)…

React Fiber 原理

React Fiber 在 React 16 之前的版本对比更新 VirtualDOM 的过程是采用 Stack 架构实现的,也就是循环加递归,这种方式的问题是一旦任务开始进行就无法被中断。 如果应用中的组件数量庞大, Virtual DOM 的层级比较深,主线程被长期占用,知道整颗 Virtual DOM 树比对更新完成…

视野修炼-技术周刊第104期 | 下一代 JavaScript 工具链

① 🐙 尤大创办公司 VoidZero ② Tauri 2.0 稳定版发布 ③ Vite 时髦的新主页 ④ qrframe - 漂亮二维码生成 ⑤ HTTP QUERY 方法提案 ⑥ TinyJS - 轻量级的创建DOM元素 ⑦ 9月 Web 平台的新功能 ⑧ ESLint 现在正式支持 Linting JSON 和 Markdown欢迎来到第 104 期的【视野修…

雅礼国庆集训 day1 T2 折射

题面 题面下载 算法 转化题意 说白了就是给了你一堆点,让你数这种折线有多少个 (严格向下走,并且横坐标之间的差越来越小)看着像一种在 y 轴方向排序的 dp 但是由于是折线, 所以需要加一维来判断转向 dp 设计 状态设计 \(dp_{i, 0/1}\) 表示第 i 个点, 是向左下还是右上 状态…

React 中的 diff 算法

React diff为什么使用虚拟 DOM ? 浏览器在处理 DOM 的时候会很慢,处理 JavaScript 会很快,页面复杂的时候,频繁操作 DOM 会有很大的性能开销(每次数据变化都会引起整个 DOM 树的重绘和重排)。 为了避免频繁操作 DOM,React 会维护两个虚拟 DOM,如果有数据更新,会借此计…

20222315 2024-2025-1 《网络与系统攻防技术》实验一实验报告

1.实验内容 1.掌握反汇编与十六进制编程器 2.能正确修改机器指令改变程序执行流程 3.能正确构造payload进行bof攻击 2.实验过程 1.直接修改程序机器指令,改变程序执行流程 将pwn1文件下载至kali中并将pwn1文件改名为pwn20222315,并将其内容复制到pwn2反汇编文件objdump -d…

多校A层冲刺NOIP2024模拟赛03

多校A层冲刺NOIP2024模拟赛03\(T1\) A. 五彩斑斓(colorful) \(90/100pts\)部分分\(20pts\) :枚举左上 \((k,h)\) 、右下端点 \((i,j)\) ,时间复杂度为 \(O(n^{2}m^{2})\) 。 \(90/100pts\) :当 \(a_{i,j} \ne a_{k,j}\) 时任意的 \(h \in [1,j]\) 都符合题意、不妨钦定 \(…