洛谷P4390 [BOI2007]Mokia 摩基亚 题解

已于 2022-05-02 18:10:13 修改
阅读量372 收藏
点赞数 1
分类专栏： OI 文章标签： 算法 数据结构
于 2022-05-02 12:41:15 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_50332374/article/details/124539417
版权

洛谷P4390 [BOI2007]Mokia 摩基亚 题解

题目链接：P4390 [BOI2007]Mokia 摩基亚

题意：摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户 C 的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。

在定位系统中，世界被认为是一个 w × w w×w w×w 的正方形区域，由 1 × 1 1\times 1 1×1 的方格组成。每个方格都有一个坐标 ( x , y ) (x,y) (x,y)， 1 ≤ x , y ≤ w 1\leq x,y\leq w 1≤x,y≤w。坐标的编号从 1 1 1 开始。对于一个 4 × 4 4\times 4 4×4 的正方形，就有 1 ≤ x ≤ 4 1\leq x\leq 4 1≤x≤4， 1 ≤ y ≤ 4 1\leq y\leq 4 1≤y≤4（如图）：

img

请帮助 Mokia 公司编写一个程序来计算在某个矩形区域内有多少名用户。

这道题的题解区一片互相抄袭的烂货，实在受不了了就来水一篇题解

可以发现这道题目就是个二维数点问题，不知道的推荐看看 这篇博客

然后这道题目，我们不能用归并排序来做

为什么呢？因为这个题目它还有一个维度：时间顺序

因此我们可以使用CDQ分治（第一维是时间，第二维是 x x x ，第三维是 y y y ）

这里我们甚至不用按时间排序，它本来就是顺序的

注意到 w ≤ 2 × 1 0 6 w\le 2\times 10^6 w≤2×106 ，有一丢丢大，所以可以不必要地离散化一下 y y y

因此时间复杂度 O ( n log ⁡ 2 n ) O(n\log^2 n) O(nlog2n)

然后就跑到了最优解第一页，嘿嘿树状数组跑的飞快🤤（逃~

代码如下

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define gc() getchar()
#define pc(a) putchar(a)
#define SIZ (int)(1e5+15)
#define lb lower_bound
char buf1[SIZ],*p1,*p2;
char readchar()
{
    if(p1==p2)p1=buf1,p2=buf1+fread(buf1,1,SIZ,stdin);
    return p1==p2?EOF:*p1++;
}
template<typename T>void read(T &k)
{
    char ch=gc();T x=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=gc();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
    k=x*f;
}
template<typename T>void write(T k)
{
    if(k<0){pc('-');k=-k;}
    static T stk[66];T top=0;
    do{stk[top++]=k%10,k/=10;}while(k);
    while(top){pc(stk[--top]+'0');}
}
#define N (int)(2e6+15)
struct node
{
    int x,y,typ,add,id,ans;
}t[N],b[N];
int w,n,Q,pos,tcnt;
int tree[N],ans[N],tmp[N];
void addQ(int a,int b,int c,int d=0,int e=0)
{
    t[++pos]={a,b,c,d,e};
}
#define lowbit(x) (x&(-x))
void add(int x,int v)
{
    while(x&&x<=n)
    {
        tree[x]+=v;
        x+=lowbit(x);
    }
}
int sum(int x)
{
    int res=0;
    while(x>=1)
    {
        res+=tree[x];
        x-=lowbit(x);
    }
    return res;
}
void init()
{
    for(int i=1; i<=pos; i++)
        tmp[i]=t[i].y;
    sort(tmp+1,tmp+1+pos);
    n=unique(tmp+1,tmp+1+pos)-tmp-1;
    for(int i=1; i<=pos; i++)
        t[i].y=lb(tmp+1,tmp+1+n,t[i].y)-tmp;
}
void cdq(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1;
    cdq(l,mid);cdq(mid+1,r);
    int i=mid+1,j=l,cnt=0;
    for(;i<=r; i++)
    {
        while(t[j].x<=t[i].x&&j<=mid)
        {
            if(t[j].typ==1)add(t[j].y,t[j].add);
            b[++cnt]=t[j++];
        }
        if(t[i].typ==2)t[i].ans+=sum(t[i].y);
        b[++cnt]=t[i];
    }
    for(int i=l; i<j; i++)
        if(t[i].typ==1)add(t[i].y,-t[i].add);
    while(j<=mid)
        b[++cnt]=t[j++];
    for(int i=1; i<=cnt; i++)
        t[l+i-1]=b[i];
}
signed main()
{
    read(w);read(w);
    int op;
    while(read(op),op!=3)
    {
        int a,b,c,d;
        if(op==1)
        {
            read(a);read(b);read(c);
            addQ(a,b,1,c);
        }else if(op==2)
        {
            read(a);read(b);read(c);read(d);++Q;
            addQ(c,d,2,1,Q);
            addQ(a-1,d,2,-1,Q);
            addQ(c,b-1,2,-1,Q);
            addQ(a-1,b-1,2,1,Q);
        }
    }
    init();
    cdq(1,pos);
    for(int i=1; i<=pos; i++)
    {
        if(t[i].typ==2)
            ans[t[i].id]+=t[i].add*t[i].ans;
    }
    for(int i=1; i<=Q; i++)
        write(ans[i]),pc('\n');
    return 0;
}

