当前位置: 首页 > news >正文

P1903 [国家集训队] 数颜色 / 维护队列 Solution

Description

给定序列 a = ( a 1 , a 2 , ⋯ , a n ) a=(a_1,a_2,\cdots,a_n) a=(a1,a2,,an),有 m m m 个操作分两种:

  • modify ⁡ ( i , x ) \operatorname{modify}(i,x) modify(i,x):执行 a i ← x a_i\gets x aix.
  • query ⁡ ( l , r ) \operatorname{query}(l,r) query(l,r):求 ∣ { a l , a l + 1 , ⋯ , a r } ∣ |\{a_l,a_{l+1},\cdots,a_r\}| {al,al+1,,ar}.

Limitations

1 ≤ n , m ≤ 133333 1\le n,m\le 133333 1n,m133333
1 ≤ a i , x ≤ 1 0 6 1 \le a_i,x \le 10^6 1ai,x106
2.5 s , 512 MB 2.5\text{s},512\text{MB} 2.5s,512MB

Solution

依然是不好直接做的询问,考虑莫队,维护 cnt i \textit{cnt}_i cnti 表示当前区间内 i i i 出现次数.
adddel 时更新 cnt \textit{cnt} cnt,若一个数新出现或消失,我们更新答案.
但这里带了修改,我们给莫队直接加上时间维,每个查询多存一个 上次修改的时间戳.
然后需要修改与回溯函数.
拿修改讲,如果 i i i 在当前区间内,那它就会影响答案,我们 del 掉原来的数,add 上新的即可.

有几点要提一下:

  • 修改完后可以将 a i a_i ai x x x 交换,下次遇到是就是回溯,可以省一个函数.
  • 不能用 map 维护 cnt \textit{cnt} cnt,直接用数组.
  • 要调块长,取 B = n 2 3 B=n^{\frac{2}{3}} B=n32 可以过,笔者试了好几发.

Code

1.79 KB , 3.93 s , 8.99 MB (in total, C++20 with O2) 1.79\text{KB},3.93\text{s},8.99\text{MB}\;\texttt{(in total, C++20 with O2)} 1.79KB,3.93s,8.99MB(in total, C++20 with O2)

#include <bits/stdc++.h>
using namespace std;using i64 = long long;
using ui64 = unsigned long long;
using i128 = __int128;
using ui128 = unsigned __int128;
using f4 = float;
using f8 = double;
using f16 = long double;template<class T>
bool chmax(T &a, const T &b){if(a < b){ a = b; return true; }return false;
}template<class T>
bool chmin(T &a, const T &b){if(a > b){ a = b; return true; }return false;
}
const int L = 1e6 + 10;
struct Modify { int pos, x; };
struct Query { int l, r, t, id; }; signed main() {ios::sync_with_stdio(0);cin.tie(0), cout.tie(0);int n, m;cin >> n >> m;vector<int> a(n);for (int i = 0; i < n; i++) cin >> a[i];vector<Modify> M;vector<Query> Q;for (int i = 0, x, y; i < m; i++) {char op;cin >> op >> x >> y, x--;if (op == 'R') M.push_back({x, y});else y--, Q.push_back({x, y, M.size() - 1, Q.size()});}const int B = pow(n, 2.0 / 3);sort(Q.begin(), Q.end(), [&](Query &a, Query &b) {if (a.l / B == b.l / B) {if (a.r / B == b.r / B) return a.t < b.t;return a.r / B < b.r / B;}return a.l / B < b.l / B;});array<int, L> cnt{};int ans = 0, l = 0, r = -1, t = -1;auto add = [&](int v) {if (cnt[v]++ == 0) ans++;};auto del = [&](int v) {if (--cnt[v] == 0) ans--;};auto upd = [&](int t, int l, int r) {auto &[pos, x] = M[t];if (l <= pos && pos <= r) del(a[pos]), add(x);swap(a[pos], x);};vector<int> res(Q.size());for (auto [ql, qr, qt, id] : Q) {while (ql < l) add(a[--l]);while (r < qr) add(a[++r]);while (l < ql) del(a[l++]);while (qr < r) del(a[r--]);while (t < qt) upd(++t, l, r);while (qt < t) upd(t--, l, r);res[id] = ans;}for (int i = 0; i < Q.size(); i++) cout << res[i] << endl;return 0;
}

相关文章:

  • 【趣谈】Cyber、Web、Network都是网络有什么区别
  • WinSW注册服务
  • Windows环境下用pyinstaller将python脚本编译为exe文件
  • 使用pip3安装软件包报错`externally-managed-environment`的几种解决方式
  • Java SE(5)——数组
  • Cesium添加WMS,WMTS,地形图图,3D Tiles数据
  • AimRT 从零到一:官方示例精讲 —— 二、HelloWorld示例.md
  • 【C++11】新的类功能、lambda
  • 存储器分类
  • 【网络】HTTP报文首部字段
  • Qt的WindowFlags窗口怎么选?
  • LeetCode 2906 统计最大元素出现至少K次的子数组(滑动窗口)
  • oracle怎样通过固化较优执行计划来优化慢sql
  • ant design pro 项目发布遇到登录页访问404
  • 【免费下载】2012-2023年全国夜间灯光数据
  • 从遍历序列构造二叉树:前序+中序与中序+后序的递归解法详解
  • JavaScript面试问题
  • tbb parallel_for 使用
  • 颜色分类,不靠“调色盘”:双指针 VS 计数排序的正面PK
  • 【Linux】服务自启动设置的方式
  • 全国人民代表大会常务委员会公告〔十四届〕第十号
  • 阿斯利康中国区一季度收入增5%,或面临最高800万美元新罚单
  • 陕西省副省长窦敬丽已任宁夏回族自治区党委常委、统战部部长
  • 广东省副省长刘红兵跨省任湖南省委常委、宣传部部长
  • 中吉乌铁路重点控制性工程开工建设,包括三座隧道
  • 交行一季度净利253.72亿元增1.54%,不良率微降