[BJOI2006][bzoj1001] 狼抓兔子 [最小割]

题面：

传送门

思路：

其实就是一道最小割的题目......

我的写法加了两个优化，常数比较小，所以过掉了

一个是当前弧，一个是若当前点并不能流出去，那么标记dep为-1

听说正解是对偶图最短路？可以找时间学一学......

Code：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#define inf 1e9
#define id(i,j) (i-1)*m+j
using namespace std;
inline int read(){
    int re=,flag=;char ch=getchar();
    while(ch<''||ch>''){
        if(ch=='-') flag=-;
        ch=getchar();
    }
    while(ch>=''&&ch<='') re=(re<<)+(re<<)+ch-'',ch=getchar();
    return re*flag;
}
int n,m,cnt=-,ans,first[],dep[],cur[];
int q[],head,tail;
struct edge{
    int to,next,cap;
}a[];
inline void add(int u,int v,int cap){
    a[++cnt]=(edge){v,first[u],cap};first[u]=cnt;
    a[++cnt]=(edge){u,first[v],cap};first[v]=cnt;
}
bool bfs(int s,int t){
    int i,u,v;
    head=,tail=;q[]=s;
    for(i=s;i<=t;i++) dep[i]=-,cur[i]=first[i];
    dep[s]=;
    while(head<tail){
        u=q[head++];
        for(i=first[u];~i;i=a[i].next){
            v=a[i].to;
            if(~dep[v]||(!a[i].cap)) continue;
            dep[v]=dep[u]+;
            q[tail++]=v;
        }
    }
    return ~dep[t];
}
int dfs(int u,int t,int limit){
    if(u==t||!limit) return limit;
    int i,v,f,flow=;
    for(i=cur[u];~i;i=a[i].next){
        v=a[i].to;cur[u]=i;//优化一
        if(dep[v]==dep[u]+&&a[i].cap>){
            f=dfs(v,t,min(limit,a[i].cap));
            if(f){
                flow+=f;limit-=f;
                a[i].cap-=f;a[i^].cap+=f;
                if(!limit) return flow;
            }
            else{
                dep[v]=-;continue;//优化二
            }
        }
    }
    return flow;
}
void dinic(int s,int t){
    int tmp;
    while(bfs(s,t)){
        while(tmp=dfs(s,t,inf)) ans+=tmp;
    }
}
int main(){
    memset(first,-,sizeof(first));
    int i,j,t1;
    n=read();m=read();
    for(i=;i<=n;i++){
        for(j=;j<m;j++){
            t1=read();
            add(id(i,j),id(i,j+),t1);
        }
    }
    for(i=;i<n;i++){
        for(j=;j<=m;j++){
            t1=read();
            add(id(i,j),id(i+,j),t1);
        }
    }
    for(i=;i<n;i++){
        for(j=;j<m;j++){
            t1=read();
            add(id(i,j),id(i+,j+),t1);
        }
    }
    dinic(,n*m);
    printf("%d",ans);
}