POJ 3020 Antenna Placement 匈牙利算法,最大流解法 难度:1

http://poj.org/problem?id=3020

 #include <cstdio>
 #include <cstring>
 #include <vector>
 using namespace std;
 char maz[][];
 int e[][];
 int d[][]={,,,-,,,-,};
 int n,m;
 vector <int > G[];
 bool vis[];
 void addedge(int from,int to){
     e[from][to]=;
     G[from].push_back(to);
 }
 int dfs(int s){
     vis[s]=true;
     if(s==n*m+){//printf("dfs %d ok\n",s);
     return ;}
     for(int i=;i<G[s].size();i++){
         if(!vis[G[s][i]]&&e[s][G[s][i]]&&dfs(G[s][i])){
             e[s][G[s][i]]=;
             e[G[s][i]][s]=;
           //  printf("dfs %d ok\n",s);
             return ;
         }
     }
     //printf("dfs %d failed \n",s);
     return ;
 }
 int maxflow(){
     int ans=,f;
     while(){
         f=dfs(n*m);
         memset(vis,,sizeof(vis));
         if(f==)break;
         ans+=f;
     }
     return ans;
 }
 void printe(){
     for(int i=;i<m*n+;i++){
         bool f=false;
         for(int j=;j<m*n+;j++){
             if(e[i][j]){printf("e[%d][%d] ",i,j);f=true;}
         }
        if(f) printf("\n");
     }
 }
 int main(){
     int t;
     scanf("%d",&t);
     while((t--)&&scanf("%d%d",&n,&m)==){
         for(int i=;i<=n*n+;i++){G[i].clear();}
         memset(e,,sizeof(e));

         for(int i=;i<n;i++){
             scanf("%s",maz[i]);
         }
         int ans=;
         int star=;
         for(int i=;i<n;i++){
             for(int j=;j<m;j++){
                 if(maz[i][j]=='*'){
                     star++;
                     if((i+j)&){
                         addedge(n*m,i*m+j);//n*n s n*n+1 t
                        // printf("addedge %d %d %d %d:%d to %d\n",n,0,i,j,n*m,i*m+j);
                         for(int k=;k<;k++){
                             int nx=i+d[k][];
                             int ny=j+d[k][];
                             if(nx>=&&nx<n&&ny>=&&ny<m&&maz[nx][ny]=='*'){
                                 addedge(i*m+j,nx*m+ny);
                          //       printf("addedge %d %d %d %d:%d to %d\n",i,j,nx,ny,i*m+j,nx*m+ny);
                             }
                         }
                     }
                     else {
                         addedge(i*m+j,n*m+);
                         //printf("addedge %d %d %d %d:%d to %d\n",i,j,n,1,i*m+j,n*m+1);
                         for(int k=;k<;k++){
                             int nx=i+d[k][];
                             int ny=j+d[k][];
                             if(nx>=&&nx<n&&ny>=&&ny<m&&maz[nx][ny]=='*'){
                                 addedge(i*m+j,nx*m+ny);
                                 e[i*m+j][nx*m+ny]=;
                           //      printf("addedge %d %d %d %d:%d to %d\n",i,j,nx,ny,i*m+j,nx*m+ny);
                             }
                         }
                     }
                 }
             }
         }

         ans=maxflow();
         printf("%d\n",star-ans);
     }
     return ;
 }