本文介绍了一种高效算法，用于解决长度为n的序列中进行mmm次操作的问题，包括对[l,r]区间加ddd的操作和查询[l,r]区间的最大公约数。通过维护原数组和差分数组，利用递推公式简化查询过程，提供C++代码实现。
题目链接：点击查看
题目大意：给出一个长度为 n n n 的序列，需要执行 m m m 次操作，每次操作分为下列两种类型：
- C l r d：将区间 [ l , r ] [l,r] [l,r] 位置的数字都加上 d d d
- Q l r：询问区间 [ l , r ] [l,r] [l,r] 的最大公约数
题目分析：首先是推广一下最大公约数的公式：
- g c d ( a , b ) = g c d ( a , b − a ) gcd(a,b)=gcd(a,b−a) gcd(a,b)=gcd(a,b−a)
- g c d ( a , b , c ) = g c d ( a , b − a , c − b ) gcd(a,b,c)=gcd(a,b-a,c-b) gcd(a,b,c)=gcd(a,b−a,c−b)
- g c d ( a 1 , a 2 , . . . a n ) = g c d ( a 1 , a 2 − a 1 , . . . , a n − a n − 1 ) gcd(a_1,a_2,...a_n)=gcd(a_1,a_2-a_1,...,a_n-a_{n-1}) gcd(a1,a2,...an)=gcd(a1,a2−a1,...,an−an−1)
如此一来我们可以同时维护一下原数组以及差分数组，设 b i = a i − a i − 1 b_i=a_i-a_{i-1} bi=ai−ai−1，对于每次修改而言，差分数组是单点修改，而原数组是区间修改。对于每次查询而言，差分数组是区间查询，原数组是单点查询。
最后说一下改如何查询，对于 [ l , r ] [l,r] [l,r] 的查询，我们只需要输出 g c d ( a [ l ] , b [ l + 1 ] , . . . , b [ r ] ) gcd(a[l],b[l+1],...,b[r]) gcd(a[l],b[l+1],...,b[r]) 就是答案了
代码：
// Problem: 区间最大公约数
// Contest: AcWing
// URL: https://www.acwing.com/problem/content/247/
// Memory Limit: 64 MB
// Time Limit: 1000 ms
//
// Powered by CP Editor (https://cpeditor.org)
// #pragma GCC optimize(2)
// #pragma GCC optimize("Ofast","inline","-ffast-math")
// #pragma GCC target("avx,sse2,sse3,sse4,mmx")
#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
#include<cassert>
#include<bitset>
#define lowbit(x) x&-x
using namespace std;
typedef long long LL;
typedef unsigned long long ull;
template<typename T>
inline void read(T &x)
{
T f=1;x=0;
char ch=getchar();
while(0==isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
x*=f;
}
template<typename T>
inline void write(T x)
{
if(x<0){x=~(x-1);putchar('-');}
if(x>9)write(x/10);
putchar(x%