Fast Matrix Operations ,uva 11992

#include <iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
const int maxn=100010<<2;
int setv[30][maxn],addv[30][maxn],sum[30][maxn],maxnum[30][maxn],minnum[30][maxn];
void PushUp(int num,int rt)
{
    int ls=rt<<1,rs=rt<<1|1;
    sum[num][rt]=sum[num][ls]+sum[num][rs];
    maxnum[num][rt]=max(maxnum[num][ls],maxnum[num][rs]);
    minnum[num][rt]=min(minnum[num][ls],minnum[num][rs]);
}
void PushDown(int num,int l,int r,int rt)
{
    int m=(l+r)>>1;
    int ls=rt<<1,rs=rt<<1|1;
    if(setv[num][rt]>0)
    {
        setv[num][ls]=setv[num][rs]=setv[num][rt];
        addv[num][ls]=addv[num][rs]=0;
        sum[num][ls]=(m-l+1)*setv[num][rt];
        sum[num][rs]=(r-m)*setv[num][rt];
        maxnum[num][ls]=maxnum[num][rs]=setv[num][rt];
        minnum[num][ls]=minnum[num][rs]=setv[num][rt];
        setv[num][rt]=-1;
    }
    if(addv[num][rt]>0)
    {
        addv[num][ls]+=addv[num][rt];
        addv[num][rs]+=addv[num][rt];
        sum[num][ls]+=(m-l+1)*addv[num][rt];
        sum[num][rs]+=(r-m)*addv[num][rt];
        maxnum[num][ls]+=addv[num][rt];
        maxnum[num][rs]+=addv[num][rt];
        minnum[num][ls]+=addv[num][rt];
        minnum[num][rs]+=addv[num][rt];
        addv[num][rt]=0;
    }
}
void build(int num,int l,int r,int rt)
{
    setv[num][rt]=-1;
    addv[num][rt]=0;
    if(l==r)
    {
        sum[num][rt]=maxnum[num][rt]=minnum[num][rt]=0;
        return;
    }
    int m=(l+r)>>1;
    build(num,l,m,rt<<1);
    build(num,m+1,r,rt<<1|1);
    PushUp(num,rt);
}
void Update(int num,int type,int value,int L,int R,int l,int r,int rt)
{
    if(l>=L&&r<=R)
    {
        if(type==1)
        {
            addv[num][rt]+=value;
            sum[num][rt]+=(r-l+1)*value;
            minnum[num][rt]+=value;
            maxnum[num][rt]+=value;
        }
        if(type==0)
        {
            addv[num][rt]=0;
            setv[num][rt]=value;
            sum[num][rt]=(r-l+1)*value;
            minnum[num][rt]=maxnum[num][rt]=value;
        }
        return;
    }
    PushDown(num,l,r,rt);
    int m=(l+r)>>1;
    if(m>=L)
      Update(num,type,value,L,R,l,m,rt<<1);
    if(m<R)
      Update(num,type,value,L,R,m+1,r,rt<<1|1);
    PushUp(num,rt);
}
void Query(int num,int L,int R,int l,int r,int rt,int& s,int& maxv,int& minv)
{
    if(l>=L&&r<=R)
    {
        s=sum[num][rt];
        maxv=maxnum[num][rt];
        minv=minnum[num][rt];
        return;
    }
    PushDown(num,l,r,rt);
    int m=(l+r)>>1;
    int tsum,tminv,tmaxv;
    s=0;maxv=-1;minv=inf;
    if(m>=L)
    {
        Query(num,L,R,l,m,rt<<1,tsum,tmaxv,tminv);
        s+=tsum;
        maxv=max(maxv,tmaxv);
        minv=min(minv,tminv);
    }
    if(m<R)
    {
        Query(num,L,R,m+1,r,rt<<1|1,tsum,tmaxv,tminv);
        s+=tsum;
        maxv=max(maxv,tmaxv);
        minv=min(minv,tminv);
    }
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    int m,n,q;
    while(cin>>m>>n>>q)
    {
        for(int i=1;i<=m;++i)
        {
            build(i,1,n,1);
        }
        int type,x1,y1,x2,y2,v;
        while(q--)
        {
            cin>>type;
            if(type==1)
            {
                cin>>x1>>y1>>x2>>y2>>v;
                for(int i=x1;i<=x2;++i)
                  Update(i,1,v,y1,y2,1,n,1);
            }
            else if(type==2)
            {
                cin>>x1>>y1>>x2>>y2>>v;
                for(int i=x1;i<=x2;++i)
                  Update(i,0,v,y1,y2,1,n,1);
            }
            else
            {
                cin>>x1>>y1>>x2>>y2;
                int s=0,minv=inf,maxv=-1;
                int ts,tmin,tmax;
                for(int i=x1;i<=x2;++i)
                {
                    Query(i,y1,y2,1,n,1,ts,tmax,tmin);
                    s+=ts;
                    maxv=max(maxv,tmax);
                    minv=min(minv,tmin);
                }
                cout<<s<<" "<<minv<<" "<<maxv<<endl;
            }
        }
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值