HDU 5726 GCD ( 线段树 + map )

Bug原产地

于 2020-09-18 15:21:47 发布

阅读量109

点赞数

分类专栏：数学文章标签： gcd

本文链接：https://blog.csdn.net/qhq889/article/details/108665326

版权

数学专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题目链接
在这里插入图片描述
大致题意：
给你长度为 $n$ 的数列， $m$ 次询问，每次询问给定两个数 $l, r$ ,求区间 $G C D = g c d (a [l], a [l + 1], . . . . ., a [r]),$ 并给出原数列有多少个区间的 $g c d = G C D$

解题报告：
求区间 $g c d$ 有方法：线段树，RMQ，ST表
该题关键是如何预处理求出所有区间的 $g c d$ 呢？？
参考网上博客一般有两种解法：
$1 .$ 枚举区间左端点不断二分右端点（本人表示不太懂参考博客）
$2 .$ 暴力预处理，巧用 map 。

对于数列而言，区间右端点往右走时，它们的 $g c d$ 单调不增，且区间 $g c d$ 的种类数也不多。
跑了 $2000 m s$ 左右
代码展示：

#include<bits/stdc++.h>
#define LL long long
#define pii pair<int,LL>
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
const int maxn=1e5+5;
const int inf=0x3f3f3f3f;
const LL MOD=998244353;
inline int read() {int s = 0, f = 1;char ch = getchar();while (!isdigit(ch)) {f = -1;ch = getchar();}while (isdigit(ch)) {s = (s << 1) + (s << 3) + ch - '0';ch = getchar();}return s * f;}
inline void write(int x) {if (x < 0) x = -x,putchar('-');if (x > 9) write(x / 10);putchar(x % 10 + '0');}

int a[maxn],GCD[maxn<<2];
map<int,LL>mp[maxn];
map<int,LL> ans;
void build(int l,int r,int rt){
    if(l==r){
        GCD[rt]=a[l];return ;
    }
    int mid=(l+r)>>1;
    build(l,mid,rt<<1);
    build(mid+1,r,rt<<1|1);
    GCD[rt]=__gcd(GCD[rt<<1],GCD[rt<<1|1]);
}
int query(int L,int R,int l,int r,int rt){
    if(L<=l&&R>=r) return GCD[rt];
    int mid=(l+r)>>1,ans=0;
    if(L<=mid) ans=__gcd(ans,query(L,R,l,mid,rt<<1));
    if(R>mid) ans=__gcd(ans,query(L,R,mid+1,r,rt<<1|1));
    return ans;
}
int main() {
    int T,n,q,l,r;
    scanf("%d",&T);
    for(int t=1;t<=T;t++){
        ans.clear();
        scanf("%d",&n);
        for(int i=1;i<=n;i++) scanf("%d",&a[i]),mp[i].clear();
        for(int i=1;i<=n;i++){
            ans[a[i]]++;mp[i][a[i]]++;
            for(auto it=mp[i-1].begin();it!=mp[i-1].end();it++){
                int tmp=__gcd(a[i],it->first);
                mp[i][tmp]+=it->second;
                ans[tmp]+=it->second;
            }
        }
        build(1,n,1);
        printf("Case #%d:\n",t);
        scanf("%d",&q);
        for(int i=1;i<=q;i++){
            scanf("%d%d",&l,&r);
            int tmp=query(l,r,1,n,1);
            printf("%d %lld\n",tmp,ans[tmp]);
        }
    }
    return 0;
}