莫队的各种杂题，博主真不要脸分明是水题

最新推荐文章于 2023-08-09 21:25:19 发布
最新推荐文章于 2023-08-09 21:25:19 发布
阅读量345 收藏 1
点赞数 1
分类专栏： 暴力算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Blue_ACAC89/article/details/79399292
版权

莫队这东西其实就是一种优化了的暴力，具体原理很好理解，不用我吹，复杂度证明也不用我吹。

NO.1：小Z的袜子

如题是一道水题，其实是莫队经典题，莫队就是出自这道题的。这道题是当年莫涛队长出的一道题，是第一道莫队题。

如果你会莫队就一定会，代码如下：

#include"cstdio"
#include"cmath"
#include"algorithm"
#include"cstring"
#include"iostream"
#define LL long long
#define Rep(i,n,x) for(i=x;i<=n;i++)
using namespace std;
struct mes
{
    LL l,r,ans1,ans2,dir;
};
mes c[50001];
LL a[50001],b[50001];
int k;
inline LL read()
{
    char c=getchar();int x=0;
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    return x;
}
bool cmp1(mes x,mes y)
{
    return ((x.l/k<y.l/k)||((x.l/k==y.l/k)&&(x.r<y.r)));
}
bool cmp2(mes x,mes y)
{
    return x.dir<y.dir;
}
inline void gcd(LL &x,LL &y)
{
    if(y==0)x=0,y=1;
    else
    { 
        int p1=x,p2=y;
        if(p1<p2)swap(p1,p2);
        while(p1!=0&&p2!=0)
        {
            p1=p1%p2;
            if(p1<p2)swap(p1,p2);
        }
        int p=p1?p1:p2;
        x/=p;y/=p;
    }	
}
int main()
{
    int n=read(),m=read();
    k=sqrt(n);
    int i;
    Rep(i,n,1)
        b[i]=read();
    LL l=1,r=0;
    Rep(i,m,1)
    {
        c[i].l=read(),c[i].r=read();
        c[i].dir=i;
    }
    sort(c+1,c+m+1,cmp1);
    memset(a,0,sizeof(a));
    LL ans=0;
    Rep(i,m,1)
    {
        while(l>c[i].l){
            l--;
            ans+=a[b[l]];
            a[b[l]]++;
        }
        while(r<c[i].r){
            r++;
            ans+=a[b[r]];
            a[b[r]]++;
        }
        while(l<c[i].l){
            a[b[l]]--;
            ans-=a[b[l]];
            l++;
        }
        while(r>c[i].r){
            a[b[r]]--;
            ans-=a[b[r]];
            r--;
        }
        c[i].ans1=ans;
        c[i].ans2=(r-l+1)*(r-l)/2;
        gcd(c[i].ans1,c[i].ans2);
    }
    sort(c+1,c+m+1,cmp2);
    Rep(i,m,1)
    printf("%lld/%lld\n",c[i].ans1,c[i].ans2);
    return 0;
}

第一次打莫队码风似乎有点骚。还有由于BZOJ很爆炸所以为了大家浏览方便些接了洛谷的链接。

NO.2：小B的询问

比小Z的袜子水，就更不用解说了。

代码如下：

#include"cstdio"
#include"algorithm"
#include"iostream"
#include"cmath"
#include"cstring"
#define LL long long
#define MAXN 50000+10
#define Rep(i,n,x) for(i=x;i<=n;i++)
using namespace std;
inline LL read()
{
    char c=getchar();LL x=0;
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    return x; 
}
struct asks
{
    LL l,r,NO,ans;
};
asks a[MAXN];
LL num[MAXN],cnt[MAXN],n,m,k,p;
bool cmp1(asks x,asks y)
{
    return ((x.l/p<y.l/p)||((x.l/p==y.l/p)&&(x.r<y.r)));
}
bool cmp2(asks x,asks y)
{
    return x.NO<y.NO;
}
int main()
{
    n=read(),m=read(),k=read();
    p=sqrt(n);
    int i;
    Rep(i,n,1)
    num[i]=read();
    Rep(i,m,1)
    {
        a[i].l=read(),a[i].r=read();
        a[i].NO=i;
    }
    sort(a+1,a+m+1,cmp1);
    memset(cnt,0,sizeof(cnt));
    LL l=1,r=0,ans=0;
    Rep(i,m,1)
    {
        while(l>a[i].l)
        {
            l--;
            ans+=cnt[num[l]]*2+1;
            cnt[num[l]]++;
        }
        while(r<a[i].r)
        {
            r++;
            ans+=cnt[num[r]]*2+1;
            cnt[num[r]]++;
        }
        while(l<a[i].l)
        {
            ans-=cnt[num[l]]*2-1;
            cnt[num[l]]--;
            l++;
        }
        while(r>a[i].r)
        {
            ans-=cnt[num[r]]*2-1;
            cnt[num[r]]--;
            r--;
        }
        a[i].ans=ans;
    }
    sort(a+1,a+m+1,cmp2);
    Rep(i,m,1)
    printf("%lld\n",a[i].ans);
    return 0;
}

