简单树型dp介绍与常用板子+经典例题

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求法 维护以每个点为根节点的子树种的最长路径和次长路径

1072. 树的最长路径 - AcWing题库

注意法一是不能维护负权值的而法二可以

代码如下:

#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:

树的中心  1073. 树的中心 - AcWing题库

一个节点距离其他点的最远的点 要么是在以自己为根的子树的距离最远点,要么是当前考虑的点向上出发到其父节点,再从父节点出发到最远的点

需要注意的是 由树的拓扑结构可知 某个节点父节点对其贡献是不能计算父节点走”回头路“的距离的,也就是从父节点出发的那条可贡献距离不能再经过该点了

下面是代码:

 

#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 的简单题

 数字转换 1075. 数字转换 - AcWing题库 
 

简单数论建图+树的直径

#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;



}

大臣的旅费1207. 大臣的旅费 - AcWing题库 
 

#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;


}

树的直径裸题(蓝桥)

旅游规划1078. 旅游规划 - AcWing题库 
 

判断点在某一条直径上吗? 最大值+次大值==直径长度(树型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中码量多一些的

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值