BZOJ4241 历史研究 莫队算法 堆

欢迎访问~原文出处——博客园-zhouzhendong&AK

去博客园看该题解

题目

Description

IOI国历史研究的第一人——JOI教授,最近获得了一份被认为是古代IOI国的住民写下的日记。JOI教授为了通过这份日记来研究古代IOI国的生活,开始着手调查日记中记载的事件。

日记中记录了连续N天发生的时间,大约每天发生一件。

事件有种类之分。第i天(1<=i<=N)发生的事件的种类用一个整数Xi表示,Xi越大,事件的规模就越大。

JOI教授决定用如下的方法分析这些日记:

1. 选择日记中连续的一些天作为分析的时间段

2. 事件种类t的重要度为t*(这段时间内重要度为t的事件数)

3. 计算出所有事件种类的重要度,输出其中的最大值

现在你被要求制作一个帮助教授分析的程序,每次给出分析的区间,你需要输出重要度的最大值。

Input

第一行两个空格分隔的整数N和Q,表示日记一共记录了N天,询问有Q次。

接下来一行N个空格分隔的整数X1...XN,Xi表示第i天发生的事件的种类

接下来Q行,第i行(1<=i<=Q)有两个空格分隔整数Ai和Bi,表示第i次询问的区间为[Ai,Bi]。

Output

输出Q行,第i行(1<=i<=Q)一个整数,表示第i次询问的最大重要度

Sample Input

5 5

9 8 7 8 9

1 2

3 4

4 4

1 4

2 4

Sample Output

9

8

8

16

16

HINT

1<=N<=10^5

1<=Q<=10^5

1<=Xi<=10^9 (1<=i<=N)

Source

JOI 2013~2014 春季training合宿 竞技1 By PoPoQQQ

题目概括

  给出一个序列,其中有n个数字。

  现在给出m个询问,每次询问格式为L  R。我们设一个值在L~R范围内的重要度为该值乘上该值在L~R范围内出现过的次数。求L~R范围内重要度最大的数值的重要度。

Sample Input

5 5

9 8 7 8 9

1 2

3 4

4 4

1 4

2 4

Sample Output

9

8

8

16

16

样例解释

对于询问1,1~2范围内,9出现了1次,重要度为9;8出现了1次,重要度为8;所以该区间内重要度最大为9;

对于询问2,3~4范围内,7出现了1次,重要度为7;8出现了1次,重要度为8;所以该区间内重要度最大为8;

对于询问3,4~4范围内,8出现了1次,重要度为8;所以该区间内重要度最大为8;

对于询问4,1~4范围内,7出现了1次,重要度为7;8出现了2次,重要度为16;9出现了1次,重要度为9;所以该区间内重要度最大为8;

对于询问5,2~4范围内,7出现了1次,重要度为7;8出现了2次,重要度为16;所以该区间内重要度最大为8;

题解

  首先,无可置疑,离散化总要做的,闭着眼睛先做了。

  假设我们对于一个数值,设置其哈希值为棋离散化后的位置。

  然后思考,在线不会做,那么离线;

  那么就可以使用莫队算法——不会,没事,现学~ 链接1 链接2

  

  然后我还是来概括一下。

  

  我们把询问按照一定的顺序排列,然后大暴力修改边界值(比如从Li~Li-1,Ri~Ri-1这种段的修改),那么时间复杂度为Σ(|Li-Li-1|+|Ri-Ri-1|)  (1<i<=m)

  我们要使得这个值最小,可以采用分块的方法。

  我们把整个区间(共n)分成每 sqrt(n) 一份,然后把询问以 L 所在的块的位置为第一关键字,把R作为第二关键字排序,然后就可以得到一个大约是 m sqrt(n) 的时间复杂度。

  

  不要把莫队之后的处理看的很麻烦,其实就是大暴力!

  

  关键在于莫队的排序。

  

  而本题要求的是最大值,所以要开一个堆来维护(要打线段树也可以,我不拦你,有可能是被卡常~)

  这个堆不简单~

  

  它要支持寻找一个数值所对应的位置。

  

  一开始打成了一边交换值,一边交换位置的,实际上是有些漏洞的。

  

  对于我的算法,我的堆要保存3个量:

  1.堆中元素的值(在L和R移动的时候,进入一个就在其哈希值的位置上加上其值,减掉的话就反一反)。

  2.堆中某一哈希值的元素所在位置。

  3.这样还是不够的,还要给每个元素配上一个它所对应的哈希值。

  至于具体怎么维护~看代码。

  这题我没写过线段树,不知道怎样。

  反正我的堆是跑了46.7秒,写线段树就自己估计一下吧……也许不是所有的线段树都可以过的(要卡常?)……

  再在提醒一下,在改变区间范围时候的大暴力时,对于L、R分别要分两种情况讨论。

  具体还是看代码呗~

