HDOJ 5726 GCD 二分+ST表

原创 2016年08月28日 18:33:12

HDOJ 5726

题意:给定区间【L,R】求出区间【L,R】内数组中的数的GCD,并且找出所有的数对L0,R0,1<=L0<R0<=n,使得区间【L0,R0】内的GCD等于前一个GCD值


看到这个题,跟codeforces 689 D似曾相识啊!所以把两个题都一起补了

题解统一写在这儿吧

首先说做法:

ST表:用O(nlogn)的离线处理好任意区间的内的GCD值,然后做到O(1)的查询

ST表的搞法是典型的二分想法

拿最小值举例,我们要求某一段区间的最小值

可以把区间二分,先求左边区间的最小值,右边区间的最小值

那么整段区间的最小值就是上述两个值的较小值

因为GCD,MAX,MIN,LCM等都有这个性质,那么我们都可以用ST表来这样进行处理

然后呢,枚举左端点,用二分找到右端点,然后判断GCD值是否相等

然后把值统计到GCD的hash之中


为什么可以这么样做?!

因为:GCD,MAX,MIN,LCM等都是不增函数


也就是说,我们可以利用这个性质,使用ST表和二分来避免时间和空间上浪费

这个题呢,因为g值并不多,所以可以提前都给预处理好,建立一个mp的映射

codeforces 689D呢

因为不好预处理,所以涉及到两次二分的问题

第一次二分,找到右端点的左值(如果有的话)

如果有,那么进行第二次二分,找到右端点的右值,那么对应该左端点,答案的增长是右端点的区间长度

如果没有,那么说明对应该左端点,没有符合条件的右端点


代码:

#include<bits/stdc++.h>
using namespace std;

const int maxn=1e5+50; 
int a[maxn][20],T,n,m;
map<int,long long> mp;

int gcd(int x,int y){
	return y==0?x:gcd(y,x%y);
}

void init(){
	for(int j=1;(1<<j)<=n;j++)
		for(int i=1;i+(1<<j)-1<=n;i++)
			a[i][j]=gcd(a[i][j-1],a[i+(1<<(j-1))][j-1]);
}

int GETGCD(int L,int R){
	int K=(int)log2((double)(1.0*(R-L+1)));
	return gcd(a[L][K],a[R-(1<<K)+1][K]);
}

void getans(){
	mp.clear();
	int l,r,g,i,j,mid;
	for(i=1;i<=n;i++){
		g=a[i][0];
		j=i;
		while(j<=n){
			l=j,r=n;
			while(l<r){
				mid=(l+r+1)>>1;
				if (GETGCD(i,mid)==g) l=mid;
				else r=mid-1;
			}
			mp[g]+=(l-j+1);
			j=l+1;
			g=GETGCD(i,j);
		}
	}
}

int main(){
	//freopen("input.txt","r",stdin);
	scanf("%d",&T);
	for(int Case=1;Case<=T;Case++){
		scanf("%d",&n);
		for(int i=1;i<=n;i++) scanf("%d",&a[i][0]);
		init();
		getans();
		printf("Case #%d:\n",Case);
		scanf("%d",&m);
		int l,r,g;
		while(m--){
			scanf("%d%d",&l,&r);
			g=GETGCD(l,r);
			printf("%d %I64d\n",g,mp[g]);
		}
	}
	return 0;
}


相关文章推荐

ST表 HDOJ5726 区间GCD

题目传送门:http://acm.hdu.edu.cn/showproblem.php?pid=5726ST表适用于静态数据的查询(不支持修改),经过一次O(nlogn)的离线预处理之后,查询O(1)...

hdu5726 GCD st表 + 二分

题目描述:给出n( 1 思路:这个问题的突破口在于从第i个数开始,连续向后取gcd,其gcd的值是逐渐变小的,且可能出现的gcd的值只有不超过log(ai)个。            那...
  • jijijix
  • jijijix
  • 2017年02月06日 11:07
  • 88

RMQ -ST ——hdu5726 GCD

题目:http://acm.hdu.edu.cn/showproblem.php?pid=5726AC代码:#include #include #include #include #inclu...

[2016ACM多校] HDU5726 st表 离线

题意给n个数{a[i]},q个询问(l,r),与[l,r]最大公约数相同的区间有多少个。思路首先是查询最大公约数,然后统计相同的区间。查询可以用各种结构,既然只用查询那就用st简单高效,主要考虑统计。...
  • SongBn
  • SongBn
  • 2016年07月19日 22:31
  • 224

HDU 5726 GCD(RMQ+二分,详解,ST算法)

题目链接:HDU 5726 题面: GCD Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Ja...

HDU 5726 求gcd=k的区间的个数 (二分+RMQ)

题意:给一个数组a,大小为n,接下来有m个询问,每次询问给出l、r,定义f[l,r]=gcd(al,al+1,...,ar),问f[l,r]的值 和 有多少对(l',r')使得f[l',r']=f[l...

HDU5726 GCD 二分查找加RMQ 多校联赛第一场

GCD Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Sub...

RMQ+二分——GCD ( HDU 5726 ) ( 2016 Multi-University Training Contest 1 1004 )

题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5726 分析: 给出一串数,然后进行区间查询,每次查询一个区间的gcd,并且要求同时输出和这个区间...
  • FeBr2
  • FeBr2
  • 2016年07月25日 17:28
  • 229

HDU 5726 GCD(RMQ+二分)(线段树也可)

hdu 5726 GCD

hdu 5726 GCD RMQ+二分枚举

GCD Time Limit: 10000/5000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Submi...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:HDOJ 5726 GCD 二分+ST表
举报原因:
原因补充:

(最多只允许输入30个字)