BZOJ2561 最小生成树 网络流

题意:初始时给定一张图,然后往图上不断加边,每次加边后询问最少需要删除多少条边,使得新加入的边可以出现在最小和最大生成树上

题解:

如果u,v之间不存在路径,直接塞进去就好,着重讨论一下之前已经存在路径的情况。(以下以最小生成树为例,最大生成树随便YY一下就出来了- -)

显然加入新的边之后会形成环(不一定只有一个,可能之前有多个u到v的路径),而只要环上存在比L小的边,那么就不会将新加入的边加入MST中,所以就要删除生成的环上尽量少的长度大于L的边,使得不存在环(也就不存在除L之外的从u到v的路径)。这样就好办了,将所有长度小于L的边建网络流,跑从u到v的最小割就可以了。

显然长度大于L的边和长度小于L的边互不干扰,跑两边Dinic把答案加起来就好。

#include <queue>
#include <vector>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <climits>
#include <iostream>
#include <algorithm>
using namespace std;

const int MAXN=600000+2;
struct HASH{
    int u,w;
    HASH *next;
    HASH(){}
    HASH(int _u,int _w,HASH *_next):u(_u),w(_w),next(_next){}
}*table[MAXN],mem[MAXN];
struct EDGE{
    int u,c;
    EDGE(){}
    EDGE(int _u,int _c):u(_u),c(_c){}
}e[MAXN];
int N,M,U,V,L,cnt,cur[MAXN],d[MAXN];
queue<int> q;
vector<int> tab[MAXN];

void Insert1(int u,int v,int w){ table[u]=&(mem[cnt++]=HASH(v,w,table[u]));}

void Insert2(int u,int v,int c){
    tab[u].push_back(cnt),e[cnt++]=EDGE(v,c);
    tab[v].push_back(cnt),e[cnt++]=EDGE(u,c);
}

bool BFS(int s,int t){
    memset(d,-1,sizeof(d));
    d[s]=0,q.push(s);

    int x;
    while(!q.empty()){
        x=q.front(),q.pop();
        for(int i=0;i<tab[x].size();i++)
            if(e[tab[x][i]].c && d[e[tab[x][i]].u]==-1)
                d[e[tab[x][i]].u]=d[x]+1,q.push(e[tab[x][i]].u);
    }
    return d[t]>0;
}

int DFS(int x,int f,int t){
    if(x==t) return f;

    int flow,used=0;
    for(int i=cur[x];i<tab[x].size();i++)
        if(e[tab[x][i]].c && d[e[tab[x][i]].u]==d[x]+1){
            flow=DFS(e[tab[x][i]].u,min(f-used,e[tab[x][i]].c),t);
            e[tab[x][i]].c-=flow,e[tab[x][i]^1].c+=flow,used+=flow;
            if(e[tab[x][i]].c) cur[x]=i;
            if(used==f) return f;
        }

    if(!used) d[x]=-1;
    return used;
}

int Dinic(int s,int t){
    int ret=0;
    while(BFS(s,t)){
        memset(cur,0,sizeof(cur));
        ret+=DFS(s,INT_MAX,t);
    }
    return ret;
}

int Query(int u,int v,int l,int r){
    cnt=0;
    for(int i=1;i<=N;i++) tab[i].clear();

    for(int i=1;i<=N;i++)
        for(HASH *p=table[i];p;p=p->next)
            if(l<p->w && p->w<r) Insert2(i,p->u,1);

    return Dinic(u,v);
}

int main(){
    scanf("%d %d",&N,&M);
    for(int i=1,u,v,w;i<=M;i++){
        scanf("%d %d %d",&u,&v,&w);
        Insert1(u,v,w);
    }
    scanf("%d %d %d",&U,&V,&L);

    printf("%d\n",Query(U,V,-1,L)+Query(U,V,L,INT_MAX));

    return 0;
}
View Code

 

转载于:https://www.cnblogs.com/WDZRMPCBIT/p/6481308.html

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值