HDU3468 Treasure Hunting(中等) [最大流](二分匹配)+最短路

题意:在一张寻宝图上，走最短路从一个汇合点到下一个汇合点，如果到汇合点上的路径有金币的话，只能拿一个，拿过后这个点就变为了可以通行的没有金币的点，问最多能拿几个金币。（汇合点按顺序编号是A~Z||a~z的点）
思路：最大流建图挺难的，参考网上的思路，用BFS预处理所有的汇合点到其他点的最短路径，然后找到在当前汇合点和下一个汇合点之间的金币，把当前的汇点与这个金币连一条边权是1的边。最后源点与每一个汇点连边，每一个金币的点与汇点连边。
跑一遍最大流即可。实现起来还是有很多的细节要注意。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<queue>
#include<cstring>
#include<string>
using namespace std;
#define pb push_back
#define cl(a,b) memset(a,b,sizeof(a))
const int maxn=105;
const int inf=1<<28;

template<int nv,int ne>
struct isap{
    int n,size;
    int head[nv];
    int dis[nv],gap[nv],cur[nv],pre[nv];
    int maxflow;

    struct edge{
        int v,w,next;
        edge(){}
        edge(int _v,int _w,int _next):v(_v),w(_w),next(_next){}
    }E[ne];

    void init(int n){
        this->n=n,size=0;
        cl(head,-1);
    }
    void insert(int u,int v,int w){
        E[size]=edge(v,w,head[u]);
        head[u]=size++;
        E[size]=edge(u,0,head[v]);
        head[v]=size++;
    }
    int maxFlow(int src,int des){
        maxflow=0;
        for(int i=0;i<=n;i++){
            dis[i]=gap[i]=0;
            cur[i]=head[i];
        }
        int u=pre[src]=src;
        int aug=0;///or aug=-1
        while(dis[src]<n){
           loop:for(int &i=cur[u];i!=-1;i=E[i].next){
                int v=E[i].v;
                if(E[i].w&&dis[u]==dis[v]+1){
                    aug=min(aug,E[i].w);
                    pre[v]=u;
                    u=v;
                    if(v==des){
                        maxflow+=aug;
                        for(u=pre[u];v!=src;v=u,u=pre[u]){
                            E[cur[u]].w-=aug;
                            E[cur[u]^1].w+=aug;
                        }
                        aug=inf;
                    }
                    goto loop;
                }
            }
            int mdis=n;
            for(int i=head[u];i!=-1;i=E[i].next){
                int v=E[i].v;
                if(E[i].w&&mdis>dis[v]){
                    cur[u]=i;
                    mdis=dis[v];
                }
            }
            if(--gap[dis[u]]==0)break;
            gap[dis[u]=mdis+1]++;
            u=pre[u];
        }
        return maxflow;
    }
};
isap<10005,200005> G;

int n,m;
char a[maxn][maxn];
int dist[60][maxn*maxn];//计算每一个汇合点到其他点的距离
bool vis[maxn][maxn];
int pos[maxn*maxn];//汇合点的位置
int gold[maxn*maxn];//有金币的位置
int dir[4][2]={0,1,1,0,-1,0,0,-1};
int charToInt(char x){
    if(x>='a'&&x<='z')return x-'a'+26;
    return x-'A';
}
void bfs(int s){
    queue<int> q;
    cl(vis,false);
    cl(dist[s],0x3f3f3f3f);

    q.push(pos[s]);
    vis[pos[s]/m][pos[s]%m]=true;
    dist[s][pos[s]]=0;

    while(!q.empty()){
        int u=q.front();q.pop();
        int x=u/m;
        int y=u%m;
        for(int i=0;i<4;i++){
            int xx=dir[i][0]+x;
            int yy=dir[i][1]+y;
            if(xx<0||xx>=n||yy<0||yy>=m)continue;
            if(vis[xx][yy]||a[xx][yy]=='#')continue;
            vis[xx][yy]=true;
            dist[s][xx*m+yy]=dist[s][x*m+y]+1;
            q.push(xx*m+yy);
        }

    }
}
int main(){
    while(~scanf("%d%d",&n,&m)){
        int cnt1=0,cnt2=0;
        for(int i=0;i<n;i++){
            scanf("%s",a[i]);
            for(int j=0;j<m;j++){
                if(isalpha(a[i][j])){
                    pos[charToInt(a[i][j])]=i*m+j;
                    cnt1++;//汇合点的个数
                }
                else if(a[i][j]=='*'){
                    gold[cnt2++]=i*m+j;//金币的个数
                }
            }
        }
        G.init(cnt1+cnt2);
        for(int i=0;i<cnt1;i++)bfs(i);
        bool flag=true;
        for(int i=1;i<cnt1;i++){
            for(int j=0;j<cnt2;j++){
                if(dist[i][gold[j]]+dist[i-1][gold[j]]==dist[i-1][pos[i]]){//金币在最短路上
                    G.insert(i,cnt1+j,1);
                }
                if(dist[i-1][pos[i]]==0x3f3f3f3f){
                    flag=false;
                }
            }
        }
        if(!flag){
            puts("-1");continue;
        }
        for(int i=1;i<cnt1;i++){
            G.insert(0,i,1);
        }
        for(int i=0;i<cnt2;i++){
            G.insert(i+cnt1,cnt1+cnt2,1);
        }
        printf("%d\n",G.maxFlow(0,cnt1+cnt2));
    }
    return 0;
}