[BZOJ 1176&COGS 1752][BOI2007]Mokia(CDQ分治)

Description

维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000.

Solution

这道在BZOJ上是权限题啊…然而我在cogs上找到了它

依旧是CDQ分治的论文题,CDQ太强辣…

把查询拆成四个前缀和的形式,在CDQ分治中用时间把操作分为两部分,保证两部分的x各自单调(只有x较小的操作才能对较大的操作产生影响),然后用树状数组维护y这一维

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#define MAXN 2000005
using namespace std;
int s,w,c[MAXN],cnt=0,sign[MAXN],tim=0;
struct Node
{
    int opt,x,y,v,id,ans;
    Node(int opt=0,int x=0,int y=0,int v=0,int id=0):opt(opt),x(x),y(y),v(v),id(id),ans(0){}
}O[200005],t[200005];
bool cmp1(Node a,Node b){return a.x<b.x;}
bool cmp2(Node a,Node b){return a.id<b.id;}
int read()
{
    int x=0,f=1;char c=getchar();
    while(c<‘0‘||c>‘9‘){
        if(c==‘-‘)f=-1;c=getchar();
    }
    while(c>=‘0‘&&c<=‘9‘){
        x=x*10+c-‘0‘;c=getchar();
    }
    return x*f;
}
int lowbit(int x){return x&-x;}
void add(int pos,int x)
{
    while(pos<=w)
    {
        if(sign[pos]!=tim)c[pos]=0,sign[pos]=tim;
        c[pos]+=x,pos+=lowbit(pos);
    }
}
int query(int pos)
{
    int res=0;
    while(pos>0)
    {
        if(sign[pos]==tim)res+=c[pos];
        pos-=lowbit(pos);
    }
    return res;
}
void merge(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1;
    int i=l,j=mid+1,k=l;
    while(i<=mid&&j<=r)
    {
        if(O[i].x<O[j].x)t[k++]=O[i++];
        else t[k++]=O[j++];
    }
    while(i<=mid)t[k++]=O[i++];
    while(j<=r)t[k++]=O[j++];
    for(int i=l;i<=r;i++)O[i]=t[i];
}
void solve(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1;
    for(int i=l,j=mid+1,k=l;k<=r;k++)
    if(O[k].id<=mid)t[i++]=O[k];
    else t[j++]=O[k];
    for(int i=l;i<=r;i++)O[i]=t[i];
    solve(l,mid);
    int j=l;tim++;
    for(int i=mid+1;i<=r;i++)
    {
        while(O[j].x<=O[i].x&&j<=mid)
        {
            if(O[j].opt==1)add(O[j].y,O[j].v);
            j++;
        }
        if(O[i].opt==2)O[i].ans+=query(O[i].y);
    }
    solve(mid+1,r);
    merge(l,r);
}
int main()
{
    s=read(),w=read();
    while(1)
    {
        int opt=read();
        if(opt==3)break;
        if(opt==1)
        {
            int x=read(),y=read(),a=read();
            O[++cnt]=Node(opt,x,y,a,cnt);
        }
        else
        {
            int x1=read(),y1=read(),x2=read(),y2=read();
            O[++cnt]=Node(opt,x1-1,y1-1,1,cnt);
            O[++cnt]=Node(opt,x1-1,y2,-1,cnt);
            O[++cnt]=Node(opt,x2,y1-1,-1,cnt);
            O[++cnt]=Node(opt,x2,y2,1,cnt);
        }
    }
    sort(O+1,O+1+cnt,cmp1);
    solve(1,cnt);
    sort(O+1,O+1+cnt,cmp2);
    for(int i=1;i<=cnt;i++)
    if(O[i].opt==2)
    {
        int res=0;
        res+=O[i].ans*O[i].v,++i;
        res+=O[i].ans*O[i].v,++i;
        res+=O[i].ans*O[i].v,++i;
        res+=O[i].ans*O[i].v;
        printf("%d\n",res);
    }
    return 0;
}
时间: 2024-10-13 19:21:07

[BZOJ 1176&COGS 1752][BOI2007]Mokia(CDQ分治)的相关文章

bzoj 1176 [Balkan2007]Mokia - CDQ分治 - 树状数组

Description 维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000. Input 第一行两个整数,S,W;其中S为矩阵初始值;W为矩阵大小 接下来每行为一下三种输入之一(不包含引号): "1 x y a" "2 x1 y1 x2 y2" "3" 输入1:你需要把(x,y)(第x行第y列)的格子权值增加a 输入

BZOJ 1176 Balkan 2007 Mokia CDQ分治

