POJ 2226 Muddy Fields（二分匹配【最大流】）

存档纪念
存档的原因是WA了很多次 然后试了很多样例终于试出来了
有很多要注意的地方 等以后详细说
参考样例在最后

#include <iostream>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <cmath>
#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;

//typedef pair<int, int> P;

/*
    二分图最大匹配最大流简化做法，复杂度O(|V||E|)
*/

const int inf = 2147483647;
const double eps = 1e-8;
const int maxn = 2500;
const int N = 55;
const ll mod = 1e9+7;

typedef pair<int,int> pa;
char map[N][N];
pa save[N][N];
int r, c;
int p, q;
vector<int> G[maxn];
int match[maxn];        //匹配 
bool used[maxn];

void addedge(int u, int v){
    G[u].push_back(v);
    G[v].push_back(u);
}

//dfs寻找增广路 
bool dfs(int v){
    used[v] = true;
    for(int i=0; i<G[v].size(); i++){
        int u = G[v][i], w = match[u];
        if(w<0 || !used[w] && dfs(w)){  //u未匹配过，或还有增广路 
            match[v] = u;
            match[u] = v;
            return true;
        }
    }
    return false;
}

//求解二分图的最大匹配 
int bipartite_matching(){
    int res = 0;
    memset(match, -1, sizeof(match));
    for(int v = 0; v< p + q; v++){
        if(match[v] < 0){
            memset(used, 0, sizeof(used));
            if(dfs(v)){
                res++;
            }
        }
    } 
    return res;
}

int main(){
    while(~scanf("%d %d", &r, &c)){
        for(int i=0; i<maxn; i++){
            G[i].clear();
            used[i] = 0;

        }
        memset(save, -1, sizeof(save));
        memset(map, 0, sizeof(map));
        for(int i = 0; i < r; i++){
            scanf("%s",map[i]);
        } 
        p = 0;
        int flag = 0;
        for(int i = 0; i < c; i++){
            flag = 0;
            for(int j = 0; j < r; j++){
                if(map[j][i] == '*'){
                    save[j][i].first = p;
                    flag = 1;
                }
                else{
                    save[j][i].first = -1;
                    if(flag) p++;
                    flag = 0;

                }
            }
            if(flag) p++;
        }

        q = 0;
        flag = 0;
        for(int i = 0; i < r; i++){
            flag = 0;
            for(int j = 0; j < c; j++){
                if(map[i][j] == '*'){
                    save[i][j].second = q;
                    flag = 1;
                }
                else{
                    save[i][j].second = -1;
                    if(flag) q++;
                    flag = 0;
                }
            }
            if(flag) q++;
        }

//      
//      for(int i = 0; i < r; i++){
//          for(int j = 0; j < c; j++){
//              cout<<save[i][j].first<<","<<save[i][j].second<<"  ";
//          }
//          puts("");
//      }
//      
        for(int i = 0; i < r; i++){
            for(int j = 0; j < c; j++){
                if(save[i][j].first!=-1&&save[i][j].second!=-1){
//                  system("pause");
//                  cout<<"i: "<<i<<" j: "<<j<<endl;
//                  cout<<save[i][j].first<<" "<<save[i][j].second<<endl;
                    addedge(save[i][j].first,p + save[i][j].second );
                }
            }   
        }

//      int ans = bipartite_matching();


//      for(int i=0; i<n; i++){
//          int m;
//          scanf("%d", &m);
//          for(int j=0; j<m; j++){
//              int b;
//              scanf("%d", &b);
//              addedge(i, n+b-1);
//          }
//      }
        printf("%d\n", bipartite_matching());
    }
    return 0;
} 
/*
3 4
.***
.***
.***
5 5
.*.*.
*.*.*
.*.*.
*.*.*
.*.*.
16 16
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
.***.***.***.***
50 50
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
.***.***.***.***.***.***.***.***.***.***.***.*****
50 50
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*
*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.*.
*/

小数据和大数据的样例是必要的！
一定要注意数组长度&用模板时总点数的范围