码迷,mamicode.com
首页 > 其他好文 > 详细

P5979 [PA2014]Druzyny dp 分治 线段树 分类讨论 启发式合并

时间:2020-07-28 22:54:01      阅读:115      评论:0      收藏:0      [点我收藏+]

标签:register   mes   end   怎么   freopen   std   details   its   lib   

LINK:Druzyny

这题研究了一下午 终于搞懂了.

\(n^2\)的dp很容易得到.

考虑优化.又有大于的限制又有小于的限制这个非常难处理.

不过可以得到在限制人数上界的情况下能转移到的最远端点 且这个数组是单调的.

而下界是随意的.

这个可以利用单调队列做 也可以暴力线段树.

然后考虑怎么优化 考虑CDQ分治 不过这里算是基于时间的分治吧.

大体上就是求出来左边让左边对右边进行贡献.

考虑这部分怎么做 左边为 \([l,mid]\)右边为\(mid+1,r\) 分别设两个决策\(j,i\)分属这两个集合.

容易得到\(left_i\leq j,C_{j+1,i}\leq i-j\)后者的限制条件还是需要暴力转移.

考虑对最大值进行分治 即对笛卡尔树上每个节点进行分治 这样C的条件就很明朗了.

设分治中心为x 那么左区间为\(l,x-1\)右区间为\(x,r\)

刚才的条件就变成了\(left_i\leq j,C_x\leq i-j\)

考虑右边区间小的话就暴力扫右边区间 算出对应左区间 然后线段树上查区间最大值.

左边小的话 暴力搜啊左边二分出\(left\)的合法端点 算出对应右区间 然后线段树上区间修改.

可以发现上述暴力扫的过程类似启发式合并 这样复杂度为\(nlog^2n\)还不足以通过此题.

正解是分类讨论了几种情况:

首先我们知道 对于\(left_i\)是单调递增的 根据这个\(left_i\)可以划分三个区间.

分别为\(left_i<l,l<=left_i<x,x=<left_i\)

第三个可以忽略 第二个考虑这是一个横跨分治中心的 且在分治区间中的.

每个i最多在一个分治区间中会发生这种情况.这个部分可以线段树暴力查区间最大值.所以复杂度为\(nlogn\)

考虑第三种情况 这种情况下.

继续分类讨论考虑左边大还是右边大.

右边大的话 左边暴力扫 分两种情况 一种是 \(i<x+C_x\)这种情况可以在暴力扫的时候更新.一种是\(i>=x+C_x\)这种可以等右边扫完线段树区间修改.

前者启发式合并\(nlogn\) 后者每个节点处至多进行一次复杂度\(nlogn\)

左边大的话 右边暴力扫 先对右边最近节点扫到的j即\(i-C_x\)这个位置之前的位置线段树求max.然后再扫不断进行更新.

前者每个节点至多进行一次 复杂度\(nlogn\) 后者还是启发式合并\(nlogn\)

还有就是当前区间分成三段的节点可以二分.

这样总复杂度为\(nlogn\) 可以通过本题.

非常妙的分类讨论.

