bzoj 3524 [Poi2014]Couriers

最新推荐文章于 2019-01-20 11:10:19 发布

原创最新推荐文章于 2019-01-20 11:10:19 发布 · 345 阅读

0 ·

CC 4.0 BY-SA版权

可持久化线段树专栏收录该内容

1 篇文章

订阅专栏

本文介绍了一道名为Couriers的算法题目（POI2014），并提供了一种使用主席树的数据结构解决方案。任务是在给定的序列中查找是否存在一个数，在指定区间内出现的频率超过区间长度的一半。

3524: [Poi2014]Couriers
Time Limit: 20 Sec Memory Limit: 256 MB
Submit: 2203 Solved: 850
[Submit][Status][Discuss]
Description

给一个长度为n的序列a。1≤a[i]≤n。
m组询问，每次询问一个区间[l,r]，是否存在一个数在[l,r]中出现的次数大于(r-l+1)/2。如果存在，输出这个数，否则输出0。

Input

第一行两个数n，m。
第二行n个数，a[i]。
接下来m行，每行两个数l,r，表示询问[l,r]这个区间。

Output

m行，每行对应一个答案。

Sample Input

7 5

1 1 3 2 3 4 3

1 3

1 4

3 7

1 7

6 6

Sample Output

HINT

【数据范围】

n,m≤500000

2016.7.9重设空间，但未重测！

Source

By Dzy

【分析】
发现性质：答案唯一
然后就可以上主席树了= =

【代码】

//bzoj 3524
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mid (l+r>>1)
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int N=500005;
const int mxn=10000005;
int n,m,tot,cnt;
int a[N],b[N],root[N];
int ls[mxn],rs[mxn],c[mxn];
inline void build(int l,int r,int u,int &v,int k)
{
    v=(++cnt),c[v]=c[u]+1;
    if(l==r) return;
    ls[v]=ls[u],rs[v]=rs[u];
    if(k<=mid) build(l,mid,ls[u],ls[v],k);
    else build(mid+1,r,rs[u],rs[v],k);
}
inline int query(int l,int r,int u,int v,int k)
{
    if(l==r) return l;
    if(c[rs[v]]-c[rs[u]]>=k) return query(mid+1,r,rs[u],rs[v],k);
    if(c[ls[v]]-c[ls[u]]>=k) return query(l,mid,ls[u],ls[v],k);
    return 0;
}
int main()
{
    int i,j,l,r;
    scanf("%d%d",&n,&m);
    fo(i,1,n) scanf("%d",&a[i]),b[i]=a[i];
    sort(b+1,b+n+1);
    tot=unique(b+1,b+n+1)-b-1;
    fo(i,1,n) a[i]=lower_bound(b+1,b+tot+1,a[i])-b;
    fo(i,1,n) build(1,tot,root[i-1],root[i],a[i]);
    while(m--)
    {
        scanf("%d%d",&l,&r);
        printf("%d\n",b[query(1,tot,root[l-1],root[r],(r-l+1)/2+1)]);
    }
    return 0;
}
/*
7 5
1 1 3 2 3 4 3
1 3
1 4
3 7
1 7
6 6
*/