P3237 [HNOI2014]米特運輸【題解】
題目描述
米特是 星球上一種非常神祕的物質,蘊含著巨大的能量。在以米特為主要能源的 星上,這種米特能源的運輸和儲存一直是一個大問題。
星上有 個城市,我們將其順序編號為 到 , 號城市為首都。這 個城市由 條單向高速通道連線起來,構成一棵以 號城市(首部)為根的樹,高速通道的方向由樹中的兒子指向父親。樹按深度分層:根結點深度為 ,屬於第 層;根結點的子節點深度為 ,屬於第 層;依此類推,深度為 的結點屬於第 層。
建好高速通道之後, 星人開始考慮如何具體地儲存和傳輸米特資源。由於發展程度不同,每個城市儲存米特的能力不盡相同,其中第 個城市建有一個容量為 的米特儲存器。這個米特儲存器除了具有儲存的功能,還具有自動收集米特的能力。
如果到了晚上六點,有某個儲存器處於未滿的狀態,它就會自動收集大氣中蘊含的米特能源,在早上六點之前就能收集滿;但是,只有在儲存器完全空的狀態下啟動自動收集程式才是安全的,未滿而又非空時啟動可能有安全隱患。
早上六點到七點間,根節點城市( 號城市)會將其儲存器裡的米特消耗殆盡。根節點不會自動蒐集米特,它只接受子節點傳輸來的米特。
早上七點,城市之間啟動米特傳輸過程,傳輸過程逐層遞進:先是第 層節點城市向第 層(根節點城市,即 號城市)傳輸,直到第 層的儲存器滿或第2層的儲存器全為空;然後是第 層向第 層傳輸,直到對於第 層的每個節點,其儲存器滿或其予節點(位於第 層)的儲存器全為空;依此類推,直到最後一層傳輸完成。傳輸過程一定會在晚上六點前完成。
由於技術原因,運輸方案需要滿足以下條件:
不能讓某個儲存器到了晚上六點傳輸結束時還處於非空但又未滿的狀態,這個時候儲存器仍然會啟動自動收集米特的程式,而給已經儲存有米特的儲存器啟動收集程式可能導致危險,也就是說要讓儲存器到了晚上六點時要麼空要麼滿;
關於首都——即 號城市的特殊情況, 每天早上六點到七點間 號城市中的米特儲存器裡的米特會自動被消耗殆盡,即運輸方案不需要考慮首都的米特怎麼運走;
除了 號城市,每個節點必須在其子節點城市向它運輸米特之前將這座城市的米特儲存器中原本存有的米特全部運出去給父節點,不允許儲存器中殘存的米特與外來的米特發生混合;
運向某一個城市的若干個來源的米特數量必須完全相同,不然,這些來源不同的米特按不同比例混合之後可能發生危險。
現在 星人已經建立好高速通道,每個城市也有了一定儲存容量的米特儲存器。為了滿足上面的限制條件,可能需要重建一些城市中的米特儲存器。你可以,也只能,將某一座城市(包括首都)中原來存在的米特儲存器摧毀,再新建一座任意容量的新的米特儲存器,其容量可以是小數(在輸入資料中,儲存器原始容量是正整數,但重建後可以是小數),不能是負數或零,使得需要被重建的米特儲存器的數目儘量少。
輸入輸出格式
輸入格式:
第一行是一個正整數
,表示城市的數目。
接下來
行,每行一個正整數,其中的第i行表示第i個城市原來存在的米特儲存器的容量。
再接下來是
行,每行兩個正整數
,
表示城市
到城市
有一條高速通道
。
輸出格式:
輸出檔案僅包含一行,一個整數,表示最少的被重建(即修改儲存器容量)的米特儲存器的數目。
輸入輸出樣例
輸入樣例#1:
5
5
4
3
2
1
1 2
1 3
2 4
2 5
輸出樣例#1:
3
說明
【樣例解釋】
一個最優解是將 改成 , 改成 , 改成 。這樣, 和 運給 的 量相等, 和 運給 的量相等,且每天晚上六點的時候, , 滿, , , 空,滿足所有限制條件。
對於 %的資料滿足 ,
題解:
題目太長,先轉述一下題意吧
求使得父節點的權值為子節點的權值之和,且子節點權值相等 的最少的修改次數
由樣例解釋可以發現,到晚上六點時葉子結點為空,其他節點滿,而且每層的容量是相同的
那麼我們有:對於樹上任一個點,其權值一旦確定,整棵樹的權值即可確定。
例如:
從根往下遞推,累計路徑上的權值,如下圖:
!:不想寫高精用 將乘法轉化為加法
注意精度損失,要判EPS
程式碼:
#include<iostream>
#include<cstdio>
#include<ctype.h>
#include<cstring>
#include<cmath>
#include<algorithm>
const double Eps=1e-8;
using namespace std;
inline int read(){
int x=0,f=0;char ch=getchar();
while(!isdigit(ch))f|=ch=='-',ch=getchar();
while(isdigit(ch))x=x*10+(ch^48),ch=getchar();
return f?-x:x;
}
int head[500007],cnt;
double a[500007];
int w[500007],out[500007];
struct Edge{int next,to;}edge[1000007];
inline void add_edge(int from,int to){
edge[++cnt].next=head[from];
edge[cnt].to=to;head[from]=cnt;
}
void dfs(int x,int fa,double sum){
a[x]=sum+log(w[x]);
for(int i=head[x];i;i=edge[i].next){
int to=edge[i].to;
if(to==fa)continue;out[x]++;//統計當前節點的初度,也就是每層的節點數
}
for(int i=head[x];i;i=edge[i].next){
int to=edge[i].to;
if(to==fa)continue;
dfs(to,x,sum+log(out[x]));
}
}
int main(){
int n=read(),lll=1,ans=1;
for(int i=1;i<=n;++i)w[i]=read();
for(int i=1;i<n;++i){
int u=read(),v=read();
add_edge(u,v);add_edge(v,u);
}
dfs(1,0,log(1));
sort(a+1,a+n+1);
for(int i=2;i<=n;++i)if(a[i]-a[i-1]<Eps)ans=max(ans,++lll);else lll=1;
printf("%d",n-ans);
return 0;
}