CF1422F Boring Queries(ST表 + 主席树)

CF1422F Boring Queries

给定一个长度为nnn的数组a,(1≤ai≤2×105)a,(1 \leq a_i \leq 2 \times 10 ^ 5)a,(1ai2×105),有mmm次询问,每次询问给定l,rl, rl,r,要我们求区间[l,r][l, r][l,r]aia_iailcmlcmlcm,强制在线。

对于质因子ppp,如果p2>2×105p ^ 2 > 2 \times 10 ^ 5p2>2×105,那么它在每个数中出现的次数只有两种可能0/10/10/1,可以如HH的项链的那题,一样对每个因子维护一个最后出现的乘积。

如果p2≤2×105p ^ 2 \leq 2 \times 10 ^ 5p22×105,我们可以考虑用RMQRMQRMQ表来维护区间最大值,容易知道最多只要维护868686RMQRMQRMQ表,所以整体复杂度O(86×nlog⁡n)O(86 \times n \log n)O(86×nlogn)

#include <bits/stdc++.h>using namespace std;const int N = 1e5 + 10, M = 2e5 + 10, mod = 1e9 + 7;int root[N], ls[N * 40], rs[N * 40], sum[N * 40], num;int prime[N], Log[N], inv[M], pre[M], fac[M], a[N], cnt, n, m;vector<int> p[87];bool st[M];struct RMQ {char f[N][18];void init() {for (int j = 1; j < 18; j++) {for (int i = 1; i + (1 << j) - 1 <= n; i++) {f[i][j] = max(f[i][j - 1], f[i + (1 << j - 1)][j - 1]);}}}int query(int l, int r) {int s = Log[r - l + 1];return max(f[l][s], f[r - (1 << s) + 1][s]);}
}rmq[87];void init() {inv[1] = 1;for (int i = 2; i < M; i++) {inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;if (!st[i]) {prime[++cnt] = i;}for (int j = 1; j <= cnt && 1ll * i * prime[j] < M; j++) {st[i * prime[j]] = 1;if (i % prime[j] == 0) {break;}}}for (int i = 2; i < N; i++) {Log[i] = Log[i / 2] + 1;}for (int j = 1; prime[j] * prime[j] < M; j++) {for (int i = 1; i <= n; i++) {while (a[i] % prime[j] == 0) {a[i] /= prime[j];rmq[j].f[i][0]++;}}rmq[j].init();for (int cur = 1; cur < M; cur *= prime[j]) {p[j].push_back(cur);}}for (int i = 87; i <= cnt; i++) {for (int j = prime[i]; j < M; j += prime[i]) {fac[j] = prime[i];}}
}void build(int &rt, int l, int r) {rt = ++num;if (l == r) {sum[rt] = 1;return ;}int mid = l + r >> 1;build(ls[rt], l, mid);build(rs[rt], mid + 1, r);sum[rt] = sum[ls[rt]] * sum[rs[rt]];
}void update(int &rt, int pre, int l, int r, int x, int v) {rt = ++num, ls[rt] = ls[pre], rs[rt] = rs[pre], sum[rt] = 1ll * sum[pre] * v % mod;if (l == r) {return ;}int mid = l + r >> 1;if (x <= mid) {update(ls[rt], ls[pre], l, mid, x, v);}else {update(rs[rt], rs[pre], mid + 1, r, x, v);}
}int query(int rt, int l, int r, int L, int R) {if (l >= L && r <= R) {return sum[rt];}int ans = 1, mid = l + r >> 1;if (L <= mid) {ans = 1ll * ans * query(ls[rt], l, mid, L, R) % mod;}if (R > mid) {ans = 1ll * ans * query(rs[rt], mid + 1, r, L, R) % mod;}return ans;
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);scanf("%d", &n);for (int i = 1; i <= n; i++) {scanf("%d", &a[i]);}init();build(root[0], 1, n);for (int i = 1; i <= n; i++) {root[i] = root[i - 1];if (!fac[a[i]]) {continue;}update(root[i], root[i], 1, n, i, fac[a[i]]);if (pre[fac[a[i]]]) {update(root[i], root[i], 1, n, pre[fac[a[i]]], inv[fac[a[i]]]);}pre[fac[a[i]]] = i;}scanf("%d", &m);int ans = 0;for (int i = 1, l, r; i <= m; i++) {scanf("%d %d", &l, &r);l = (ans + l) % n + 1, r = (ans + r) % n + 1;if (l > r) {swap(l, r);}ans = 1;for (int j = 1; j <= 86; j++) {ans = 1ll * ans * p[j][rmq[j].query(l, r)] % mod;}ans = 1ll * ans * query(root[r], 1, n, l, r) % mod;printf("%d\n", ans);}return 0;
}

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

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

