[IOI2007] sails 船帆

显然答案与杆的顺序无关(按每个高度考虑)。

从高到低考虑杆,设此时的状态为(S)(S[i])是高度(i)的帆的数目,初始全为(0),当前杆的高度为(h),杆上需要放置的帆的数目为(k),贪心地,假设两个高度的同等宜选,优先选择更高的;帆尽量放置在(S[i]=0(ile h))的高度上,若还有剩余,设剩下(t(tle k))个,则放置在除去以选择的高度以外,(S[i](ile h))(t)小的位置。

整理一下,每次选出(S[1cdots h])中前(k)小的(相同大小选下标较大的)高度放置帆。

{5,3} 0 0 1 1 1  +0
{4,3} 1 1 2 1 1  +1
{4,1} 1 2 2 1 1  +1
{3,2} 2 2 3 1 1  +3
{3,2} 3 3 3 1 1  +4
{2,1} 3 4 3 1 1  +3

但是此时(S)似乎不好维护。。考虑将过程倒过来,从低到高考虑杆,((S)一开始为空集,其余定义相同),假设两个高度同等宜选,优先选择更低的,其余大致相同。即每次选出(S[1cdots h])中前(k)小的(相同大小选下标较小的)高度放置帆。

{2,1} 1 0        +0
{3,2} 1 1 1      +0
{3,2} 2 2 1      +2
{4,1} 2 2 1 1    +0
{4,3} 3 2 2 2    +4
{5,3} 3 3 3 2 1  +4

这样用平衡树维护(S)(初始大小为0,每次长度变化时补0节点),每次查询全局前(k)小然后整体加一即可。(省去了下标范围的约束)。进一步可发现,连下标都不用维护了


巧妙的分割线(之前的splay已经弃坑了)

观察各个时态的(S),发现它总是一个单调不增的序列,把某时态的(S)中相同且相邻分成一块,如图,蓝色方框表示将选出前(k)小的位置,注意被完全覆盖的块的高度可以直接区间+1,而割开的块区间+1的范围是反过来的,显然这样的块最多一块,于是可以上线段树来维护。

pic

答案也不用每步累加,设(S[i])表示最终状态上高度为(i)的帆的数目,显然总答案是(sum_ifrac{S[i](S[i]-1)}2),这与因为从本文第一句照应。

#include <bits/stdc++.h>
#define LL long long 
const int N=1e5+10;

namespace sgt {
	struct sgtNode {
		int mx,mn,add;
	} t[N<<2];
#define ls (x<<1)
#define rs (x<<1|1)
	void update(int x) {
		t[x].mx=std::max(t[ls].mx,t[rs].mx);
		t[x].mn=std::min(t[ls].mn,t[rs].mn);
	}
	void pushr(int x,int add) {t[x].mn+=add,t[x].mx+=add,t[x].add+=add;}
	void pushdown(int x) {
		if(t[x].add) pushr(ls,t[x].add),pushr(rs,t[x].add),t[x].add=0;
	}
	LL calc(int x,int l,int r) {
		if(l==r) return 1LL*t[x].mx*(t[x].mx-1)/2;
		int mid=(l+r)>>1; pushdown(x);
		return calc(ls,l,mid)+calc(rs,mid+1,r);
	}
	void build(int x,int l,int r) {
		t[x].mn=+1e9,t[x].mx=-1e9;
		if(l==r) return; int mid=(l+r)>>1;
		build(ls,l,mid); build(rs,mid+1,r);
	}
	void insert(int x,int l,int r,int p) {
		if(l==r) return void(t[x].mx=t[x].mn=0);
		int mid=(l+r)>>1; pushdown(x);
		if(p<=mid) insert(ls,l,mid,p);
		else insert(rs,mid+1,r,p);
		update(x);
	}
	void modify(int x,int l,int r,int L,int R) {
		if(L>R) return;
		if(L<=l&&r<=R) return pushr(x,1);
		int mid=(l+r)>>1; pushdown(x);
		if(L<=mid) modify(ls,l,mid,L,R);
		if(mid<R) modify(rs,mid+1,r,L,R);
		update(x);
	}
	int getVal(int x,int l,int r,int p) {
		if(t[x].mn==t[x].mx) return t[x].mn; 
		int mid=(l+r)>>1; pushdown(x);
		if(p<=mid) return getVal(ls,l,mid,p);
		else return getVal(rs,mid+1,r,p); 
	}
	int getRangeL(int x,int l,int r,int w) {
		if(l==r) return l;
		int mid=(l+r)>>1; pushdown(x);
		if(t[ls].mn<=w) return getRangeL(ls,l,mid,w);
		else return getRangeL(rs,mid+1,r,w);
	}
	int getRangeR(int x,int l,int r,int w) {
		if(l==r) return l;
		int mid=(l+r)>>1; pushdown(x);
		if(t[rs].mx>=w) return getRangeR(rs,mid+1,r,w);
		else return getRangeR(ls,l,mid,w);
	}
}

int n,m;
std::pair<int,int> a[N];

int main() {
	scanf("%d",&n);
	for(int i=1; i<=n; ++i) {
		scanf("%d%d",&a[i].first,&a[i].second);
		m=std::max(m,a[i].first);
	}
	std::sort(a+1,a+n+1);
	sgt::build(1,1,m);
	for(int i=1,R=0; i<=n; ++i) {
		while(R<a[i].first) sgt::insert(1,1,m,++R);
		int pos=R-a[i].second+1; //被割开的位置 
		int val=sgt::getVal(1,1,m,pos);
		int rgl=sgt::getRangeL(1,1,m,val);
		int rgr=sgt::getRangeR(1,1,m,val);
		
//		std::cout<<val<<' '<<rgl<<' '<<rgr<<std::endl;
		
		sgt::modify(1,1,m,rgr+1,R);
		sgt::modify(1,1,m,rgl,rgl+a[i].second-(R-rgr)-1);
	}
	printf("%lld
",sgt::calc(1,1,m));
	return 0; 
} 
原文地址:https://www.cnblogs.com/nosta/p/10557737.html