必备知识

​ 树链剖分,最大权独立集(即没有上司的舞会(树上DP)),矩阵乘法;

D-DP

模版简述

​ 模板

​ 关于动态DP,其实是关于一类动态修改点权的问题,但是很难去处理;

​ 我们平常的DP经常是离线DP,而当在线时,就会出现事故;

​ D-DP是关于求最大权独立集的,支持动态修改点值,其思想即DP的通项递推式改为矩阵乘的形式进行递推,而树上问题就可以使用树链剖分提前处理出区间矩阵乘;

​ 大家肯定都会 \(O(nm)\) 做法,即暴力修改点值进行DP.

​ 求最大权独立集的DP式子为 \(f[x][1]=\sum f[y][0] , f[x][0]=\sum max(f[y][1],f[y][0])\) ,如果我们将其改成矩阵乘的形式,就变成了下面这样;

\[[^{f[y][0]}_{f[y][1]}] * [^{0,0}_{v[x],-inf}]=[^{f[x][0]}_{f[x][1]}]
\]

​ 这个式子并不完整,这只是从一个点的矩阵乘,这个式子只是类似于矩阵乘,我们可以重新定义一下矩阵乘的定义,即将加法改为 \(max\) ,将乘法改为加法,那么就可以通项递推了;

​ 但是这个式子只限于从一个点递推过来,那么其他点怎么办?

​ 没有办法...但是这是一个思路,这个式子用线段树维护,可以让我们可以快速处理出链上信息.

​ 既然在树上,我们不如将重链看成这个链,而将链上的轻儿子的信息提前处理出来,在重链上的点的信息通过进行矩阵速推.

\[[^{f[y][0]}_{f[y][1]}]*[^{f'[x][0],f'[x][0]}_{f'[x][1],-inf}]=[^{f[x][0]}_{f[x][1]}]
\]

​ (注: \(f'[x][0]\) 是点 \(x\) 只处理轻儿子和自己信息时的值)

​ 那么如何修改?

​ 同一条重链上点的修改对于这条链上的点是没有影响的,因为重链每个点只保存其轻儿子的信息,有影响的是对于这条链 \(top\) 的父亲节点,那么我们的思路就出来了,修改点时,修改链 \(top\) 父亲的矩阵信息即可;

​ 总结一下思路: 首先将DP式子化为通项矩阵乘的形式,先处理出每个节点轻儿子的信息,在树上用树链剖分维护链上矩阵乘信息,修改时,每次修改链首父亲的矩阵值,修改至树根重链.修改细节即修改点信息,之后从链尾到链首矩阵乘,与之前的值作差修改.最后求答案时直接处理链上信息即可;

code
#include<bits/stdc++.h>

#define le(p) p<<1

#define re(p) p<<1|1

#define ll long long

using namespace std;

const int maxn=1e5+7,inf=1e9;

int n,m,head[maxn],son[maxn],top[maxn],bot[maxn],fl[maxn];

int fa[maxn],id[maxn],cent,tot,a[maxn],sz[maxn],f[maxn][2];

struct node{

	int next,to;

}edge[maxn<<1];

struct matrix{

	ll g[2][2];

	matrix operator *(const matrix x) const{

		matrix ans;

		ans.g[0][0]=max(g[0][0]+x.g[0][0],g[0][1]+x.g[1][0]);

		ans.g[0][1]=max(g[0][0]+x.g[0][1],g[0][1]+x.g[1][1]);

		ans.g[1][0]=max(g[1][0]+x.g[0][0],g[1][1]+x.g[1][0]);

		ans.g[1][1]=max(g[1][0]+x.g[0][1],g[1][1]+x.g[1][1]);

		return ans;

	}//矩阵乘

}tr[maxn<<2];

matrix ori[maxn];

template<typename type_of_scan>

inline void scan(type_of_scan &x){

	type_of_scan f=1;x=0;char s=getchar();

	while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}

	while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}

	x*=f;

}

template<typename top,typename... tops>

inline void scan(top &x,tops&... X){

	scan(x),scan(X...);

}

inline void add(int u,int v){

	edge[++cent]=(node){head[u],v};head[u]=cent;

	edge[++cent]=(node){head[v],u};head[v]=cent;

}

//树链剖分预处理

void dfs1(int x){

	sz[x]=1;

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]) continue;

		fa[y]=x,dfs1(y);

		sz[x]+=sz[y];

		if(sz[y]>sz[son[x]]) son[x]=y;

	}

}