相关文章

高性能最终一致性框架Ray之基本概念原理

一、Actor介绍Actor是一种并发模型&#xff0c;是共享内存并发模型的替代方案。共享内存模型的缺点&#xff1a;共享内存模型使用各种各样的锁来解决状态竞争问题&#xff0c;性能低下且让编码变得复杂和容易出错。共享内存受限于单节点的服务器资源限制。Actor模型的优点&…

B.The Tortoise and the Hare 长春

B. The Tortoise and the Hare 给定一个长度为nnn的数组a,(1≤ai<m)a, (1 \leq a_i < m)a,(1≤ai​<m)&#xff0c;mmm是一个给定的数(1≤m≤109)(1 \leq m \leq 10 ^ 9)(1≤m≤109)&#xff0c;有QQQ次操作&#xff0c;分为两类&#xff1a; 给定u,v,(1≤u≤n,1≤v…

漫画:程序员带娃日常(1)

我是一名IT界资深人士做过乙方苦哈哈做过甲方做过项目、搞点管理做过培训、也上台演讲不过现在最大挑战、交期最长的项目是。。。带 娃娃子啊&#xff0c;爸爸也是第一次当爸爸&#xff0c;要请你多多指教啦&#xff01;1“小刘&#xff0c;把这个需求细化一下。”“兄弟们&…

L. Coordinate Paper(CCPC 长春)构造

L. Coordinate Paper 构造一个长度为nnn的序列aaa&#xff0c;满足ai≥0a_i \geq 0ai​≥0&#xff0c;∑i1nais\sum\limits_{i 1} ^{n} a_i si1∑n​ai​s&#xff0c;对于任意的i∈[1,n−1]i \in [1, n - 1]i∈[1,n−1]&#xff0c;都有ai−ai1korai1−ai1a_i - a_{i 1} …

servlet生命周期

1.servlet的生命周期 主要有三个方法&#xff1a; init()初始化阶段service()处理客户端请求阶段destroy()终止阶段 初始化阶段&#xff1a; Servlet容器加载Servlet&#xff0c;加载完成后&#xff0c;Servlet容器会创建一个Servlet实例并调用init()方法&#xff0c;init(…

表达式树练习实践:入门基础

什么是表达式树来自微软官方文档的定义&#xff1a;表达式树以树形数据结构表示代码。它能干什么呢&#xff1f;你可以对表达式树中的代码进行编辑和运算。这样能够动态修改可执行代码、在不同数据库中执行 LINQ 查询以及创建动态查询。好不好玩&#xff1f;表达式树还能用于动…

Servlet与线程安全

Servlet与线程安全 先说结论&#xff0c;Servlet本身是单例的&#xff0c;线程安全的。但是如果引入共享变量&#xff0c;则可能会变得线程不安全。 1. 什么是线程安全 首先说明一下对线程安全的讨论&#xff0c;哪种情况我们可以称作线程安全&#xff1f; 《Java并发编程实…

P3591 [POI2015]ODW(分块)

P3591 [POI2015]ODW 给定一颗有nnn个节点的树&#xff0c;点有点权&#xff0c;给定一个长度为nnn的排列ppp&#xff0c;给定一个长度为n−1n - 1n−1的数组ccc&#xff0c; 我们会在树上进行n−1n - 1n−1次行走&#xff0c;第iii次我们会从p[i]p[i]p[i]走到p[i1]p[i 1]p[i…

尝鲜体验 VS Code Python 原生 Jupyter Notebook 支持

9 月 21 日&#xff0c;PyCon China 2019 在上海举办。微软的 VS Code 发布了原生的Juypter Notebook 支持。 虽然要等到月底 Marketplace 里面才会正式更新&#xff0c;但 GitHub 上的 repo 已经包含了相应的 commit&#xff0c;所以还是可以通过手动安装的方式尝鲜。现有 VSC…

Java Web之filter、listener、Interceptor

