• 【HAOI2015】 T1


    为什么感觉越来越迷了X. X

    原题:

    有一棵点数为 N 的树,树边有边权。给你一个在 0~N 之内的正整数 K,你要在这棵树中选择 K 个点,将其染成黑色,并将其他的N-K 个点染成白色。将所有点染色后,你会获得黑点两两之间的距离加上白点两两之间的距离的和的受益。问受益最大值是多少。

    n<=2000

    一眼树形DP,然而状态似乎不太好转移啊

    我最开始是f[i][j]表示第i个点为根的子树中j个黑点到i的距离和,然后发现还有白点的距离和要计算,再开一个数组又可能出现两种方案不一样的情况

    思考无果,只能膜拜别人的题解,最后看到萌帝的

    正解是f[i][j]表示第i个点为根的子树中有j个黑点的最优答案,转移的时候因为不论是白点还是黑点想要和这个子树外的白点/黑点连接一定会经过i和某子节点之间的边,所以"(外面黑点个数*里面黑点个数+外面白点个数*里面白点个数)*这条边的权值"就是对答案的贡献

    然后枚举子树中有k个黑点,子树选k个黑点的最优答案+这条边对答案的贡献就是以i为根的子树中的最优答案,酱紫就可以往上转移了

    还有一个问题就是f不能在一开始就初始化成0,而是要初始化成-oo并在每次dfs开始的时候f[x][0]=f[x][1]=0,因为这是一个类似背包的东西,而且必须保证这个包装满

    感觉好迷啊,关于这题的具体思路和上面提到的背包的问题如何往更广泛的应用去扩展感觉很困难啊

    还要再多想/看

    (也许是我状态不好X. X

    代码:

     1 #include<iostream>
     2 #include<cstdio>
     3 #include<algorithm>
     4 #include<cstring>
     5 #include<cmath>
     6 using namespace std;
     7 #define ll long long
     8 ll rd(){ll z=0,mk=1;  char ch=getchar();
     9     while(ch<'0'||ch>'9'){if(ch=='-')mk=-1;  ch=getchar();}
    10     while(ch>='0'&&ch<='9'){z=(z<<3)+(z<<1)+ch-'0';  ch=getchar();}
    11     return z*mk;
    12 }
    13 struct ddd{int nxt,y;  ll v;}e[4100];  int lk[2100],ltp=0;
    14 inline void ist(int x,int y,ll z){  e[++ltp].nxt=lk[x],lk[x]=ltp,e[ltp].y=y,e[ltp].v=z;}
    15 int n,m;
    16 ll f[2100][2100],sz[2100];
    17 void dfs(int x,int y){
    18     sz[x]=1,f[x][0]=f[x][1]=0;
    19     ll bwl=0;
    20     for(int i=lk[x];i;i=e[i].nxt)if(e[i].y!=y){
    21         dfs(e[i].y,x);
    22         sz[x]+=sz[e[i].y];
    23         for(int j=sz[x];j>=0;--j)
    24             for(int k=0;k<=sz[e[i].y] && k<=j;++k){
    25                  bwl=(ll)k*(m-k)+(ll)(sz[e[i].y]-k)*(n-m-sz[e[i].y]+k);
    26                  bwl=bwl*e[i].v+f[e[i].y][k];
    27                  f[x][j]=max(f[x][j],f[x][j-k]+bwl);
    28             }
    29     }
    30 }
    31 int main(){//freopen("ddd.in","r",stdin);
    32     memset(f,-10,sizeof(f));
    33     cin>>n>>m;
    34     ll l,r,v;
    35     for(int i=1;i<n;++i){
    36         l=rd(),r=rd(),v=rd();
    37         ist(l,r,v),ist(r,l,v);
    38     }
    39     dfs(1,0);
    40     cout<<f[1][m]<<endl;
    41     return 0;
    42 }
    View Code
  • 相关阅读:
    201771010111 李瑞红 《第十二周学习总结》
    201771010111李瑞红《第十一周学习总结》
    李瑞红201771010111第十周 学习总结
    李瑞红201771010111《第九周学习总结》
    201771010111李瑞红《面向对象的程序设计》第八周实验总结
    2017710101111李瑞红《第七周学习总结》
    李瑞红201771010111《第六周学习总结》
    李瑞红201771010111《面向对象程序设计(java)》第四周学习总结
    201771010113 李婷华《面向对象程序设计(Java)》第十二周总结
    201771010113 李婷华 《面向对象程序设计(Java)》第十一周总结
  • 原文地址:https://www.cnblogs.com/JSL2018/p/6472193.html
Copyright © 2020-2023  润新知