题面:

传送门:http://codeforces.com/problemset/problem/475/D

Given a sequence of integers a1, …, an and q queries x1, …, xq on it. For each query xi you have to count the number of pairs (l, r) such that 1 ≤ l ≤ r ≤ n and gcd(al, al + 1, …, ar) = xi.

题目大意:

对于n个数的序列,给出q个询问,问有多少个区间满足区间最大公约数为x

分析:

求区间最大公约数很简单,可以利用Sparse-Table算法O(nlog2n)" role="presentation" style="position: relative;">O(nlog2n)O(nlog2n)预处理,O(1)" role="presentation" style="position: relative;">O(1)O(1)查询

我们只要把标准的ST表的递推公式改一下就可以了

st[i][j]=gcd(st[i][j−1],st[i+2j−1][j−1]" role="presentation" style="position: relative;">st[i][j]=gcd(st[i][j−1],st[i+2j−1][j−1]st[i][j]=gcd(st[i][j−1],st[i+2j−1][j−1]

最暴力的方法是直接枚举所有区间gcd值,再求区间个数,然后将x对应的区间个数保存在一个哈希表中,这种算法时间复杂度为O(n2)" role="presentation" style="position: relative;">O(n2)O(n2),实际上考虑到STL的map使用红黑树实现,查询为O(log2n)" role="presentation" style="position: relative;">O(log2n)O(log2n),时间复杂度还会更高

从枚举的过程中我们可以发现,固定区间左端点l,枚举右端点r时,很多区间的gcd值是和上一个区间相同的,导致这一次枚举是重复且无用的.因此,我们应该用更高效的方法去枚举区间.

因此,我们用二分查找去枚举gcd每次变化的位置,假设变化的位置为k(即gcd([l,k])!=gcd(l,k+1),则值为gcd([l,k])的区间个数会增加k-l+1个.这样重复枚举的情况会被消除.

由数论知识得左端点固定,随着右端点右移,gcd会变化log2n" role="presentation" style="position: relative;">log2nlog2n次!

(区间gcd的值必须是ai的因数,而且每次减少至少要除2(少了一个公因数,最小是2))

所以总时间复杂度为O(nlog2n)" role="presentation" style="position: relative;">O(nlog2n)O(nlog2n)

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#define maxn 100005
#define maxlog 32
#define INF 10000005
using namespace std;
int a[maxn];
inline int gcd(int a,int b) {
return b==0?a:gcd(b,a%b);
}
struct sparse_table {//ST表
int st[maxn][maxlog];
void init(int *a,int n) {
for(int i=1; i<=n; i++) st[i][0]=a[i];
for(int j=1; (1<<j)<=n; j++) {
for(int i=1; i+(1<<j)-1<=n; i++) {
st[i][j]=gcd(st[i][j-1],st[i+(1<<(j-1))][j-1]);
}
}
}
int query(int l,int r) {
if(l==r) return st[l][0];
int k=log2(r-l+1);
return gcd(st[l][k],st[r-(1<<k)+1][k]);
}
};
int n,q;
sparse_table gcda;
map<int,long long>ans;//预处理答案
int bin_search(int L,int R,int sta,int v){//二分查找区间gcd变化的位置的前一位
int l=L,r=R;
int ans=INF;
while(l<=r){
int mid=(l+r)>>1;
if(gcda.query(sta,mid)==v){
l=mid+1;
}else{
r=mid-1;
ans=min(ans,mid-1);
}
}
if(ans==INF) return r;
else return ans;
}
int main() {
// freopen("input.txt","r",stdin);
scanf("%d",&n);
for(int i=1; i<=n; i++) scanf("%d",&a[i]);
gcda.init(a,n);
ans.clear();
int now_gcd,now_pos;
for(int i=1; i<=n; i++) {
now_gcd=a[i];
now_pos=i;
while(1) {
int pre_pos=now_pos;//原来的位置
now_pos=bin_search(pre_pos,n,i,now_gcd);//变化的位置,下一个数gcd不同
now_gcd=gcda.query(i,now_pos);//原来的gcd
ans[now_gcd]+=now_pos-pre_pos+1;//更新答案
if(now_pos<n) {
now_pos++;//左端点移动到下一个数,开始枚举新区间
now_gcd=gcda.query(i,now_pos);//更新新的gcd
} else break; }
}
scanf("%d",&q);
int x;
while(q--) {
scanf("%d",&x);
printf("%I64d\n",ans[x]);
}
}

Codeforces 475D 题解(二分查找+ST表)的更多相关文章

  1. P7599-[APIO2021]雨林跳跃【二分,倍增,ST表】

    正题 题目链接:https://www.luogu.com.cn/problem/P7599 题目大意 \(n\)棵树,在某棵树上时可以选择向左右两边第一棵比它高的树跳,现在\(q\)次询问从\([A ...

  2. bzoj3277 串 (后缀数组+二分答案+ST表)

    常见操作:先把所有串都连到一起,但中间加上一个特殊的符号(不能在原串中/出现过)作为分割 由于全部的子串就等于所有后缀的所有前缀,那我们对于每一个后缀,去求一个最长的前缀,来满足这个前缀在至少K个原串 ...

  3. 【CSP膜你赛】柠檬的密码(manacher 二分 单调性 st表)

    题目描述 Lemon觉得他需要一个复杂的密码来保证他的帐号的安全.他经过多日思考,决定使用一个长度为奇数的回文串来作为他的密码.  但是这个回文串太长了,Lemon记不住,于是Lemon决定把它记在本 ...

  4. Codeforces 873E Awards For Contestants ST表

    原文链接https://www.cnblogs.com/zhouzhendong/p/9255885.html 题目传送门 - CF873E 题意 现在要给 $n(n\leq 3000)$ 个学生颁奖 ...

  5. [BZOJ3277/BZOJ3473] 串 - 后缀数组,二分,双指针,ST表,均摊分析

    [BZOJ3277] 串 Description 现在给定你n个字符串,询问每个字符串有多少子串(不包括空串)是所有n个字符串中至少k个字符串的子串(注意包括本身). Solution 首先将所有串连 ...

  6. Codeforces 803G Periodic RMQ Problem ST表+动态开节点线段树

    思路: (我也不知道这是不是正解) ST表预处理出来原数列的两点之间的min 再搞一个动态开节点线段树 节点记录ans 和标记 lazy=-1 当前节点的ans可用  lazy=0 没被覆盖过 els ...

  7. 2019.03.04 bzoj5308: [Zjoi2018]胖(二分答案+st表)

    传送门 想题5分钟调题两小时系列 其实还是我tcl 读完题之后自然会知道一个关键点能够更新的点是一段连续的区间,于是我们对于每个点能到的左右区间二分答案,用ststst表维护一下查询即可. 代码: # ...

  8. 2018.10.14 NOIP训练 直线(二分答案+st表+切比雪夫距离转化)

    传送门 二分答案好题. 这已经是当年普及组模拟时挖的坑了233. 这道题还是很不错的. 考虑把坐标系转个45度再操作. 为了不爆精度可以直接转切比雪夫距离. 然后就直接二分答案. 其中竖线就按二分的答 ...

  9. 说一说ST表 讲一讲水题

    ST表 一.算法介绍 如何快速求解RMQ问题呢?暴力复杂度O(n),线段树复杂度O(n)~O(logn),要是数据规模达到10^7或者更高呢?我们需要一种可以做到O(1)查询的算法,这时就可以用到ST ...

随机推荐

  1. django开发环境搭建(参考流程)

    django开发环境搭建(参考流程) 2013-08-08 01:09:06 分类: LINUX 原文地址:django开发环境搭建(参考流程) 作者:bailiangcn 对于一个初学者,在实际的开 ...

  2. VUE+Ionic,项目搭建&打包成APK

    安装Vue&创建Vue工程 1.安装Vue CLI: npm install -g vue-cli 2.创建新的Vue项目,创建最后一步会提醒是否使用npm install 自动安装,如果选择 ...

  3. 封装插件并发布到npm的历程

    1.封装插件 2.使用webpack打包的时候发生了很多问题 3.过程很反复 4.但最终还是发布成功了

  4. 【微信小程序】setData的使用以及注意事项

    Page.prototype.setData(Object data, Function callback) setData 函数用于将数据从逻辑层发送到视图层(异步),同时改变对应的 this.da ...

  5. android 7.0适配(总结)

    file_paths.xml <?xml version="1.0" encoding="utf-8"?><paths xmlns:andro ...

  6. jQuery 遍历 - 同胞(siblings)

    同胞拥有相同的父元素. 通过 jQuery,您能够在 DOM 树中遍历元素的同胞元素. 在 DOM 树中水平遍历 有许多有用的方法让我们在 DOM 树进行水平遍历: siblings() next() ...

  7. Acvitivi网关(十一)

    1排他网关 1.1 什么是排他网关 排他网关(也叫异或(XOR)网关,或叫基于数据的排他网关),用来在流程中实现决策. 当流程执行到这个网关,所有分支都会判断条件是否为 true,如果为 true 则 ...

  8. HDU - 6578 Blank DP + 滚动数组

    HDU - 6578 Blank 题意 给你\(\{0,1,2, 3\}\)四个数,分别填入长度为\(n\)的数列中,有\(m\)个限制条件,\(l_{i}, r_{i}, x_{i}\)表示在\([ ...

  9. 三、angular7登录请求和路由带参传递

    在 app.module.ts 中引入 HttpClientModule 并注入 import {HttpClientModule} from '@angular/common/http'; impo ...

  10. SpringBoot:使用IDEA快速构建项目

    西部开源-秦疆老师:基于SpringBoot 2.1.6 的博客教程 秦老师交流Q群号: 664386224 未授权禁止转载!编辑不易 , 转发请注明出处!防君子不防小人,共勉! SpringBoot ...