BZOJ[5070]危险的迷宫 费用流

传送门ber~

同Codevs[1227]方格取数2

每个点拆成入和出，入向出连容量为1费用为代价的边
相邻(可以互相走到)的点出向入连容量为1费用为0的边
源向起点的入连容量为1费用为0的边
出口得出向汇连容量为1费用为0的边

跑最小费用最大流

判断是否有解看是否满流就可以了

代码如下：

#include<algorithm>
#include<cstring>
#include<ctype.h>
#include<cstdio>
#include<queue>
#define INF 2147483647
using namespace std;
inline int read(){
    int x=0,f=1;char c;
    do c=getchar(),f=c=='-'?-1:f; while(!isdigit(c));
    do x=(x<<3)+(x<<1)+c-'0',c=getchar(); while(isdigit(c));
    return x*f;
}
const int S=501;
const int T=502;
queue<int>q;
int n,m,k,x,y,x1,y1,cnt,flow,cost,top=1;
int fir[505],d[505];
bool b[505];
struct Point{
    int s,t;
}a[15][15];
struct Edge{
    int to,nex,k,v;
    Edge(int _to=0,int _nex=0,int _k=0,int _v=0):to(_to),nex(_nex),k(_k),v(_v){}
}nex[5020];
inline void add(int x,int y,int k,int v){
    nex[++top]=Edge(y,fir[x],k,v);
    fir[x]=top;
}
inline bool spfa(){
    for(int i=1;i<=502;i++)
        d[i]=INF,b[i]=false;
    d[S]=0;
    q.push(S);
    while(!q.empty()){
        int x=q.front();q.pop();
        b[x]=false;
        for(int i=fir[x];i;i=nex[i].nex)
            if(nex[i].k && d[nex[i].to]>d[x]+nex[i].v){
                d[nex[i].to]=d[x]+nex[i].v;
                if(!b[nex[i].to]) q.push(nex[i].to),b[nex[i].to]=true;
            }
    }
    return d[T]!=INF;
}
int dfs(int x,int v){
    if(!v || x==T){
        cost+=v*d[T];
        return v;
    }
    int tmp=0;
    b[x]=true;
    for(int i=fir[x];i;i=nex[i].nex)
        if(nex[i].k && d[nex[i].to]==d[x]+nex[i].v && !b[nex[i].to]){
            int f=dfs(nex[i].to,min(v,nex[i].k));
            nex[i].k-=f;nex[i^1].k+=f;tmp+=f;v-=f;
            if(!v) break;
        }
    if(!tmp) d[x]=INF;
    return tmp;
}
inline void dinic(int &f){
    while(spfa()) f=f+dfs(S,INF);
    return;
}
int main(){
    n=read();m=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            a[i][j].s=++cnt;a[i][j].t=++cnt;
            k=read();
            add(a[i][j].s,a[i][j].t,1,k);
            add(a[i][j].t,a[i][j].s,0,-k);
        }
    k=read();
    for(int i=1;i<=k;i++){
        x=read();y=read();x1=read();y1=read();
        add(a[x][y].t,a[x1][y1].s,1,0);add(a[x1][y1].s,a[x][y].t,0,0);
        add(a[x1][y1].t,a[x][y].s,1,0);add(a[x][y].s,a[x1][y1].t,0,0);
    }
    k=read();
    for(int i=1;i<=k;i++){
        x=read();y=read();
        add(S,a[x][y].s,1,0);
        add(a[x][y].s,S,0,0);
    }
    for(int i=1;i<=k;i++){
        x=read();y=read();
        add(a[x][y].t,T,1,0);
        add(T,a[x][y].t,0,0);
    }
    dinic(flow);
    printf("%d",flow==k?cost:-1);
return 0;
}