void dfs2(int x,int tp){

	id[x]=++tot,fl[tot]=x;top[x]=tp;

	if(son[x]) dfs2(son[x],tp),bot[x]=bot[son[x]];

	else return bot[x]=x,void();

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]||y==son[x]) continue;

		dfs2(y,y);

	}

}

//树上DP

void dfs3(int x){

	f[x][0]=0;f[x][1]=a[x];

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]) continue;

		dfs3(y);

		f[x][0]+=max(f[y][0],f[y][1]);

		f[x][1]+=f[y][0];

	}

}

void push_up(int p){

	tr[p]=tr[le(p)]*tr[re(p)];

}

//建树时每个点只保存自己和轻儿子的信息

void build(int l,int r,int p){

	if(l==r){

		tr[p].g[0][0]=tr[p].g[0][1]=0;tr[p].g[1][1]=-inf;

		tr[p].g[1][0]=a[fl[l]];int x=fl[l];

		for(int i=head[x];i;i=edge[i].next){

			int y=edge[i].to;

			if(y==fa[x]||y==son[x]) continue;

			tr[p].g[0][0]+=max(f[y][0],f[y][1]);

			tr[p].g[1][0]+=f[y][0];

		}

		tr[p].g[0][1]=tr[p].g[0][0];

		ori[fl[l]]=tr[p];

		return ;

	}

	int mid=l+r>>1;

	build(l,mid,le(p));build(mid+1,r,re(p));

	push_up(p);

}

void add(int x,int l,int r,int p){

	if(l==r) return tr[p]=ori[fl[l]],void();

	int mid=l+r>>1;

	if(x<=mid) add(x,l,mid,le(p));

	else add(x,mid+1,r,re(p));

	push_up(p);

}

matrix query(int nl,int nr,int l,int r,int p){

	if(nl<=l&&r<=nr) return tr[p];

	int mid=l+r>>1;

	if(nl<=mid&&mid<nr) return query(nl,nr,l,mid,le(p))*query(nl,nr,mid+1,r,re(p));

	else if(nl<=mid) return query(nl,nr,l,mid,le(p));

	else return query(nl,nr,mid+1,r,re(p));

}

matrix query(int x){

	return query(id[top[x]],id[bot[x]],1,n,1);

}

matrix solve(int x,int k){

	ori[x].g[1][0]+=k-a[x];a[x]=k;

	while(x){

		matrix old,news;

		old=query(top[x]);//没修改的

		add(id[x],1,n,1);//change

		news=query(top[x]);//修改过的

		x=fa[top[x]];if(!x) break;

		ori[x].g[0][0]+=max(news.g[0][0],news.g[1][0])-max(old.g[0][0],old.g[1][0]);

		ori[x].g[0][1]=ori[x].g[0][0];

		ori[x].g[1][0]+=news.g[0][0]-old.g[0][0];

	}

	return query(1);

}

int main(){

	scan(n,m);

	for(int i=1;i<=n;i++) scan(a[i]);

	for(int i=1,u,v;i<=n-1;i++)

		scan(u),scan(v),add(u,v);

	dfs1(1),dfs2(1,1),dfs3(1);build(1,n,1);

	while(m--){

		int x,k;

		scan(x,k);

		matrix ans=solve(x,k);

		printf("%d\n",max(ans.g[0][0],ans.g[1][0]));

	}

	return 0;

}