** 1、Servlet ** Servlet 是服务端的 Java 应用程序&#xff0c;用于处理HTTP请求&#xff0c;做出相应的响应。 当客户端向服务器发出HTTP请求时&#xff0c;首先会由服务器中的 Web 容器&#xff08;如Tomcat&#xff09;对请求进行路由&#xff0c;交给该URL对应的 Servl…

L - Lookup Performance(主席树)

L - Lookup Performance 问对于一颗二叉搜索树来说&#xff0c;如果我们要找一个值域区间的值有多少个&#xff0c;他会向下递归查找几次&#xff0c; 设&#xff0c;第iii个节点所代表的最大最小值为li,ril_i, r_ili​,ri​&#xff0c;此时我们要查询L,RL, RL,R之间的值有多…

abp vNext微服务框架分析

abp vNext新框架的热度一直都很高&#xff0c;于是最近上手将vNext的微服务Demo做了一番研究。我的体验是&#xff0c;vNext的微服务架构确实比较成熟&#xff0c;但是十分难以上手&#xff0c;对于没有微服务开发经验的.net人员来说几乎是看不懂的&#xff0c;所以研究一番后再…

Java web之web.xml配置详解

什么是web.xml web.xml是web项目的配置文件&#xff0c;一般的web工程都会用到web.xml来配置&#xff0c;方便大型开发。web.xml主要用来配置Filter&#xff0c;Listener&#xff0c;Servlet等。但是web.xml并不是必须的&#xff0c;一个web工程可以没有web.xml文件。 web工程…

.NET Conf 2019日程(北京时间)

一年一度的 .NET Conf马上就要开始了,我将日程简易的翻译了一下,并且时间全部转换为北京时间,以方便国内.NETer.第1天 (北京时间9月24日).NET Conf 2019 基调 - Scott Hunter Mads Torgersen James Montemagno Olia Gavrysh Daniel Roth Glenn Condron Bri Achtman欢迎来到 .NE…

C - Swaps 2(树状数组,思维)

C - Swaps 2 给定两个长度为nnn的数组A,BA, BA,B&#xff0c;我们可以进行若干次如下操作&#xff0c;使AAA变成BBB&#xff0c; 选定i<ni < ni<n&#xff0c;将aia_iai​减小111&#xff0c;将ai1a_{i 1}ai1​增加111。交换ai,ai1a_i, a_{i 1}ai​,ai1​。 问我…

Hibernate基础

** 一、什么是Hibernate&#xff1f; **   Hibernate是一个轻量级的ORMapping框架   ORMapping原理&#xff08;Object Relational Mapping&#xff09; ORMapping基本对应规则&#xff1a; 1&#xff1a;类跟表相对应2&#xff1a;类的属性跟表的字段相对应3&#xff1…

进击的.NET 在云原生时代的蜕变

你一定看过这篇文章 《进击的 Java &#xff0c;云原生时代的蜕变》, 本篇文章的灵感来自于这篇文章。北京时间9.24 就将正式发布.NET Core 3.0, 所以写下这篇文章让大家全面认识.NET Core。.NET 生态系统是一个不断变化的生态圈&#xff0c;我相信它正在朝着一个伟大的方向发…

F. Strange Array(Codeforces Round #727 (Div. 2))(主席树)

F. Strange Array 给定一个长度为nnn的数组aaa&#xff0c;1≤ai≤n1 \leq a_i \leq n1≤ai​≤n&#xff0c;对于每个aia_iai​&#xff0c;我们要找到一个l≤i,r≥il \leq i, r \geq il≤i,r≥i&#xff0c; 使得&#xff0c;我们对区间[l,r][l, r][l,r]升序后&#xff0c;…

Hibernate懒加载问题的5种解决方案

** 1、Hibernate基础 ** Hibernate基础&#xff0c;传送门 ** 2、什么是Hibernate懒加载 ** 当我们查询一个对象的时候&#xff0c;在默认情况下&#xff0c;返回的只是该对象的代理对象&#xff0c;当用户去使用该对象的属性时&#xff0c;才会向数据库再一次发出查询语…

程序员过关斩将--要想获取我的用户信息,就得按照规矩来

菜菜君&#xff0c;我又来啦又有什么事吗&#xff1f;我按照你上篇文章写的JWT的方式已经把网站认证写完了&#xff0c;而且效果还不错那恭喜你呀&#xff0c;下次面试又多了一项技能不过&#xff0c;现在又有一个问题&#xff0c;我做的系统有一个合作商想要利用我们的用户信息…