CodeForces - 438D The Child and Sequence 线段树取模

链接：

D. The Child and Sequence

题意：

给定一个长度为n的序列，要求能够执行m次下列操作：

1.查询区间[l,r]的和

2.将区间[l,r]的每一个数%=mod

3.修改第x个数为y

思路：

1、3分别是普通线段树的区间查询和单点修改，只要对线段树稍作修改，使得支持操作2就可以。

如何对一个区间内每个数取模？

首先明确，取模的顺序不能更换，(4%3)%2不能换成(4%2)%3，因此用类似lazy标记的思想想要把模的数记录下来是不可以的。

那么只能硬着头皮每个数取模，可是，这样不会TLE吗？复杂度O(n*m)肯定会T的，这时候就要用到一个性质：

设x%y = z，

①当y<=x/2时，z<x/2，

②当x/2 < y <= x时，z = x-y < x/2，

③当y>x时，z = x，此时没必要取模。

因此，对于一个确定的数a[i]，实际取模的次数不会超过log2(a[i]),当a[i] = 1e9时，取模次数≈30

因此，只需要维护一个区间的最大值，当区间最大值<mod时跳过，否则继续细分，直到为单点时执行取模操作。

复杂度O(mlogn)

代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
static const int maxn = 100010;
static const int INF = 0x3f3f3f3f;
static const double eps = 1e-6;
static const double pi = acos(-1);

void redirect(){
    #ifdef LOCAL
        freopen("test.txt","r",stdin);
    #endif
}
struct node{
    int l,r,mv;
    ll w;
}tree[4*maxn];
void build(int l,int r,int k){
    tree[k].l = l;tree[k].r = r;
    if(l == r){
        scanf("%I64d",&tree[k].w);
        tree[k].mv = tree[k].w;
        return; 
    }
    int mid = (l+r)>>1;
    build(l,mid,k<<1);
    build(mid+1,r,k<<1|1);
    tree[k].w = tree[k<<1].w + tree[k<<1|1].w;
    tree[k].mv = max(tree[k<<1].mv,tree[k<<1|1].mv);
}
int x,y,md;
inline void modify(int k){
    if(tree[k].mv < md)return;
    if(tree[k].l == tree[k].r){
        tree[k].w %= md,tree[k].mv = tree[k].w;
        return;
    }
    int mid = (tree[k].l+tree[k].r)>>1;
    if(x <= mid)modify(k<<1);
    if(y > mid)modify(k<<1|1);
    tree[k].w = tree[k<<1].w + tree[k<<1|1].w;
    tree[k].mv = max(tree[k<<1].mv,tree[k<<1|1].mv);
}
inline void change(int k){
    if(tree[k].l == tree[k].r){
        tree[k].w = y;
        tree[k].mv = tree[k].w;
        return;
    }
    int m = (tree[k].l + tree[k].r)>>1;
    if(x <= m)change(k<<1);
    else change(k<<1|1);
    tree[k].w = tree[k<<1].w + tree[k<<1|1].w;
    tree[k].mv = max(tree[k<<1].mv,tree[k<<1|1].mv);
}
ll ans;
inline void sum(int k){
    if(tree[k].l >= x && tree[k].r <= y){
        ans += tree[k].w;
        return;
    }
    int m=(tree[k].l+tree[k].r)>>1;
    if(x <= m)sum(k<<1);
    if(y > m)sum(k<<1|1);
}
int main(){
    redirect();
    int n,m,op,mod;
    scanf("%d %d",&n,&m);
    build(1,n,1);
    while(m--){
        scanf("%d %d %d",&op,&x,&y);
        if(op == 1){
            ans = 0;
            sum(1);
            printf("%I64d\n",ans);
        }
        else if(op == 2){
            scanf("%d",&md);
            modify(1);
        }
        else change(1);
    }
    return 0;
}