hdu 1533 KM或费用流

以前用KM写过，现在再用费用流写。

 #include <iostream>
 #include <cstdio>
 #include <cstring>
 #include <queue>
 #include <vector>
 #include <utility>
 #define abs(a) ((a)<0?-(a):(a))
 #define maxn 210
 #define oo 0x3f3f3f3f
 using namespace std;

 struct Edge {
     int u, v, c, f;
     Edge( int u, int v, int c, int f ):u(u),v(v),c(c),f(f){}
 };
 struct Mcmf {
     int n, src, dst;
     vector<Edge> edge;
     vector<int> g[maxn];
     int dis[maxn], pth[maxn], ext[maxn];

     void init( int n, int src, int dst ) {
         this->n = n;
         this->src = src;
         this->dst = dst;
         for( int i=; i<=n; i++ )
             g[i].clear();
         edge.clear();
     }
     void add_edge( int u, int v, int c, int f ) {
         g[u].push_back( edge.size() );
         edge.push_back( Edge(u,v,c,f) );
         g[v].push_back( edge.size() );
         edge.push_back( Edge(v,u,-c,) );
     }
     bool spfa( int &flow, int &cost ) {
         queue<int> qu;
         memset( dis, 0x3f, sizeof(dis) );
         qu.push( src );
         ext[src] = true;
         dis[src] = ;
         pth[src] = -;
         while( !qu.empty() ) {
             int u=qu.front();
             qu.pop();
             ext[u] = false;
             for( int t=; t<g[u].size(); t++ ) {
                 Edge &e = edge[g[u][t]];
                 if( e.f && dis[e.v]>dis[e.u]+e.c ) {
                     dis[e.v] = dis[e.u]+e.c;
                     pth[e.v] = g[u][t];
                     if( !ext[e.v] ) {
                         ext[e.v] = true;
                         qu.push( e.v );
                     }
                 }
             }
         }
         if( dis[dst]==oo ) return false;
         int flw = oo;
         for( int eid=pth[dst]; eid!=-; eid=pth[edge[eid].u] )
             flw = min( flw, edge[eid].f );
         for( int eid=pth[dst]; eid!=-; eid=pth[edge[eid].u] ) {
             edge[eid].f -= flw;
             edge[eid^].f += flw;
         }
         flow += flw;
         cost += flw*dis[dst];
         return true;
     }
     void mcmf( int &flow, int &cost ) {
         flow = cost = ;
         while( spfa(flow,cost) );
     }
 };

 int n, m;
 int aa[maxn][], bb[maxn][], ta, tb;
 Mcmf M;

 int main() {
     ios::sync_with_stdio( false );
     while() {
         cin>>n>>m;
         if( n== && m== ) return ;
         ta = tb = ;
         for( int i=; i<=n; i++ )
             for( int j=; j<=m; j++ ) {
                 char ch;
                 cin>>ch;
                 if( ch=='m' ) {
                     ta++;
                     aa[ta][] = i;
                     aa[ta][] = j;
                 } else if( ch=='H' ) {
                     tb++;
                     bb[tb][] = i;
                     bb[tb][] = j;
                 }
             }
         M.init( ta+tb+, ta+tb+, ta+tb+ );
         for( int i=; i<=ta; i++ )
             M.add_edge( M.src, i, ,  );
         for( int j=ta+; j<=ta+tb; j++ )
             M.add_edge( j, M.dst, ,  );
         for( int i=; i<=ta; i++ )
             for( int j=ta+; j<=ta+tb; j++ ) {
                 int dis = abs(aa[i][]-bb[j-ta][])+abs(aa[i][]-bb[j-ta][]);
                 M.add_edge( i, j, dis,  );
             }
         int flow, cost;
         M.mcmf( flow, cost );
         cout<<cost<<endl;
     }
 }