例题

​ [NOIP2018]保卫王国;

​ 有两个点,强制选点或强制不选,求最小覆盖集;

​ 最小覆盖集=全集-最大权独立集;

​ 强制选点或不选可以通过将其改为 \(-inf\) 和 \(inf\) ,之后就是D-DP模板了;

code
#include<bits/stdc++.h>

#define ll long long

#define le(x) x<<1

#define re(x) x<<1|1

using namespace std;

const int maxn=1e5+7;

const ll inf=1e9;

int n,m,fa[maxn],id[maxn],fl[maxn],sz[maxn],head[maxn];

int son[maxn],top[maxn],cent,tot,bot[maxn];

ll f[maxn][2],a[maxn],sum,ans;

vector<int>to[maxn];

char type[2];

struct node{

	int next,to;

}edge[maxn<<1];

struct matrix{

	ll g[2][2];

	matrix operator *(const matrix x) const{

		matrix ans;

		ans.g[0][0]=max(g[0][0]+x.g[0][0],g[0][1]+x.g[1][0]);

		ans.g[0][1]=max(g[0][0]+x.g[0][1],g[0][1]+x.g[1][1]);

		ans.g[1][0]=max(g[1][0]+x.g[0][0],g[1][1]+x.g[1][0]);

		ans.g[1][1]=max(g[1][0]+x.g[0][1],g[1][1]+x.g[1][1]);

		return ans;

	}

}tr[maxn<<3];

matrix ori[maxn];

template<typename type_of_scan>

inline void scan(type_of_scan &x){

	type_of_scan f=1;x=0;char s=getchar();

	while(s<'0'||s>'9'){if(s=='-')f=-1;s=getchar();}

	while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}

	x*=f;

}

template<typename top,typename... tops>

inline void scan(top &x,tops&... X){

	scan(x),scan(X...);

}

inline void add(int u,int v){

	edge[++cent]=(node){head[u],v};head[u]=cent;

	edge[++cent]=(node){head[v],u};head[v]=cent;

	to[u].push_back(v),to[v].push_back(u);

}

void dfs1(int x){

	sz[x]=1;

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]) continue;

		fa[y]=x;dfs1(y);

		sz[x]+=sz[y];

		if(sz[y]>sz[son[x]]) son[x]=y;

	}

}

void dfs2(int x,int tp){

	id[x]=++tot,fl[tot]=x,top[x]=tp;

	if(son[x]) dfs2(son[x],tp),bot[x]=bot[son[x]];

	else return bot[x]=x,void();

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]||y==son[x]) continue;

		dfs2(y,y);

	}

}

void dfs3(int x){

	f[x][1]=a[x],f[x][0]=0;

	for(int i=head[x];i;i=edge[i].next){

		int y=edge[i].to;

		if(y==fa[x]) continue;

		dfs3(y);

		f[x][0]+=max(f[y][0],f[y][1]);

		f[x][1]+=f[y][0];

	}

}

void push_up(int p){

	tr[p]=tr[le(p)]*tr[re(p)];

}

void build(int l,int r,int p){

	if(l==r){

		tr[p].g[0][0]=tr[p].g[0][1]=0;int x=fl[l];

		tr[p].g[1][0]=a[fl[l]];tr[p].g[1][1]=-inf;

		for(int i=head[x];i;i=edge[i].next){

			int y=edge[i].to;

			if(y==fa[x]||y==son[x]) continue;

			tr[p].g[0][0]+=max(f[y][0],f[y][1]);

			tr[p].g[1][0]+=f[y][0];

		}

		tr[p].g[0][1]=tr[p].g[0][0];

		ori[fl[l]]=tr[p];

		return ;

	}

	int mid=l+r>>1;

	build(l,mid,le(p)),build(mid+1,r,re(p));

	push_up(p);

}

