zoukankan      html  css  js  c++  java
  • [cf516D]Drazil and Morning Exercise

    令$f_{k}$为离$k$最远的点到$k$的距离,任取树的一条直径$(x,y)$,有$f_{k}=max(dis(k,x),dis(k,y))$

    更进一步的,取直径中点$mid$(这里定义为$f_{mid}$最小的点,有多个任取一点)并以其为根建树,则所有节点儿子的$f$不小于父亲的$f$

    根据这个性质,枚举这个连通块的根$k$(也可以说是连通块内$f$最小的点),即统计其子树内权值在$[f_{k},f_{k}+l]$中的点数量,但可持久化线段树去维护会TLE

    事实上,这个还有一种比较巧妙地方式计算:

    将所有$f$从大到小排序,依次枚举$k$,去计算$k$子树内权值在$[f_{k},f_{k}+l]$的点

    通过单调性,在左边删除、右边加入可以维护出在排序后的对应区间,现在将这个区间内的所有点的导出子图用并查集去维护,那么所求的也就是$k$子树大小

    关于如何维护并查集,也就是要支持插入和删除,插入将其树上所有儿子都与其合并,删除的节点必然是导出子图上的叶子,直接删除即可

      1 #include<bits/stdc++.h>
      2 using namespace std;
      3 #define N 100005
      4 #define ll long long
      5 #define fi first
      6 #define se second
      7 struct Edge{
      8     int nex,to,len;
      9 }edge[N<<1];
     10 int V,E,n,m,x,y,z,head[N],sh[N],f[N][21],fa[N],sz[N];
     11 ll l,dep[N];
     12 pair<ll,int>ans[N];
     13 int find(int k){
     14     if (k==fa[k])return k;
     15     return fa[k]=find(fa[k]);
     16 }
     17 void merge(int x,int y){
     18     x=find(x),y=find(y);
     19     if (x!=y){
     20         fa[x]=y;
     21         sz[y]+=sz[x];
     22     }
     23 }
     24 void add(int x,int y,int z){
     25     edge[E].nex=head[x];
     26     edge[E].to=y;
     27     edge[E].len=z;
     28     head[x]=E++;
     29 }
     30 int lca(int x,int y){
     31     if (sh[x]<sh[y])swap(x,y);
     32     for(int i=20;i>=0;i--)
     33         if (sh[f[x][i]]>=sh[y])x=f[x][i];
     34     if (x==y)return x;
     35     for(int i=20;i>=0;i--)
     36         if (f[x][i]!=f[y][i]){
     37             x=f[x][i];
     38             y=f[y][i];
     39         }
     40     return f[x][0];
     41 }
     42 ll dis(int x,int y){
     43     return dep[x]+dep[y]-2*dep[lca(x,y)];
     44 }
     45 void pre_diam(int k,int fa,ll s){
     46     dep[k]=s;
     47     for(int i=head[k];i!=-1;i=edge[i].nex)
     48         if (edge[i].to!=fa)pre_diam(edge[i].to,k,s+edge[i].len);
     49 }
     50 void pre_lca(int k,int fa,ll s,int ss){
     51     sh[k]=ss;
     52     dep[k]=s;
     53     f[k][0]=fa;
     54     for(int i=1;i<=20;i++)f[k][i]=f[f[k][i-1]][i-1];
     55     for(int i=head[k];i!=-1;i=edge[i].nex)
     56         if (edge[i].to!=fa)pre_lca(edge[i].to,k,s+edge[i].len,ss+1);
     57 }
     58 void pre_son(int k,int fa){
     59     f[k][0]=fa;
     60     for(int i=head[k];i!=-1;i=edge[i].nex)
     61         if (edge[i].to!=fa)pre_son(edge[i].to,k);
     62 }
     63 int main(){
     64     scanf("%d",&n);
     65     memset(head,-1,sizeof(head));
     66     for(int i=1;i<n;i++){
     67         scanf("%d%d%d",&x,&y,&z);
     68         add(x,y,z);
     69         add(y,x,z);
     70     }
     71     pre_diam(1,0,0);
     72     x=1;
     73     for(int i=2;i<=n;i++)
     74         if (dep[i]>dep[x])x=i;
     75     pre_diam(x,0,0);
     76     y=1;
     77     for(int i=2;i<=n;i++)
     78         if (dep[i]>dep[y])y=i;
     79     pre_lca(1,1,0,0);
     80     for(int i=1;i<=n;i++)ans[i]=make_pair(max(dis(x,i),dis(y,i)),i);
     81     sort(ans+1,ans+n+1);
     82     pre_son(ans[1].se,0);
     83     scanf("%d",&m);
     84     for(int i=1;i<=m;i++){
     85         scanf("%lld",&l);
     86         int max_ans=0;
     87         for(int j=1;j<=n;j++){
     88             fa[j]=j;
     89             sz[j]=1;
     90         }
     91         for(int j=n,k=n;j;j--){
     92             for(int t=head[ans[j].se];t!=-1;t=edge[t].nex)
     93                 if (edge[t].to!=f[ans[j].se][0])merge(ans[j].se,edge[t].to);
     94             while (ans[k].fi>ans[j].fi+l){
     95                 sz[find(ans[k].se)]--;
     96                 k--;
     97             }
     98             max_ans=max(max_ans,sz[find(ans[j].se)]);
     99         }
    100         printf("%d
    ",max_ans);
    101     }
    102 }
    View Code
  • 相关阅读:
    怎样检测数据类型
    怎样理解函数参数的传递
    怎样在微信H5中点击直接跳转到公众号
    怎样启动和关闭nginx服务器
    怎样测试nginx.conf配置文件的正确性
    怎样查找/搜索文件
    怎样创建用户组并添加用户进用户组
    怎样查看Nginx版本号
    怎样写一个Hello World!
    怎样理解基本类型(原始类型)的数据和引用类型的数据
  • 原文地址:https://www.cnblogs.com/PYWBKTDA/p/14521924.html
Copyright © 2011-2022 走看看