BZOJ 1758

mengbing。。。。。

大家都说这是一道卡常题 、名不虚传、233333

树分治 再加上迭代答案 而不是二分 就可以过了

迭代的次数 平均下来 大概就只有 2.8次的样子

需要单调队列优化更新答案的过程

  1 #include <bits/stdc++.h>
  2 #define N 100010
  3 #define inf 1e10
  4 using namespace std;
  5 inline int read()
  6 {
  7     int x=0,f=1; char ch=getchar();
  8     while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1; ch=getchar();}
  9     while(ch>=‘0‘&&ch<=‘9‘){x=x*10+ch-‘0‘; ch=getchar();}
 10     return x*f;
 11 }
 12 struct edge
 13 {
 14     int v,next,w;
 15 }vs[N<<1];
 16 int n,L,R,q[N],vis[N],sum,Rtt;
 17 int st[N],ee,Rt,size[N],tp[N],dep[N];
 18 int stk[N],top,fa[N],mxd;
 19 double prestans=0.0,ans,vopt[N],,f[N];
 20 vector <int> Gp[N];
 21 inline void addedge(int u,int v,int w)
 22 {
 23     vs[++ee].v=v;vs[ee].w=w;
 24     vs[ee].next=st[u];st[u]=ee;
 25 }
 26 void dfs1(int rt,int pr)
 27 {
 28     size[rt]=1;
 29     for(int i=st[rt];i;i=vs[i].next)
 30     {
 31         if(vs[i].v==pr||vis[vs[i].v]) continue;
 32         dfs1(vs[i].v,rt);
 33         size[rt]+=size[vs[i].v];
 34     }
 35 }
 36 void dfs2(int rt,int pr)
 37 {
 38     tp[rt]=0;
 39     for(int i=st[rt];i;i=vs[i].next)
 40     {
 41         if(vs[i].v==pr||vis[vs[i].v]) continue;
 42         dfs2(vs[i].v,rt);
 43         tp[rt]=max(tp[rt],size[vs[i].v]);
 44     }
 45     tp[rt]=max(tp[rt],sum-size[rt]);
 46     if(tp[rt]<tp[Rtt]) Rtt=rt;
 47 }
 48 void predeal(int rt)
 49 {
 50     vis[rt]=1;
 51     for(int i=st[rt];i;i=vs[i].next)
 52     {
 53         if(vis[vs[i].v]) continue;
 54         Rtt=0; dfs1(vs[i].v,rt);
 55         sum=size[vs[i].v];
 56         dfs2(vs[i].v,rt);
 57         Gp[rt].push_back(Rtt);
 58         predeal(Rtt);
 59     }
 60 }
 61 void getans(int rt,int val)
 62 {
 63     stk[top=1]=rt; fa[1]=0; f[1]=val; dep[1]=1;
 64     for(int i=1;i<=top;i++)
 65     {
 66         for(int j=st[stk[i]];j;j=vs[j].next)
 67         {
 68             if(vis[vs[j].v]||vs[j].v==fa[i]) continue;
 69             stk[++top]=vs[j].v; dep[top]=dep[i]+1;
 70             fa[top]=stk[i]; f[top]=f[i]+vs[j].w;
 71         }
 72     }
 73     int l=1,r=0,pre=mxd;
 74     for(int i=1,j;i<=top;i=j)
 75     {
 76         double mx=-inf;
 77         for(j=i;dep[j]==dep[i]&&j<=top;j++)
 78             mx=max(mx,f[j]);
 79         while(pre>=0&&pre>=L-dep[i])
 80         {
 81             while(vopt[q[r]]<vopt[pre]&&r>=l) r--;
 82             q[++r]=pre; pre--;
 83         }
 84
 85         while(q[l]+dep[i]>R&&l<=r) l++;
 86         if(l<=r) ans=max(ans,(vopt[q[l]]+mx+prestans*q[l])/(q[l]+dep[i]));
 87     }
 88     for(int i=1;i<=top;i++)
 89         vopt[dep[i]]=max(vopt[dep[i]],f[i]-dep[i]*prestans);
 90     mxd=max(mxd,dep[top]);
 91 }
 92 void solve(int rt)
 93 {
 94     vis[rt]=1; mxd=0; vopt[0]=0;
 95     for(int i=st[rt];i;i=vs[i].next)
 96     {
 97         if(vis[vs[i].v]) continue;
 98         getans(vs[i].v,vs[i].w);
 99     }
100     for(int i=0;i<=mxd;i++) vopt[i]=-inf;
101     for(unsigned i=0;i<Gp[rt].size();i++)
102         solve(Gp[rt][i]);
103 }
104 double check()
105 {
106     memset(vis,0,sizeof vis);
107     ans=-inf;
108     solve(Rt);
109     return ans;
110 }
111 int main()
112 {
113     //freopen("read.in","r",stdin);
114     n=read();L=read();R=read();
115     for(int i=1;i<n;i++)
116     {
117         int x=read(),y=read(),w=read();
118         addedge(x,y,w);addedge(y,x,w);
119     }
120     dfs1(1,0);
121     Rtt=0; tp[0]=n+1; sum=n;
122     dfs2(1,0);
123     Rt=Rtt; predeal(Rt);
124     for(int i=1;i<=n;i++) vopt[i]=-inf;
125         for(int i=1;i<=3;i++)
126         prestans=check();
127     printf("%.3lf",prestans);
128     return 0;
129 }

↓↓↓↓ dmk ↓↓↓↓

 1 #include <bits/stdc++.h>
 2 #define N 100000
 3 using namespace std;
 4
 5 struct edge
 6 {
 7     int u,v,next,w;
 8 }vs[N<<1];
 9 int st[N+100],ee,dep[N+100],mx1,mx2;
10 inline void addedge(int u,int v,int w)
11 {
12     vs[++ee].v=v;vs[ee].next=st[u];
13     vs[ee].u=u;vs[ee].w=w;st[u]=ee;
14 }
15 int main()
16 {
17     srand(time(0));
18     freopen("read.in","w",stdout);
19     printf("%d\n",N);
20     for(int i=2;i<=N;i++)
21     {
22         int u=rand()*rand()%(i-1)+1;
23         addedge(u,i,rand()*rand()%1000000);
24         dep[i]=dep[u]+1;
25         if(dep[i]>=mx1) mx2=mx1,mx1=dep[i];
26         else if(dep[i]>mx2) mx2=dep[i];
27     }
28     int L=rand()*rand()%(mx1+mx2)+1,R=rand()*rand()%(mx1+mx2)+1;
29     if(L>R) swap(L,R);
30 /*     sort(dep+1,dep+1+N);
31     for(int i=1;i<=N;i++)
32         printf("%d ",dep[i]); printf("\n"); */
33     printf("%d %d\n",L,R);
34     for(int i=1;i<=N;i++)
35         printf("%d %d %d\n",vs[i].u,vs[i].v,vs[i].w);
36     return 0;
37 }

好好讲道理 应该求出直径 所以这个dmk 有问题2333 但总体还是能拍得出来 逃