题目链接##

Mychael vs Kid

题解##

先说说这题的由来及前身

前身##

首先有一个很经典的题目:

维护区间加,查询区间\(gcd\)

如果强行用线段树维护的话,区间加之后就没法直接确定当前区间的\(gcd\),不可直接维护

这个时候就用到了\(gcd\)的一个性质:

\[(a,b) = (a - b,b)
\]

三个的\(gcd\)也是符合的:

\[(a,b,c) = (a,b,c - b) = (a,b - a,c - b)
\]

同样可以推广出\(n\)个的情况

\[gcd\{a_i\} = gcd(a_1,gcd\{a_j - a_{j - 1}\}) \qquad i \in [1,n] \quad j \in [2,n]
\]

也就是说,区间差分了之后,我们要求原区间\(gcd[l,r]\),就相当于求\(a_l\)和差分后区间\([l + 1,r]\)的\(gcd\)

所以我们只需维护差分区间和区间每个位置的值即可

维护区间值用线段树一点问题都没有

在加法下维护差分区间,一个区间加上一个数,仅影响区间端点的值,所以单点修改即可

所以我们开两棵线段树,一棵维护权值,一棵维护差值的\(gcd\),即可\(O(nlog^2n)\)解决这道题

本题##

本题多加入了一个区间乘法操作

本来想\(yy\)分块的,可后来发现线段树依旧可写并且暴艹分块

思考一下发现,区间乘法后,区间差分值也乘上那个数,对应\(gcd\)也乘上一个数,可以使用线段树维护

至于端点的差值改变,只需取出端点的值计算出差值后单点加法修改即可

复杂度依旧是\(O(nlog^2n)\)

部分分##

前\(5\%\),咳,,,

对于\(n,m \le 300\),暴力求即可

对于没有操作\(2\)的,就是原版题目

对于\(n,m \le 3 \times 10^4\)的,可以考虑分块

