zoukankan      html  css  js  c++  java
  • Max Flow

    Max Flow

    题目描述

    Farmer John has installed a new system of N−1 pipes to transport milk between the N stalls in his barn (2≤N≤50,000), conveniently numbered 1…N. Each pipe connects a pair of stalls, and all stalls are connected to each-other via paths of pipes.

    FJ is pumping milk between K pairs of stalls (1≤K≤100,000). For the ith such pair, you are told two stalls si and ti, endpoints of a path along which milk is being pumped at a unit rate. FJ is concerned that some stalls might end up overwhelmed with all the milk being pumped through them, since a stall can serve as a waypoint along many of the K paths along which milk is being pumped. Please help him determine the maximum amount of milk being pumped through any stall. If milk is being pumped along a path from si to ti, then it counts as being pumped through the endpoint stalls si and ti, as well as through every stall along the path between them.

    输入

    The first line of the input contains N and K.

    The next N−1 lines each contain two integers x and y (x≠y) describing a pipe between stalls x and y.

    The next K lines each contain two integers s and t describing the endpoint stalls of a path through which milk is being pumped.

    输出

     An integer specifying the maximum amount of milk pumped through any stall in the barn.

    样例输入

    5 10
    3 4
    1 5
    4 2
    5 4
    5 4
    5 4
    3 5
    4 3
    4 3
    1 3
    3 5
    5 4
    1 5
    3 4
    

    样例输出

    9
    分析:Tarjan+差分思想;
       难点在于怎么处理树上的区间加减问题;
       把数列差分思想用到树上,区间[i~j]加上k,等价于a[i]+=k,a[j+1]-=k,ans[t]=Σa[p](1=<p<=t);
       树上的s,t区间分成两个,设s,t公共祖先为p,一个是s到p的儿子,另一个是t到p,两个区间都差分一下;
       最后dfs递归处理类似前缀和的答案;
    代码:
    #include <iostream>
    #include <cstdio>
    #include <cstdlib>
    #include <cmath>
    #include <algorithm>
    #include <climits>
    #include <cstring>
    #include <string>
    #include <set>
    #include <map>
    #include <queue>
    #include <stack>
    #include <vector>
    #include <list>
    #define rep(i,m,n) for(i=m;i<=n;i++)
    #define rsp(it,s) for(set<int>::iterator it=s.begin();it!=s.end();it++)
    #define mod 1000000007
    #define inf 0x3f3f3f3f
    #define vi vector<int>
    #define pb push_back
    #define mp make_pair
    #define fi first
    #define se second
    #define ll long long
    #define pi acos(-1.0)
    #define pii pair<int,int>
    const int maxn=1e5+10;
    const int dis[4][2]={{0,1},{-1,0},{0,-1},{1,0}};
    using namespace std;
    ll gcd(ll p,ll q){return q==0?p:gcd(q,p%q);}
    ll qpow(ll p,ll q){ll f=1;while(q){if(q&1)f=f*p%mod;p=p*p%mod;q>>=1;}return f;}
    int n,m,k,t,fa[maxn],ans[maxn],vis[maxn],faa[maxn],ma;
    vi a[maxn],query[maxn];
    int find(int x)
    {
        return faa[x]==x?x:faa[x]=find(faa[x]);
    }
    void dfs(int now,int pre)
    {
        for(int x:a[now])
        {
            if(x!=pre)
            {
                dfs(x,now);
                fa[x]=now;
            }
        }
    }
    void dfs1(int now)
    {
        vis[now]=1;
        for(int x:query[now])
        {
            if(vis[x])
            {
                int p=find(x);
                ans[p]--,ans[fa[p]]--;
            }
        }
        for(int x:a[now])
        {
            if(!vis[x])
            {
                dfs1(x);
                faa[x]=now;
            }
        }
    }
    void dfs2(int now,int pre)
    {
        for(int x:a[now])
        {
            if(x!=pre)
            {
                dfs2(x,now);
                ans[now]+=ans[x];
            }
        }
        ma=max(ma,ans[now]);
    }
    int main()
    {
        int i,j;
        scanf("%d%d",&n,&k);
        rep(i,1,n)faa[i]=i;
        rep(i,1,n-1)
        {
            scanf("%d%d",&j,&t);
            a[j].pb(t),a[t].pb(j);
        }
        rep(i,1,k)
        {
            scanf("%d%d",&j,&t);
            ans[j]++,ans[t]++;
            query[j].pb(t);
            query[t].pb(j);
        }
        dfs(1,-1);
        dfs1(1);
        dfs2(1,-1);
        printf("%d
    ",ma);
        //system("pause");
        return 0;
    }
  • 相关阅读:
    圣杯布局总结
    新手用pyCharm编辑器创建项目不知道如何选择(适用于遇到同样问题的新手)
    支付宝小程序日期选择组件datePicker封装
    关于 WebView 的一些笔记
    ES6学习之箭头函数
    MongoDB--在windows下的安装过程及基本配置
    使用HTML5的canvas做图片剪裁
    git常见操作和常见错误
    算法——找出缺失的整数
    入坑涨姿势
  • 原文地址:https://www.cnblogs.com/dyzll/p/5770722.html
Copyright © 2011-2022 走看看