LOJ 6029 市场--线段树区间更新之分块除法

最新推荐文章于 2021-10-21 15:39:56 发布
最新推荐文章于 2021-10-21 15:39:56 发布
阅读量486 收藏 1
点赞数
分类专栏： 树 ————数据结构————— 分块 文章标签： ACM
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/zjyang12345/article/details/90299871
版权

题目链接

参考链接

一、题意：

四种区间操作：加、除（向下取整）、询问区间和、询问最小值

 

二、思路：

第一想到的是将除d变成乘1/d，但是由于不一定是整除，所以除法和加法运算顺序不同，会影响结果。比如 3和3，先除2后加结果是2，先加后除2答案是3

第二既然区间更新不好做，那可以暴力单点更新？显然这种第51号元素的做法妥妥的TLE~

所以我们要想怎么保持区间更新，然后尽量把除法变为减法，

首先有个玄学玩意叫除法分块，对于n/i而言，总存在连续的区间值是相等的（更多应用链接）于是我们可以在所求区间基础上，去维护这个变化值相同（即max和min变化一样）的小区间，然后用区间减法维护（式子：a=a-(a-a/d)）

简单的连续区间更新见poj1540

 

三、前置技能

我们先看一下乘法和加法共存怎么维护，由于可以用结合律拆分（a+b）*c=a*c+b*c，处理乘法时，同时将乘加两种懒惰标记乘c即可

  模板P3373

HDU 4578 加上维护p方

#include<iostream>
#include<cstdio>
#include<cstring>
#include<climits>
#include<iomanip>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
#define ll long long

