BZOJ1431 : MLand

Claris發表於2018-11-17

考慮任意一棵生成樹,它的代價是一個一次函式。

因此所有生成樹的最小值隨著時間變化呈現出的是一個上凸殼。

三分查詢最大值即可。

時間複雜度$O(m\log m\log w)$。

 

#include<cstdio>
#include<algorithm>
int n,m,t1,t2,f[205],i;
struct E{
  int u,v,a,b;double w;
  void set(double t){w=a*t+b;}
}e[1205];
inline bool cmp(const E&a,const E&b){return a.w<b.w;}
int F(int x){return f[x]==x?x:f[x]=F(f[x]);}
double cal(double t){
  int i;
  for(i=0;i<m;i++)e[i].set(t);
  std::sort(e,e+m,cmp);
  for(i=0;i<n;i++)f[i]=i;
  double ret=0;
  for(i=0;i<m;i++)if(F(e[i].u)!=F(e[i].v))ret+=e[i].w,f[f[e[i].u]]=f[e[i].v];
  return ret;
}
int main(){
  scanf("%d%d%d%d",&n,&m,&t1,&t2);
  for(i=0;i<m;i++)scanf("%d%d%d%d",&e[i].u,&e[i].v,&e[i].a,&e[i].b);
  double l=t1,r=t2;
  for(i=100;i--;){
    double len=(r-l)/3,m1=l+len,m2=r-len;
    double f1=cal(m1),f2=cal(m2);
    if(f1>f2)r=m2;else l=m1;
  }
  return printf("%.3f",cal(l)),0;
}