HDU 5869 Different GCD Subarray Query rmq+离线+数状数组
题目链接:
http://acm.hdu.edu.cn/showproblem.php?pid=5869
Different GCD Subarray Query
Time Limit: 6000/3000 MS (Java/Others)Memory Limit: 65536/65536 K (Java/Others)
#### 问题描述
> This is a simple problem. The teacher gives Bob a list of problems about GCD (Greatest Common Divisor). After studying some of them, Bob thinks that GCD is so interesting. One day, he comes up with a new problem about GCD. Easy as it looks, Bob cannot figure it out himself. Now he turns to you for help, and here is the problem:
>
> Given an array a of N positive integers a1,a2,⋯aN−1,aN; a subarray of a is defined as a continuous interval between a1 and aN. In other words, ai,ai+1,⋯,aj−1,aj is a subarray of a, for 1≤i≤j≤N. For a query in the form (L,R), tell the number of different GCDs contributed by all subarrays of the interval [L,R].
输入
There are several tests, process till the end of input.
For each test, the first line consists of two integers N and Q, denoting the length of the array and the number of queries, respectively. N positive integers are listed in the second line, followed by Q lines each containing two integers L,R for a query.
You can assume that
1≤N,Q≤100000
1≤ai≤1000000
输出
For each query, output the answer in one line.
样例输入
5 3
1 3 4 6 9
3 5
2 5
1 5
样例输出
6
6
6
题意
给你n个数排成一行,每次询问求一段区间内的所有子串的不同的gcd的种数有多少。
题解
首先,要处理出所有的子区间是比较困难的,而且事实上很多子区间的gcd都是重复出现的,我们把区间右端点固定,那么随着左区间从右往左移,区间gcd的值成倍递减,所以我们取到的不同的gcd的值只有logAi种,处理出来之后,我们就可以用离线的方式,用线段树来求区间不同的值有几个,每个事件(处理出来的sigma(logA[i])个区间gcd)以左端点为准插入树状数组中维护一下(离线查询的线段树/树状数组可以参考[这个])。
代码
#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<stack>
#include<ctime>
#include<vector>
#include<cstdio>
#include<string>
#include<bitset>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<functional>
using namespace std;
#define X first
#define Y second
#define mkp make_pair
#define lson (o<<1)
#define rson ((o<<1)|1)
//#define mid (l+(r-l)/2)
#define sz() size()
#define pb(v) push_back(v)
#define all(o) (o).begin(),(o).end()
#define clr(a,v) memset(a,v,sizeof(a))
#define bug(a) cout<<#a<<" = "<<a<<endl
#define rep(i,a,b) for(int i=a;i<(b);i++)
#define scf scanf
#define prf printf
typedef long long LL;
typedef vector<int> VI;
typedef pair<int,int> PII;
typedef vector<pair<int,int> > VPII;
const int INF=0x3f3f3f3f;
const LL INFL=0x3f3f3f3f3f3f3f3fLL;
const double eps=1e-8;
const double PI = acos(-1.0);
//start----------------------------------------------------------------------
const int maxn=1e5+10;
const int maxa=1e6+10;
int gcd(int a,int b) {
return b==0?a:gcd(b,a%b);
}
struct Node {
int l,r,v;
Node(int l,int r,int v):l(l),r(r),v(v) {}
};
bool cmp1(const Node& t1,const Node& t2) {
return t1.r<t2.r;
}
int n,m;
int arr[maxn],ans[maxn];
int dp[maxn][20];
void rmq_init() {
for(int i=1; i<=n; i++) dp[i][0]=arr[i];
for(int j=1; (1<<j)<=n; j++) {
for(int i=1; i+(1<<j)-1<=n; i++) {
dp[i][j]=gcd(dp[i][j-1],dp[i+(1<<(j-1))][j-1]);
}
}
}
int query(int l,int r) {
int k=0;
while(l+(1<<k)-1<=r) k++;
k--;
return gcd(dp[l][k],dp[r-(1<<k)+1][k]);
}
int sumv[maxn],mp[maxa];
void add(int x,int v) {
while(x<maxn) {
sumv[x]+=v;
x+=(-x)&x;
}
}
int sum(int x) {
int ret=0;
while(x>0) {
ret+=sumv[x];
x-=x&(-x);
}
return ret;
}
void init(){
clr(mp,0);
clr(sumv,0);
}
int main() {
while(scf("%d%d",&n,&m)==2&&n) {
init();
for(int i=1; i<=n; i++) scanf("%d",&arr[i]);
rmq_init();
//v表示区间的gcd.
//处理出代表性的事件
vector<Node> events;
for(int i=1; i<=n; i++) {
int r=i;
while(r>0) {
int l=0,v=query(r,i);
events.pb(Node(r,i,v));
while(l+1<r) {
int mid=l+(r-l)/2;
int tmp=query(mid,i);
if(tmp==v) r=mid;
else l=mid;
}
r--;
}
}
//v表示查询的id
vector<Node> que;
for(int i=0; i<m; i++) {
int l,r;
scf("%d%d",&l,&r);
que.pb(Node(l,r,i));
}
sort(all(que),cmp1);
sort(all(events),cmp1);
//离线处理
int j=0;
rep(i,0,que.sz()) {
while(j<events.sz()&&events[j].r<=que[i].r){
Node& e=events[j];
if(mp[e.v]){
add(mp[e.v],-1);
}
add(e.l,1);
mp[e.v]=e.l;
j++;
}
ans[que[i].v]=sum(que[i].r)-sum(que[i].l-1);
}
rep(i,0,m) prf("%d\n",ans[i]);
}
return 0;
}
//end-----------------------------------------------------------------------
HDU 5869 Different GCD Subarray Query rmq+离线+数状数组的更多相关文章
- HDU 5869 Different GCD Subarray Query 离线+树状数组
Different GCD Subarray Query Problem Description This is a simple problem. The teacher gives Bob a ...
- hdu 5869 Different GCD Subarray Query BIT+GCD 2016ICPC 大连网络赛
Different GCD Subarray Query Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/65536 K ( ...
- HDU 5869 Different GCD Subarray Query 树状数组+离线
Problem Description This is a simple problem. The teacher gives Bob a list of problems about GCD (Gr ...
- HDU 5869 Different GCD Subarray Query 树状数组 + 一些数学背景
http://acm.hdu.edu.cn/showproblem.php?pid=5869 题意:给定一个数组,然后给出若干个询问,询问[L, R]中,有多少个子数组的gcd是不同的. 就是[L, ...
- HDU 5869 Different GCD Subarray Query (GCD种类预处理+树状数组维护)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5869 问你l~r之间的连续序列的gcd种类. 首先固定右端点,预处理gcd不同尽量靠右的位置(此时gc ...
- HDU 5869 Different GCD Subarray Query
离线操作,树状数组,$RMQ$. 这个题的本质和$HDU$ $3333$是一样的,$HDU$ $3333$要求计算区间内不同的数字有几个. 这题稍微变了一下,相当于原来扫描到$i$的之后是更新$a[i ...
- HDU 5869 Different GCD Subarray Query(2016大连网络赛 B 树状数组+技巧)
还是想不到,真的觉得难,思路太巧妙 题意:给你一串数和一些区间,对于每个区间求出区间内每段连续值的不同gcd个数(该区间任一点可做起点,此点及之后的点都可做终点) 首先我们可以知道每次添加一个值时gc ...
- 【刷题】HDU 5869 Different GCD Subarray Query
Problem Description This is a simple problem. The teacher gives Bob a list of problems about GCD (Gr ...
- HDU 5869.Different GCD Subarray Query-区间gcd+树状数组 (神奇的标记右移操作) (2016年ICPC大连网络赛)
树状数组... Different GCD Subarray Query Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 65536/6 ...
随机推荐
- MySQL+Service+Servlet+Jsp实现Table表格分页展示数据
下面以一个示例讲解如何使用MySQL+Service+Servlet+Jsp实现Table表格分页展示数据: eg:请假管理系统 要求如下: 一.打开首页页面, 访问查询请假记录的 servlet , ...
- HTML标签速记整理W3C
标题 <h1>段落<p>链接< href="">图像<img src="">自关闭元素,不需要结束标记换行标志& ...
- php添加数据库转义特殊字符串
addslashes()
- BigData:值得了解的十大数据发展趋势
当今,世界无时无刻不在发生着变化.对于技术领域而言,普遍存在的一个巨大变化就是为大数据(Big data)打开了大门,并应用大数据技相关技术来改善各行业的业务并促进经济的发展.目前,大数据的作用已经上 ...
- STM32(5)——通用定时器基本定时器
1.STM32的Timer简介 STM32中一共有11个定时器,其中2个高级控制定时器,4个普通定时器和2个基本定时器,以及2个看门狗定时器和1个系统嘀嗒定时器. 其中系统嘀嗒定时器是前文中所描述的S ...
- python 模块路径查找 及 添加
**** python 模块路径查找: 通过模块的__file__属性来确定: **** 模块路径添加: 方法一:函数添加1 import sys2 查看sys.path3 添加sys.path.ap ...
- 20155220 2016-2017-2《java程序设计》第三周学习总结
教材学习内容总结 1.注意java与c语言的区别,在java中,数组是一个对象. 2.了解java.util.Scanner和java.math.BigDecimal这两个标准类 3.System.a ...
- 【LG2257】YY的GCD
[LG2257]YY的GCD 题面 洛谷 题解 题目大意: 给定\(n,m\)求\(\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)为质数]\). 我们设\(f(x)=[x为 ...
- MySQL入门篇(三)之my.cnf配置文件详解【转】
转自:https://www.cnblogs.com/panwenbin-logs/p/8360703.html #*** client options 相关选项 ***# #以下选项会被MySQL客 ...
- 【redis的搭建】centos6.4下搭建redis
说明:本文内容参考自一些资料,如有雷同,还请见谅. 部分参考: http://blog.csdn.net/su377486/article/details/51803616 http://blog.c ...