P5979 [PA2014]Druzyny dp 分治 线段树 分类讨论 启发式合并
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\) 可以通过本题.
非常妙的分类讨论.
需要注意的是 这题卡空间 不动态开点 利用id函数即可使得儿子之间有序标号.空间为(2n).
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 分治 线段树 分类讨论 启发式合并的更多相关文章
- BZOJ-1067 降雨量 线段树+分类讨论
这道B题,刚的不行,各种碎点及其容易忽略,受不鸟了直接 1067: [SCOI2007]降雨量 Time Limit: 1 Sec Memory Limit: 162 MB Submit: 2859 ...
- bzoj 3779 重组病毒 好题 LCT+dfn序+线段树分类讨论
题目大意 1.将x到当前根路径上的所有点染成一种新的颜色: 2.将x到当前根路径上的所有点染成一种新的颜色,并且把这个点设为新的根: 3.查询以x为根的子树中所有点权值的平均值. 分析 原题codec ...
- SPOJ 2916 Can you answer these queries V(线段树-分类讨论)
题目链接:http://www.spoj.com/problems/GSS5/ 题意:给出一个数列.每次查询最大子段和Sum[i,j],其中i和j满足x1<=i<=y1,x2<=j& ...
- 美团CodeM初赛B轮 合并字符串的价值 (线段树,分类讨论)
输入两个字符串a和b,合并成一个串c,属于a或b的字符在c中顺序保持不变.如"ACG"和"UT"可以被组合成"AUCTG"或"AC ...
- 2020牛客暑假多校训练营 第二场 H Happy Triangle set 线段树 分类讨论
LINK:Happy Triangle 这道题很容易. 容易想到 a+b<x a<x<b x<a<b 其中等于的情况在第一个和第三个之中判一下即可. 前面两个容易想到se ...
- UVALive 7148 LRIP【树分治+线段树】
题意就是要求一棵树上的最长不下降序列,同时不下降序列的最小值与最大值不超过D. 做法是树分治+线段树,假设树根是x,y是其当前需要处理的子树,对于子树y,需要处理出两个数组MN,MX,MN[i]表示以 ...
- 【loj6145】「2017 山东三轮集训 Day7」Easy 动态点分治+线段树
题目描述 给你一棵 $n$ 个点的树,边有边权.$m$ 次询问,每次给出 $l$ .$r$ .$x$ ,求 $\text{Min}_{i=l}^r\text{dis}(i,x)$ . $n,m\le ...
- 【BZOJ4372】烁烁的游戏 动态树分治+线段树
[BZOJ4372]烁烁的游戏 Description 背景:烁烁很喜欢爬树,这吓坏了树上的皮皮鼠.题意:给定一颗n个节点的树,边权均为1,初始树上没有皮皮鼠.烁烁他每次会跳到一个节点u,把周围与他距 ...
- 【bzoj4372】烁烁的游戏 动态点分治+线段树
题目描述 给一颗n个节点的树,边权均为1,初始点权均为0,m次操作:Q x:询问x的点权.M x d w:将树上与节点x距离不超过d的节点的点权均加上w. 输入 第一行两个正整数:n,m接下来的n-1 ...
随机推荐
- 揭秘JAVA JVM内幕
在之前的文章 一步步解析java执行内幕 中,比较详细分析了java代码是如何一步一步在jvm中执行的,然而设计的的jvm核心技术点,并未做深入分析,本篇文章将重点分析jvm,涉及到的内容包括jvm内 ...
- sql语句-CASE WHEN用法详解
当我们需要从数据源上 直接判断数据显示代表的含义的时候 ,就可以在SQL语句中使用 Case When这个函数了. Case具有两种格式.简单Case函数和Case搜索函数. 第一种 格式 : 简单C ...
- 大致掌握django的功能
目录 静态文件配置 request对象方法初识 pycharm链接数据库(mysql) django链接数据库(mysql) django orm 字段的增删查改 数据的增删查改 数据的查,改,删 d ...
- 我终于弄懂了Python的装饰器(四)
此系列文档: 1. 我终于弄懂了Python的装饰器(一) 2. 我终于弄懂了Python的装饰器(二) 3. 我终于弄懂了Python的装饰器(三) 4. 我终于弄懂了Python的装饰器(四) 四 ...
- Centos7 python 安装 Ignoring ensurepip failure: pip 9.0.1 requires SSL/TLS
安装python时出现Ignoring ensurepip failure:pip required SSL/TLS 因为没有安装OpenSSL 使用yum install openssl-devel
- 关于springboot中过滤器和拦截器
在解决跨域问题中,发现拦截器和过滤器用得不是熟练.就参考了下一下两个作者的文档.希望大家也可以汲取精华 文档1 https://blog.csdn.net/moonpure/article/det ...
- 【RPA Starter第二课】Introduction to the UiPath Enterprise Platform UiPath企业平台简介
Introduction to the UiPath Enterprise Platform UiPath 企业平台简介 课程目标: 了解UiPath实现RPA的步骤 描述每个UiPath解决方案的关 ...
- linq介绍及工作中应用两例——左联与内联,linq循环方法
目录 1 linq介绍 1.1 linq产生背景 1.2 linq使用范围 1.3 linq核心程序集 1.4 linq架构图 1.5 linq使用形式对比 1.5.1 linq To Objects ...
- Zookeeper ----- 系统模型
数据模型 Zookeeper的数据模型与文件系统非常相似,唯一不同的它的每个节点(ZNode)都可以存放数据,无论父节点还是子节点. 事务ID 即前面提到的ZXID.对每个事务请求,Zookeeper ...
- Flutter 实现酷炫的3D效果
老孟导读:此文讲解3个酷炫的3D动画效果. 下面是要实现的效果: Flutter 中3D效果是通过 Transform 组件实现的,没有变换效果的实现: class TransformDemo ext ...