树形dp专栏

前言

自己树形dp太菜了,要重点搞

219D Choosing Capital for Treeland

终于自己做了一道不算那么毒瘤的换根dp

\(f[u]\) 表示以 \(u\) 为根,子树内总共需要交换的边数, \(up[u]\) 表示以 \(u\) 为根,子树外总共需要交换的边数。

Dfs1 求出 \(f[u]\) ,就有:

\[f[u]=\sum_{v∈son[u]} f[v] + (edge[u->v] == 1)\]

edge[u->v] 表示 u->v 这条边的方向是不是 u->v

Dfs2 求出 \(up[v]\)注意,是从u点求u的儿子点v),容斥一下,就有:

\[up[v]=f[u]-f[v]+up[u]+(+1 / -1)\]

(+1 / -1) 是看 edge[u->v]是否等于 1,是的话就有多一条边交换方向,不是的话就要-1,因为多算了一条边

Code

#include<bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while(ch<'0' || ch>'9') { if(ch=='-') f=-1; ch=getchar(); }
    while(ch>='0'&&ch<='9') { x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); }
    return x * f;
}
const int N = 2e5+7;
int n,cnt;
int head[N],f[N],up[N];
struct Edge {
    int next,to,flag;
}edge[N<<1];
inline void add(int u,int v,int flag) {
    edge[++cnt] = (Edge)<%head[u],v,flag%>;
    head[u] = cnt;
}
void Dfs1(int u,int fa) {
    for(int i=head[u];i;i=edge[i].next) {
        int v = edge[i].to;
        if(v != fa) {
            Dfs1(v,u);
            f[u] += f[v] + (edge[i].flag==0); //反边 
        }
    }
}
void Dfs2(int u,int fa) {
    for(int i=head[u];i;i=edge[i].next) {
        int v = edge[i].to;
        if(v != fa) {
            up[v] = f[u] - f[v] + up[u];
            if(edge[i].flag == 1) up[v]++;
            else up[v]--;
            Dfs2(v,u);
        }
    }
}
int main()
{
    n = read();
    for(int i=1,u,v;i<=n-1;++i) {
        u = read(), v = read();
        add(u,v,1), add(v,u,0);
    }
    Dfs1(1,0);
    Dfs2(1,0);
    int ans = INF;
    for(int i=1;i<=n;++i)
        ans = min(ans,f[i]+up[i]);
    printf("%d\n",ans);
    for(int i=1;i<=n;++i)
        if(f[i]+up[i] == ans) printf("%d ",i);
    return 0;
}

533B Work Group

这题都看了题解才会(虽然说想到了题解这个状态,但不会转移)。。。dp功力还不行啊qwq

\(f[u][0/1]\) 表示子树总和为 偶数/奇数 的最大价值(且包括根,但是根可以选或不选)

有人也许会说,奇数怎么可呢,不是说一定要偶数吗?(其实这个自己推推数据在纸上画画就差不多知道了)

像下面这张图

黑色代表选了。不选根就可以选奇数的儿子呗。转移

\[f[u][0]=\max_{v∈son[u]}\{f[v][0]+f[u][0],f[v][1]+f[v][1]\}\]

\[f[u][1]=\max_{v∈son[u]}\{f[v][0]+f[u][1],f[v][1]+f[u][0]\}\]

最后 \(f[u][1]=\max\{f[u][1],f[u][0]+p[u]\}\)

(1代表奇,0代表偶)1+1=0,0+0=0; 1+0=0+1=1; 这个很好理解。

可是这个 \(f[u][0]\) 来更新 \(f[u][0]\) 怎么理解呢?

其实就是前面这个 \(f[u][0]\) 是我们待更新的,后面这个 \(f[u][0]\) 是之前遍历的子树里的总最优解,意义有所不同。遍历过程中的 \(f[u][0]\) 稍别与\(f[u][0]\)的定义的,只有所有更新结束后它才是u的子树选偶数个的最大值qwq。(这个不是很简单的东西吗,你怎么想了这么久啊我确实想了这么久

Code

#include<bits/stdc++.h>
#define INF 1e18
#define int long long
using namespace std;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while(ch<'0' || ch>'9') { if(ch=='-') f=-1; ch=getchar(); }
    while(ch>='0'&&ch<='9') { x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); }
    return x * f;
}
const int N = 2e5+7;
int n,cnt;
int head[N],p[N];
int f[N][2];    //f[i,0/1] 表示以i为根 子树总数是偶数/是奇数的  
struct Edge {
    int next,to;
}edge[N<<1];
inline void add(int u,int v) {
    edge[++cnt] = (Edge)<%head[u],v%>;
    head[u] = cnt;
}
void Dfs1(int u,int fa) {
    //printf("QLL:: %d %d\n",u,fa);
    f[u][1] = -INF; //f[u,1] 开始不能是奇数 
    for(int i=head[u];i;i=edge[i].next) {
        int v = edge[i].to;
        if(v != fa) {
            Dfs1(v,u);
            int x0 = f[u][0], x1 = f[u][1];
            f[u][0] = max(f[u][0],max(x0+f[v][0],x1+f[v][1]));
            f[u][1] = max(f[u][1],max(x0+f[v][1],x1+f[v][0]));
        }
    }
    f[u][1] = max(f[u][1],f[u][0]+p[u]);
    //printf("ELL :: %d %d %d\n",u,f[u][0],f[u][1]);
}
signed main()
{
    n = read();
    for(int i=1,u;i<=n;++i) {
        u = read(); p[i] = read();
        if(u!=-1) add(u,i), add(i,u);
    }
    Dfs1(1,0);
    printf("%lld\n",max(f[1][0],f[1][1]));
    return 0;
}

