NC15557 连续区间的最大公约数
题目
题目描述
给一个数列共n(n<=100,000)个数,a1,a2,...,an.(0<=ai<=1000,000,000).有q(q<=100,000)个询问。每个询问为l,r(1<=l<=r<=n).求gcd(al,al+1,...,ar).
再求区间[l,r]的子区间中(l<=l'<=r'<=r)满足gcd(al,al+1,...,ar) = gcd(al',al'+1,...ar')的子区间个数.
输入描述
第一行一个数T表示数据组数
第二行一个数n
接下来一行n个数,a1,a2,...,an
接下一行一个数q
接下来一行2个数l和r。
输出描述
首先输出一行“Case #:t”,t代表当前是第几组数据。
接下来q行,每行输出2个数,第一个是gcd(al,al+1,...,ar),
第二个是区间[l,r]的子区间中(l<=l'<=r'<=r)满足gcd(al,al+1,...,ar) = gcd(al',al'+1,...ar')的子区间个数。
示例1
输入
2
5
1 2 4 6 7
4
1 5
2 4
3 4
4 4
5
1 2 4 2 1
6
1 1
1 3
2 2
2 3
2 4
3 3
输出
Case #1:
1 8
2 4
2 1
6 1
Case #2:
1 1
1 3
2 1
2 2
2 5
4 1
题解
知识点:线段树,GCD与LCM,双指针,枚举。
区间信息需要维护区间 \(\gcd\) \(GCD\)、区间答案 \(ans\) 。但是,仅通过孩子区间 \(ans\) 是无法获得区间 \(ans\) 的,因为符合条件的子区间可能还包括区间中凭空产生的跨越两个孩子区间的子区间,即和孩子区间 \(ans\) 中的子区间没有任何关系的子区间。因此,我们需要维护额外的信息。
首先,考虑如何从区间中间开始,枚举出跨越两个孩子区间的符合条件的子区间。最朴素的想法是 \(O(n^2)\) 枚举每个端点,但是显然是不行的,我们考虑如何优化。
注意到,子区间 \(\gcd\) 随着区间变大一定是不增的,而区间 \(\gcd\) 一定是最小的,因此我们可以采用尺取法。初始时,左右端点 \(l,r\) 分别在左右子区间的左端点,随后枚举 \(l\) 扩展 \(r\) ,每次 \(r\) 扩展使得 \([l,r]\) 的 \(\gcd\) 变小直到等于区间 \(gcd\) ,那么 \(r\) 到区间右端点的所有点都是当前 \(l\) 的合法右端点(因此还需要维护区间长度 \(len\) ),最后 \(l\) 向右挪一个点,直到 \(l,r\) 到左右子区间的右端点结束。但是,这样的复杂度是 \(O(n)\) 的,显然还是不够的,我们需要继续优化。
我们发现,尺取的过程中 \(\gcd\) 变化频率通常是很小的。实际上, \(l,r\) 的移动有很长一段距离, \([l,r]\) 的 \(\gcd\) 是不变的。这是因为对于一个区间,其 \(gcd\) 的变化种类最多只有其区间最大值的质因数个种类,在 \(10^9\) 内不超过 \(10\) 个。因此,我们考虑维护区间从左端点开始的前缀 \(\gcd\) 数组 \(L\) ,从右端点开始的后缀 \(\gcd\) 数组 \(R\) 。例如, 区间有 \(7\) 个数 \(\{ 24,12,16,6,30,18,9\}\) ,那么其 \(L = \{(24,1),(12,1),(4,1),(2,3),(1,1) \}\) , \(R=\{ (9,2),(3,2),(1,3) \}\) (其中第一个数表示 \(\gcd\) ,第二个数表示持续长度) 。于是,尺取法的端点枚举这两个数字即可,最多移动 \(10^2\) 次。
合并时, \(len,GCD\) 直接求, \(ans\) 等于尺取法的答案再加上 \(GCD\) 和区间 \(GCD\) 一致的左右子区间 \(ans\) 。对于 \(L,R\) ,我们不妨考虑 \(L\) ,首先继承左子区间的 \(L\) ,然后逐步加上右子区间的 \(R\) ,如果要加上的一块是当前 \(L\) 末尾 \(\gcd\) 的倍数则直接加在 \(L\) 末尾的长度上,否则 push_back
一个新元素, \(R\) 同理。
时间复杂度 \(O((n+q) \log n)\)
空间复杂度 \(O(n)\)
代码
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
template<class T>
class SegmentTree {
int n;
vector<T> node;
T query(int rt, int l, int r, int x, int y) {
if (r < x || y < l) return T();
if (x <= l && r <= y) return node[rt];
int mid = l + r >> 1;
return query(rt << 1, l, mid, x, y) + query(rt << 1 | 1, mid + 1, r, x, y);
}
public:
SegmentTree(int _n = 0) { init(_n); }
SegmentTree(const vector<T> &src) { init(src); }
void init(int _n = 0) {
n = _n;
node.assign(n << 2, T());
}
void init(const vector<T> &src) {
assert(src.size() >= 2);
init(src.size() - 1);
function<void(int, int, int)> build = [&](int rt, int l, int r) {
if (l == r) return node[rt] = src[l], void();
int mid = l + r >> 1;
build(rt << 1, l, mid);
build(rt << 1 | 1, mid + 1, r);
node[rt] = node[rt << 1] + node[rt << 1 | 1];
};
build(1, 1, n);
}
T query(int x, int y) { return query(1, 1, n, x, y); }
};
struct T {
int len = 0;
int GCD = 0;
vector<pair<int, int>> L = {}, R = {};
ll ans = 0;
friend T operator+(const T &a, const T &b) {
if (!a.len) return b;
if (!b.len) return a;
T x = T();
x.len = a.len + b.len;
x.GCD = gcd(a.GCD, b.GCD);
if (a.GCD == x.GCD) x.ans += a.ans;
if (b.GCD == x.GCD) x.ans += b.ans;
int l = a.R.size() - 1, r = 0, sum = 0;
while (l >= 0) {
while (r < b.L.size()) {
if (gcd(a.R[l].first, b.L[r].first) == x.GCD) break;
sum += b.L[r].second;
r++;
}
x.ans += 1LL * a.R[l].second * (b.len - sum);
l--;
}
x.L = a.L;
for (int i = 0;i < b.L.size();i++) {
if (b.L[i].first % x.L.back().first) x.L.push_back({ gcd(x.L.back().first,b.L[i].first),b.L[i].second });
else x.L.back().second += b.L[i].second;
}
x.R = b.R;
for (int i = 0;i < a.R.size();i++) {
if (a.R[i].first % x.R.back().first) x.R.push_back({ gcd(a.R[i].first,x.R.back().first),a.R[i].second });
else x.R.back().second += a.R[i].second;
}
return x;
}
};
bool solve() {
int n;
cin >> n;
vector<T> a(n + 1);
for (int i = 1;i <= n;i++) {
int x;
cin >> x;
a[i] = { 1,x,{{x,1}},{{x,1}},1 };
}
SegmentTree<T> sgt(a);
int q;
cin >> q;
while (q--) {
int l, r;
cin >> l >> r;
auto tar = sgt.query(l, r);
cout << tar.GCD << ' ' << tar.ans << '\n';
}
return true;
}
int main() {
std::ios::sync_with_stdio(0), cin.tie(0), cout.tie(0);
int t = 1;
cin >> t;
for (int i = 1;i <= t;i++) {
cout << "Case #" << i << ":" << '\n';
if (!solve()) cout << -1 << '\n';
}
return 0;
}
NC15557 连续区间的最大公约数的更多相关文章
- 2018年长沙理工大学第十三届程序设计竞赛 I 连续区间的最大公约数
连续区间的最大公约数 思路:参照BZOJ 4488: [Jsoi2015]最大公约数脑补出的一个\(map\)套\(vector\)的写法,写起来比线段树短,运行时间比线段树快. 代码: #pragm ...
- 长沙理工校赛I题题解-连续区间的最大公约数
题目来源https://www.nowcoder.com/acm/contest/96/I 解题前们需要先知道几个结论: 首先,gcd是有区单调性的: gcd(L,R)>=gcd(L,R+d) ...
- CERC 2013 Magical GCD
题目大意如下:给定一个序列,每个序列有值xi,现给定t个数列,对于每个长n的数列,求一段[l,r]使 [r-l+1]*gcd(l,r)最大,gcd(l,r)指的是该连续区间的最大公约数. 不难想到n^ ...
- C语言辗转相除法求2个数的最小公约数
辗转相除法最大的用途就是用来求两个数的最大公约数. 用(a,b)来表示a和b的最大公约数. 有定理: 已知a,b,c为正整数,若a除以b余c,则(a,b)=(b,c). (证明过程请参考其它资料) 例 ...
- 求两个数字的最大公约数-Python实现,三种方法效率比较,包含质数打印质数的方法
今天面试,遇到面试官询求最大公约数.小学就学过的奥数题,居然忘了!只好回答分解质因数再求解! 回来果断复习下,常用方法辗转相除法和更相减损法,小学奥数都学过,很简单,就不细说了,忘了的话可以百度:ht ...
- BZOJ4488: [Jsoi2015]最大公约数
Description 给定一个长度为 N 的正整数序列Ai对于其任意一个连续的子序列{Al,Al+1...Ar},我们定义其权值W(L,R )为其长度与序列中所有元素的最大公约数的乘积,即W(L,R ...
- 求N个数的最大公约数和最小公倍数(转)
除了分解质因数,还有另一种适用于求几个较小数的最大公约数.最小公倍数的方法 下面是数学证明及算法实现 令[a1,a2,..,an] 表示a1,a2,..,an的最小公倍数,(a1,a2,..,an)表 ...
- 辗转相除法求最大公约数,非goto
#include<iostream> using namespace std; //不推荐用goto,当然用它更快 //辗转相除法求两数的最大公约数 int gcd(long int a, ...
- ZOJ Problem Set - 1337 Pi 最大公约数
这道题目的关键在于怎么求两个整数的最大公约数,这里正好复习一下以前的知识,如下: 1.设整数a和b 2.如果a和b都为0,则二者的最大公约数不存在 3.如果a或b等于0,则二者的最大公约数为非0的一个 ...
- Euclid求最大公约数
Euclid求最大公约数算法 #include <stdio.h> int gcd(int x,int y){ while(x!=y){ if(x>y) x=x-y; else y= ...
随机推荐
- Jinfo 查看 jvm 配置及使用 Jstat 查看堆内存使用与垃圾回收
本文为博主远传,未经允许不得转载: 1. Jinfo 查看正在运行的Java应用程序的扩展参数: 包含 JVM 参数与 java 系统参数 命令: jinfo pid 2. 使用 jstat 查看堆 ...
- 基于Java 的商城网站系统设计与实现(8000字论文)
摘要 随着我国经济活力的不断提升和互联网的快速发展,信息的重要性正在显现出来.电子商务作为经济发展的重要一环取得了突飞猛进的发展.由于具有高效便捷的优点,网上购物已经成为一种不可或缺的新型生活方式,近 ...
- Ubuntu Linux下的PDF阅读器推荐——Okular
安装方法 在Ubuntu下直接使用sudo apt-get install okular即可,如果中间遇到依赖项的问题,可以通过运行sudo apt --fix-broken install来自动修复 ...
- android应用申请加入电池优化白名单
首先,在 AndroidManifest.xml 文件中配置一下权限: 1 <uses-permission android:name="android.permission.REQU ...
- Go-连接redis
- 海思Hi35xx 实现本地和远程升级程序的方法
前言 嵌入式linux设备要进行软件升级有很种多方式方法,总的来说可以分为本地升级和远程升级. 本地升级包括升级工具升级,存储介质升级等,远程升级是指通过网络进行程序升级. 这里介绍一种同时至此本地和 ...
- 鲲鹏920上面 Docker 部署 clickhouse 的方式方法
鲲鹏920上面 Docker 部署 clickhouse 的方式方法 背景 最近有一套鲲鹏920的测试环境, 研发同事想纯Dcoker部署一套环境. 其中就包括了 Clickhouse 之前发现Cli ...
- [转帖]金仓数据库KingbaseES数据目录结构
KingbaseES数据库结构 [kingbase@postgresV8]$ tree -LP2data/ . ├── data │ ├── base # 存储用户创建的数据库文件及隶属于用户数据库的 ...
- [转帖]Sosreport:收集系统日志和诊断信息的工具
https://zhuanlan.zhihu.com/p/39259107 如果你是 RHEL 管理员,你可能肯定听说过 Sosreport :一个可扩展.可移植的支持数据收集工具.它是一个从类 Un ...
- Find -mtime 的图解