inline int read(){
    int x=0; bool f=1; char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(; isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+(c^'0');
    if(f) return x;
    return 0-x;
}
int n,m;
ll mod;
ll a[100005],sum[400005],min1[400005],max1[400005],mul[400005],laz[400005];
void pushup(int i){sum[i]=(sum[(i<<1)]+sum[(i<<1)|1]);min1[i]=min(min1[(i<<1)],min1[(i<<1)|1]);max1[i]=max(max1[(i<<1)],max1[(i<<1)|1]);}
void pushdown(int i,int s,int t){
    int l=(i<<1),r=(i<<1)|1,mid=(s+t)>>1;
    if(mul[i]!=1){
    mul[l]=mul[i]*mul[l]%mod;mul[r]=mul[i]*mul[r]%mod;
    laz[l]=mul[i]*laz[l]%mod;laz[r]=laz[r]*mul[i]%mod;
    sum[l]=sum[l]*mul[i]%mod;sum[r]=sum[r]*mul[i]%mod;
    mul[i]=1;
    }
    if(laz[i]){
    sum[l]+=laz[i]*(mid-s+1);
    sum[l]%=mod;
    sum[r]+=laz[i]*(t-mid);
    sum[r]%=mod;
//    min1[l]+=laz[i];
//    min1[r]+=laz[i];
//    
//    max1[l]+=laz[i];
//    max1[r]+=laz[i];
    laz[l]+=laz[i];
    laz[l]%=mod;
    laz[r]+=laz[i];
    laz[r]%=mod;
    laz[i]=0;
    }
}
void build(int l,int r,int i){
     mul[i]=1;
    if(l==r){sum[i]=min1[i]=max1[i]=a[l];return;}
    int mid=(l+r)>>1;
    build(l,mid,i<<1);build(mid+1,r,(i<<1)|1);pushup(i);
}
void cheng(int l,int r,int s,int t,int i,ll z){
    int mid=(s+t)>>1;
    if(l<=s&&t<=r){
    	//简单的乘加维护 
    	mul[i]=(mul[i]*z)%mod;laz[i]=(laz[i]*z)%mod;
        sum[i]=sum[i]*z%mod;return;
    }
    pushdown(i,s,t);
    if(mid>=l)cheng(l,r,s,mid,(i<<1),z);
    if(mid+1<=r)cheng(l,r,mid+1,t,(i<<1)|1,z);
    pushup(i);
}
void add(int l,int r,int s,int t,int i,ll z){
    int mid=(s+t)>>1;
    if(l<=s&&t<=r){
        sum[i]+=z*(t-s+1)%mod;laz[i]=(laz[i]+z)%mod;
        max1[i]+=z,max1[i]%=mod,min1[i]+=z,min1[i]%=mod;
        return;
    }
    pushdown(i,s,t);
    if(mid>=l)add(l,r,s,mid,(i<<1),z);
    if(mid+1<=r)add(l,r,mid+1,t,(i<<1)|1,z);
    pushup(i);
}
ll query(int l,int r,int s,int t,int i){
    int mid=(s+t)>>1;ll tot=0;
    if(l<=s&&t<=r){return sum[i];}
    pushdown(i,s,t);
    if(mid>=l)tot+=query(l,r,s,mid,(i<<1));
    tot%=mod;
    if(mid+1<=r)tot+=query(l,r,mid+1,t,(i<<1)|1);
    pushup(i);
    return tot%mod;
}
ll querymin(int l,int r,int s,int t,int i){
	  int mid=(s+t)>>1;
	  ll tot=(ll)1e12;
    if(l<=s&&t<=r){return min1[i];}
    pushdown(i,s,t);
    if(mid>=l)tot=min(tot,querymin(l,r,s,mid,(i<<1)));
    
    if(mid+1<=r)tot=min(tot,querymin(l,r,mid+1,t,(i<<1)|1));
    pushup(i);
    return tot;
}
int main()
{
    int i,j,x,y,q;
    ll z;
    cin>>n>>m>>mod;//n=read();m=read();
    for(i=1;i<=n;i++)a[i]=read();
    build(1,n,1);
    for(i=1;i<=m;i++){
        q=read();
        if(q==1){
            x=read();y=read();scanf("%lld",&z);
            
            cheng(x,y,1,n,1,z);
        }
        else if(q==2){
            x=read();y=read();scanf("%lld",&z);
            
            add(x,y,1,n,1,z);
        }
        else if(q==3){
            x=read();y=read();
            
            printf("%lld\n",query(x,y,1,n,1));
        }
//        else if(q==3){
//            x=read();y=read();
//            ++x,++y;
//            printf("%lld\n",querymin(x,y,1,n,1));
//        }
    }
       return 0;
}

 

四、本题代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<climits>
#include<iomanip>
#include<cmath>
#include<vector>
#include<algorithm>
using namespace std;
#define ll long long
inline int read(){
    int x=0; bool f=1; char c=getchar();
    for(;!isdigit(c);c=getchar()) if(c=='-') f=0;
    for(; isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+(c^'0');
    if(f) return x;
    return 0-x;
}
int n,m;
ll mod;
ll a[100005],sum[400005],min1[400005],max1[400005],mul[400005],laz[400005];
void pushup(int i){sum[i]=(sum[(i<<1)]+sum[(i<<1)|1]);min1[i]=min(min1[(i<<1)],min1[(i<<1)|1]);max1[i]=max(max1[(i<<1)],max1[(i<<1)|1]);}
void pushdown(int i,int s,int t){
    int l=(i<<1),r=(i<<1)|1,mid=(s+t)>>1;
//    if(mul[i]!=1){
//    mul[l]*=mul[i];mul[r]*=mul[i];
//    laz[l]*=mul[i];laz[r]*=mul[i];
//    sum[l]/=mul[i];sum[r]/=mul[i];
//    mul[i]=1;
//    }
    if(laz[i]){
    sum[l]+=laz[i]*(mid-s+1);
    sum[r]+=laz[i]*(t-mid);
    min1[l]+=laz[i];
    min1[r]+=laz[i];
    max1[l]+=laz[i];
    max1[r]+=laz[i];
    laz[l]+=laz[i];
    laz[r]+=laz[i];
    laz[i]=0;
    }
}
void build(int s,int t,int i){
//    mul[i]=1;
    if(s==t){sum[i]=min1[i]=max1[i]=a[s];return;}
    int mid=(s+t)>>1;
    build(s,mid,i<<1);build(mid+1,t,(i<<1)|1);pushup(i);
}
void chu(int l,int r,int s,int t,int i,ll z){
    int mid=(s+t)>>1;
    if(l<=s&&t<=r){
    	//分块除法，不能用简单的乘加维护 
    	ll c1=max1[i]-(long long)floor((double)max1[i]/z);
		ll c2=min1[i]-(long long)floor((double)min1[i]/z);
        if(c1==c2){

            laz[i]-=c1, sum[i]-=(t-s+1)*c1, max1[i]-=c1, min1[i]-=c1;
            return;
        }
//        mul[i]*=z;laz[i]/=z;
//        sum[i]/=z;return;
    }
    pushdown(i,s,t);
    if(mid>=l)chu(l,r,s,mid,(i<<1),z);
    if(mid+1<=r)chu(l,r,mid+1,t,(i<<1)|1,z);
    pushup(i);
}
void add(int l,int r,int s,int t,int i,ll z){
    int mid=(s+t)>>1;
    if(l<=s&&t<=r){
        sum[i]+=z*(t-s+1);laz[i]+=z;
        max1[i]+=z,min1[i]+=z;
        return;
    }
    pushdown(i,s,t);
    if(mid>=l)add(l,r,s,mid,(i<<1),z);
    if(mid+1<=r)add(l,r,mid+1,t,(i<<1)|1,z);
    pushup(i);
}
ll query(int l,int r,int s,int t,int i){
    int mid=(s+t)>>1;ll tot=0;
    if(l<=s&&t<=r){return sum[i];}
    pushdown(i,s,t);
    if(mid>=l)tot+=query(l,r,s,mid,(i<<1));
    
    if(mid+1<=r)tot+=query(l,r,mid+1,t,(i<<1)|1);
    pushup(i);
    return tot;
}
ll querymin(int l,int r,int s,int t,int i){
	  int mid=(s+t)>>1;
	  ll tot=(ll)1e12;
    if(l<=s&&t<=r){return min1[i];}
    pushdown(i,s,t);
    if(mid>=l)tot=min(tot,querymin(l,r,s,mid,(i<<1)));
    
    if(mid+1<=r)tot=min(tot,querymin(l,r,mid+1,t,(i<<1)|1));
    pushup(i);
    return tot;
}
int main()
{
    int i,j,x,y,q;
    ll z;
    cin>>n>>m;//n=read();m=read();
    for(i=1;i<=n;i++)a[i]=read();
    build(1,n,1);
    for(i=1;i<=m;i++){
        q=read();
        if(q==2){
            x=read();y=read();scanf("%lld",&z);
            ++x,++y;
            chu(x,y,1,n,1,z);
        }
        else if(q==1){
            x=read();y=read();scanf("%lld",&z);
            ++x,++y;
            add(x,y,1,n,1,z);
        }
        else if(q==4){
            x=read();y=read();
            ++x,++y;
            printf("%lld\n",query(x,y,1,n,1));
        }
        else if(q==3){
            x=read();y=read();
            ++x,++y;
            printf("%lld\n",querymin(x,y,1,n,1));
        }
    }
       return 0;
}

 

zjyang12345
关注
专栏目录
【题解】sdoj2790【雅礼集训 2017 Day1】市场 线段树区间整除
翻过这座山，他们就会听到你的故事
08-09 549
把除变成减这种操作惊了……写到后面写昏头，一个sb错一直没查出来…… 问题描述 输入 第一行为两个空格隔开的整数 n,q 分别表示商贩个数和政令 + 询问个数。 第二行包含 n 个由空格隔开的整数 a0∼an−1 接下来 q 行，每行表示一个操作，第一个数表示操作编号 1∼4 ，接下来的输入和问题描述一致。 输出 对于每个 3、4 操作，输出询问答案。 样例输入 10 1...
LOJ---6285：数列分块入门 9【区间众数】
KobeDuu
03-29 881
题目： https://loj.ac/problem/6285 分析： 如果只查询众数的个数，完全可以莫队，加数时容易维护众数的数量，删除数时，众数的数量要么减1，要么不变，只需再开一个标记数组维护众数的数量即可 根据陈立杰---《区间众数解题报告》实现了下面两种解法 解法一： 块数分成sqrt(n)超时了，150可以过（分块真毒瘤） #pragma GCC optimize(2)...
「雅礼集训 2017 Day1」市场（线段树维护区间除）
Stargazer的博客
01-04 536
输入 第一行为两个空格隔开的整数 n,q 分别表示商贩个数和政令 + 询问个数。 第二行包含 n 个由空格隔开的整数 a0∼an−1 接下来 q 行，每行表示一个操作，第一个数表示操作编号 1∼4 ，接下来的输入和问题描述一致。 输出 对于每个 3、4 操作，输出询问答案。 样例输入 10 10 -5 -4 -3 -2 -1 0 1 2 3 4 1 0 4 1 1 5 9 1 2 0 9 3 ...
CCF 201709-5 除法（线段树）
aoy29758的博客
09-15 516
操作1：给[l, r]中v的倍数除v 操作2：查询[l, r]的和 思路：类似势能线段树（例题：HDU4027）的思想，首先忽略v=1的操作，然后，1e6数据范围内的数就算每次都除以2，也不用太多次就能变为1，而变为1之后就不用再处理了，所以遍历区间[l, r]找到需要修改的数进行单点修改就行了 代码： 1 #include<cstdio> 2 #inclu...
线段树详解 （原理，实现与应用）
热门推荐
岩之痕
09-09 14万+
线段树详解 By 岩之痕 一：综述 线段树是一种可以快速进行区间修改和区间查询的数据结构。点修改，区间修改和区间查询的复杂度都是O(log2(n)) 并且，线段树可以维护很多种类的信息。说到线段树就不得不提一下树状数组，树状数组维护的信息必须满足区间减法，因为树状数组只是维护1到X的信息，通过相减来得到区间信息，代码简单而且速度快，占用空间又少于线段树，所以一般情况下可以用树状数组的地方就
hdu 6315 Naive Operations（线段树区间整除区间）
Endeavor_G的博客
11-08 330
题意: 给出一个长度为nnn初值为0的a数组，以及长度为nnn的bbb数组，b数组是开始给出的固定的。 接下来qqq次操作，add(l,r)add(l,r)add(l,r) 使得数组a区间[l,r][l,r][l,r]所有元素+1，或者查询∑i=lra[i]b[i]\sum^{r}_{i=l}{\frac{a[i]}{b[i]}}∑i=lr​b[i]a[i]​ 思路: 区间操作一般就用线段树来维护...
LOJ-6279-数列分块入门3(分块, 二分)
最新发布
05-24
(分块) $O(n\sqrt n)$ 我们可以预处理出每一个块的信息，每个块又分为若干小块，每个小块中的数都小于等于块中的最大值，可以用二分查找确定小块中满足条件的数的个数，然后用根号分治处理。 时间复杂度 每次查询...
线段树 3（区间最值与历史最值）【线段树】
//(^ o ^)//
10-21 531
>Link luogu P6242 >Description >解题思路 这个大佬写的好好%% 对于每个节点，我们维护：当前最大值，当前次大值，历史最大值，历史次大值（这道题好像不用用到这个），当前最大值的个数，区间和，和四个懒标记 四个懒标记分别为： lazy1：当前最大值要加上的值 lazy2：历史最大值要加上的值 lazy3：当前非最大值要加上的值 lazy4：历史非最大值要加上的值 加操作：全部都加上 取min操作：考虑要取的min值valvalval的不同情况。valv
[Loj.ac-java] Binary Search Tree - 二分搜索树
qq_28033719的博客
08-03 209
问题： 说明： https://oi-wiki.org/ds/bst/ 来 写 一 个 BST。 问题链接：https://loj.ac/problem/104 提交记录： https://loj.ac/submission/889961 输入案例： 传了一份测试数据：https://download.csdn.net/download/qq_28033719/12609664（假装骗积分） 我的代码： 准备知识 1、以往知识： 左旋转 右旋转前驱 后继https://blog...
【LOJ】#2537. 「PKUWC2018」Minimax-线段树合并
远行客
06-23 433
题解 我们从叶子节点逐层递归到根节点。 我们可以建动态开点权值线段树，每个结点上建一颗，递归时考虑合并左右子树的信息。 可以得到，当合并左右儿子到某父结点上时，可以这样转移： 设f[i][j]是在i结点上，权值为j(令权值为j只能从左儿子中转移上来)的概率，则f[i][j]=f[lc[i]][j]∗((1−p)∗sum1[rc[i]](所有值大于i的节点上的概率之和)+p∗sum2[rc[...
ccf除法-线段树模板
imotolove的博客
06-09 793
这题的要求是对区间处理，查询区间和，一看就是线段树模板题，下面是代码 #include<iostream> using namespace std; int pp[100001]; struct q{ int l,r; long long sum;//记录区间和 }p[2000001]; void build(int l,int r,long long x)//建树函数 { ...
线段树维护区间开方/除法
weixin_30340353的博客
06-18 317
今天考试考了一些神仙数据结构 T1 线段树维护区间加，区间开方，区间和 （数据范围：5e5) T2 线段树维护区间加，区间除，区间和，区间最值 对于这些题目，就像是之前考的区间与，区间或一样，除法，开方的操作会让各个数字之间越来越相近，最后变成一串一串连续的数字都是一样的，所以对于这一部分的操作我们一定程度上使用暴力，而如果一段都相等就相当于直接进行区间剪发的操作 那么我们来看如何判断区间一段...
sdoj2790 市场（线段树区间整除）
常胜将军的博客
07-17 618
题意： 分析：区间整除和区间开根号有点像。有人也叫势能线段树，由于每次区间整除之后区间的数都至少要减少到原来的一半，logai次之后区间的数就变为 1了。加上区间加操作，操作的次数也不会多太多。维护一个区间 max、min和sum，当max-max/k=min-min/k时，就相当于区间每个数都减去max-max/k。（最大值整除后作差与最小值整除后作差的差值相同，则可以把差值打标记，即由除变...
WHU 1464(线段树减，除，区间求和)
Feynman1999的博客
08-22 634
problemN个数，M个操作，操作有三类：　　1. D,A,B,C 区间[a,b]中 所有数除以C，计算机整除　　2. M,A,B,C 区间[a,b]所有数减去C　　3. S,A,B 区间[a,b]求和，输出结果思路注意 1.除法不能打lazy（计算机整除问题）而是更新到底，复杂度可以保证（log级别），因为很快降到0 2.注意减法与除法的融合代码示例#include<cst
CCF(除法)：线段树区间修改（50分）+线段树点修改（100分）+线段树（100分）
09-04 194
除法 201709-5 这道题有很多种方法来做，最常用的就是线段树和树状数组。 如果使用线段树来做，就会想到区间修改的update函数。但是这里可能会涉及到v是1或者a[j]是0的情况，所以用这种方法会超时，最多50分。 可以修改一下代码，使用点修改来做这道题。在main函数里面增加一个循环，用来判断。 当然，还有一种方法就是树状数组，这种方法和上面这种方法运行时间相差无几，但是代码量...
线段树 从入门到进阶
dfssbqyrq96616258的博客
09-21 324
线段树是什么？？线段树怎么写？？ 如果你在考提高组前一天还在问这个问题，那么你会与一等奖失之交臂；如果你还在冲击普及组一等奖，那么这篇博客会浪费你人生中宝贵的5~20分钟。 上面两句话显而易见，线段树这个数据结构是一个从萌新到正式OI选手的过渡，是一个非常重要的算法，也是一个对于萌新来说较难的算法。不得不说，我学习了这个算法5遍左右才有勇气写的这篇博客。 但是，对于OI正式选手来说...
[题解]LOJ_6029市场（线段树区间加/除
weixin_38170853的博客
09-24 326
https://www.cnblogs.com/scx2015noip-as-php/p/loj6029.html 一开始以为裸的区间除，维护这一段是否全为0，后来发现还有加法，这样就会被卡了 这题除数达到了1e9，题解说在数很小除数很大的时候区间除会变成区间减，如一段序列由$x,x-1$组成，除一个$d$相当于减去了同样的变化量 所以我们需要维护区间最大最小值，每次判断最大值和最...
CCF201709-5 除法(树状数组 / 线段树)
坤斤拷的博客
09-08 403
传送门：CCF201709-5 除法 同样是17年10月刚接触编程一个月的蒟蒻的我做过的题，当初是暴力做的，结果必然是超时的QAQ。 这道题想不到会是第五题，直接上树状数组就好了，但还是建了一棵树练了练手，下面会有树状数组和线段树版本的代码。 虽然题目要求的更新是对部分数进行除法，实际上，我们可以对原数组进行预处理，找出满足条件的数，得到商，然后在树上单点修改，就变得非常容易了。但仍有一些需...
「雅礼集训 2017 Day1」市场 (线段树除法，区间最小，区间查询)
lisuier的博客
08-09 578
老师说，你们暴力求除法也整不了多少次就归一了，暴力就好了（应该只有log（n）次） 于是暴力啊暴力，结果我归天了。 好吧，在各种题解的摧残下，我终于出了一篇巨好看（chou lou）代码（很多结构体党嫌丑） 那么具体除法怎么实现就是关键了 对于单个点或者区间内的数完全相同的区间，可以做成区间减法 因为除法会使数变小，而相同的数减小的量是相同的， 那么怎么判断区间内的数是否完全相同呢...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值