HDU - 4630 No Pain No Game 线段树离线操作 求区间内任意两个数的最大公约数

HDU - 4630 No Pain No Game

题意：

有n个数，m个区间询问，每次输出区间内任意两个数的最大公约数。

思路：

线段树离线+点更新+区间查询

我们维护一颗线段树，线段树维护的值是区间内最大的gcd。首先我们预处理所有数字的因子，然后我们将排序过后的区间（按照右端点从小到大排序）从左到右依次处理（为了加强理解，设定k<i<j），每扫描一个值a[i],我们查看a[i]的所有因子，对于某一个因子x，如果之前有某个值a[k]也存在这样的因子x，就更新上一次 x 出现的位置（k）。接着我们考虑数组pre[x]代表x这个因子上一次（最近一次）出现的位置（即某个数值包含因子x），如果没有出现过就标记为-1，如果之前处理了所有的pre[x]，那么我们枚举每个数值的因子，就可以根据pre数组判断能否更新pre[x]对应的点（pre[x]的初始值为 -1，因子 x 第一次出现时改变pre[x]的值，但不更新。因为只有出现两次或两次以上 x 才为两个数或者更多数的公约数）。那么对于按右端点排序好的查询，如果在不断update的过程中遇到了查询的右端点，那么我们就可以做查询即可。

4
8 2 4 9
3
2 3
1 3
2 4
i = 1时， 线段树的叶子节点[0, 0, 0, 0];
i = 2时， 线段树的叶子节点[2, 0, 0, 0];
i = 3时， 线段树的叶子节点[4, 2, 0, 0]; 遇到两个询问， 输出区间的最大值 Query(2, 3)  = 2  Query(1, 3) = 4;
i = 4时， 线段树的叶子节点[4, 2, 1, 0]; 遇到一个询问， Query(2, 4) = 2;

 

#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<sstream>
#include<cstring>
#include<bitset>
#include<cstdio>
#include<string>
#include<deque>
#include<stack>
#include<cmath>
#include<queue>
#include<set>
#include<map>
#define mod 1000000007
#define ls rt<<1
#define rs rt<<1|1
using namespace std;
typedef long long ll;
const int maxn = 5e4+2;
int a[maxn],sum[maxn<<2],pre[maxn],res[maxn];
vector<int>v[maxn];
struct node
{
    int l,r;
    int id;
}q[maxn];
void init()
{
    for(int i=1;i<maxn;i++)
    {
        for(int j=i;j<maxn;j+=i)
            v[j].push_back(i);
    }
}
bool cmp(node a,node b)
{
    return a.r<b.r;
}
void pushup(int rt)
{
    sum[rt] = max(sum[ls],sum[rs]);
}
void update(int x,int val,int l,int r,int rt)  //点更新（更新x点最值）
{

    if(l==r)
    {
        sum[rt] = max(sum[rt],val);
        return ;
    }
    int mid = (l+r)>>1;
    if(x<=mid)
        update(x,val,l,mid,ls);
    else
        update(x,val,mid+1,r,rs);
    pushup(rt);
}
int query(int L,int R,int l,int r,int rt)  //  叶子结点存储的即为最大公约数
{
    if(L<=l && r<=R)
        return sum[rt];
    int mid = (l+r)>>1;
    int ans = 0;
    if(L<=mid)
        ans = max(ans,query(L,R,l,mid,ls));
    if(R>mid)
        ans = max(ans,query(L,R,mid+1,r,rs));
    return ans;
}
int main()
{
    int t;
    int n,m;
    init();
    scanf("%d",&t);
    while(t--)
    {
        memset(sum,0,sizeof(sum));
        memset(pre,-1,sizeof(pre));
        scanf("%d",&n);
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]);
        scanf("%d",&m);
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&q[i].l,&q[i].r);
            q[i].id = i;
        }
        sort(q,q+m,cmp);
        for(int i=1,j=0;i<=n && j<m;i++)
        {
            for(int k=0;k<v[a[i]].size();k++)
            {
                int x = v[a[i]][k];
                if(pre[x]!=-1)
                    update(pre[x],x,1,n,1);  //更新上一次出现的位置

                pre[x] = i;
            }
            while(j<m && q[j].r==i)
            {
                res[q[j].id] = query(q[j].l,q[j].r,1,n,1);
                j++;
            }
        }
        for(int i=0;i<m;i++)
            printf("%d\n",res[i]);
    }
    return 0;
}

参考博客：https://blog.csdn.net/u010033217/article/details/38156507