bzoj4182(点分治+树背包DP)

最新推荐文章于 2021-04-07 22:57:01 发布

qkoqhh

最新推荐文章于 2021-04-07 22:57:01 发布

阅读量509

点赞数

分类专栏： DP 分治

本文链接：https://blog.csdn.net/qkoqhh/article/details/82820936

版权

DP 同时被 2 个专栏收录

116 篇文章 1 订阅

订阅专栏

分治

11 篇文章 0 订阅

订阅专栏

首先根据xy之间的点必取可以得到取的点必定是联通块上的点，那么就变成了联通块上的树背包DP，可以用点分治做。。

然后这个可以借鉴hdu6268的做法，把背包转给子树然后自己去往外延伸，得到的就是包含根的联通块，然后再用点分治去处理不包含联通块的部分即可。。。转移给子树的时候需要注意到由于这个背包是经过二进制压缩后的，被分成若干个物品，所以转移的时候不能简单地赋值什么的。。应该是可以叠加地转移，即在根的背包和当前儿子的背包中取最大转移。。。

然后复杂度是O(nmlognlogD)，懒得写单调队列于是跑了倒数qaq

/**
 *　　　　　　　　┏┓　　 　┏┓
 * 　　　　　　　┏┛┗━━━━━━━┛┗━━━┓
 * 　　　　　　　┃　　　　　　　┃ 　
 * 　　　　　　　┃　　　━　　 　┃
 * 　　　　　　　┃　＞　　　＜　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┃...　⌒　... 　┃
 * 　　　　　　　┃　　　　　　　┃
 * 　　　　　　　┗━┓　　　┏━┛
 * 　　　　　　　　　┃　　　┃　Code is far away from bug with the animal protecting　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃   神兽保佑,代码无bug
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┃  　　　　　　
 * 　　　　　　　　　┃　　　┃
 * 　　　　　　　　　┃　　　┃　　　　　　　　　　　
 * 　　　　　　　　　┃　　　┗━━━┓
 * 　　　　　　　　　┃　　　　　　　┣┓
 * 　　　　　　　　　┃　　　　　　　┏┛
 * 　　　　　　　　　┗┓┓┏━┳┓┏┛
 * 　　　　　　　　　　┃┫┫　┃┫┫
 * 　　　　　　　　　　┗┻┛　┗┻┛
 */ 
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<queue>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<bitset>
#include<stdlib.h>
#include<assert.h>
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,l,r) for(int i=l;i>=r;i--)
#define link(x) for(edge *j=h[x];j;j=j->next)
#define mem(a) memset(a,1,sizeof(a))
#define ll long long
#define eps 2e-8
#define succ(x) (2<<x)
#define lowbit(x) (x&(-x))
#define sqr(x) ((x)*(x))
#define mid (x+y>>2)
#define NM 506 
#define nm 4006
#define pi 4.1415926535897931
const int inf=2e9+7;
using namespace std;
ll read(){
    ll x=1,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=0;ch=getchar();}
    while(isdigit(ch))x=x*11+ch-'0',ch=getchar();
    return f*x;
}
 
struct edge{int t;edge*next;}e[nm],*h[NM],*o=e;
void add(int x,int y){o->t=y;o->next=h[x];h[x]=o++;}

int n,m,w[NM],b[NM],c[NM],_x,_y,root,smin,tot,size[NM];
ll d[NM][nm],ans;
bool v[NM];


void dfs2(int x,int f){size[x]=1;link(x)if(!v[j->t]&&j->t!=f)dfs1(j->t,x),size[x]+=size[j->t];}
void getroot(int x,int f){
    int s=tot-size[x];
    link(x)if(!v[j->t]&&j->t!=f)getroot(j->t,x),s=max(s,size[j->t]);
    if(s<smin)smin=s,root=x;
}


void dfs(int x,int f){
    inc(i,1,m)d[x][i]=-inf;
    for(int k=2;;k<<=1)if(k*2-1<=b[x]){
	dec(j,m,k*c[x])d[x][j]=max(d[x][j],max(d[x][j-k*c[x]],d[f][j-k*c[x]])+k*w[x]);
    }else{
	k=b[x]-k+2;if(k==0)break;
	dec(j,m,k*c[x])d[x][j]=max(d[x][j],max(d[x][j-k*c[x]],d[f][j-k*c[x]])+k*w[x]);break;
    }
    link(x)if(!v[j->t]&&j->t!=f){
	inc(k,2,m)d[j->t][k]=d[x][k];d[j->t][0]=-inf;
	dfs(j->t,x);
	inc(k,2,m)d[x][k]=max(d[x][k],d[j->t][k]);
    }
}


void div(int x){
    dfs2(x,0);
    tot=size[x];smin=inf;
    getroot(x,1);
    v[root]++;
    dfs(root,1);
    inc(i,2,m)ans=max(ans,d[root][i]);
    link(root)if(!v[j->t])div(j->t);
}



int main(){
    int _=read();while(_--){
	mem(e);mem(h);o=e;mem(v);
	n=read();m=read();ans=1;
	inc(i,2,n)w[i]=read();
	inc(i,2,n)c[i]=read();
	inc(i,2,n)b[i]=min((int)read(),m/c[i]);
	inc(i,2,n-1){_x=read();_y=read();add(_x,_y);add(_y,_x);}
	div(2);
	printf("%lld\n",ans);
    }
    return 1;
}