给定一个长度为 N 的数列 A,以及 M 条指令,每条指令可能是以下两种之一:
C l r d
,表示把 A[l],A[l+1],…,A[r] 都加上 d。Q l r
,表示询问 A[l],A[l+1],…,A[r] 的最大公约数(GCD)。
对于每个询问,输出一个整数表示答案。
输入格式
第一行两个整数 N,M。
第二行 N 个整数 A[i]。
接下来 M 行表示 M 条指令,每条指令的格式如题目描述所示。
输出格式
对于每个询问,输出一个整数表示答案。
每个答案占一行。
数据范围
N≤500000,M≤100000
1≤A[i]≤1018
|d|≤1018
保证数据在计算过程中不会超过 long long 范围。
输入样例:
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
解析:
gcd部分:
一般而言要证明 ,要证明 a=b,可以转化为证明 a>=b,b>=a.
证明 gcd(a1,a2,...,an)<=gcd (a1,a2-a1,a3-a2,...,an-an-1);
令 d = gcd(a1,a2,…,an),
则 a1%d=0,a2%d=0,…,an%d=0
则 a1%d=0,(a2-a1)%d=0,(a3-a2)%d=0,…,(an-an-1)%d=0
所以 gcd(a1,a2,...,an)<=gcd(a1,a2-a1,a3-a2,...,an-an-1);(子集包含关系)
类似的,可以证明>= 关系。
详细解释参考蓝书
#include<iostream>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<utility>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<math.h>
#include<map>
#include<sstream>
#include<deque>
#include<unordered_map>
#include<unordered_set>
#include<bitset>
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int, int> PII;
typedef pair<double, double> PDD;
const int INF = 0x3f3f3f3f;
const LL Mod = 1e9;
const int N = 5e5 + 10, M = 250 + 10, P = 110;
int n, m;
LL w[N];
struct Node {
int l, r;
LL sum, d;
}tr[N*4];
LL gcd(LL a, LL b) {
if (b == 0)return a;
return gcd(b, a % b);
}
void pushup(Node& u, Node& l, Node& r) {
u.sum = l.sum + r.sum;
u.d = gcd(l.d, r.d);
}
void pushup(int u) {
pushup(tr[u], tr[u << 1], tr[u << 1 | 1]);
}
void build(int u, int l, int r) {
if (l == r) {
LL b = w[l] - w[l - 1];
tr[u] = { l,r,b,b};
}
else {
tr[u].l = l, tr[u].r = r;
int mid = l + r >> 1;
build(u<<1, l, mid), build(u<<1|1, mid + 1, r);
pushup(u);
}
}
Node query(int u, int l, int r) {
if (tr[u].l >= l && tr[u].r <= r) {
return tr[u];
}
else {
int mid = tr[u].l + tr[u].r >> 1;
if (r <= mid)return query(u << 1, l, r);
else if (l > mid)return query(u << 1 | 1, l, r);
else {
Node ret;
Node left = query(u << 1, l, r);
Node right = query(u << 1 | 1, l, r);
pushup(ret, left, right);
return ret;
}
}
}
void modify(int u, int x, LL c) {
if (tr[u].l == x && tr[u].r == x) {
LL d = tr[u].sum + c;
tr[u] = { x,x,d,d };
}
else {
int mid = tr[u].l + tr[u].r >> 1;
if (x <= mid)modify(u << 1, x, c);
else modify(u << 1 | 1, x, c);
pushup(u);
}
}
int main() {
cin >> n >> m;
for (int i = 1; i <= n; i++) {
scanf("%lld", &w[i]);
}
build(1, 1, n);
char op[2];
int l, r;
LL d;
while (m--) {
scanf("%s%d%d", &op, &l, &r);
if (op[0] == 'Q') {
auto left = query(1, 1, l);
//printf("+____________%d\n", left.sum);
Node right({ 0, 0, 0, 0 });
if (l + 1 <= r) right = query(1, l + 1, r);
printf("%lld\n", abs(gcd(left.sum, right.d)));
}
else {
scanf("%lld", &d);
modify(1, l, d);
if(r+1<=n)modify(1, r + 1, -d);
}
}
return 0;
}