BZOJ P1001 [BJOI2006] 狼抓兔子【最大流】【挖坑】

很久没有写网络流…现在写这道最小割板题居然调了两小时…很难受啊

连边然后直接上网络流就可以了….没什么思维难度和代码难度(?你明明写两小时)…

看了其他题解说这道题可以用什么叫做对偶图然后跑最短路?这道题最开始确实有写最短路的想法,但是完全不知道最短路应该如何连边。看了网上的其他题解,说有一个东西叫做对偶图,感觉很厉害,所以挖一个坑,找个时间去了解一下对偶图,顺便用最短路做一下。

关键代码:

bool BFS(){
    int I,J,K;
    memset(Dis,0x3f3f3f3f,sizeof(Dis));
    queue<int>Q;Q.push(Start);Dis[Start]=1;int Inf=Dis[End];
    while(Q.size()!=0){
        int X=Q.front();Q.pop();
        for(I=Head[X];I;I=Next[I]){
            if(Edge[I]&&Dis[To[I]]==Inf){
                Dis[To[I]]=Dis[X]+1;Q.push(To[I]);
            }
        }
        if(Dis[End]!=Inf){
            return true;
        }
    }
    return false;
}
int MaxFlow(int Cur,int Tmp){
    int I,J,K,Flow=0;
    if(Cur==End||Tmp==0){
        return Tmp;
    }
    for(I=Head[Cur];I;I=Next[I]){
        if(Dis[To[I]]==Dis[Cur]+1&&Edge[I]){
            int _Flow=MaxFlow(To[I],min(Tmp,Edge[I]));
            if(_Flow){
                Edge[I]-=_Flow;Edge[I^1]+=_Flow;
                Flow+=_Flow;Tmp-=_Flow;
                if(Tmp==0){
                    return Flow;
                } 
            } else {
                Dis[To[I]]=-1;
            }
        }
    }
    return Flow;
}
int main(){
    int I,J,K;
    scanf("%d%d",&N,&M);
    for(I=1;I<=N;I++){
        for(J=1;J<M;J++){
            int X=(I-1)*M+J,Y=X+1,Z;scanf("%d",&Z);Insert(X,Y,Z);Insert(Y,X,Z);
        }
    }
    for(I=1;I<N;I++){
        for(J=1;J<=M;J++){
            int X=(I-1)*M+J,Y=X+M,Z;scanf("%d",&Z);Insert(X,Y,Z);Insert(Y,X,Z);
        }
    }
    for(I=1;I<N;I++){
        for(J=1;J<M;J++){
            int X=(I-1)*M+J,Y=X+M+1,Z;scanf("%d",&Z);Insert(X,Y,Z);Insert(Y,X,Z);
        }
    }
    ...
    printf("%d",Ans);
    return 0;
}

猜你喜欢

转载自blog.csdn.net/yanzhenhuai/article/details/81160832