NO.3：数颜色

一道待修改的莫队，其实就是比普通莫队多加一维，这道题难度比上道高些。 
#include"cstdio"
#include"algorithm"
#include"cmath"
#include"iostream"
#define Rep(i,n,x) for(i=x;i<=n;i++)
#define MAXN 10000+10
using namespace std;
struct asks
{
    int l,r,time,ans,dir;
};
asks a[MAXN];
int num[MAXN],anses[MAXN],n,m,p;
int color[1000001];
struct change
{
    int place,next;
};
change b[1002];
inline int read()
{
    char c=getchar();int x=0;
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    return x;
}
bool cmp(asks x,asks y)
{
    return ((x.l/p<y.l/p)||((x.l/p==y.l/p)&&(x.r/p<y.r/p))||((x.l/p==y.l/p)&&(x.r/p==y.r/p)&&(x.time<y.time)));
}
int turn(int now,int l,int r,int &ans)
{
    if(b[now].place>=l&&b[now].place<=r)
    {
        color[num[b[now].place]]--;
        if(color[num[b[now].place]]==0)
        ans--;
        color[b[now].next]++;
        if(color[b[now].next]==1)
        ans++;
    }
    swap(num[b[now].place],b[now].next);
}
int main()
{
    n=read(),m=read();
    int i,j;
    Rep(i,n,1)
        num[i]=read();
    int v=1,k=0;
    Rep(i,m,1)
    {
        char p1=getchar();
        while(p1!='Q'&&p1!='R')p1=getchar();
        if(p1=='Q')
        {
            a[++k].l=read();
            a[k].r=read();
            a[k].time=v;
            a[k].dir=k;
        }
        else
        {
            v++;
            b[v].place=read();
            b[v].next=read();
        }
    }
    p=sqrt(n);
    sort(a+1,a+k+1,cmp);
    int l=1,r=0,t=1,ans=0;
    Rep(i,k,1)
    {
        while(a[i].l<l){l--;color[num[l]]++;if(color[num[l]]==1)ans++;}
        while(r<a[i].r){r++;color[num[r]]++;if(color[num[r]]==1)ans++;}
        while(l<a[i].l){color[num[l]]--;if(color[num[l]]==0)ans--;l++;}
        while(r>a[i].r){color[num[r]]--;if(color[num[r]]==0)ans--;r--;}
        while(t<a[i].time){t++;turn(t,l,r,ans);}
        while(t>a[i].time){turn(t,l,r,ans);t--;}
        anses[a[i].dir]=ans;
    }
    Rep(i,k,1)
    printf("%d\n",anses[i]);
    return 0;
}

NO.4:SDOI2009 HH的项链

