• 【洛谷P3128】最大流


    本题当然可以用树剖解决,而且是树剖的模板题。但是对于本题来说,有一种更巧妙的办法

    树上差分

    类似于普通的差分,我们对于题目中的每一条路线i->j,把它拆成i->lca(i,j)->j.然后将i点、j点的权值加一,将lca(i,j)的权值减二

    乍看起来没什么毛病,但是仔细想想,对lca(i,j)的子树进行前缀和,则lca(i,j)点的权值加二又减二,而题目要求lca(i,j)权值应当+1

    所以我们将lca(i,j)的权值减一,将lca(i,j)的父亲的权值减一,这样就能符合题意了

    差分之前预处理每个点的祖先,LCA+差分即可

     1 #include <iostream>
     2 #include <cstdio>
     3 #include <cstring>
     4 typedef long long ll;
     5 inline int read() {
     6     int ret=0,f=1;
     7     char c=getchar();
     8     while(c<'0'||c>'9') {if(c=='-') f=-1;c=getchar();}
     9     while(c<='9'&&c>='0') ret=ret*10+c-'0',c=getchar();
    10     return ret*f;
    11 }
    12 using namespace std;
    13 int n,k;
    14 struct edge {
    15     int next,to;
    16 }a[50010<<1];
    17 int num,head[50010<<1];
    18 int sum[50010],ans;
    19 int d[50010],fa[50010][31];
    20 inline void add(int from,int to) {
    21     a[++num].next=head[from]; a[num].to=to; head[from]=num;
    22     swap(from,to);
    23     a[++num].next=head[from]; a[num].to=to; head[from]=num;
    24 }
    25 void find(int u,int f) {
    26     d[u]=d[f]+1;
    27     fa[u][0]=f;
    28     for(int i=0;i<30;i++) fa[u][i+1]=fa[fa[u][i]][i];
    29     for(int i=head[u];i;i=a[i].next) {
    30         int v=a[i].to;
    31         if(v==f) continue ;
    32         find(v,u);
    33     }
    34 }
    35 int lca(int x,int y) {
    36     if(d[x]<d[y]) swap(x,y);
    37     for(int i=30;i>=0;i--) {
    38         if(d[fa[x][i]]>=d[y]) x=fa[x][i];
    39         if(x==y) return x;
    40     }
    41     for(int i=30;i>=0;i--)    
    42         if(fa[x][i]!=fa[y][i]) {
    43             x=fa[x][i];
    44             y=fa[y][i];
    45         }
    46     return fa[x][0];
    47 }
    48 void dfs(int u,int f) {
    49     for(int i=head[u];i;i=a[i].next) {
    50         int v=a[i].to;
    51         if(v==f) continue ;
    52         dfs(v,u);
    53         sum[u]+=sum[v];
    54     }
    55     ans=max(ans,sum[u]);
    56 }
    57 int main() {
    58     n=read(); k=read();
    59     for(int i=1;i<n;i++) add(read(),read());
    60     find(1,0);
    61     for(int i=1;i<=k;i++) {
    62         int x=read(),y=read();
    63         sum[x]++;
    64         sum[y]++;
    65         int ret=lca(x,y);
    66         sum[ret]--;
    67         sum[fa[ret][0]]--; 
    68     }
    69     dfs(1,0);
    70     printf("%d
    ",ans);
    71     return 0;
    72 }
    AC Code
  • 相关阅读:
    git配置公钥,私钥
    vscode之vue文件代码格式化代码无效解决办法
    [python 并行3]进程
    [spring 并行6]分布式
    [spring 并行5]GPU
    [python 并行2]线程
    [spring 并行4]异步
    [python 并行1]简介
    [flask] flask api + vue 跨域问题
    [spring学习4] MVC
  • 原文地址:https://www.cnblogs.com/shl-blog/p/10544014.html
Copyright © 2020-2023  润新知