题意

题目链接

Sol

正经做法不会,听lxl讲了一种很神奇的方法

我们考虑如果满足条件,那么需要具备什么条件

设mx为询问区间最大值,mn为询问区间最小值

  1. mx - mn = (r - l) * k

  2. 区间和 = mn * len + \(\frac{n * (n - 1)}{2} k\)

  3. \(\text{立方和} = \sum_{i = 1}^{len} (mn + (i - 1)k) ^2\)

第三条后面的可以直接推式子推出来(\(\sum_{i = 1}^n i^2 = \frac{n(n+1)(2n+1)}{6}\))

最后的/6可以直接乘一下然后ull自然溢出。

#include<bits/stdc++.h>
#define ull unsigned long long
#define LL long long
#define Fin(x) {freopen(#x".in","r",stdin);}
#define Fout(x) {freopen(#x".out","w",stdout);}
using namespace std;
const int MAXN = 4e6 + 10;
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, M;
#define ls k << 1
#define rs k << 1 | 1
struct {
int l, r, mn, mx;
ull s, s2;
}T[MAXN];
void update(int k) {
T[k].mn = min(T[ls].mn, T[rs].mn);
T[k].mx = max(T[ls].mx, T[rs].mx);
T[k].s = T[ls].s + T[rs].s;
T[k].s2 = T[ls].s2 + T[rs].s2;
}
void Build(int k, int ll, int rr) {
T[k].l = ll; T[k].r = rr;
if(ll == rr) {T[k].mn = T[k].mx = T[k].s = read(); T[k].s2 = T[k].s * T[k].s ; return ;}
int mid = ll + rr >> 1;
Build(ls, ll, mid); Build(rs, mid + 1, rr);
update(k);
}
void Modify(int k, int p, int v) {
if(T[k].l == T[k].r) {T[k].mn = T[k].mx = T[k].s = v; T[k].s2 = T[k].s * T[k].s; return ;}
int mid = T[k].l + T[k].r >> 1;
if(p <= mid) Modify(ls, p, v);
if(p > mid) Modify(rs, p, v);
update(k);
}
int QueryMn(int k, int ql, int qr) {
if(ql <= T[k].l && T[k].r <= qr) return T[k].mn;
int mid = T[k].l + T[k].r >> 1;
if(ql > mid) return QueryMn(rs, ql, qr);
else if(qr <= mid) return QueryMn(ls, ql, qr);
else return min(QueryMn(ls, ql, qr), QueryMn(rs, ql, qr));
}
int QueryMx(int k, int ql, int qr) {
if(ql <= T[k].l && T[k].r <= qr) return T[k].mx;
int mid = T[k].l + T[k].r >> 1;
if(ql > mid) return QueryMx(rs, ql, qr);
else if(qr <= mid) return QueryMx(ls, ql, qr);
else return max(QueryMx(ls, ql, qr), QueryMx(rs, ql, qr));
}
ull QuerySum(int k, int ql, int qr) {
if(ql <= T[k].l && T[k].r <= qr) return T[k].s;
int mid = T[k].l + T[k].r >> 1;
if(ql > mid) return QuerySum(rs, ql, qr);
else if(qr <= mid) return QuerySum(ls, ql, qr);
else return QuerySum(ls, ql, qr) + QuerySum(rs, ql, qr);
}
ull QuerySum2(int k, int ql, int qr) {
if(ql <= T[k].l && T[k].r <= qr) return T[k].s2;
int mid = T[k].l + T[k].r >> 1;
if(ql > mid) return QuerySum2(rs, ql, qr);
else if(qr <= mid) return QuerySum2(ls, ql, qr);
else return QuerySum2(ls, ql, qr) + QuerySum2(rs, ql, qr);
}
signed main() {
N = read(); M = read();
Build(1, 1, N);
int GG = 0;
while(M--) {
int opt = read();
if(opt == 1) {
int x = read() ^ GG, y = read() ^ GG;
Modify(1, x, y);
} else {
int l = read() ^ GG, r = read() ^ GG; ull k = read() ^ GG;
ull n = r - l + 1, mn = QueryMn(1, l, r), mx = QueryMx(1, l, r);
ull s = QuerySum(1, l, r), s2 = QuerySum2(1, l, r), ns = mn * n + n * (n - 1) * k / 2;
ull gg = (6 * mn * mn * n) + (6 * mn * k * n * (n - 1)) + (k * k * (n - 1) * n * (2 * (n - 1) + 1)),
gg2 = 6 * s2;
if((mx - mn == (r - l) * k) &&
(mn * n + n * (n - 1) / 2 * k == s) &&
(gg == gg2)) puts("Yes"), GG++;
else puts("No");
}
}
return 0;
}
/*
5 3
1 3 2 5 6
2 2 2 23333
1 5 4
2 1 5 1
*/