void add(int x,int l,int r,int p){

	if(l==r) return tr[p]=ori[fl[l]],void();

	int mid=l+r>>1;

	if(x<=mid) add(x,l,mid,le(p));

	else add(x,mid+1,r,re(p));

	push_up(p);

}

matrix query(int nl,int nr,int l,int r,int p){

	if(nl<=l&&r<=nr) return tr[p];

	int mid=l+r>>1;

	if(nl<=mid&&mid<nr) return query(nl,nr,l,mid,le(p))*query(nl,nr,mid+1,r,re(p));

	else if(nl<=mid) return query(nl,nr,l,mid,le(p));

	else return query(nl,nr,mid+1,r,re(p));

}

matrix query(int x){

	return query(id[top[x]],id[bot[x]],1,n,1);

}

void modify(int x,ll k){

	ori[x].g[1][0]+=k-a[x];a[x]=k;

	while(x){

		matrix old,news;

		old=query(top[x]);add(id[x],1,n,1);

		news=query(top[x]);

		x=fa[top[x]];if(!x) break;

		ori[x].g[0][0]+=max(news.g[0][0],news.g[1][0])-max(old.g[0][0],old.g[1][0]);

		ori[x].g[0][1]=ori[x].g[0][0];

		ori[x].g[1][0]+=news.g[0][0]-old.g[0][0];

	}

	return ;

}

int main(){

	scan(n,m);cin>>type;

	for(int i=1;i<=n;i++) scan(a[i]),sum+=a[i];

	for(int i=1,u,v;i<=n-1;i++)

		scan(u),scan(v),add(u,v);

	for(int i=1;i<=n;i++) sort(to[i].begin(),to[i].end());

	dfs1(1),dfs2(1,1),dfs3(1),build(1,n,1);

	while(m--){

		int f1,f2,g1,g2,t1,t2;

		scan(f1,f2,g1,g2);

		if(f2==g2&&!f2){

			vector<int>::iterator pos=lower_bound(to[f1].begin(),to[f1].end(),g1);

			if(pos!=to[f1].end()&&(*pos)==g1){

				puts("-1");continue;

			}

		}

		t1=a[f1],t2=a[g1];

		modify(f1,f2?-inf:inf);

		modify(g1,g2?-inf:inf);

		matrix ol=query(1);

		ans=max(ol.g[0][0],ol.g[1][0]);

		if(!f2) ans=ans-inf+t1;

		if(!g2) ans=ans-inf+t2;

		printf("%lld\n",sum-ans);

		modify(f1,t1),modify(g1,t2);

	}

	return 0;

}

​ 当然还有倍增做法,但这里不再讲述,D-DP中一个重要思想就是讲DP式子变成矩阵乘的形式,用数据结构进行维护,其高效并支持拓展;