这里还有非离散化版本的

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define gc() getchar()
#define pc(a) putchar(a)
#define SIZ (int)(1e5+15)
#define lb lower_bound
char buf1[SIZ],*p1,*p2;
char readchar()
{
    if(p1==p2)p1=buf1,p2=buf1+fread(buf1,1,SIZ,stdin);
    return p1==p2?EOF:*p1++;
}
template<typename T>void read(T &k)
{
    char ch=gc();T x=0,f=1;
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=gc();}
    while(isdigit(ch)){x=(x<<1)+(x<<3)+(ch^48);ch=gc();}
    k=x*f;
}
template<typename T>void write(T k)
{
    if(k<0){pc('-');k=-k;}
    static T stk[66];T top=0;
    do{stk[top++]=k%10,k/=10;}while(k);
    while(top){pc(stk[--top]+'0');}
}
#define N (int)(2e6+15)
struct node
{
    int x,y,typ,add,id,ans;
}t[N],b[N];
int w,n,Q,pos,tcnt;
int tree[N],ans[N],tmp[N];
void addQ(int a,int b,int c,int d=0,int e=0)
{
    t[++pos]={a,b,c,d,e};
}
#define lowbit(x) (x&(-x))
void add(int x,int v)
{
    while(x&&x<=w)
    {
        tree[x]+=v;
        x+=lowbit(x);
    }
}
int sum(int x)
{
    int res=0;
    while(x>=1)
    {
        res+=tree[x];
        x-=lowbit(x);
    }
    return res;
}
void cdq(int l,int r)
{
    if(l==r)return;
    int mid=(l+r)>>1;
    cdq(l,mid);cdq(mid+1,r);
    int i=mid+1,j=l,cnt=0;
    for(;i<=r; i++)
    {
        while(t[j].x<=t[i].x&&j<=mid)
        {
            if(t[j].typ==1)add(t[j].y,t[j].add);
            b[++cnt]=t[j++];
        }
        if(t[i].typ==2)t[i].ans+=sum(t[i].y);
        b[++cnt]=t[i];
    }
    for(int i=l; i<j; i++)
        if(t[i].typ==1)add(t[i].y,-t[i].add);
    while(j<=mid)
        b[++cnt]=t[j++];
    for(int i=1; i<=cnt; i++)
        t[l+i-1]=b[i];
}
signed main()
{
    read(w);read(w);
    int op;
    while(read(op),op!=3)
    {
        int a,b,c,d;
        if(op==1)
        {
            read(a);read(b);read(c);
            addQ(a,b,1,c);
        }else if(op==2)
        {
            read(a);read(b);read(c);read(d);++Q;
            addQ(c,d,2,1,Q);
            addQ(a-1,d,2,-1,Q);
            addQ(c,b-1,2,-1,Q);
            addQ(a-1,b-1,2,1,Q);
        }
    }
    cdq(1,pos);
    for(int i=1; i<=pos; i++)
    {
        if(t[i].typ==2)
            ans[t[i].id]+=t[i].add*t[i].ans;
    }
    for(int i=1; i<=Q; i++)
        write(ans[i]),pc('\n');
    return 0;
}