BZOJ4373: 算术天才⑨与等差数列(线段树 hash?)的更多相关文章

  1. [BZOJ4373]算术天才⑨与等差数列(线段树)

    [l,r]中所有数排序后能构成公差为k的等差数列,当且仅当: 1.区间中最大数-最小数=k*(r-l) 2.k能整除区间中任意两个相邻数之差,即k | gcd(a[l+1]-a[l],a[l+2]-a ...

  2. bzoj4373 算术天才⑨与等差数列——线段树+set

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4373 一个区间有以 k 为公差的数列,有3个条件: 1.区间 mx - mn = (r-l) ...

  3. 【BZOJ4373】算术天才⑨与等差数列 线段树+set

    [BZOJ4373]算术天才⑨与等差数列 Description 算术天才⑨非常喜欢和等差数列玩耍.有一天,他给了你一个长度为n的序列,其中第i个数为a[i].他想考考你,每次他会给出询问l,r,k, ...

  4. 【BZOJ4373】算术天才⑨与等差数列 [线段树]

    算术天才⑨与等差数列 Time Limit: 10 Sec  Memory Limit: 128 MB[Submit][Status][Discuss] Description 算术天才⑨非常喜欢和等 ...

  5. BZOJ 4373 算术天才⑨与等差数列 线段树+set(恶心死我了)

    mdzz,这道题重构了4遍,花了一个晚上... 满足等差数列的条件: 1. 假设min是区间最小值,max是区间最大值,那么 max-min+k(r−l) 2. 区间相邻两个数之差的绝对值的gcd=k ...

  6. BZOJ 4373算术天才⑨与等差数列(线段树)

    题意:给你一个长度为n的序列,有m个操作,写一个程序支持以下两个操作: 1. 修改一个值 2. 给出三个数l,r,k, 询问:如果把区间[l,r]的数从小到大排序,能否形成公差为k的等差数列. n,m ...

  7. BZOJ 4373: 算术天才⑨与等差数列 线段树

    Description 算术天才⑨非常喜欢和等差数列玩耍. 有一天,他给了你一个长度为n的序列,其中第i个数为a[i]. 他想考考你,每次他会给出询问l,r,k,问区间[l,r]内的数从小到大排序后能 ...

  8. bzoj 4373 算术天才⑨与等差数列——线段树+set

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4373 能形成公差为k的等差数列的条件:mx-mn=k*(r-l) && 差分 ...

  9. BZOJ4373 算术天才⑨与等差数列 【线段树】*

    BZOJ4373 算术天才⑨与等差数列 Description 算术天才⑨非常喜欢和等差数列玩耍. 有一天,他给了你一个长度为n的序列,其中第i个数为a[i]. 他想考考你,每次他会给出询问l,r,k ...

随机推荐

  1. js获取n分钟(或n小时或n个月)后(或前)的时间(日期)

    标题有点绕,其实意思就是根据系统当前时间,获取n分钟或n小时或n个月后的时间. 例如:当前时间下,获取10分钟后的时间. var date=new Date(); //1. js获取当前时间 var ...

  2. 整理版jq 复习贴子

    1绝对定位(abs)与相对定位(relative) 区别是相对定位参照自己的位置进行移动(当然需要设置top left这些生效)并且原来的位置保留着 偏移后会把 其它的层遮罩住 绝对定位就是的参照位置 ...

  3. CENTOS安装ElasticSearch(转)

    From: https://my.oschina.net/topeagle/blog/591451?fromerr=mzOr2qzZ CENTOS安装ElasticSearch ElasticSear ...

  4. Ruby:Mechanize的使用教程

    小技巧 puts Mechanize::AGENT_ALIASES 可以打印出所有可用的user_agent puts Mechanize.instance_methods(false) 输出Mech ...

  5. 全网最详细的Sublime Text 3的插件官方网站(图文详解)

    不多说,直接上干货! 全网最详细的Windows里下载与安装Sublime Text *(图文详解) 全网最详细的Sublime Text 3的激活(图文详解) 全网最详细的Sublime Text ...

  6. MySql 高可用架构Atlas

    Atlas是由 Qihoo 360公司Web平台部基础架构团队开发维护的一个基于MySQL协议的数据中间层项目.它在MySQL官方推出的MySQL-Proxy 0.8.2版本的基础上,修改了大量bug ...

  7. KVM:日常管理常用命令

    1.查看.编辑及备份KVM 虚拟机配置文件 以及查看KVM 状态: 1.1.KVM 虚拟机默认的配置文件在 /etc/libvirt/qemu 目录下,默认是以虚拟机名称命名的.xml 文件,如下,: ...

  8. html中img图片居中

    直接看代码 style="display:block; margin:0 auto;" 可以看到,蓝色的代码是让 块状元素居中 而红色的代码,是让img转成块状元素 ···原来im ...

  9. MYSQLI_USE_RESULT or MYSQLI_STORE_RESULT

    之前都是使用同事封装好的mysql类,今天做性能测试时自己手动编写了查询mysql的操作.偶然发现mysqli::query(或者mysqli_query)有一个参数$resultmode取值为MYS ...

  10. setup&hold

    setup time:建立时间,也就是在时钟上升沿到来前,数据需要稳定的时间.hold time:保持时间,指的是在时钟上升沿到来后,数据还需要保持的时间.实际上设置setup time和hold t ...