【CH 4302】Interval GCD

【题目】

传送门

题目描述：

给定一个长度为 $n$ 的数列 $a$，以及 $m$ 条指令 $(n\le 5\times 10^5,m\le 10^5)$，每条指令可能是以下两种之一：

$“Clrd”$，表示把 $a_l,a_{l+1},\dots ,a_r$ 都加上 $d$。
$“Qlr”$，表示询问 $a_l,a_{l+1},\dots ,a_r$ 的最大公约数 $(gcd)$。

输入格式：

第一行两个整数 $n,m$，第二行 $n$ 个整数 $a_i$，接下来 $m$ 行每条指令的格式如题目描述所示。

输出格式：

对于每个询问，输出一个整数表示答案。

样例数据：

输入
5 5
1 3 5 7 9
Q 1 5
C 1 5 1
Q 1 5
C 3 3 6
Q 2 4

输出
1
2
4

数据范围与约定：

$n,m\le 2\times 10^5$，$l\le r$，数据保证任何时刻序列中的数都是不超过 $2^{62}-1$的正整数。

【分析】

数据范围为什么都不统一啊。。。

做这道题之前，要先知道一个东西，叫做辗转相减法

有一个关于它的推论，即 $gcd(x,y,z)=gcd(x,y-x,z-y)$，推广到 $n$ 个数也成立

这里还是简单证一下（如果有错请提出来）
假设 $gcd(x,y,z)=a$，那么必有互质的三个数 $m_1,m_2,m_3$ 使 $x=m_{1}a,y=m_{2}a,z=m_{3}a$
那么 $gcd(x,y-x,z-y)=gcd(m_{1}a,(m_2-m_1)a,(m_3-m_2)a)$
由于 $m_1,m_2,m_3$ 互质，所以 $m_1,m_2-m_1,m_3-m_2$ 必然也互质，所以 $gcd(x,y-x,z-y)=a=gcd(x,y,z)$

那么，知道了这个之后又有什么用呢

我们可以用差分的思想，把原数组变为 $a_1,a_2-a_1,a_3-a_2,\dots ,a_n-a_{n-1}$（用 $su{m}_i$ 表示 $a_i-a_{i-1}$）

这样的话，修改就好办了，直接在 $su{m}_l$ 加 $d$，在 $su{m}_{r+1}$ 减 $d$（中间的数加的 $d$ 相互抵消了）

查询的时候，根据辗转相减法，先找出 $gcd(su{m}_{l+1},su{m}_{l+2},\dots ,su{m}_r)$，由于此时 $su{m}_l$ 记录的是 $a_l-a_{l-1}$，要单独处理一下，而不难看出 $su{m}_1+su{m}_2+\cdots +su{m}_l=a_l$，所以答案为

$$ans=gcd(\sum _{i=1}^{l}su{m}_i,gcd(su{m}_{l+1},su{m}_{l+2},\dots ,su{m}_r))$$

因此，用一个树状数组记录 $sum$ 的前缀和，用一个线段树记录 $sum$ 的区间 $gcd$ 就可以了

【代码】

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 500005
#define ll long long
#define lowbit(x) x&-x
using namespace std;
int n,m;
ll a[N],bit[N],Gcd[N<<2];
void add(int i,ll x)
{
	while(i<=n)
	{
		bit[i]+=x;
		i+=lowbit(i);
	}
}
ll sum(int i)
{
	ll ans=0;
	while(i)
	{
		ans+=bit[i];
		i-=lowbit(i);
	}
	return ans;
}
void build(int root,int l,int r)
{
	if(l==r)
	{
		Gcd[root]=a[l]-a[l-1];
		return;
	}
	int mid=(l+r)>>1;
	build(root<<1,l,mid);
	build(root<<1|1,mid+1,r);
	Gcd[root]=__gcd(Gcd[root<<1],Gcd[root<<1|1]);
}
void modify(int root,int l,int r,int x,ll k)
{
	if(l==r)
	{
		Gcd[root]+=k;
		return;
	}
	int mid=(l+r)>>1;
	if(x<=mid)  modify(root<<1,l,mid,x,k);
	else  modify(root<<1|1,mid+1,r,x,k);
	Gcd[root]=__gcd(Gcd[root<<1],Gcd[root<<1|1]);
}
ll query(int root,int l,int r,int x,int y)
{
	if(l>=x&&r<=y)
	  return Gcd[root];
	int mid=(l+r)>>1;
	if(y<=mid)  return query(root<<1,l,mid,x,y);
	if(x>mid)  return query(root<<1|1,mid+1,r,x,y);
	return __gcd(query(root<<1,l,mid,x,y),query(root<<1|1,mid+1,r,x,y));
}
int main()
{
	int x,y,i;  ll k;
	scanf("%d%d",&n,&m);
	for(i=1;i<=n;++i)
	  scanf("%lld",&a[i]);
	build(1,1,n);
	for(i=1;i<=m;++i)
	{
		char c;  cin>>c;
		scanf("%d%d",&x,&y);
		if(c=='C')
		{
			scanf("%lld",&k);
			add(x,k),modify(1,1,n,x,k);
			if(y<n)add(y+1,-k),modify(1,1,n,y+1,-k);
		}
		else  printf("%lld\n",abs(__gcd(a[x]+sum(x),query(1,1,n,x+1,y))));
	}
	return 0;
}