bzoj3573 [Hnoi2014]米特运输（贪心+数学）

我们发现确定了一个点的权值之后其他点的权值就也都确定了。
于是我们算出每个点不改的话根的权值是多少。
然后贪心地取相同权值最多的那个方案即可。
因为数太大，于是我们可以取个对数变成加。
然后就完了qaq

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 500010
#define eps 1e-5
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,h[N],num=0,du[N],ans=0;
double a[N],d[N];
struct edge{
    int to,next;
}data[N<<1];
inline void dfs(int x,int Fa){
    for(int i=h[x];i;i=data[i].next){
        int y=data[i].to;if(y==Fa) continue;
        d[y]=d[x]+log(du[x]);dfs(y,x);
    }
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i) a[i]=log(read());
    for(int i=1;i<n;++i){
        int x=read(),y=read();du[x]++;du[y]++;
        data[++num].to=y;data[num].next=h[x];h[x]=num;
        data[++num].to=x;data[num].next=h[y];h[y]=num;
    }for(int i=2;i<=n;++i) du[i]--;
    d[1]=log(1);dfs(1,0);
    for(int i=1;i<=n;++i) a[i]+=d[i];
    sort(a+1,a+n+1);int tmp=1;
    for(int i=2;i<=n;++i){
        if(fabs(a[i]-a[i-1])<eps) ++tmp;
        else ans=max(ans,tmp),tmp=1;
    }ans=max(ans,tmp);
    printf("%d\n",n-ans);
}