std

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<map>
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define mp(a,b) make_pair<int,int>(a,b)
#define cls(s) memset(s,0,sizeof(s))
#define cp pair<int,int>
#define LL long long int
#define ls (u << 1)
#define rs (u << 1 | 1)
using namespace std;
const int maxn = 100005,maxm = 100005,INF = 1000000000;
inline int read(){
int out = 0,flag = 1; char c = getchar();
while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
return out * flag;
}
LL n,m,A[maxn],D[maxn];
LL gcd(LL a,LL b){
if (a < 0) a = -a;
if (b < 0) b = -b;
return b ? gcd(b,a % b) : a;
}
struct Seg_Gcd{
LL Gcd[maxn << 2],tag[maxn << 2];
void upd(int u){
Gcd[u] = gcd(Gcd[ls],Gcd[rs]);
}
void pd(int u){
if (tag[u] != 1){
Gcd[ls] *= tag[u]; tag[ls] *= tag[u];
Gcd[rs] *= tag[u]; tag[rs] *= tag[u];
tag[u] = 1;
}
}
void build(int u,int l,int r){
tag[u] = 1;
if (l == r){
Gcd[u] = D[l];
return;
}
int mid = l + r >> 1;
build(ls,l,mid);
build(rs,mid + 1,r);
upd(u);
}
void Add(int u,int l,int r,int pos,int v){
if (l == r){Gcd[u] += v; return;}
pd(u);
int mid = l + r >> 1;
if (mid >= pos) Add(ls,l,mid,pos,v);
else Add(rs,mid + 1,r,pos,v);
upd(u);
}
void Mult(int u,int l,int r,int L,int R,int v){
if (l >= L && r <= R){Gcd[u] *= v; tag[u] *= v; return;}
pd(u);
int mid = l + r >> 1;
if (mid >= L) Mult(ls,l,mid,L,R,v);
if (mid < R) Mult(rs,mid + 1,r,L,R,v);
upd(u);
}
LL query(int u,int l,int r,int L,int R){
if (l >= L && r <= R) return Gcd[u];
pd(u);
int mid = l + r >> 1;
if (mid >= R) return query(ls,l,mid,L,R);
if (mid < L) return query(rs,mid + 1,r,L,R);
return gcd(query(ls,l,mid,L,R),query(rs,mid + 1,r,L,R));
}
}T2;
struct Seg{
LL val[maxn << 2],mult[maxn << 2],add[maxn << 2];
void pd(int u){
if (mult[u] != 1){
val[ls] *= mult[u]; mult[ls] *= mult[u]; add[ls] *= mult[u];
val[rs] *= mult[u]; mult[rs] *= mult[u]; add[rs] *= mult[u];
mult[u] = 1;
}
if (add[u]){
val[ls] += add[u]; add[ls] += add[u];
val[rs] += add[u]; add[rs] += add[u];
add[u] = 0;
}
}
void build(int u,int l,int r){
add[u] = 0; mult[u] = 1;
if (l == r){
val[u] = A[l];
return;
}
int mid = l + r >> 1;
build(ls,l,mid);
build(rs,mid + 1,r);
}
void Add(int u,int l,int r,int L,int R,int v){
if (l >= L && r <= R){val[u] += v; add[u] += v; return;}
pd(u);
int mid = l + r >> 1;
if (mid >= L) Add(ls,l,mid,L,R,v);
if (mid < R) Add(rs,mid + 1,r,L,R,v);
}
void Mult(int u,int l,int r,int L,int R,int v){
if (l >= L && r <= R){val[u] *= v; add[u] *= v; mult[u] *= v; return;}
pd(u);
int mid = l + r >> 1;
if (mid >= L) Mult(ls,l,mid,L,R,v);
if (mid < R) Mult(rs,mid + 1,r,L,R,v);
}
LL query(int u,int l,int r,int pos){
if (l == r) return val[u];
pd(u);
int mid = l + r >> 1;
if (mid >= pos) return query(ls,l,mid,pos);
return query(rs,mid + 1,r,pos);
}
}T1;
int main(){
n = read(); m = read();
for (int i = 1; i <= n; i++){
A[i] = read();
D[i] = A[i] - A[i - 1];
}
T1.build(1,1,n);
T2.build(1,1,n);
LL opt,l,r,v,t1,t2;
while (m--){
opt = read(); l = read(); r = read();
if (opt == 1){
v = read();
T1.Add(1,1,n,l,r,v);
T2.Add(1,1,n,l,v);
if (r + 1 <= n) T2.Add(1,1,n,r + 1,-v);
}
else if (opt == 2){
v = read();
t1 = T1.query(1,1,n,l);
if (r + 1 <= n) t2 = T1.query(1,1,n,r);
if (l < r) T2.Mult(1,1,n,l + 1,r,v);
T2.Add(1,1,n,l,t1 * v - t1);
if (r + 1 <= n) T2.Add(1,1,n,r + 1,t2 - t2 * v);
T1.Mult(1,1,n,l,r,v);
}
else {
if (l < r) printf("%lld\n",gcd(T1.query(1,1,n,l),T2.query(1,1,n,l + 1,r)));
else printf("%lld\n",T1.query(1,1,n,l));
}
}
return 0;
}

