BZOJ 3524 [Poi2014]Couriers(二分+蒙特卡罗)

 

【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=3524

 

【题目大意】

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

 

【题解】

  我们可以在[l,r]中随机一个位置检验这个位置上数是不是答案,
  检测方法可以在数组中保存每个数在序列中的不同位置,二分即可
  选中答案的概率为1/2,我们做k次蒙特卡罗,正确率就为1-(1/2)^k,
  复杂度为O(kmlogn)。

  主席树做法见 LINK

 

【代码】

#include 
#include 
#include 
using namespace std;
const int N=500010;
vector pos[N];
int a[N],l,r,n,m;
int main(){
    while(~scanf("%d%d",&n,&m)){
        for(int i=1;i<=n;i++)pos[i].clear();
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);   
            pos[a[i]].push_back(i);
        }
        while(m--){
            scanf("%d%d",&l,&r);
            int len=r-l+1,find=0;
            for(int k=1;k<=20;k++){
                int u=a[l+rand()%len];
                int st=lower_bound(pos[u].begin(),pos[u].end(),l)-pos[u].begin();
                int en=upper_bound(pos[u].begin(),pos[u].end(),r)-pos[u].begin();
                if((en-st)*2>len){find=u;break;}
            }if(find)printf("%d\n",find);
            else puts("0");
        }
    }return 0;
}

转载于:https://www.cnblogs.com/forever97/p/bzoj3524solution2.html

你可能感兴趣的:(BZOJ 3524 [Poi2014]Couriers(二分+蒙特卡罗))