代码

#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstdlib>
using namespace std;
typedef long long LL;
const int N=100000+5;
int n,m,bag_size;
int pos[N],hs,b[N],pv[N];
LL a[N],hash[N],v[N];
struct Query{
    int L,R,bh;
    LL ans;
}q[N];
bool cmp(Query a,Query b){
    int La=a.L/bag_size,Lb=b.L/bag_size;
    if (La==Lb)
        return a.R<b.R;
    return La<Lb;
}
bool cmpbh(Query a,Query b){
    return a.bh<b.bh;
}
int find(LL x){
    int le=1,ri=hs,mid;
    while (le<=ri){
        mid=(le+ri)>>1;
        if (hash[mid]==x)
            return mid;
        if (hash[mid]<x)
            le=mid+1;
        else
            ri=mid-1;
    }
}
void up_sift(int x){
    int i=x,j=i>>1;
    while (i>1&&v[i]>v[j]){
        pos[pv[i]]=j,pos[pv[j]]=i;
        swap(pv[i],pv[j]);
        swap(v[i],v[j]);
        i=j,j=i>>1;
    }
}
void down_sift(int x){
    int i=x,j=i<<1;
    while (j<=hs){
        if (j<hs&&v[j]<v[j+1])
            j++;
        if (v[i]>=v[j])
            break;
        pos[pv[i]]=j,pos[pv[j]]=i;
        swap(pv[i],pv[j]);
        swap(v[i],v[j]);
        i=j,j=i<<1;
    }
}
int main(){
//    freopen("mode.in","r",stdin);
//    freopen("mode.out","w",stdout);
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        scanf("%lld",&a[i]),hash[i]=a[i];
    sort(hash+1,hash+n+1);
    hs=1;
    for (int i=2;i<=n;i++)
        if (hash[i-1]!=hash[i])
            hash[++hs]=hash[i];
    for (int i=1;i<=n;i++)
        b[i]=find(a[i]);
    for (int i=1;i<=hs;i++)
        pos[i]=i,v[i]=0,pv[i]=i;
    for (int i=1;i<=m;i++){
        scanf("%d%d",&q[i].L,&q[i].R);
        q[i].bh=i;
    }
    bag_size=sqrt((double)n)+0.5;
    sort(q+1,q+m+1,cmp);
    int L=q[1].L,R=q[1].L-1;
    for (int i=1;i<=m;i++){
        while (R<q[i].R){
            R++;
            v[pos[b[R]]]+=a[R];
            up_sift(pos[b[R]]);
        }
        while (R>q[i].R){
            v[pos[b[R]]]-=a[R];
            down_sift(pos[b[R]]);
            R--;
        }
        while (L<q[i].L){
            v[pos[b[L]]]-=a[L];
            down_sift(pos[b[L]]);
            L++;
        }
        while (L>q[i].L){
            L--;
            v[pos[b[L]]]+=a[L];
            up_sift(pos[b[L]]);
        }
        q[i].ans=v[1];
    }
    sort(q+1,q+m+1,cmpbh);
    for (int i=1;i<=m;i++)
        printf("%lld\n",q[i].ans);
    return 0;
}
时间: 2024-10-24 16:54:53

BZOJ4241 历史研究 莫队算法 堆的相关文章

【bzoj3207】花神的嘲讽计划Ⅰ Hash+STL-map+莫队算法

题目描述 背景 花神是神,一大癖好就是嘲讽大J,举例如下: “哎你傻不傻的![hqz:大笨J]” “这道题又被J屎过了!!” “J这程序怎么跑这么快!J要逆袭了!” …… 描述 这一天DJ在给吾等众蒟蒻讲题,花神在一边做题无聊,就跑到了一边跟吾等众蒟蒻一起听.以下是部分摘录: 1. “J你在讲什么!” “我在讲XXX!” “哎你傻不傻的!这么麻烦,直接XXX再XXX就好了!” “……” 2. “J你XXX讲过了没?” “……” “那个都不讲你就讲这个了?哎你傻不傻的!” “……” DJ对这种情景

莫队算法

Beautiful Girl 题意 给定一个长度为 n 的序列 a[1], a[2], ..., a[n] . m 组询问 (l, r, K) , 求区间 [l, r] 去除重复的数之后的第 K 小. n, m <= 100000 . 分析 莫队算法 + 值域分块. 1 #include <cstdio> 2 #include <cstring> 3 #include <cstdlib> 4 #include <cctype> 5 #include &

