zoukankan      html  css  js  c++  java
  • Codeforces 618D Hamiltonian Spanning Tree(树的最小路径覆盖)

    题意:给出一张完全图,所有的边的边权都是 y,现在给出图的一个生成树,将生成树上的边的边权改为 x,求一条距离最短的哈密顿路径。

    先考虑x>=y的情况,那么应该尽量不走生成树上的边,如果生成树上有一个点的度数是n-1,那么必然需要走一条生成树上的边,此时答案为x+y*(n-2).

      否则可以不走生成树上的边,则答案为y*(n-1).

    再考虑x<y的情况,那么应该尽量走生成树上的边,由于树上没有环,于是我们每一次需要走树的一条路,然后需要从非生成树上的边跳到树的另一个点上去,

      显然跳的越少越好,于是我们只需要找到树的最小路径覆盖,跳路径覆盖数-1次就可以了。

      对于有向图的最小路径覆盖,一般是使用二分图匹配或者最大流来解决的。

      而对于树的最小路径覆盖,可以用树形DP来解决。

      令dp[x][0]表示x不与x的父亲构成路径的最小路径覆盖数,dp[x][1]表示x与x的父亲构成路径的最小路径覆盖数。

      那么则有:

        x没有儿子,dp[x][0]=dp[x][1]=1.

        x只有一个儿子,dp[x][0]=dp[x][1]=dp[son[x]][1];

        x有两个或者更多儿子,dp[x][0]=min(dp[son[x][i]][1]+dp[son[x][j]][1]+dp[son[x]][0])-1. dp[x][1]=min(dp[son[x][i]][1]+dp[son[x]][0]);

    # include <cstdio>
    # include <cstring>
    # include <cstdlib>
    # include <iostream>
    # include <vector>
    # include <queue>
    # include <stack>
    # include <map>
    # include <bitset>
    # include <set>
    # include <cmath>
    # include <algorithm>
    using namespace std;
    # define lowbit(x) ((x)&(-x))
    # define pi acos(-1.0)
    # define eps 1e-8
    # define MOD 1000000007
    # define INF 1000000000
    # define mem(a,b) memset(a,b,sizeof(a))
    # define FOR(i,a,n) for(int i=a; i<=n; ++i)
    # define FO(i,a,n) for(int i=a; i<n; ++i)
    # define bug puts("H");
    # define lch p<<1,l,mid
    # define rch p<<1|1,mid+1,r
    # define mp make_pair
    # define pb push_back
    typedef pair<int,int> PII;
    typedef vector<int> VI;
    # pragma comment(linker, "/STACK:1024000000,1024000000")
    typedef long long LL;
    int Scan() {
        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*10+ch-'0';ch=getchar();}
        return x*f;
    }
    const int N=200005;
    //Code begin...
    
    struct Edge{int p, next;}edge[N<<1];
    int head[N], cnt=1;
    int dee[N], sum, dp[N][2];
    
    void add_edge(int u, int v){edge[cnt].p=v; edge[cnt].next=head[u]; head[u]=cnt++;}
    void dfs(int x, int fa){
        int siz=0, sum=0, f=-INF, s=-INF;
        for (int i=head[x]; i; i=edge[i].next) {
            int v=edge[i].p;
            if (v==fa) continue;
            dfs(v,x); ++siz; sum+=dp[v][0];
            if (dp[v][0]-dp[v][1]>f) s=f, f=dp[v][0]-dp[v][1];
            else if (dp[v][0]-dp[v][1]>s) s=dp[v][0]-dp[v][1];
        }
        if (siz==0) dp[x][0]=dp[x][1]=1;
        else {
            if (siz==1) dp[x][0]=sum-f, dp[x][1]=sum-f;
            else dp[x][0]=sum-f-s-1, dp[x][1]=sum-f;
        }
    }
    int main ()
    {
        int n, x, y, u, v;
        scanf("%d%d%d",&n,&x,&y);
        FO(i,1,n) scanf("%d%d",&u,&v), add_edge(u,v), add_edge(v,u), ++dee[u], ++dee[v];
        if (x>=y) {
            bool flag=false;
            FOR(i,1,n) if (dee[i]==n-1) flag=true;
            if (flag) printf("%lld
    ",(LL)(n-2)*y+x);
            else printf("%lld
    ",(LL)(n-1)*y);
        }
        else {
            dfs(1,0);
            printf("%lld
    ",(LL)(dp[1][0]-1)*y+(LL)(n-dp[1][0])*x);
        }
        return 0;
    }
    View Code
  • 相关阅读:
    libcurl库进行http通讯-开篇
    hdu4059The Boss on Mars 容斥原理
    “最大子序列和”算法 java
    Android 的事件传递机制,详解
    POJ
    hdu 5652 India and China Origins
    解决Swap file &quot;.ceshi.c.swp&quot; already exists!问题
    3DMax脚本插件--改动材质&amp;贴图名称
    ASP.NET MVC+Bootstrap 实现短信验证
    [从头学数学] 第176节 概率
  • 原文地址:https://www.cnblogs.com/lishiyao/p/6974144.html
Copyright © 2011-2022 走看看