Myhchael原创题系列 Mychael vs Kid 【题解】的更多相关文章

  1. Newtonsoft.Json C# Json序列化和反序列化工具的使用、类型方法大全 C# 算法题系列(二) 各位相加、整数反转、回文数、罗马数字转整数 C# 算法题系列(一) 两数之和、无重复字符的最长子串 DateTime Tips c#发送邮件,可发送多个附件 MVC图片上传详解

    Newtonsoft.Json C# Json序列化和反序列化工具的使用.类型方法大全   Newtonsoft.Json Newtonsoft.Json 是.Net平台操作Json的工具,他的介绍就 ...

  2. NOI题库 1768最大子矩阵 题解

    NOI题库 1768最大子矩阵  题解     总时间限制: 1000ms 内存限制: 65536kB   描述   已知矩阵的大小定义为矩阵中所有元素的和.给定一个矩阵,你的任务是找到最大的非空(大 ...

  3. hdu5017:补题系列之西安网络赛1011

    补题系列之西安网络赛1011 题目大意:给定一个椭球: 求它到原点的最短距离. 思路: 对于一个椭球的标准方程 x^2/a^2 + y^2/b^2 +z^2/c^2=1 来说,它到原点的最短距离即为m ...

  4. nyoj 1208——水题系列——————【dp】

    水题系列 时间限制:1000 ms  |  内存限制:65535 KB 难度:2   描述     给你一个有向图,每条边都有一定的权值,现在让你从图中的任意一点出发,每次走的边的权值必须必上一次的权 ...

  5. Mychael原创题 洛谷T23923 Mychaelの水题 【题解】

    原题链接 题目大意: 有来自三个地区的人各a,b,c位,他们排成了一排.请问有多少种不同类型的排法,使得相邻的人都来自不同的地区 \(a,b,c<=200\) 答案取模 题解 弱弱的标程解法 设 ...

  6. NOI题库刷题日志 (贪心篇题解)

    这段时间在NOI题库上刷了刷题,来写点心得和题解 一.寻找平面上的极大点 2704:寻找平面上的极大点 总时间限制:  1000ms  内存限制:  65536kB 描述 在一个平面上,如果有两个点( ...

  7. C#版 - PAT乙级(Basic Level)真题 之 1021.个位数统计 - 题解

    版权声明: 本文为博主Bravo Yeung(知乎UserName同名)的原创文章,欲转载请先私信获博主允许,转载时请附上网址 http://blog.csdn.net/lzuacm. C#版 - P ...

  8. LeetCode算法扫题系列19

    原创作品,可以转载,但是请标注出处地址:https://www.cnblogs.com/V1haoge/p/9104677.html LeetCode算法第19题(难度:中等) 题目:给定一个链表,删 ...

  9. LeetCode算法扫题系列83

    原创作品,可以转载,但是请标注出处地址:https://www.cnblogs.com/V1haoge/p/9104582.html LeetCode算法第83题(难度:简单) 题目:给定一个排序链表 ...

随机推荐

  1. WeTest功能优化第3期:业内首创,有声音的云真机

    第3期功能优化目录 [云真机远程调试]音频同步传输实现测试有声 [兼容性测试报告]新增视频助力动态定位问题 [云真机远程调试]菜单栏优化助力机型选择 本期介绍的新功能,秉承创造用户需求的理念,在云真机 ...

  2. java 二叉树的创建 遍历

    本来说复习一下BFS和DFS,辗转就来到了二叉树...本文包括二叉树的创建和遍历 概念 数据:1 2 3 4 5 6 7生成一颗二叉树 上面的数是数据,不是位置,要区别一下数据和位置 红色的代表位置, ...

  3. (python)leetcode刷题笔记 02 Add Two Numbers

    2. Add Two Numbers You are given two non-empty linked lists representing two non-negative integers. ...

  4. 【MFC】学习与问题整合

    需要源码联系邮件:kangxlchn@163.com 1.新建一个MFC工程(基于对话框) 环境:vs2017 统统NEXT 新建完成后打开MFCPrj.cpp文件 打开类试图 每创建一个MFC项目, ...

  5. TensorFlow | ReluGrad input is not finite. Tensor had NaN values

    问题的出现 Question 这个问题是我基于TensorFlow使用CNN训练MNIST数据集的时候遇到的.关键的相关代码是以下这部分: cross_entropy = -tf.reduce_sum ...

  6. UVa 10082 - WERTYU 解题报告 - C语言

    1.题目大意: 输入一个错位的字符串(字母全为大写),输出原本想打出的句子. 2.思路: 如果将每个输入字符所对应的应输出字符一一使用if或者switch,则过于繁琐.因此考虑使用常量数组实现. 3. ...

  7. 局部加权回归(LWR) Matlab模板

    将百度文库上一份局部加权回归的代码,将其改为模板以便复用. q2x,q2y为数据集,是n*1的矩阵: r是波长参数,就是对于距离的惩罚力度: q_x是要拟合的数据横坐标,是1*n的矩阵: 得到的q_y ...

  8. apache 多端口访问 配置

    使用本地ip:端口号,或者修改hosts文件+域名的方法来进行本地多站点web调试.     注意这里是用apache 不是iis 1: 安装好AppServ2.5.9软件,官网是:[url]http ...

  9. ACM 第十天

    动态规划2 1.树形DP 2.概率DP 3.区间DP 模板 ; len < n; len++) { //操作区间的长度 , j = len; j <= n; i++, j++) { //始 ...

  10. lintcode-152-组合

    152-组合 组给出两个整数n和k,返回从1......n中选出的k个数的组合. 样例 例如 n = 4 且 k = 2 返回的解为: [[2,4],[3,4],[2,3],[1,2],[1,3],[ ...