分块大法好

用分快来水“xjr考考你数据结构”真是再爽不过了

xjr考考你数据结构一:虽然有点慢(6724ms)

#include<cstdio>
#include<cmath>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read()
{
    char ch=getchar();int sig=1,x=0;
    while(!isdigit(ch)) {if(ch==‘-‘) sig=-1;ch=getchar();}
    while(isdigit(ch)) x=x*10+ch-‘0‘,ch=getchar();
    return x*sig;
}
const int maxn=1000010;
const int INF=1000000000;
int n,SIZE,A[maxn];
int bl[maxn],st[maxn],en[maxn];
int maxv[maxn],minv[maxn],sumv[maxn];
void recal(int x)
{
    maxv[x]=-INF;minv[x]=INF;sumv[x]=0;
    for(int i=st[x];i<=en[x];i++)
    {
        maxv[x]=max(maxv[x],A[i]);
        minv[x]=min(minv[x],A[i]);
        sumv[x]+=A[i];
    }
}
int main()
{
    n=read();SIZE=(int)sqrt(n);
    for(int i=1;i<=n;i++)
    {
        A[i]=read();
        bl[i]=(i-1)/SIZE+1;
        en[bl[i]]=i;
        if(!st[bl[i]]) st[bl[i]]=i;
    }
    for(int i=1;i<=bl[n];i++) recal(i);
    int q=read(),l,r,x,v;
    while(q--)
    {
        char tp=getchar();
        while(!isalpha(tp)) tp=getchar();
        if(tp==‘Q‘)
        {
            l=read();r=read();
            int mx=-INF,mn=INF,sum=0;
            if(bl[l]==bl[r]) for(int i=l;i<=r;i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
            else
            {
                for(int i=l;i<=en[bl[l]];i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
                for(int i=st[bl[r]];i<=r;i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
                for(int i=bl[l]+1;i<bl[r];i++) mx=max(mx,maxv[i]),mn=min(mn,minv[i]),sum+=sumv[i];
            }
            printf("MaxNum: %d, MinNum: %d, Sum: %d\n",mx,mn,sum);
        }
        else
        {
            x=read();v=read();
            A[x]=v;recal(bl[x]);
        }
    }
    return 0;
}

xjr考考你数据结构二:加个add标记(1342ms)

#include<cstdio>
#include<cmath>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read()
{
    char ch=getchar();int sig=1,x=0;
    while(!isdigit(ch)) {if(ch==‘-‘) sig=-1;ch=getchar();}
    while(isdigit(ch)) x=x*10+ch-‘0‘,ch=getchar();
    return x*sig;
}
const int maxn=100010;
const int INF=1000000000;
int n,SIZE;
int bl[maxn],st[maxn],en[maxn];
long long A[maxn],maxv[maxn],minv[maxn],sumv[maxn],add[maxn];
void recal(int x)
{
    maxv[x]=-INF;minv[x]=INF;sumv[x]=0;
    for(int i=st[x];i<=en[x];i++)
    {
        maxv[x]=max(maxv[x],A[i]);
        minv[x]=min(minv[x],A[i]);
        sumv[x]+=A[i];
    }
}
int main()
{
    n=read();SIZE=(int)sqrt(n);
    for(int i=1;i<=n;i++)
    {
        A[i]=read();
        bl[i]=(i-1)/SIZE+1;
        en[bl[i]]=i;
        if(!st[bl[i]]) st[bl[i]]=i;
    }
    for(int i=1;i<=bl[n];i++) recal(i);
    int q=read(),l,r,v;
    while(q--)
    {
        char tp=getchar();
        while(!isalpha(tp)) tp=getchar();
        if(tp==‘Q‘)
        {
            l=read();r=read();if(l>r) swap(l,r);
            long long mx=-INF,mn=INF,sum=0;
            if(bl[l]==bl[r]) for(int i=l;i<=r;i++) mx=max(mx,A[i]+add[bl[r]]),mn=min(mn,A[i]+add[bl[r]]),sum+=A[i]+add[bl[r]];
            else
            {
                for(int i=l;i<=en[bl[l]];i++) mx=max(mx,A[i]+add[bl[l]]),mn=min(mn,A[i]+add[bl[l]]),sum+=A[i]+add[bl[l]];
                for(int i=st[bl[r]];i<=r;i++) mx=max(mx,A[i]+add[bl[r]]),mn=min(mn,A[i]+add[bl[r]]),sum+=A[i]+add[bl[r]];
                for(int i=bl[l]+1;i<bl[r];i++) mx=max(mx,maxv[i]+add[i]),mn=min(mn,minv[i]+add[i]),sum+=sumv[i]+add[i]*(en[i]-st[i]+1);
            }
            printf("MaxNum: %lld, MinNum: %lld, Sum: %lld\n",mx,mn,sum);
        }
        else
        {
            l=read();r=read();v=read();if(l>r) swap(l,r);
            if(bl[l]==bl[r])
            {
                for(int i=l;i<=r;i++) A[i]+=v;
                recal(bl[l]);
            }
            else
            {
                for(int i=l;i<=en[bl[l]];i++) A[i]+=v;
                for(int i=st[bl[r]];i<=r;i++) A[i]+=v;
                for(int i=bl[l]+1;i<bl[r];i++) add[i]+=v;
                recal(bl[l]);recal(bl[r]);
            }
        }
    }
    return 0;
}

xjr考考你数据结构三:换成set标记(偷懒将维护信息和传递信息写一块了),这时间(999ms)

#include<cstdio>
#include<cmath>
#include<cctype>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read()
{
    char ch=getchar();int sig=1,x=0;
    while(!isdigit(ch)) {if(ch==‘-‘) sig=-1;ch=getchar();}
    while(isdigit(ch)) x=x*10+ch-‘0‘,ch=getchar();
    return x*sig;
}
const int maxn=100010;
const int INF=1000000000;
int n,SIZE;
int bl[maxn],st[maxn],en[maxn];
long long A[maxn],maxv[maxn],minv[maxn],sumv[maxn],set[maxn];
void recal(int x)
{
    if(set[x]) for(int i=st[x];i<=en[x];i++) A[i]=set[x];
    maxv[x]=-INF;minv[x]=INF;sumv[x]=0;
    for(int i=st[x];i<=en[x];i++)
    {
        maxv[x]=max(maxv[x],A[i]);
        minv[x]=min(minv[x],A[i]);
        sumv[x]+=A[i];
    }
    set[x]=0;
}
int main()
{
    n=read();SIZE=(int)sqrt(n);
    for(int i=1;i<=n;i++)
    {
        A[i]=read();
        bl[i]=(i-1)/SIZE+1;
        en[bl[i]]=i;
        if(!st[bl[i]]) st[bl[i]]=i;
    }
    for(int i=1;i<=bl[n];i++) recal(i);
    int q=read(),l,r,v;
    while(q--)
    {
        char tp=getchar();
        while(!isalpha(tp)) tp=getchar();
        if(tp==‘Q‘)
        {
            l=read();r=read();if(l>r) swap(l,r);
            long long mx=-INF,mn=INF,sum=0;
            if(bl[l]==bl[r])
            {
                if(set[bl[l]]) mx=mn=set[bl[l]],sum=set[bl[l]]*(r-l+1);
                else for(int i=l;i<=r;i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
            }
            else
            {
                if(set[bl[l]]) mx=mn=set[bl[l]],sum=set[bl[l]]*(en[bl[l]]-l+1);
                else for(int i=l;i<=en[bl[l]];i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
                if(set[bl[r]]) mx=max(mx,set[bl[r]]),mn=min(mn,set[bl[r]]),sum+=set[bl[r]]*(r-st[bl[r]]+1);
                else for(int i=st[bl[r]];i<=r;i++) mx=max(mx,A[i]),mn=min(mn,A[i]),sum+=A[i];
                for(int i=bl[l]+1;i<bl[r];i++) mx=max(mx,maxv[i]),mn=min(mn,minv[i]),sum+=sumv[i];
            }
            printf("MaxNum: %lld, MinNum: %lld, Sum: %lld\n",mx,mn,sum);
        }
        else
        {
            l=read();r=read();v=read();if(l>r) swap(l,r);
            if(bl[l]==bl[r])
            {
                recal(bl[l]);
                for(int i=l;i<=r;i++) A[i]=v;
                recal(bl[l]);
            }
            else
            {
                recal(bl[l]);recal(bl[r]);
                for(int i=l;i<=en[bl[l]];i++) A[i]=v;
                for(int i=st[bl[r]];i<=r;i++) A[i]=v;
                for(int i=bl[l]+1;i<bl[r];i++) set[i]=v,maxv[i]=minv[i]=set[i],sumv[i]=set[i]*(en[i]-st[i]+1);
                recal(bl[l]);recal(bl[r]);
            }
        }
    }
    return 0;
}

时间: 2024-11-02 11:39:49

分块大法好的相关文章

[bzoj2141][排队] (分块大法好)

Description 排排坐,吃果果,生果甜嗦嗦,大家笑呵呵.你一个,我一个,大的分给你,小的留给我,吃完果果唱支歌,大家乐和和.红星幼儿园的小朋友们排起了长长地队伍,准备吃果果.不过因为小朋友们的身高有所区别,排成的队伍高低错乱,极不美观.设第i个小朋友的身高为hi,我们定义一个序列的杂乱程度为:满足ihj的(i,j)数量.幼儿园阿姨每次会选出两个小朋友,交换他们的位置,请你帮忙计算出每次交换后,序列的杂乱程度.为方便幼儿园阿姨统计,在未进行任何交换操作时,你也应该输出该序列的杂乱程度. I

分块大法吼2

莫对大法,O(nsqrt(n))解决无修改的区间询问问题: czy的后宫3 [题目描述] 上次czy在机房妥善安排了他的后宫之后,他发现可以将他的妹子分为c种,他经常会考虑这样一个问题:在[l,r]的妹子中间,能挑选出多少不同类型的妹子呢? 注意:由于czy非常丧尸,所以他要求在所挑选的妹子类型在[l,r]中出现次数为正偶数,你懂得. 问题简述:n个数,m次询问,每次问[l,r]区间有多少个数恰好出现正偶数次 [输入格式] 第一行3个整数,表示n,c,m 第二行n个数,每个数Ai在[1,c]之间

骗分大法之-----分块||迷之线段树例题a

什么是分块呢? 就是一种可以帮你骗到不少分的神奇的算法. 分块的写法有几种,我所知道的有①预处理②不预处理 不预处理的代码我看得一脸懵逼 所以我在这里就谈一下预处理的版本www 首先看一道题: 给定一个包含n个数的序列,初值全为0,现对这个序列有两种操作:操作1:把 给定 第k1 个数改为k2;操作2:查询 从第k1个数到第k2个数得最大值.(k1<=k2<=n) 所有的数都 <=100000 好的,如果我们搞遍历肯定超时到爆炸. 那么就要用到分块大法了 把这n个数分成若干块,然后每个块

BZOJ 2738 矩阵乘法 分块

题目大意:给定一个矩阵,多次求一个子矩阵中的第k小 正解:CDQ分治 不会 二维莫队? 不会 于是果断分块大法好(又是 我们将这n*n个数排序 分n次插入 每次插入n个 每次插入后 去链表上处理尚未出解的询问(我懒得写链表写了并查集) 如果当前询问的子矩阵内已经插入大于等于k个数 那么答案一定在当次插入的n个数中 暴力查找即可 时间复杂度O(n^3+nq) 好卡-- #include<cstdio> #include<cstring> #include<iostream>

分块感想

分块大法好 数据结构+算法=程序,这句话是个大神讲的.可见数据结构几乎占了算法题中的半壁江山.   数据结构的基本思想就是分块:树状数组是按其二进制表达来分块的,线段树是根据其区间长度分块--(将原序列分解为若干连续不想相交的子序列)   分块就被誉为"万能数据结构",其思想个人认为比较好理解. 就是讲数据打包,对整个包的数据的询问可以直接预处理得出.修改包中的元素则需维护整个包的信息.(这是分块我认为最重要的两个性质(优缺点)):但分块的大小不同于树状数组,线段树--(其块的大小不变

BZOJ 2038: [2009国家集训队]小Z的袜子(hose) 分块

分块大法好 2038: [2009国家集训队]小Z的袜子(hose) Time Limit: 20 Sec  Memory Limit: 259 MB Submit: 2938  Solved: 1303 [Submit][Status] Description 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命-- 具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只

【分块】bzoj1901 Zju2112 Dynamic Rankings

区间k大,分块大法好,每个区间内存储一个有序表. 二分答案,统计在区间内小于二分到的答案的值的个数,在每个整块内二分.零散的暴力即可. 还是说∵有二分操作,∴每个块的大小定为sqrt(n*log2(n))比较快呢. 1 #include<cstdio> 2 #include<algorithm> 3 #include<cstring> 4 #include<cmath> 5 using namespace std; 6 int n,a[10001],num[

codevs 1082 线段树练习 3 --分块练习

时间限制: 3 s 空间限制: 128000 KB 题目等级 : 大师 Master 题目描述 Description 给你N个数,有两种操作: 1:给区间[a,b]的所有数增加X 2:询问区间[a,b]的数的和. 输入描述 Input Description 第一行一个正整数n,接下来n行n个整数, 再接下来一个正整数Q,每行表示操作的个数, 如果第一个数是1,后接3个正整数, 表示在区间[a,b]内每个数增加X,如果是2, 表示操作2询问区间[a,b]的和是多少. pascal选手请不要使用

HDU5057(分块)

Argestes and Sequence Time Limit: 5000/2500 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 1279    Accepted Submission(s): 373 Problem Description Argestes has a lot of hobbies and likes solving query problems espec