$poj$

$Description$

一颗树有 $n$ 个节点,这些节点被标号为:$1,2,3…n,$每个节点 $i$ 都有一个权值 $A[i]$。

现在要把这棵树的节点全部染色,染色的规则是:

根节点R可以随时被染色;对于其他节点,在被染色之前它的父亲节点必须已经染上了色.

每次染色的代价为$T*A[i]$,其中$T$代表当前是第几次染色.

求把这棵树染色的最小总代价。

$Sol$

贪心:树中权值最大的点,一定会在它的父结点染色后立即染色

可以这样想,假设没有树的约束,只是一个序列,那么显然是按照从大到小排序的顺序染色为最优,现在有树的约束条件,也应该让权值最大的尽量早地染色.

因为权值最大的点与它的父结点的染色是接连进行的,所以我们可以把它们合并起来.

现在有权值为$x,y,z$的三个点,已知$x$和$y$的染色是接连着进行的,其中$x$是$y$的父结点.那么有两种决策.

1.先染$x,y$,再染$z$ ,$x+2y+3z=y+(x+y)+3z$

2.先染$z$,再染$x,y, z+2x+3y=z+y+2*(x+y)$

发现无论是上面的哪种情况有一种做法是通用的:先把$y$累加进答案,然后把$y$结点与$x$结点合并,具体来说,$x$的权值改为$(x+y)/2$,删除$y$结点.为什么是把权值改成$(x+y)/2$呢?

有一个十分重要的问题,就是如何判断这两种做法的优劣.设合并后的结点为$i$,这个结点所包含的结点数(它自己和合并进去的)为$s[i]$,包含的所有的权值和为$a[i]$.

1.$a[i]+(s[i]+1)*z.$

2.$z+a[i]*2.$

发现应该把$1$式中的$z$的系数变成$2$,把两个式子同时加上$(s[i]-1)*z$,再同除$s[i]$,变成

1.$a[i]/s[i]+2*z$

2.$z+a[i]/s[i]$

可以看成权值为$a[i]/s[i]$的结点和$z$的两个结点,根据最开始所说的贪心可知比较它们的大小就能决定先给谁染色了,于是问题就解决了.

还是要强调的是上面所赋予的新权值只能用于与别的结点比较大小从而决定染色顺序,至于答案的累加,就在对于最初始的两个式子的分析中"把$y$累加进答案"这句话.具体来说$j$结点累加进$i$结点,$ans+=a[j]*s[i]$.

$Code$

#include<iostream>
#include<cstdio>
#define il inline
#define Rg register
#define go(i,a,b) for(Rg int i=a;i<=b;i++)
#define yes(i,a,b) for(Rg int i=a;i>=b;i++)
#define db double
#define ll long long
using namespace std;
il int read()
{
int x=,y=;char c=getchar();
while(c<''||c>''){if(c=='-')y=-;c=getchar();}
while(c>=''&&c<=''){x=(x<<)+(x<<)+c-'';c=getchar();}
return x*y;
}
int n,rt,a[],fa[],s[];
ll as;
int main()
{
//while(1)
{
n=read(),rt=read();as=;
if(!n && !rt)return ;
go(i,,n)a[i]=read(),s[i]=;
go(i,,n-){int x=read(),y=read();fa[y]=x;}
go(T,,n-)
{
db maxs=;int pos;
go(i,,n)
if(i!=rt && 1.0*a[i]/s[i]>=maxs)maxs=1.0*a[i]/s[i],pos=i;
go(i,,n)
if(fa[i]==pos)fa[i]=fa[pos];
as+=a[pos]*s[fa[pos]];
s[fa[pos]]+=s[pos];
a[fa[pos]]+=a[pos];
a[pos]=;
}
as+=a[rt];
printf("%lld\n",as);
}
return ;
}

随机推荐

  1. 《C语言深度解剖》学习笔记之内存管理

    第5章 内存管理 1.野指针 定义指针变量的同时最好初始化为NULL,用完指针后也将变量的值设置为NULL.也就是说除了使用时,别的时间都把它设置为NULL 2.堆,栈和静态区 堆:由malloc系列 ...

  2. Android 高仿微信(QQ)滑动弹出编辑、删除菜单效果,增加下拉刷新功能

    不可否认,微信.QQ列表的滑动删除.编辑功能着实很经典(从IOS那边模仿过来的),然.Android这边,对列表的操作,其实大多还停留上下文菜单来实现. Android如何实现list item的滑动 ...

  3. hdu 4114 Disney's FastPass(最短路+状态压缩)

    Disney's FastPass Time Limit: 20000/10000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Othe ...

  4. win10提示 磁盘包含不是“PARTITION_BASIC_DATA_GUID"类型的分区

    在win10创建新的磁盘分区的时候,有时候会提示 磁盘包含不是"PARTITION_BASIC_DATA_GUID"类型的分区 如果你试了其他方法都不凑效,那么看看你已经有几个盘( ...

  5. Codeforces Round #577 (Div 2)

    A. Important Exam 水题 #include<iostream> #include<string.h> #include<algorithm> #in ...

  6. idea乱码问题(全)

    中文乱码问题分类: 编码普通中文乱码 properties文件中文乱码 console控制台中文乱码 搜索框中文乱码 svn注释中文乱码 问题截图: 2.properties文件中文乱码 4,.搜索框 ...

  7. python中break和continue的区别

    python中break和continue的区别   break和continue 1.break 意思为结束循环   例: i = 0 while i<10:     i+=1     if ...

  8. 函数的渐近的界&阶的比较

    一.函数的渐近的界   我们在研究算法性能的时候,往往会在意算法的运行时间,而运行时间又与算法输入的规模相关,对于一个算法,我们可以求出运行时间和输入规模的函数,当输入规模足够大时,站在极限的角度看, ...

  9. supersocket Silverlight 策略服务器

    <?xml version="1.0"?> <configuration> <configSections> <section name= ...

  10. 【t001】不高兴的津津

    Time Limit: 1 second Memory Limit: 50 MB [问题描述] 津津上初中了.妈妈认为津津应该更加用功学习,所以津津除了上学之外,还要参加妈妈为她报名的各科复习班.另外 ...