题目大意:有一些操作,给一个坐标代表的点加上一个数,和求出一个矩形中的所有数的和. 思路:一眼题,二位树状数组水过. ... .. . 哪里不对?W<=2000000.逗我?这n^2能开下? 这个时候CDQ神牛又来帮助我们了. 这个题应该算是CDQ分治的模板题了吧,简单分析一下,其实不难. 写这个题之前建议写一下BZOJ 1935 SHOI 2007 Tree 园丁的烦恼 树状数组这个题,是本题的简化版. 按照正常的解法,我们应该建立一个二位的数据结构,然后分别维护两维的信息.如果用动态开点的线

BZOJ 1176 [Balkan2007]Mokia CDQ分治

题目大意: 维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数Q<=10000,W<=2000000. POJ1195的加强版 没记错的话上午这题还没有中文题目描述的说0.0 好迅速 首先这题看到W就知道二维树状数组挂了 看到M就发现离散化了也搞不了 0.0 这题似乎是CDQ分治被发现之后第二个解决的题目...不过只有会员才知道的世界,今天反应过来刷刷... 修改和询问放在一起分治,一个询问拆分成4个,树状数组处

BZOJ 2726: [SDOI2012]任务安排( dp + cdq分治 )

考虑每批任务对后面任务都有贡献, dp(i) = min( dp(j) + F(i) * (T(i) - T(j) + S) ) (i < j <= N)  F, T均为后缀和. 与j有关的量只有t = dp(j) - F(i) * T(j) , 我们要最小化它. dp(j)->y, T(j)->x, 那么y = F(i) * x + t, 就是给一些点和一个斜率...然后最小化截距, 显然维护下凸包就可以了. 然后因为无比坑爹的出题人....时间可以为负数, 所以要用平衡树维护(

BZOJ 2716: [Violet 3]天使玩偶 [CDQ分治]

传送门 题意: 维护二维点集P,支持以下两个操作(1)插入点(x,y)(2)给定询问(x,y),求点集中离询问点最近的点距离定义为曼哈顿距离Dis(P1,P2)=|x1-x2|+|y1-y2|n,m<=500000x,y<=1000000 时间,$x$,$y$ $CDQ$分治里需要四个象限分类讨论,树状数组维护最大值 然后有两个象限是后缀和 然后跟$PoPoQQQ$学了一个神奇的技巧,树状数组加一个时间戳,就可以不用每次清空之前的操作了 #include <iostream> #i

[BZOJ 1492][NOI2007]货币兑换Cash(CDQ分治+斜率优化Dp)

Description 小Y最近在一家金券交易所工作.该金券交易所只发行交易两种金券:A纪念券(以下简称A券)和 B纪念券(以下 简称B券).每个持有金券的顾客都有一个自己的帐户.金券的数目可以是一个实数.每天随着市场的起伏波动, 两种金券都有自己当时的价值,即每一单位金券当天可以兑换的人民币数目.我们记录第 K 天中 A券 和 B券 的 价值分别为 AK 和 BK(元/单位金券).为了方便顾客,金券交易所提供了一种非常方便的交易方式:比例交易法 .比例交易法分为两个方面:(a)卖出金券:顾客提

BZOJ 1492 NOI 2007 货币兑换Cash CDQ分治+斜率优化DP

题目大意:有两种金券,A和B.每一天有一个rate值,表示购入的比例:还有每一天AB金券的售价.现在给出初始的钱数,问最后能够获得多少钱. 思路:这算是神题了吧,啃论文啃别人代码将近一天才算有点明白. 首先题目中说的可以买一部分或者卖一部分是扯淡的,因为为了最大获利一定要全部买入,全部卖出.朴素的DP方程就好弄了. 设f[i]为第i天最多的B券的数量.那么f[i] = (rate[j] * f[j] * a[i] + f[j] * b[i]) / (rate[i] * a[i] + b[i])

[bzoj] 1176 Mokia || CDQ分治

原题 给出W×W的矩阵(S没有用,题目有误),给出无限次操作,每次操作的含义为: 输入1:你需要把(x,y)(第x行第y列)的格子权值增加a 输入2:你需要求出以左下角为(x1,y1),右上角为(x2,y2)的矩阵内所有格子的权值和,并输出 输入3:表示输入结束 因为修改之间相互独立,所以可以用CDQ. 三个维度分别为时间,x轴,y轴. 简单的三维偏序即可. #include<cstdio> #include<algorithm> #define N 100010 #define

【BZOJ 4170】 4170: 极光 (CDQ分治)

4170: 极光 Time Limit: 30 Sec  Memory Limit: 512 MBSubmit: 121  Solved: 64 Description "若是万一琪露诺(俗称rhl)进行攻击,什么都好,冷静地回答她的问题来吸引她.对方表现出兴趣的话,那就慢 慢地反问.在她考虑答案的时候,趁机逃吧.就算是很简单的问题,她一定也答不上来." --<上古之魔书> 天空中出现了许多的北极光,这些北极光组成了一个长度为n的正整数数列a[i],远古之魔书上记载到:2个