bzoj 1758 [Wc2010]重建计划 01分数规划 点分治 单调队列

直接01分数规划，然后把树分治查找所有链，然后合并时就是查一个定长区间的最大值，宽搜离线下来然后单调队列搞就行了。

不过我想说这题卡常数。。。卡常数。。。常数。。数。。
不过好像三个月之前做就没有这事了。。。

树分治可以预处理快了两倍，01分数规划从二分改成迭代快了四倍+，然而bz上还是跑了20+s。。。

#include <bits/stdc++.h>
using namespace std;
#define N 110000
const double inf=1e15;
int n,L,R,tot;
int head[N],nex[N<<1],to[N<<1],val[N<<1];
int vis[N],f[N],size[N],root,rt,sum,mxd,top;
int deep[N],fa[N],st[N],q[N];
double now,sv[N],rec[N],ans;
vector<int>sr[N];
char getc()
{
    static const int LEN = 4096;
    static char buf[LEN],*S=buf,*T=buf;
    if(S == T)
    {
        T = (S=buf)+fread(buf,1,LEN,stdin);
        if(S == T)return EOF;
    }
    return *S++;
}
int read()
{
    static char ch;
    static int D;
    while(!isdigit(ch=getc()));
    for(D=ch-'0'; isdigit(ch=getc());)
        D=(D<<3)+(D<<1)+(ch-'0');
    return D;
}
void add(int x,int y,int z)
{
    tot++;
    nex[tot]=head[x];head[x]=tot;
    to[tot]=y;val[tot]=z;
}
void init()
{
    ans=-inf;
    memset(vis,0,sizeof(vis));
}
void dfs1(int x,int y)
{
    size[x]=1;
    for(int i=head[x];i;i=nex[i])
        if(to[i]!=y&&!vis[to[i]])
        {
            dfs1(to[i],x);
            size[x]+=size[to[i]];
        }
}
void dfs2(int x,int y)
{
    f[x]=0;
    for(int i=head[x];i;i=nex[i])
        if(to[i]!=y&&!vis[to[i]])
        {
            dfs2(to[i],x);
            f[x]=max(f[x],size[to[i]]);
        }
    f[x]=max(f[x],sum-size[x]);
    root=f[x]<f[root] ? x:root;
}
void get(int x,int v)
{
    st[top=1]=x;sv[1]=v;
    deep[1]=1;fa[1]=0;
    for(int i=1,t;i<=top;i++)
    {
        for(int j=head[t=st[i]];j;j=nex[j])
            if(!vis[to[j]]&&to[j]!=fa[i])
            {
                st[++top]=to[j];
                deep[top]=deep[i]+1;
                fa[top]=t;sv[top]=sv[i]+val[j];
            }
    }
    int h=1,r=0,pos=mxd;
    for(int i=1,j,t;i<=top;i=j)
    {
        double mx=-inf;t=deep[i];
        for(j=i;deep[j]==t&&j<=top;j++)
            mx=max(mx,sv[j]);
        while(pos>=0&&pos>=L-t)
        {
            while(h<=r&&rec[q[r]]<rec[pos])r--;
            q[++r]=pos;pos--;
        }
        while(h<=r&&q[h]>R-t)h++;
        if(h<=r)ans=max(ans,(rec[q[h]]+q[h]*now+mx)/(t+q[h]));
    }
    for(int i=1;i<=top;i++)
        rec[deep[i]]=max(rec[deep[i]],sv[i]-now*deep[i]);
    mxd=max(mxd,deep[top]);
}
void cal(int x)
{
    vis[x]=1;mxd=0;rec[0]=0;
    for(int i=head[x];i;i=nex[i])
        if(!vis[to[i]])
            get(to[i],val[i]);
    for(int i=0;i<=mxd;i++)rec[i]=-inf;
    for(int i=0;i<sr[x].size();i++)
        cal(sr[x][i]);
}
double check()
{
    init();
    cal(rt);
    return ans;
}
void pre(int x)
{
    vis[x]=1;
    for(int i=head[x];i;i=nex[i])
    if(!vis[to[i]])
    {
        dfs1(to[i],x);
        root=0;sum=size[to[i]];
        dfs2(to[i],x);
        sr[x].push_back(root);
        pre(root);
    }
}
int main()
{
    n=read();L=read();R=read();
    for(int i=1,x,y,z;i<n;i++)
    {
        x=read();y=read();z=read();
        add(x,y,z);add(y,x,z);
    }
    dfs1(1,0);
    f[0]=n+1;root=0;sum=n;
    dfs2(1,0);
    rt=root;pre(root);
    for(int i=0;i<=n;i++)rec[i]=-inf;
    double l=0,r=1e6;now=0;
    while(r-l>=5e-5)
    {
        l=now;
        now=check();
        r=now;
    }
    printf("%.3lf\n",l);
    return 0;
}