POJ2104 K-th Number（线段树，二分，vector）
题意
不带修改区间第k小。（n<=100000)
题解
建立线段数和vector数组（vector为当前区间排列之后的序列）（归并）
然后对于每一个询问二分答案。
问题就转化为区间有多少数小于等于二分值。
对于我们每一个遍历的区间（线段数的节点）。
若与询问区间不相交return0。
若完全包含于询问区间则在此区间的vector上二分查找有多少数小于二分值（因为已经排好序，所以很好做）
若有相交部分则继续遍历子树。
1 #include<cstdio> 2 #include<cstring> 3 #include<cmath> 4 #include<iostream> 5 #include<algorithm> 6 #include<vector> 7 using namespace std; 8 const int N=100010; 9 vector<int> d[N*5]; 10 struct tree{ 11 int l,r; 12 }tr[N*8]; 13 const int INF=1000000007; 14 int a[N*2]; 15 int n,m; 16 void build(int l,int r,int now){ 17 tr[now].l=l;tr[now].r=r; 18 if(l==r){ 19 d[now].push_back(-INF); 20 d[now].push_back(a[l]); 21 return; 22 } 23 int mid=(l+r)>>1; 24 build(l,mid,now*2); 25 build(mid+1,r,now*2+1); 26 int size1=d[now*2].size()-1; 27 int size2=d[now*2+1].size()-1; 28 int i=1,j=1; 29 d[now].push_back(-INF); 30 while(i<=size1&&j<=size2){ 31 if(d[now*2][i]>d[now*2+1][j])d[now].push_back(d[now*2+1][j++]); 32 else d[now].push_back(d[now*2][i++]); 33 } 34 while(i<=size1)d[now].push_back(d[now*2][i++]); 35 while(j<=size2)d[now].push_back(d[now*2+1][j++]); 36 return; 37 } 38 int getnum(int now,int K){ 39 int x=1;int y=d[now].size()-1; 40 int tmp=0; 41 while(x<=y){ 42 int mid=(x+y)>>1; 43 // cout<<x<<" "<<y<<endl; 44 if(d[now][mid]<=K){ 45 tmp=mid; 46 x=mid+1; 47 } 48 else y=mid-1; 49 } 50 return tmp; 51 } 52 int check(int l,int r,int now,int K){ 53 // cout<<l<<" "<<r<<endl; 54 if(tr[now].l>r||tr[now].r<l)return 0; 55 if(l<=tr[now].l&&tr[now].r<=r){ 56 // cout<<l<<" "<<r<<" "<<getnum(now,K)<<endl; 57 return getnum(now,K); 58 } 59 int mid=(tr[now].l+tr[now].r)>>1; 60 int tmp1=check(l,r,now*2,K); 61 int tmp2=check(l,r,now*2+1,K); 62 return tmp1+tmp2; 63 } 64 int search(int l,int r,int k){ 65 int tmp; 66 int x=-INF; 67 int y=INF; 68 while(x<=y){ 69 int mid=(x+y)>>1; 70 int num=check(l,r,1,mid); 71 // cout<<mid<<" "<<num<<endl; 72 if(num>=k){ 73 tmp=mid; 74 y=mid-1; 75 } 76 else x=mid+1; 77 } 78 return tmp; 79 } 80 int main(){ 81 scanf("%d%d",&n,&m); 82 for(int i=1;i<=n;i++){ 83 scanf("%d",&a[i]); 84 } 85 build(1,n,1); 86 for(int i=1;i<=m;i++){ 87 int l,r,k; 88 scanf("%d%d%d",&l,&r,&k); 89 int ans=search(l,r,k); 90 printf("%d\n",ans); 91 } 92 return 0; 93 }