naive，还可以用主席树来做，虽然很想用主席树来做但我还是忍住了（大雾）。
代码如下： 
#include"cstdio"
#include"cmath"
#include"algorithm"
#define MAXN 1000001
#define Rep(i,n,x) for(i=x;i<=n;i++)
using namespace std;
struct ask
{
    int l,r,dis;
};
inline int read()
{
    char c=getchar();int x=0;
    while(c<'0'||c>'9')c=getchar();
    while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
    return x;
}
ask asks[200001];
int color[MAXN];
int num[50001],anses[200001],p;
bool cmp(ask x,ask y)
{
    return ((x.l/p<y.l/p)||(x.l/p==y.l/p&&x.r<y.r));
}
int main()
{
    int n=read(),m;
    p=sqrt(n);
    int i;
    Rep(i,n,1)
    num[i]=read();
    m=read();
    Rep(i,m,1)
    asks[i].l=read(),asks[i].r=read(),asks[i].dis=i;
    sort(asks+1,asks+m+1,cmp);
    int l=1,r=0,ans=0;
    Rep(i,m,1)
    {
        while(l>asks[i].l){l--;color[num[l]]++;if(color[num[l]]==1)ans++;}
        while(r<asks[i].r){r++;color[num[r]]++;if(color[num[r]]==1)ans++;}
        while(l<asks[i].l){color[num[l]]--;if(color[num[l]]==0)ans--;l++;}
        while(r>asks[i].r){color[num[r]]--;if(color[num[r]]==0)ans--;r--;}
        anses[asks[i].dis]=ans;
    }
    Rep(i,m,1)
    printf("%d\n",anses[i]);
    return 0;
}


