辗转相减法的扩展 $gcd(x, y, z) = gcd(x, y - x, z - y)$ 当有n个数时也成立

所以构造$a_{i}$的差分数组$b_{i} = a_{i} - a_{i - 1}$,用一个线段树来维护b数组的gcd,这样每次区间修改相当于两次单点修改

考虑到询问的时候$ans = gcd(a_{l}, query(l +1, r))$所以我们再维护原数组a的值,直接差分之后用一个树状数组就好了

注意判断边界情况。

Code:

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll; const int N = 5e5 + ; int n, qn;
ll a[N], b[N]; ll gcd(ll x, ll y) {
return (!y) ? x : gcd(y, x % y);
} template <typename T>
inline void read(T &X) {
X = ;
char ch = ;
T op = ;
for(; ch > '' || ch < ''; ch = getchar())
if(ch == '-') op = -;
for(; ch >= '' && ch <= ''; ch = getchar())
X = (X << ) + (X << ) + ch - ;
X *= op;
} struct BinaryIndexTree {
ll s[N]; #define lowbit(x) ((x) & (-x)) inline void add(int x, ll v) {
for(; x <= n; x += lowbit(x))
s[x] += v;
} inline ll query(int x) {
ll res = ;
for(; x > ; x -= lowbit(x))
res += s[x];
return res;
} } B; struct SegT {
ll s[N << ]; #define lc p << 1
#define rc p << 1 | 1
#define mid ((l + r) >> 1) inline void up(int p) {
if(p) s[p] = gcd(s[lc], s[rc]);
} void build(int p, int l, int r) {
if(l == r) {
s[p] = b[l];
return;
} build(lc, l, mid);
build(rc, mid + , r);
up(p);
} void modify(int p, int l, int r, int x, ll v) {
if(x == l && x == r) {
s[p] += v;
return;
} if(x <= mid) modify(lc, l, mid, x, v);
else modify(rc, mid + , r, x, v);
up(p);
} ll query(int p, int l, int r, int x, int y) {
if(x > y) return 1LL;
if(x <= l && y >= r) return s[p];
ll res;
if(x <= mid && y > mid)
res = gcd(query(lc, l, mid, x, y), query(rc, mid + , r, x, y));
else if(y <= mid) res = query(lc, l, mid, x, y);
else if(x > mid) res = query(rc, mid + , r, x, y);
return res;
} } A; inline ll abs(ll x) {
return x > ? x : -x;
} int main() {
read(n), read(qn);
for(int i = ; i <= n; i++) read(a[i]); for(int i = ; i <= n; i++) b[i] = a[i] - a[i - ];
for(int i = ; i <= n; i++) B.add(i, b[i]);
A.build(, , n); char op[];
for(int x, y; qn--; ) {
scanf("%s", op);
read(x), read(y);
if(op[] == 'C') {
ll v;
read(v);
B.add(x, v), A.modify(, , n, x, v);
if(y < n) B.add(y + , -v), A.modify(, , n, y + , -v);
} else {
if(x < y) printf("%lld\n", gcd(B.query(x), abs(A.query(, , n, x + , y))));
else printf("%lld\n", B.query(x));
}
} return ;
}

最新文章

  1. Tesseract-OCR字符识别简介
  2. jdk环境变量的配置并检测是否配置成功
  3. CF# 334 Alternative Thinking
  4. 44. Decode Ways &amp;&amp; Gray Code
  5. vs2013 error c4996: &#39;fopen&#39;: This function or varia
  6. DOCTYPE、指定语言、字符集
  7. Visual Assist X破解方法
  8. pdf文件流生成pdf文件
  9. html5标签兼容ie6,7,8
  10. 全面理解Unity加载和内存管理
  11. 【HDOJ】1494 跑跑卡丁车
  12. 快速搭建vsftp 服务器并配置指定目录
  13. Linux本地网络脚本配置(内网与外网)
  14. 学习Spring Boot:(一)入门
  15. C#分解质因数
  16. 梅尔罗斯百度云在线观看迅雷下载Patrick Melrose磁力BT下载
  17. 解决python2.x文件读写编码问题
  18. Unity学习笔记草稿篇(一)为unity配置添加VS智能感知
  19. 爬虫之动态HTML处理(Selenium与PhantomJS )动态页面模拟点击
  20. vue的实例

热门文章

  1. winform中的状态栏,以及在状态栏目上显示时间
  2. 关于 freetds pymssql 的安装部署
  3. Mybatis_generator自动化生成代码
  4. BeetleX高性能通讯开源组件
  5. 从内存中直接运行PE程序
  6. CF 622F The Sum of the k-th Powers——拉格朗日插值
  7. HDU4135(容斥原理)
  8. JAX-WS注解
  9. 机器学习:逻辑回归(OvR 与 OvO)
  10. C# 保存文件如有重名在原名后加(*)