数列分块入门 7

新技能：乘法标记
题目链接

#include<map>
#include<set>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long int ll;
const int N=100000+10;
const int MOD=10007;
int v[N],tag[N],mu[N],pos[N];
int n,m;
void reset(int x)
{
    for(int i=(x-1)*m+1;i<=min(x*m,n);i++)
        v[i]=(v[i]*mu[x]+tag[x])%MOD;
    tag[x]=0,mu[x]=1;
}
void add(int l,int r,int c)
{
    reset(pos[l]);
    for(int i=l;i<=min(pos[l]*m,r);i++)
        v[i]+=c,v[i]%=MOD;
    if(pos[l]!=pos[r])
    {
        reset(pos[r]);
        for(int i=(pos[r]-1)*m+1;i<=r;i++)
        v[i]+=c,v[i]%=MOD;
    }
    for(int i=pos[l]+1;i<=pos[r]-1;i++)
        tag[i]+=c,tag[i]%=MOD;
}
void muil(int l,int r,int c)
{
    reset(pos[l]);
    for(int i=l;i<=min(pos[l]*m,r);i++)
        v[i]*=c,v[i]%=MOD;
    if(pos[l]!=pos[r])
    {
        reset(pos[r]);
        for(int i=(pos[r]-1)*m+1;i<=r;i++)
        v[i]*=c,v[i]%=MOD;
    }
    for(int i=pos[l]+1;i<=pos[r]-1;i++)
    {
        tag[i]=(tag[i]*c);
        tag[i]%=MOD;
        mu[i]=(mu[i]*c);
        mu[i]%=MOD;
    }
}
int query(int x)
{
    return (v[x]*mu[pos[x]]+tag[pos[x]])%MOD;
}
int main()
{
  //  int n;
    scanf("%d",&n);
    m=sqrt(n);
    for(int i=1;i<=n;i++)scanf("%d",v+i);
    for(int i=1;i<=n;i++)
        pos[i]=(i-1)/m+1;
    for(int i=1;i<=pos[n];i++)mu[i]=1;
    for(int i=1;i<=n;i++)
    {
        int f,l,r,c;
        scanf("%d%d%d%d",&f,&l,&r,&c);
        if(f==0)add(l,r,c);
        if(f==1)muil(l,r,c);
        if(f==2)printf("%d\n",query(r));

    }
    return 0;
}