【最小割】Gym - 101982 - E - Cops And Roobers

题目链接<http://codeforces.com/gym/240837>

---

题意：

给你一张地图，里面有一个小偷。要求建造一些围栏把小偷圈住，有些地方不可造围栏，有些地方可造，但代价不同。问最少的代价是多少。

---

题解：

想到最小割就好做了。对于每个点拆成两个点，中间连一条边，使得点权变成边权，这样跑一遍就出来了。

注意n和m是反的。

---

#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<iostream>
using namespace std;
typedef long long ll;
const ll N=3e3+7;
const ll inf=1<<30;
struct Edge{
    ll u,v,w,nxt;
    Edge(ll u=0,ll v=0,ll w=0,ll nxt=0):u(u),v(v),w(w),nxt(nxt){}
}edge[N*4];
ll edn,sp,tp;
ll p[N],d[N],c[N];
void add(ll u,ll v,ll w){
    edge[++edn]=Edge(u,v,w,p[u]);p[u]=edn;
    edge[++edn]=Edge(v,u,0,p[v]);p[v]=edn;
}
bool bfs(){
    memset(d,-1,sizeof(d));d[sp]=0;
    queue<ll>q;q.push(sp);
    while(!q.empty()){
        ll u=q.front();
        q.pop();
        for(ll i=p[u];~i;i=edge[i].nxt){
            ll v=edge[i].v;
            if(d[v]==-1&&edge[i].w){
                d[v]=d[u]+1;
                q.push(v);
                if(v==tp) return true;
            }
        }
    }
    return ~d[tp];
}
ll dfs(ll u,ll b){
    if(u==tp) return b;
    ll r=0;
    for(ll i=c[u];~i;i=edge[i].nxt){
        ll v=edge[i].v;
        if(edge[i].w&&d[v]==d[u]+1){
            ll x=min(edge[i].w,b-r);
            c[u]=i;
            x=dfs(v,x);
            r+=x;
            edge[i].w-=x;
            edge[i^1].w+=x;
            if(r==b) break;
        }
    }
    if(!r)d[u]=-2;
    return r;
}
ll dinic(){
    ll total=0,t;
    while(bfs()){
        memcpy(c,p,sizeof(p));
        while(t=dfs(sp,inf))
        total+=t;
    }
    return total;
}
ll n,m,cc;
char s[50][50];
ll cs[50];
ll dir[4][2]={1,0,0,1,-1,0,0,-1};
int main()
{
    memset(p,-1,sizeof(p));edn=-1;
    scanf("%lld%lld%lld",&m,&n,&cc);
    for(ll i=0;i<n;i++){
        scanf("%s",s[i]);
    }
    for(ll i=0;i<cc;i++){
        scanf("%lld",&cs[i]);
    }
    sp=2*n*m+100,tp=sp+1;
    for(ll i=0;i<n;i++){
        for(ll j=0;j<m;j++){
            ll u=i*m+j;
            if(s[i][j]=='.'||s[i][j]=='B') add(u,u+n*m,inf);
            else add(u,u+n*m,cs[s[i][j]-'a']);
            for(ll k=0;k<4;k++){
                ll di=dir[k][0]+i;
                ll dj=dir[k][1]+j;
                if(di<0||di>=n||dj<0||dj>=m) add(u+n*m,tp,inf);
                else{
                    ll v=di*m+dj;
                    add(u+n*m,v,inf);
                }
            }
            if(s[i][j]=='B') add(sp,u,inf);
        }
    }
    ll ans=dinic();
    if(ans>=inf) printf("-1");
    else printf("%lld\n",ans);
    return 0;
}