codechef Chef and Churu

最新推荐文章于 2021-08-13 11:08:29 发布

L_0_Forever_LF

最新推荐文章于 2021-08-13 11:08:29 发布

阅读量499

点赞数

分类专栏： codechef 分块大法好

本文链接：https://blog.csdn.net/L_0_Forever_LF/article/details/52942604

版权

分块大法好同时被 2 个专栏收录

14 篇文章 0 订阅

订阅专栏

codechef

5 篇文章 0 订阅

订阅专栏

传送门
将n个函数分块，预处理出每个位置在每个函数里用到多少次，修改时维护每个块内所有函数的和，询问时被完整覆盖的块直接累加和，被部分覆盖的块用树状数组暴力求和
要用unsigned long long

code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll unsigned long long
using namespace std;

const int maxn = 110000;
const int maxnn = 350;
int n,m,N;
int lowbit(int x){return x&(-x);}
void update(ll *s,int x,ll c){while(x<=n){s[x]+=c;x+=lowbit(x);}}
ll query(ll *s,int x)
{
    ll ret=0;
    while(x){ret+=s[x];x-=lowbit(x);}
    return ret;
}
int ti[maxn],tt[maxnn][maxn],L[maxn],R[maxn];
int id[maxn],st[maxn];
ll s[maxn],f[maxn],a[maxn];

int main()
{
    scanf("%d",&n); N=sqrt(n);
    for(int i=1;i<=n;i++)
    {
        scanf("%llu",&a[i]);
        update(s,i,a[i]);
    }
    for(int i=1;i<=n;i++)
    {
        id[i]=(i-1)/N+1;
        scanf("%d%d",&L[i],&R[i]);
    }
    for(int i=1;i<=id[n];i++) st[i]=(i-1)*N+1;
    st[id[n]+1]=n+1;
    for(int i=1;i<=id[n];i++)
    {
        for(int j=st[i];j<st[i+1];j++) 
        {
            f[i]+=query(s,R[j])-query(s,L[j]-1);
            ti[L[j]]++; ti[R[j]+1]--;
        }
        for(int j=1;j<=n;j++) tt[i][j]=tt[i][j-1]+ti[j];
        for(int j=st[i];j<st[i+1];j++)
        {
            ti[L[j]]--; ti[R[j]+1]++;
        }
    }

    scanf("%d",&m);
    while(m--)
    {
        ll x,y,c; scanf("%llu%llu%llu",&c,&x,&y);
        if(c==1)
        {
            for(int i=1;i<=id[n];i++)
                f[i]+=tt[i][x]*(y-a[x]);
            update(s,x,y-a[x]); a[x]=y;
        }
        else
        {
            int t1=id[x],t2=id[y];
            ll ret=0;
            if(t1==t2)
            {
                for(int i=x;i<=y;i++)
                    ret+=query(s,R[i])-query(s,L[i]-1);
            }
            else
            {
                for(int i=x;i<st[t1+1];i++)
                    ret+=query(s,R[i])-query(s,L[i]-1);
                for(int i=t1+1;i<t2;i++)
                    ret+=f[i];
                for(int i=st[t2];i<=y;i++)
                    ret+=query(s,R[i])-query(s,L[i]-1);
                /*for(int i=x;i<=y;i++)
                    ret+=query(s,R[i])-query(s,L[i]-1);*/
            }
            printf("%llu\n",ret);
        }
    }

    return 0;
}