700B Connecting Universities

挺思维的一题。从点与点的配对完全没有办法下手,从整体的考虑,一条边可以有几条路径经过,这题就迎刃而解了

假如一条边连的两个点 \((x,y)\)\(x\)这边这一团的大学有 \(f[x]\) 座,\(y\)这边的这一团大学有\(f[y]\),我们一定要让 \(\min(f[x],f[y])\) 座大学经过这条边与另一端的大学配对。为什么?

\(f[x]<f[y]\) 如果 \(x\) 这边这\(f[x]\)个大学不去经过这条边 \((x,y)\) 向另一端配对,在 \(x\) 这边这一端自己给自己配对,答案就不是最优,自己给自己配对没有发展,最后:

\[ans=\sum\min(f[u],2*K-f[u])\]

Code

#include<bits/stdc++.h>
#define int long long
using namespace std;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while(ch<'0' || ch>'9') { if(ch=='-') f=-1; ch=getchar(); }
    while(ch>='0'&&ch<='9') { x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); }
    return x * f;
}
const int N = 200007;
int n,m,cnt,K,ans;
int head[N],f[N];
struct Edge {
    int next,from,to;
}edge[N<<1];
inline void add(int u,int v) {
    edge[++cnt] = (Edge)<%head[u],u,v%>;
    head[u] = cnt;
}
void Dfs(int u,int fa) {
    for(int i=head[u];i;i=edge[i].next) {
        int v = edge[i].to;
        if(v != fa) {
            Dfs(v,u); f[u] += f[v];
        }
    }
    ans += min(f[u],K-f[u]);
}
signed main()
{
    n = read(), K = read(); K <<= 1;
    for(int i=1,x;i<=K;++i)
        x = read(), f[x]++;
    for(int i=1,u,v;i<=n-1;++i) {
        u = read(), v = read();
        add(u,v), add(v,u);
    }
    Dfs(1,0);
    printf("%lld\n",ans);
    return 0;
}

Anton and Tree

又是一道思维难度很高的题。。。

作为提高组选手应该都能想到把一团相同颜色的点缩成一团。这个图就变成了黑白相间的一棵树。

首先这是样例里的树

如我所说,把同色点缩成一个点就是这个样子

接下来就是推推结论了,(作为提高组选手应该会觉得很好推

先给结论: 最少点击次数=(缩点后树的直径+1)/2。为什么呢?

我们不妨把直径拎出来,假设现在的直径就是下面这个图:

我们最优策略是对直径中间的一点点击一下,这样它周围的两个点就和他缩在一起了,就相当于这条链的长度-2,比如下面这张图;

因此对把直径单独拎出来最后缩成一个点的次数可以算出是 \(\frac{d+1}{2}\) (当然这里的直径d是指边数,结果向下取整)

为什么这棵树要操作的次数就是直径要操作的次数呢?

对此,我们可以把其他点看成直径上一些点的分支,如下图:

自己手推一下直径的缩点过程,发现缩点的同时,分支也缩了一些点进去,而且,缩了一圈。进而发现这些分支会先缩完,为什么?

每缩完一个点,这个点的周围就会缩小一圈,那树上最长的一条链是什么啊?直径呗。所以比直径小的会在缩直径的同时一起缩完(这样讲能理解吧,因为我很菜,不会严格的证明,再有问题就自己手推吧

Code

#include<bits/stdc++.h>
using namespace std;
inline int read() {
    int x=0,f=1; char ch=getchar();
    while(ch<'0' || ch>'9') { if(ch=='-') f=-1; ch=getchar(); }
    while(ch>='0'&&ch<='9') { x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); }
    return x * f;
}
const int N = 200007;
int n,m,cnt,mxc,ans;
int head[N];
bool col[N];
struct Edge {
    int next,to;
}edge[N<<1];
inline void add(int u,int v) {
    edge[++cnt] = (Edge)<%head[u],v%>;
    head[u] = cnt;
}
void Dfs(int u,int fa,int dep) {
    if(dep > ans) {
        mxc = u; ans = dep;
    }
    for(int i=head[u];i;i=edge[i].next) {
        int v = edge[i].to;
        if(v != fa) {
            if(col[u]==col[v]) Dfs(v,u,dep);
            else Dfs(v,u,dep+1);
        }
    }
}
int main()
{
    n = read();
    for(int i=1;i<=n;++i) col[i] = read();
    for(int i=1,u,v;i<=n-1;++i) {
        u = read(), v = read();
        add(u,v), add(v,u);
    }
    Dfs(1,0,0);
    Dfs(mxc,0,0);
    printf("%d\n",(ans+1)>>1);
    return 0;
}
全部评论

相关推荐

11-09 01:22
已编辑
东南大学 Java
高级特工穿山甲:羡慕,我秋招有家企业在茶馆组织线下面试,约我过去“喝茶详谈”😢结果我去了发现原来是人家喝茶我看着
点赞 评论 收藏
分享
评论
点赞
收藏
分享
牛客网
牛客企业服务