6月18日练习赛总结

最新推荐文章于 2024-05-21 23:16:53 发布

xiwang233

最新推荐文章于 2024-05-21 23:16:53 发布

阅读量44

点赞数

分类专栏：总结文章标签：算法 c++ 动态规划

本文链接：https://blog.csdn.net/xiwang233/article/details/131276288

版权

总结专栏收录该内容

4 篇文章 0 订阅

订阅专栏

T1 [模拟赛20221101] 子集

偶数很好想，奇数先按偶数的搞了中间留三列找规律（乱搞）

T2 最大公约数

gcd的递减是log的段数，所以可以二分每一段统计答案，st表比线段树要快很多

也可以不用数据结构，分成log段，直接从左右端点往里统计答案

#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=1e5+5;
int n,q,a[N],l,r;
map<int,int>cnt;
inline int gcd(int a,int b){
	if(b==0){
		return a;
	}else{
		return gcd(b,a%b);
	}
}
int _log[N],f[N][20];
void init(){
	for(int i=1;i<=n;i++){
		_log[i]=log2(i);
	}
	for(int i=1;i<=n;i++){
		f[i][0]=a[i];
	}
	for(int j=1;j<=_log[n];j++){
		for(int i=1;i<=n-(1<<j)+1;i++){
			f[i][j]=gcd(f[i][j-1],f[i+(1<<(j-1))][j-1]);
		}
	}
}
int ask(int x,int y){
	int k=_log[y-x+1];
	return gcd(f[x][k],f[y-(1<<k)+1][k]);
}
int find(int fr,int l,int r,int x){
	while(l<r){
		int mid=(l+r+1)>>1;
		if(ask(fr,mid)==x){
			l=mid;
		}else{
			r=mid-1;
		}
	}
	return l;
}
void solve(){
	for(int i=1;i<=n;i++){
		int j=i;
		while(j<=n){
			int jj=find(i,j,n,ask(i,j));
			cnt[ask(i,j)]+=jj-j+1;
			j=jj+1;
		}
	}
}
signed main(){
	scanf("%lld%lld",&n,&q);
	for(int i=1;i<=n;i++){
		scanf("%lld",&a[i]);
	}
	init();
	solve();
	int l,r;
	for(int i=1;i<=q;i++){
		scanf("%lld%lld",&l,&r);
		int ans=ask(l,r);
		printf("%lld %lld\n",ans,cnt[ans]);
	}
	return 0;
}