CodeForces - 86D 莫队算法

http://codeforces.com/problemset/problem/86/D 莫队算法就是调整查询的顺序,然后暴力求解. 每回可以通过现有区间解ans(l,r)得到区间(l+1,r),(l-1,r),(l,r+1),(l,r-1)的区间解. 调整方式http://blog.csdn.net/bossup/article/details/39236275 这题比那个还要简单,查询的是K^2*Z,很清楚就是莫队算法,然而做的时候没有学过,回来补题补到 关键是我一直没明白为什么重载小于号

codeforces 617 E. XOR and Favorite Number(莫队算法)

题目链接:http://codeforces.com/problemset/problem/617/E 题目: 给你a1 a2 a3 ··· an 个数,m次询问:在[L, R] 里面又多少中 [l, r] 使得 al xor al+1 xor ··· ar 为 k. 题解: 本题只有区间查询没有区间修改,而且数据量不大(10w),所以可以用离线的方法解决. 使用莫队算法来解决,就需要O(1)的修改[L, R+1] .[L, R-1].[L+1, R].[L-1, R]. 详细的莫队可以百度学一

(普通的)莫队算法简单介绍

莫队算法(由莫涛发明的)是一种离线的暴力算法(至少我这么认为).使用莫队算法的条件是,知道一个区间[l, r]的结果,那么也可以快速知道[l + 1, r],[l - 1, r], [l, r - 1], [l, r + 1]这四个区间的结果.于是可以想到,直接通过这样转移来解决一些问题.当然有些出题人机智,故意卡这种暴力,让你从头跑到尾然后从尾跑到头,于是时间复杂度高达O(n2) 而莫队算法就是通过改变处理询问的顺序来降低时间复杂度. 比如说现在知道一个区间[l1, r1],又要转移到[l2,

莫队算法良心讲解

问题:有n个数组成一个序列,有m个形如询问L, R的询问,每次询问需要回答区间内至少出现2次的数有哪些. 朴素的解法需要读取O(nm)次数.如果使用STL的Map来保存出现的次数,每次需要O(nmlogn)的复杂度.有没有更快的方法呢? 注意到询问并没有强制在线,因此我们可以使用离线方法.注意到一点,如果我们有计算完[L, R]的map,那么[L - 1, R].[L + 1, R].[L, R - 1].[L, R + 1]都能够在O(logn)的复杂度得出.是否能安排适当的询问顺序,使得每次

清橙A1206 小Z的袜子(莫队算法)

A1206. 小Z的袜子 时间限制:1.0s   内存限制:512.0MB 总提交次数:744   AC次数:210   平均分:44.44 将本题分享到: 查看未格式化的试题   提交   试题讨论 试题来源 2010中国国家集训队命题答辩 问题描述 作为一个生活散漫的人,小Z每天早上都要耗费很久从一堆五颜六色的袜子中找出一双来穿.终于有一天,小Z再也无法忍受这恼人的找袜子过程,于是他决定听天由命-- 具体来说,小Z把这N只袜子从1到N编号,然后从编号L到R(L 尽管小Z并不在意两只袜子是不是

HYSBZ 2038 小Z的袜子(hose) (莫队算法)

题意:中文题. 析:很著名的莫队算法,先把这个求概率的式子表达出来,应该是分子:C(x1, 2) + C(x2, 2) + C(x3, 2) + ... + C(xn, 2)  分母:C(n, 2),然后化成分数的表达形式,[x1(x1-1)+x2(x2-1)+...+xn(xn-1)] / (n*(n-1))  然后再化简得到 (sigma(xi*xi)  - n) / (n*(n-1)) ,然后就是对每个区间进行运算,离线,把所以的序列分成sqrt(n)块,然后用两个指针,进行对数据的计算.

(预处理+莫队算法)HDU - 5381 The sum of gcd

题意: 一个长度为n的数列,m次查询L到R之间所有连续子序列的gcd之和. 分析: 很明显的莫队算法. 很明显发现了gcd是单调递减的,并且最多存在32个的性质. 想了很久,脑补了许多种方法来拉伸L和R,但是都有漏洞. 实际上,这道题还是比较复杂的.. 在思考的过程中,我没有充分利用gcd的递减性质. 这题其实这题有共通之处,至少在我的做法上是这样的. 可以发现,在R向右拉伸的过程中,增加的和只是从L到R+1中的每一个后缀的和. 向左则为减,L的移动同理. 那么我们只要提前预处理每个位置的前缀所