D-DP的更多相关文章

  1. BZOJ 1911: [Apio2010]特别行动队 [斜率优化DP]

    1911: [Apio2010]特别行动队 Time Limit: 4 Sec  Memory Limit: 64 MBSubmit: 4142  Solved: 1964[Submit][Statu ...

  2. 2013 Asia Changsha Regional Contest---Josephina and RPG(DP)

    题目链接 http://acm.hdu.edu.cn/showproblem.php?pid=4800 Problem Description A role-playing game (RPG and ...

  3. AEAI DP V3.7.0 发布，开源综合应用开发平台

    1  升级说明 AEAI DP 3.7版本是AEAI DP一个里程碑版本,基于JDK1.7开发,在本版本中新增支持Rest服务开发机制(默认支持WebService服务开发机制),且支持WS服务.RS ...

  4. AEAI DP V3.6.0 升级说明，开源综合应用开发平台

    AEAI DP综合应用开发平台是一款扩展开发工具,专门用于开发MIS类的Java Web应用,本次发版的AEAI DP_v3.6.0版本为AEAI DP _v3.5.0版本的升级版本,该产品现已开源并 ...

  5. BZOJ 1597: [Usaco2008 Mar]土地购买 [斜率优化DP]

    1597: [Usaco2008 Mar]土地购买 Time Limit: 10 Sec  Memory Limit: 162 MBSubmit: 4026  Solved: 1473[Submit] ...

  6. [斜率优化DP]【学习笔记】【更新中】

    参考资料: 1.元旦集训的课件已经很好了 http://files.cnblogs.com/files/candy99/dp.pdf 2.http://www.cnblogs.com/MashiroS ...

  7. BZOJ 1010: [HNOI2008]玩具装箱toy [DP 斜率优化]

    1010: [HNOI2008]玩具装箱toy Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 9812  Solved: 3978[Submit][St ...

  8. px、dp和sp，这些单位有什么区别？

    DP 这个是最常用但也最难理解的尺寸单位.它与“像素密度”密切相关,所以 首先我们解释一下什么是像素密度.假设有一部手机,屏幕的物理尺寸为1.5英寸x2英寸,屏幕分辨率为240x320,则我们可以计算 ...

  9. android px转换为dip/dp

    /** * 根据手机的分辨率从 dp 的单位 转成为 px(像素) */ public int dipTopx(Context context, float dpValue) { final floa ...

  10. POJ 3254. Corn Fields 状态压缩DP (入门级)

    Corn Fields Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 9806   Accepted: 5185 Descr ...

随机推荐

  1. SQL中的主键，候选键，外键，主码，外码

    1.码=超键:能够唯一标识一条记录的属性或属性集. 标识性:一个数据表的所有记录都具有不同的超键 非空性:不能为空 有些时候也把码称作"键" 2.候选键=候选码:能够唯一标识一条记 ...

  2. 【Jboss】应用中缺少宋体怎么办

    环境jboss4.2.2 系统CentOS7.2 1.新搭建的环境,但是没有字符集,在windows上的电脑上复制了一份宋体,打成zip包 将zip包上传到服务器中,解压 2.在/usr/share/ ...

  3. Unsafe Fileupload - Pikachu

    概述: 文件上传功能在web应用系统很常见,比如很多网站注册的时候需要上传头像.上传附件等等.当用户点击上传按钮后,后台会对上传的文件进行判断 比如是否是指定的类型.后缀名.大小等等,然后将其按照设计 ...

  4. 07--Docker安装Redis

    1.拉取redis:3.2 docker pull redis:3.2 2.创建redis容器 docker run -p 6379:6379 -v /zhengcj/myredis/data:/da ...

  5. springmvc 字符串转日期格式

    http://www.mamicode.com/info-detail-2485490.html

  6. 性能测试工具locust简单应用

    简介 Locust是一种易于使用的分布式用户负载测试工具.可用于对网站(或系统)负载测试,并依据响应数据计算出系统支持的并发用户数. 安装及调试(以下操作在windows环境下进行) Locust基于 ...

  7. 记一次 RocketMQ broker 因内存不足导致的启动失败

    原创:西狩 编写日期 / 修订日期:2020-01-12 / 2020-01-12 版权声明:本文为博主原创文章,遵循 CC BY-SA-4.0 版权协议,转载请附上原文出处链接和本声明. 背景 该小 ...

  8. kubernetes 核心技术-Controller 控制器

    一.什么是Controller? Controller是在集群上管理和运行容器的对象,Controller是实际存在的,Pod是抽象的,主要创建管理pod 二.Pod和Controller的关系 Po ...

  9. 处理 K8S Orphaned pod found - but volume paths are still present on disk 孤儿pod

    问题概述 查看kubelet或/var/log/messages日志一直包错,发现是孤儿pod,是由于其pod被删除后存储路径还保存在磁盘. 报错如下 [root@node5 ~]# journalc ...

  10. 3、wait和waitpid

    1. 函数介绍 wait函数:调用该函数使进程阻塞,直到任意一个子进程结束,或者该进程接收到了一个信号为止,如果该进程没有子进程或该进程的子进程已经结束,wait函数立即返回. waitpid函数:与 ...