zoukankan      html  css  js  c++  java
  • [POJ2104]K-th Number

    1. K-th Number
    Time Limit: 20000MS   Memory Limit: 65536K
    Total Submissions: 34048   Accepted: 10810
    Case Time Limit: 2000MS

    Description

    You are working for Macrohard company in data structures department. After failing your previous task about key insertion you were asked to write a new data structure that would be able to return quickly k-th order statistics in the array segment. 
    That is, given an array a[1...n] of different integer numbers, your program must answer a series of questions Q(i, j, k) in the form: "What would be the k-th number in a[i...j] segment, if this segment was sorted?" 
    For example, consider the array a = (1, 5, 2, 6, 3, 7, 4). Let the question be Q(2, 5, 3). The segment a[2...5] is (5, 2, 6, 3). If we sort this segment, we get (2, 3, 5, 6), the third number is 5, and therefore the answer to the question is 5.

    Input

    The first line of the input file contains n --- the size of the array, and m --- the number of questions to answer (1 <= n <= 100 000, 1 <= m <= 5 000). 
    The second line contains n different integer numbers not exceeding 109 by their absolute values --- the array for which the answers should be given. 
    The following m lines contain question descriptions, each description consists of three numbers: i, j, and k (1 <= i <= j <= n, 1 <= k <= j - i + 1) and represents the question Q(i, j, k).

    Output

    For each question output the answer to it --- the k-th number in sorted a[i...j] segment.

    Sample Input

    7 3
    1 5 2 6 3 7 4
    2 5 3
    4 4 1
    1 7 3

    Sample Output

    5
    6
    3

    Hint

    This problem has huge input,so please use c-style input(scanf,printf),or you may got time limit exceed.

    Source

    Northeastern Europe 2004, Northern Subregion
    一句话题意:题目意思很简单,多次询问区间k小值
    分析:没什么好分析的,模板题,原本想用划分树做,然后去贴吧问……结果被喷“时代的默泪”,ORZ……然后默默的去看主席树了
    (主席树资料网上很多,大家百度一下)
    这里说下主要思想,就是对1~n每个位置建立一颗线段树,表示从开头到此位置每个数字在各个区间出现次数(当然要先离散……),当然这样空间不行,但是可以发现所有线段树形状一样只是sum域不同,且第i颗和第i-1颗相比,只是在第i-1颗的基础上插入了第i个点(即只修改了第i-1颗线段树的一条从上到下的路径,其他都没改变),因此只要建立一个主线段树,没变的指一条边过去就可以了。
    然后询问[a,b]可以发现线段树可以减,即[a,b]的线段树相当于[1,b]线段树(对应第b颗线段树)-[1,a-1]线段树(对应第a-1颗),然后判断下左子树大小和k的关系就行了
    code:
     1 #include<cstdio>
     2 #include<algorithm>
     3 using namespace std;
     4 const int maxn=100010;
     5 struct wjmzbmr
     6 {
     7  int l,r,ls,rs,sum;
     8 }f[maxn*20];
     9 int tot,root[maxn+50],q,b[maxn+50],sortb[maxn+50];
    10 int build(int l,int r)
    11 {
    12  int k=++tot;
    13  f[k]={l,r,0,0,0};
    14  if(l==r) return tot;
    15  int mid=(l+r)>>1;
    16  f[k].ls=build(l,mid);
    17  f[k].rs=build(mid+1,r);
    18  return k;
    19 }
    20 int change(int o,int x,int v)
    21 {
    22  int k=++tot;
    23  f[k]=f[o];f[k].sum+=v;
    24  if(f[o].l==x&&f[o].r==x) return k;
    25  int mid=(f[o].l+f[o].r)>>1;
    26  if(x<=mid) f[k].ls=change(f[o].ls,x,v);else f[k].rs=change(f[o].rs,x,v);
    27  return k; 
    28 }
    29 int ask(int a,int b,int k)
    30 {
    31  if(f[b].l==f[b].r) return f[b].l;
    32  int mid=f[f[b].ls].sum-f[f[a].ls].sum;
    33  if(k<=mid) return ask(f[a].ls,f[b].ls,k);else return ask(f[a].rs,f[b].rs,k-mid);
    34 }
    35 int main()
    36 {
    37  freopen("ce.in","r",stdin);
    38  freopen("ce.out","w",stdout);
    39  int n,m;
    40  while(scanf("%d%d",&n,&m)!=EOF)
    41  {
    42   for(int i=1;i<=n;++i) scanf("%d",&b[i]),sortb[i]=b[i];
    43   sort(sortb+1,sortb+1+n);
    44   q=1;tot=0;
    45   for(int i=2;i<=n;++i) if(sortb[q]!=sortb[i]) sortb[++q]=sortb[i];
    46   root[0]=build(1,q);
    47   for(int i=1;i<=n;++i)
    48   {   
    49    int p=lower_bound(sortb+1,sortb+n+1,b[i])-sortb;
    50    root[i]=change(root[i-1],p,1);
    51   }
    52   for(int i=1;i<=m;++i)
    53   {
    54    int a,b,k;
    55    scanf("%d%d%d",&a,&b,&k);
    56    printf("%d
    ",sortb[ask(root[a-1],root[b],k)]);
    57   }
    58  }
    59  return 0;
    60 }
    View Code
     
  • 相关阅读:
    qml+opencv(二),实现人脸检测
    CPU满格的元凶,这回是由于QTimer引起的(默认interval是0,太猛)
    Qt文档系统分析(解释Qt文档的生成工具与过程)
    2015 最受开发者欢迎的开发工具和服务(摘录对自己有用的)
    RAD Studio 10 自带Demo代码汇总说明
    ACE 容器之三 ACE_Unbounded_Queue的使用
    开发快很重要——如果只看法语或者产品结果C++似乎很强大,但是参与这个C++的开发过程,就会感觉到这种痛苦(Google也是这个看法)
    我是觉得创业永远不晚 但是绝对要谨慎,不是破釜沉舟就一定能赢的(不用拿工资啊 吃喝拉撒在公司报账就好了)
    人生的成功与否,真的是以天来计算啊——时间永远都是最宝贵的!
    thiscall仅仅应用于“C++”成员函数(this指针存放于CX寄存器,参数从右到左压)
  • 原文地址:https://www.cnblogs.com/wmrv587/p/3475971.html
Copyright © 2011-2022 走看看