转载请说明出处

q779
关注
专栏目录
[BOI2007]名次排序问题(sorting)
farandnear的博客
06-09 1053
题目描述 已知参赛选手的得分，你的任务是按照得分从高到底给出选手的排名。遗憾的是，保存选手信息的数据结构只支持一种操作，即将一个选手从位置i移动到位置j，该移动不改变其他选手的相对位置，即如果i > j，位置j和位置i-1之间的选手的位置都比原来加1，相反如果 i < j，则位置i+1和位置j之间的选手的位置都比原来减一。上述移动的操作的代价定义为i+j，这里，位置编号从1开始。请你编程确定一个
洛谷p4393（BOI2007）Sequence 序列问题
hipamp的博客
08-05 199
发现以后发题目只需要截图就可以了hhh。 感觉这种题好难想啊。看看数据范围，复杂度应该是O(n)O(n)O(n)，这就要我们舍弃脑中各种奇奇怪怪的想法了（当然有大佬可以用强大的数据结构+卡常等特技）。O(n)O(n)O(n)，我们读入就已经O(n)O(n)O(n)了，这就要求我们单独考虑每一个数对答案的贡献了。 我们来看某一个数 aaa 的一边，如果左边那个数 bbb 比他小，那么 aaa 一定...
【BOI2007】Mokia 摩基亚
Stargazer的博客
01-15 264
传送门 CDQ分治板子题 其实和陌上花开差不多 把时间、l、rl、rl、r看成三个维度 像二维树状数组统计一样，把每一个询问拆成四块前缀和相减 然后统计答案的时候容斥一下 那么现在需要考虑的就只有l,rl,rl,r 由于时间自然有序，就只需要考虑二维了 直接上CDQCDQCDQ 注意树状数组空间开到2e62e62e6 #include&lt;bits/stdc++.h&gt; using name...
[BOI2007]Mokia 摩基亚
青烟绕指柔的博客
09-18 258
题目描述 摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。 在定位系统中，世界被认为是一个W×W的正方形区域，由1×1的方格组成。每个方格都有一个坐标(x,y)，1<=x,y<=W。坐标的编号从1开始。对...
BOI 2007 mokia CDQ分治 or 二维树状数组 分析+代码
初学者
07-12 1277
#include<stdio.h> #include<string.h> #include<stdlib.h> #include<algorithm> #include<iostream> using namespace std; int N,S[5000][5000]; int main() { int sign; while(scanf("%d",&sign)!=EOF) { ...
洛谷 P1892 [BOI2003]团伙 （并查集）
qq_38232157的博客
06-21 181
并查集 本题要点： 1、fa数组开 2 * MaxN ， 对于 1 <= i <= n, 用 i + n位置存放的是 i的敌人。 2、敌人的敌人是朋友，那么如果 x 和y是敌人，则 x 和 y + n 是朋友（因为 x 和 y + n 都是y的敌人）, merge(y + n, x) y 和 x + n 是朋友（因为 y 和 x + n 都是x的敌人）, merge(x + n, y) 3、 如果仅仅是 x 和y 是朋友，则只需要 merge(x, y) #include <cstdio&
洛谷P1892团伙 （并查集-反集）
最新发布
m0_72789360的博客
11-04 239
现在有nn个人，他们之间有两种关系：朋友和敌人。我们知道： 一个人的朋友的朋友是朋友 一个人的敌人的敌人是朋友 现在要对这些人进行组团。两个人在一个团体内当且仅当这两个人是朋友。请求出这些人中最多可能有的团体数。
Mokia【题解】
JH_2002的博客
08-05 276
前言 刚开始开心的二维树状数组(真裸)，后面发现他在故意卡树状数组(W有200w)离散化都不行. 所以~~开心~~CDQ。 题面 戳一下就好了。 吐槽 特意上网搜了一下mokia，发现是盗版nokia。这年头，洛基亚还有盗版，盗版比正版还火，都用来出题 sol CDQ除了能用来解决Ｎ维偏序以外，好像就没什么了，所以转换成三位偏序就好了。 具体来说就是，对于加人这种操作，记一下加的...
CDQ分治&整体二分练习1——P4390 [BOI2007]Mokia 摩基亚
Love_mona的博客
07-30 251
蒟蒻的垂死挣扎 很裸的一道CDQ分治题，这题对于每个询问处理为二维前缀的形式，那么就可以分裂成四个查询，然后以时间、x、y三维做CDQ分治，大概怎么搞都可以吧。 #include&lt;iostream&gt; #include&lt;cstdio&gt; #include&lt;cstdlib&gt; #include&lt;cmath&gt; #include&lt;cstring&gt...
cogs 1752. [BOI2007]摩基亚Mokia
樱花满地集于我心，楪舞纷飞祈愿相随
07-05 333
题目链接：http://www.cogs.pro/cogs/problem/problem.php?pid=1752 cdq分治模版题。就是处理前缀和的时候这里变成了矩阵前缀和。我们采用扫描线的方法，用树状数组来维护一个矩阵的前缀和。 然后对于一个询问，我们采用矩阵前缀和的形式来处理这个询问，转化为4个操作。也就是对4个矩阵的处理。 然后看代码就可以了。 #include #in
Bzoj1176：Mokia&Cogs1752：[BOI2007]摩基亚Mokia
Cyhlnj
01-28 265
题目 Cogs 没有Bzoj的权限号 Sol 离线，CDQCDQCDQ分治，把询问拆成444个，变成每次求二维前缀和 那么只要一个修改操作(关键字为时间，xxx，yyy)都在这个询问前，就可以累计答案 那么就成了偏序问题了，直接CDQCDQCDQ 注意当xxx相等时要把修改丢在前面 # include &lt;bits/stdc++.h&gt; # define IL inli...
cogs [BOI2007]摩基亚Mokia CDQ 分治
weixin_30588907的博客
07-15 112
【题目描述】 摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。 在定位系统中，世界被认为是一个W×W的正方形区域，由1×1的方格组成。每个方格都有一个坐标(x,y)...
1752. [BOI2007]摩基亚Mokia (cdq分治模板题)
ok_again的专栏
12-17 1632
题目：1752. [BOI2007]摩基亚Mokia cdq分治论文中的题目，通过cdq分治可以做到降维，然后用树状数组维护就可以了。 #include #include #include #include #define CLR(a, b) memset(a, b, sizeof(a)) #define LL long long using namespace std; const i
[BOI2007]摩基亚Mokia （cdq分治）
Miracle_ma的专栏
07-01 994
摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。 在定位系统中，世界被认为是一个W×W的正方形区域，由1×1的方格组成。每个方格都有一个坐标(x,y)，1<=x,y<=W。坐标的编号从1开始。对于一个4×4的正方形，
cogs 1752 [BOI2007]摩基亚Mokia(cdq分治+树状数组)
johsnows的博客
07-03 451
数据范围过大，所以没有办法二维树状数组和二维线段树什么的。 听说kdtree可做？然而还不会。 这时候cdq分治就派上用场了，首先它不需要二维的空间，所有空间上是没问题的，然后时间复杂度也可以，具体时间复杂度分析下面再说。 这个问题其实可以看成一个三维偏序问题，（时间，横坐标，纵坐标）。 三维貌似看起来有些头疼，假如是二维呢，比如 问题是动态加点，询问一条线段里包含的权值和， 那样就是（
【COGS】1752. [BOI2007]摩基亚Mokia cdq分治模板题
poursoul
12-02 1160
传送门：【COGS】1752. [BOI2007]摩基亚Mokia 题目分析：cdq分治的基本运用。 代码如下： #include #include #include using namespace std ; typedef long long LL ; #define rep( i , a , b ) for ( int i = ( a ) ;
BOI2007 Mokia 摩基亚
a6t2007的博客
12-19 137
题目描述 题解： 容斥，将询问变成4个加权询问。 然后就是cdq了。 代码： #include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define N 200050 #define ll long long inline int ...
Smithsonian Mathematical Formulas and Elliptic Functions by G. Petit Boi
提及的另一本书《Smithsonian Mathematical Formulæ and Tables of Elliptic Functions》是由Edwin P. Adams编写的，该书聚焦于椭圆函数的数学公式和表格。椭圆函数是复分析和数学物理中的重要工具，它们在研究周期...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值