1
树的直径的两种求法
看到这个模板想到暑假集训的第一课bwd学长讲的树型dp,当时啥都不会,只是跑去抄了一个两次dfs求树的直径的方法就跑了,现在回想当时bwd学长也讲了树型dp求树的直径只是自己没当回事...
希望自己能好好利用这个寒假~ 队友南京拿了铜牌 棒棒哒 替他们高兴 我也抓紧时间提高一下自己
不然下学期教练都够呛让我继续打了....加油吧
方法一: 两边dfs
这里是利用了图论的简单知识,从树上任意一个点出发离他最远的一个点一定是直径的某个端点
代码如下:这里是这道hdu典例的代码Problem - 2196 (hdu.edu.cn)
#include<iostream>
#include<cstring>
using namespace std;
const int N=2e5+10;
int e[N],ne[N],w[N],h[N],idx;
int n;
int f1[N],f2[N];
int root,len;
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void init()
{
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
int b,c;
scanf("%d%d",&b,&c);
add(i,b,c),add(b,i,c);
}
}
void dfs(int u,int father,int f[],int x)
{
f[u]=x;
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs(j,u,f,x+w[i]);
}
if(x>len)
{
len=x;
root=u;
}
}
void solve()
{
//找直径的一个端点
len=-1;
dfs(1,-1,f1,0);
//找另一个端点 并更新f1
len=-1;
dfs(root,-1,f1,0);
//更新f2
len=-1;
dfs(root,-1,f2,0);
for(int i=1;i<=n;i++)printf("%d\n",max(f1[i],f2[i]));
}
int main()
{
while(cin>>n)
{
init();
solve();
}
}
2.树型dp求法 维护以每个点为根节点的子树种的最长路径和次长路径
注意法一是不能维护负权值的而法二可以
代码如下:
#include<iostream>
#include<cstring>
using namespace std;
const int N=2e5+10;
int e[N],ne[N],w[N],h[N],idx;
int n;
int f1[N],f2[N];
int root,len;
int res=-0x3f3f3f3f;
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void init()
{
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
add(a,b,c),add(b,a,c);
}
}
void dfs(int u,int father)
{
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs(j,u);
if(f1[u]<=f1[j]+w[i])f2[u]=f1[u],f1[u]=f1[j]+w[i];//相等也要换 使得次长也变大
else if(f2[u]<f1[j]+w[i])f2[u]=f1[j]+w[i];
}
res=max(res,f1[u]+f2[u]);
}
void solve()
{
dfs(1,-1);
cout<<res;
}
int main()
{
while(cin>>n)
{
init();
solve();
}
}
2
树型dp简单模板题2:
一个节点距离其他点的最远的点 要么是在以自己为根的子树的距离最远点,要么是当前考虑的点向上出发到其父节点,再从父节点出发到最远的点
需要注意的是 由树的拓扑结构可知 某个节点父节点对其贡献是不能计算父节点走”回头路“的距离的,也就是从父节点出发的那条可贡献距离不能再经过该点了
下面是代码:
#include<bits/stdc++.h>
using namespace std;
const int N=2e5+10;
int e[N],ne[N],w[N],h[N],idx;
int n;
int f1[N],f2[N],s1[N],s2[N];
int up[N];
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void dfs1(int u,int father)
{
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs1(j,u);
if(f1[u]<=f1[j]+w[i])
{
f2[u]=f1[u],f1[u]=f1[j]+w[i];
s2[u]=s1[u],s1[u]=j;
}
else if(f2[u]<f1[j]+w[i])
{
f2[u]=f1[j]+w[i];
s2[u]=j;
}
}
}
void dfs2(int u,int father)
{
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
if(s1[u]==j)up[j]=max(up[u],f2[u])+w[i];
else up[j]=w[i]+max(up[u],f1[u]);
dfs2(j,u);
}
}
void init()
{
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
add(a,b,c),add(b,a,c);
}
}
void solve()
{
dfs1(1,-1);//维护以当前节点为根节点的子树的最长与次长距离
dfs2(1,-1);//注意和dfs1的不同之处在于从祖宗向孩子更新
int res=0x3f3f3f3f;
for(int i=1;i<=n;i++)res=min(res,max(f1[i],up[i]));
cout<<res;
}
int main()
{
cin>>n;
init();
solve();
}
了解了基础的模板之后 解水题就比较简单了
下面分享几道树型dp 的简单题
简单数论建图+树的直径
#include<bits/stdc++.h>
using namespace std;
const int N=5e5+10,M=1e6+10;
int n;
bool st[N];
int primes[N];
int num;
int h[N],ne[M],w[M],idx,e[N];
int f1[N],f2[N];
int res;
void init()
{
for(int i=2;i<N;i++)
{
if(!st[i])primes[num++]=i;
for(int j=0;primes[j]<=N/i;j++)
{
st[primes[j]*i]=1;
if(i%primes[j]==0)break;
}
}
}
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void dfs(int u,int father)
{
f1[u]=f2[u]=0;
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs(j,u);
if(f1[u]<=f1[j]+w[i])f2[u]=f1[u],f1[u]=w[i]+f1[j];
else if(f2[u]<f1[j]+w[i])f2[u]=f1[j]+w[i];
}
res=max(res,f1[u]+f2[u]);
}
int main()
{
cin>>n;
init();
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
if(!st[i])
{
add(i,1,1),add(1,i,1);
//printf("xixi%d\n",i);
}
else{
int sum=1;
for(int j=2;j<=i/j;j++)
{
if(i%j==0)
{
sum+=j;
if(j!=i/j)sum+=i/j;
}
if(sum>=i)break;
}
if(sum<i)
{
add(sum,i,1),add(i,sum,1);
//printf("gaga%d %d\n",i,sum);
}
}
}
//树的直径
dfs(1,-1);
cout<<res;
}
#include<iostream>
#include<cstring>
using namespace std;
const int N=2e6+10;
int dist1[N],dist2[N];
int fa,len;
int n;
int ne[N],w[N],e[N],idx,h[N];
void dfs(int u,int father,int dist[],int x)
{
dist[u]=x;
for(int i=h[u];~i;i=ne[i]){
int j=e[i];
if(father==j)continue;
dfs(j,u,dist,x+w[i]);
}
if(x>len){
len=x;
fa=u;
}
}
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
int main()
{
memset(h,-1,sizeof h);
cin>>n;
for(int i=2;i<=n;i++){
int a,b,c;
scanf("%d%d%d",&a,&b,&c);
add(a,b,c),add(b,a,c);
}
len=-1;
dfs(1,-1,dist1,0);
len=-1;
dfs(fa,-1,dist1,0);
len=-1;
dfs(fa,-1,dist2,0);
int res=0;
for(int i=1;i<=n;i++)res=max(res,max(dist1[i],dist2[i]));
long long ans=0;
int st=10;
for(int i=1;i<=res;i++)
{
int tem=st+i;
ans+=tem;
}
cout<<ans;
}
树的直径裸题(蓝桥)
判断点在某一条直径上吗? 最大值+次大值==直径长度(树型dp模板2裸题)
#include<bits/stdc++.h>
using namespace std;
const int N=4e5+10;
int e[N],ne[N],w[N],h[N],idx;
int n;
int f1[N],f2[N],s1[N],s2[N];
int up[N];
int res;
void add(int a,int b,int c)
{
e[idx]=b,ne[idx]=h[a],w[idx]=c,h[a]=idx++;
}
void dfs1(int u,int father)
{
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs1(j,u);
if(f1[u]<=f1[j]+w[i])
{
f2[u]=f1[u],f1[u]=f1[j]+w[i];
s2[u]=s1[u],s1[u]=j;
}
else if(f2[u]<f1[j]+w[i])
{
f2[u]=f1[j]+w[i];
s2[u]=j;
}
}
res=max(res,f1[u]+f2[u]);
}
void dfs2(int u,int father)
{
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
if(s1[u]==j)up[j]=max(up[u],f2[u])+w[i];
else up[j]=w[i]+max(up[u],f1[u]);
dfs2(j,u);
}
}
void init()
{
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
int a,b,c;
scanf("%d%d",&a,&b);
add(a,b,1),add(b,a,1);
}
}
void solve()
{
dfs1(0,-1);//维护以当前节点为根节点的子树的最长与次长距离
dfs2(0,-1);//注意和dfs1的不同之处在于从祖宗向孩子更新
for(int i=0;i<n;i++)//最大值加次大值等于直径
if(res==(max(f1[i],up[i])+f1[i]+up[i]+f2[i]-max(f1[i],up[i])-min(f1[i],min(up[i],f2[i]))))
printf("%d\n",i);
}
int main()
{
cin>>n;
init();
solve();
}
没有上司的舞会285. 没有上司的舞会 - AcWing题库
树型dp套简单状态机模型(基础课经典题)
#include<iostream>
#include<cstring>
using namespace std;
const int N=1e5+10;
int f[N][2];
int e[N],ne[N],w[N],idx,h[N];
int root=1;
int n;
int father[N];
void add(int a,int b)
{
e[idx]=b,ne[idx]=h[a],h[a]=idx++;
}
void dfs(int u,int father)
{
f[u][1]=w[u];
for(int i=h[u];~i;i=ne[i])
{
int j=e[i];
if(j==father)continue;
dfs(j,u);
f[u][0]+=max(f[j][0],f[j][1]);
f[u][1]+=f[j][0];
}
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++)cin>>w[i];
memset(h,-1,sizeof h);
for(int i=2;i<=n;i++)
{
int a,b;
cin>>a>>b;
add(a,b),add(b,a);
father[a]=b;
}
while(father[root])root++;
dfs(root,-1);
printf("%d",max(f[root][0],f[root][1]));
}
小结:
树型dp个人认为和图论联系紧密一些 常常结合树的直径来考察 比状压好像一点 因为偏向图论的缘故属于dp中码量多一些的