bluewzl
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
莫队算法完整总结（普通莫队、带修莫队、树上莫队、回滚莫队）
Gene_I must wait for the sunrise
09-05 2687
普通莫队 带修改莫队 树上带修改莫队 回滚莫队
小B的询问(题解)(莫队)
eternal风度
09-17 158
小B的询问(题解)(莫队) Junlier良心莫队 题目 luoguP2709 小B的询问 code #include&lt;bits/stdc++.h&gt; #define lst long long #define ldb double #define N 50050 using namespace std; const int Inf=1e9; int read() { ...
莫队经典例题
liaoxiyan123的博客
09-19 283
传送门:[AHOI2013]作业 本题做法:莫队+分块(辅助). 1用莫队统计区间内每个数的个数以及所在块中数的个数和种类数. 2.分块按照值域分,整块直接加上处理后的区间的相关数据,碎块暴力枚举并统计在范围内的数. 具体看代码: #include<bits/stdc++.h> #define N 100005 #define in read() #define re register using namespace std; int n,m,l,r,block; int a[N],zhi,f
普通莫队题集
zzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz
08-16 2440
普通莫队： P2709 小B的询问 题目描述 小B有一个序列，包含N个1~K之间的整数。他一共有M个询问，每个询问给定一个区间[L…R]，求Sigma(c(i)^2)的值,其中i的值从1到K，其中c(i)表示数字i在[L…R]中的重复次数。小B请你帮助他回答询问。 输入格式 第一行，三个整数N、M、K。 第二行，N个整数，表示小B的序列。 接下来的M行，每行两个整数L、R。 输出格式 M行，每行一...
莫队专题
倘若天下着雨。。。你还在
03-07 223
莫队的使用范围 莫队是一种暴力的区间查询的方法。复杂度为n*sqrt(n); 以下面一题为例   装饰机房 题目描述 YCJS机房中挂上了好多画！ 机房后面的墙上整齐地挂上了一排n副画，画有许多种类，分别用1,2,3...表示。 但是有许多画都重复了，XiaoC想要知道在一些区间中有多少副不同种类的画，以便他重新装饰这面墙。 输入 第一行包括
博主关于该题的解答源码
06-16
题目完整源码
Scratch天体运动 青少年编程 蓝桥杯Scratch编程 选拔赛真题源码
最新发布
05-17
Scratch天体运动 蓝桥杯Scratch编程选拔赛真题源码 此案例难度系数3；...期待小朋友们相互交流学习，有什么问题，建议或者意见可以直接给博主留言，或者私下，博主看到后会第一时间给到您相应的回复
紫色系旅游美妆博主个人博客网站html模板
02-24
紫色系旅游美妆博主个人博客网站html模板
ins风美食博主宣传推广PPT模板
08-23
ins风美食博主宣传推广PPT模板适用于商务办公设计应用。由【素材】原创
莫队，带修莫队（总结+例题）
不二君
07-28 1101
学习了莫队，感觉这个算法是一个很玄学的的算法，分块的大小和排序方式就能决定这个算法的效率，现在就对莫队，和带修莫队做一个简单的总结。 适用环境： 莫队算法使用分块的思想，可以解决一类离线区间询问问题，在使用莫队算法的时候就需要将区间询问离线（在线也可以，但是那样复杂度很高，就没有必要使用这个算法了）。 算法核心： 莫队算法离线首先是按左端点分块，然后在同一个块当中区间右端点要递增（...
【模拟试题】czy的后宫3 莫队
a15063149271的博客
03-30 229
题目描述 　　上次czy在机房妥善安排了他的后宫之后，他发现可以将他的妹子分为c种，他经常会考虑这样一个问题：在[l,r]的妹子中间，能挑选出多少不同类型的妹子呢？ 注意：由于czy非常丧尸，所以他要求在所挑选的妹子类型在[l,r]中出现次数为正偶数，你懂得。 题目大意 :n个数，m次询问，每次问[l,r]区间有多少个数恰好出现正偶数...
专题训练3 莫队算法题目清单
qq_39898877的博客
09-20 520
bzoj 2048(入门) 题解：https://blog.csdn.net/qq_39898877/article/details/82807192 hdu 4467 题解：https://blog.csdn.net/qq_39898877/article/details/82840862 hdu 4638 题解：https://blog.csdn.net/qq_39898877/ar...
分块/莫队（模板&&习题）
ojzha的博客
08-04 176
分块是一种思想，把一个整体划分为若干个小块，对整块整体处理，零散块单独处理。本文主要介绍块状数组——利用分块思想处理区间问题的一种数据结构。 块状数组把一个长度为nnn 的数组划分为aaa 块，每块长度为 n/an/an/a 。对于一次区间操作，对区间內部的整块进行整体的操作，对区间边缘的零散块单独暴力处理。（所以分块被称为“优雅的暴力”） 这里，块数既不能太少也不能太多。如果太少，区间中整块的数量会很少，我们要花费大量时间处理零散块；如果太多，又会让块的长度太短，失去整体处理的意义。一般来说，我们取块数为
【BZOJ】【P3236】【Ahoi2013】【作业】【题解】【莫队】
退役狗的专栏
11-17 1352
传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3236
NOIP模拟 - 莫队
weixin_30780649的博客
09-16 75
题目描述 给定一个元素个数为 n 的整数数组 a 和 Q 个问题，每个问题有 x，y 两个参数，请统计共有多少个整数 K 满足 K 在 a[x]…a[y] 中出现了恰好 K 次。 输入格式 第一行两个整数 n，Q，表示数组 a 的元素个数和询问数； 接下来一行 n 给整数，描述数组 a ； 接下来 Q 行，每行两个数 xi,yi（1≤xi≤yi≤n），表示询问的左右边界。 输出格式 输出 Q 行，...
莫队算法详解
热门推荐
hnshhslsh的专栏
01-26 1万+
莫队算法详解 本文翻译自MO’s Algorithm (Query square root decomposition)，作者anudeep2011，发表日期为2014-12-28。由于最近碰到一些莫队算法的题目，找到的相关中文资料都比较简略，而这篇英语文章则讲解的比较详细，故翻译成中文与大家分享。由于本人水平有限，错误在所难免，请谅解。
莫队好题选讲
wljoi的博客
01-20 220
前置知识 莫队 例1：小Z的袜子 记cnt[i]cnt[i]cnt[i]为颜色iii的出现次数，那么我们的答案为 ∑cnt[i](cnt[i]−1)(R−L+1)(R−L)\frac{\sum cnt[i](cnt[i]-1)}{(R-L+1)(R-L)}(R−L+1)(R−L)∑cnt[i](cnt[i]−1)​ =(cnt[1]2+cnt[2]2+...+cnt[n]2)−(cnt[1]+cn...
【算法笔记】莫队算法入门
luhaoren2009的博客
08-09 124
本来想学完回滚莫队、树上莫队、二离莫队之后一起写一个博客，但是一直学不会/kk，只好把已会的普通莫队和带修莫队写了（以后会补上的）
gym-101138D(后缀和,莫队算法,容斥原理,好题)
fouzhe的博客
02-08 548
题目链接 D. Strange Queries time limit per test 3 seconds memory limit per test 256 megabytes input standard input output standard output You are given an array with n in
东方博宜oj题库1170
02-08
东方博宜oj题库1170是一个在线编程题库，包含了大量的编程题目，供程序员进行练习和挑战。这个题库中的题目涵盖了各种不同的难度级别和编程语言，可以帮助程序员提升编程能力和解决问题的能力。 具体关于东方博宜oj...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值