code
//#include<bits/stdc++.h>
#include<iostream>
#include<cstdio>
#include<ctime>
#include<cctype>
#include<queue>
#include<deque>
#include<stack>
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<cctype>
#include<cstdlib>
#include<queue>
#include<deque>
#include<stack>
#include<vector>
#include<algorithm>
#include<utility>
#include<bitset>
#include<set>
#include<map>
#define ll long long
#define db double
#define INF 1000000000
#define inf 100000000000000000ll
#define ldb long double
#define pb push_back
#define put_(x) printf("%d ",x);
#define get(x) x=read()
#define gt(x) scanf("%d",&x)
#define gi(x) scanf("%lf",&x)
#define put(x) printf("%d\n",x)
#define putl(x) printf("%lld\n",x)
#define rep(p,n,i) for(RE int i=p;i<=n;++i)
#define go(x) for(int i=lin[x],tn=ver[i];i;tn=ver[i=nex[i]])
#define fep(n,p,i) for(RE int i=n;i>=p;--i)
#define vep(p,n,i) for(RE int i=p;i<n;++i)
#define pii pair<int,int>
#define mk make_pair
#define RE register
#define P 1000000007ll
#define gf(x) scanf("%lf",&x)
#define pf(x) ((x)*(x))
#define uint unsigned long long
#define ui unsigned
#define EPS 1e-10
#define sq sqrt
#define S second
#define F first
#define mod 1000000007
#define id(i,j) ((i+j)|(i!=j))
#define max(x,y) ((x)<(y)?y:x)
using namespace std;
char *fs,*ft,buf[1<<15];
inline char gc()
{
	return (fs==ft&&(ft=(fs=buf)+fread(buf,1,1<<15,stdin),fs==ft))?0:*fs++;
}
inline int read()
{
	RE int x=0,f=1;RE char ch=gc();
	while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=gc();}
	while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘;ch=gc();}
	return x*f;
}
const int MAXN=1000010;
int n,top,rt;
int c[MAXN],d[MAXN];
int le[MAXN],s[MAXN];
int ls[MAXN],rs[MAXN];
struct wy
{
	int v,cnt;
	wy(){v=-INF;cnt=0;}
	wy(int _v,int _cnt){v=_v;cnt=_cnt;}
	inline friend wy operator +(wy a,wy b)
	{
		if(a.v==-INF)return b;
		if(b.v==-INF)return a;
		if(a.v>b.v)return a;
		if(b.v>a.v)return b;
		return wy(a.v,a.cnt+b.cnt>=mod?a.cnt+b.cnt-mod:a.cnt+b.cnt);
	}
}cc;
inline wy calc(wy b)	
{
	return b.v==-INF?b:wy(b.v+1,b.cnt);
}
wy tag[MAXN<<1],mx[MAXN<<1],f[MAXN];
inline void build()
{
	top=0;
	rep(1,n,i)
	{
		int now=0;
		while(top&&c[s[top]]<=c[i])rs[s[top]]=now,now=s[top--];
		s[++top]=i;ls[i]=now;
	}
	rt=s[1];
	while(top>1)rs[s[top-1]]=s[top],--top;
}
inline void pushdown(int l,int r)
{
	int mid=(l+r)>>1;
	int ww=id(l,r),w1=id(l,mid),w2=id(mid+1,r);
	tag[w1]=tag[w1]+tag[ww];
	mx[w1]=mx[w1]+tag[ww];
	tag[w2]=tag[w2]+tag[ww];
	mx[w2]=mx[w2]+tag[ww];
	tag[ww]=f[n+1];
}
inline wy ask(int l,int r,int L,int R)
{
	if(L>R)return f[n+1];
	if(L<=l&&R>=r)return mx[id(l,r)];
	int mid=(l+r)>>1;
	if(tag[id(l,r)].v>=0)pushdown(l,r);
	if(R<=mid)return ask(l,mid,L,R);
	if(L>mid)return ask(mid+1,r,L,R);
	return ask(l,mid,L,R)+ask(mid+1,r,L,R);
}
inline void modify(int l,int r,int L,int R)
{
	if(L>R)return;
	if(L<=l&&R>=r)
	{
		int ww=id(l,r);
		mx[ww]=mx[ww]+cc;
		tag[ww]=tag[ww]+cc;
		return;
	}
	int mid=(l+r)>>1;
	if(tag[id(l,r)].v>=0)pushdown(l,r);
	if(R<=mid)modify(l,mid,L,R);
	else 
	{
		if(L>mid)modify(mid+1,r,L,R);
		else modify(l,mid,L,R),modify(mid+1,r,L,R);
	}
	mx[id(l,r)]=mx[id(l,mid)]+mx[id(mid+1,r)];
}
inline void solve(int l,int r,int x)
{
	if(l==r)
	{
		cc=f[l];modify(0,n,l,r);
		f[l]=ask(0,n,l,r);
		return;
	}
	solve(l,x-1,ls[x]);
	int p=lower_bound(le+x,le+r+1,l)-le;
	int q=lower_bound(le+x,le+r+1,x)-le;
	if(x-l<r-x+1)
	{
		wy ww;
		vep(l,x,i)
		{
			ww=ww+f[i];
			if(i+c[x]>=x&&i+c[x]<p)
				f[i+c[x]]=f[i+c[x]]+calc(ww);
		}
		cc=calc(ww);modify(0,n,x+c[x],p-1);
	}
	else
	{
		wy ww=ask(0,n,l,x-c[x]-1);
		vep(x,p,i)
		{
			if(i-c[x]>=l&&i-c[x]<x)ww=ww+f[i-c[x]];
			f[i]=f[i]+calc(ww);
		}
	}
	vep(p,q,i)f[i]=f[i]+calc(ask(0,n,le[i],min(x-1,i-c[x])));
	solve(x,r,rs[x]);
}
int main()
{
	//freopen("1.in","r",stdin);
	get(n);
	rep(1,n,i)get(c[i]),get(d[i]);
	int h=1,t=0;
	rep(1,n,i)
	{
		le[i]=le[i-1];
		while(h<=t&&d[i]<=d[s[t]])--t;
		s[++t]=i;
		while(i-le[i]>d[s[h]])
		{
			++le[i];
			if(le[i]>=s[h])++h;
		}
		//cout<<le[i]<<endl;
	}
	build();f[0]=wy(0,1);
	solve(0,n,rt);
	if(f[n].v==-INF)puts("NIE");
	else printf("%d %d\n",f[n].v,f[n].cnt);
	return 0;
}

P5979 [PA2014]Druzyny dp 分治 线段树 分类讨论 启发式合并

标签:register   mes   end   怎么   freopen   std   details   its   lib   

原文地址:https://www.cnblogs.com/chdy/p/13393782.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!