第九周周任务——可持久化线段树(主席树)

前言
之前了解过一点点，但是没有深入学习，所以趁这次周任务的机会系统的学一下，之后应该会和线段树一起写一篇博客。
周任务如下

目录
Minimum Inversion Number
K-th Number
2019牛客暑期多校训练营（第九场）Cutting Bamboos

代码

Minimum Inversion Number

//http://acm.hdu.edu.cn/showproblem.php?pid=1394
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=5010;
struct Node{
	int l,r,sum;	
}node[maxn];
void build(int k,int l,int r){
	node[k].l=l;
	node[k].r=r;
	node[k].sum=0;
	if(l==r){
		return ;
	}
	int mid=(l+r)>>1;
	build(k<<1,l,mid);
	build(k<<1|1,mid+1,r);
}
int query(int k,int l,int r){
	int ans=0;
	if(l==node[k].l && r==node[k].r) 
		return node[k].sum;
	int mid=node[k].l+node[k].r>>1;
	if(l<=mid) ans+=query(k<<1,l,min(mid,r));
	if(r>mid) ans+=query(k<<1|1,max(l,mid+1),r);
	return ans;
}
void update(int k,int x){
	if(node[k].l==node[k].r){
		node[k].sum=1;
		return;
	}
	int mid=node[k].l+node[k].r>>1;
	if(x<=mid) update(k<<1,x);
	else update(k<<1|1,x);
	node[k].sum=node[k<<1].sum+node[k<<1|1].sum;
}
int main(){
	int n;
	while(~scanf("%d",&n)){
		int a[5019],sum=0;
		build(1,0,n);
		for(int i=1;i<=n;i++){
			scanf("%d",&a[i]);
			update(1,a[i]);
			sum+=query(1,a[i]+1,n);
		}
		int ans=sum;
		for(int i=1;i<=n;i++){	
			sum+=n-2*a[i]-1;
			ans=min(sum,ans);
		}
		printf("%d\n",ans);
	} 
	return 0;
} 

K-th Number

#include<stdio.h>
#include<algorithm>
using namespace std;
const int maxn=100010;
int n,a[maxn],b[maxn],rt[maxn],tot=0;
struct Node{
	int r,l,val;
}node[maxn*20]; 
void update(int l,int r,int x,int &y,int pos){
	node[++tot]=node[x];
	node[tot].val++;
	y=tot;
	if(l==r) return ;
	int mid=r+l>>1;
	if(mid>=pos) update(l,mid,node[x].l,node[y].l,pos);
	else update(mid+1,r,node[x].r,node[y].r,pos);
}
int query(int l,int r,int x,int y,int k){
	if(l==r) return l;
	int mid=r+l>>1;
	int sum=node[node[y].l].val-node[node[x].l].val;
	if(sum>=k)query(l,mid,node[x].l,node[y].l,k);
	else query(mid+1,r,node[x].r,node[y].r,k-sum);
} 
int main(){
	int m;
	scanf("%d%d",&n,&m);
	for(int i=0;i<n;i++){
		scanf("%d",&a[i]);
		b[i]=a[i];
	}
	sort(b,b+n);
	unique(b,b+n);
	for(int i=1;i<=n;i++){
		int q=lower_bound(b,b+n,a[i-1])-b;
		update(1,n,rt[i-1],rt[i],q+1);
	}
	for(int i=0;i<m;i++){
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		printf("%d\n",b[query(1,n,rt[x-1],rt[y],z)-1]);
	}
	return 0;
} 

2019牛客暑期多校训练营（第九场）Cutting Bamboos

理清思路中

总结

这周虽然放假，但是那个出数据真的难到我了，搞了几天没搞好，难受的一匹，下周就要上交了，软考和四级口语要到了，感觉时间要抓紧了，不能混日子了。