HDOJ 3469 - Treasure Hunting BFS+二分图最大匹配...深入理解二分图最大匹配..

                   题意:

                             iSea和他朋友一起来到一个迷宫中...每次他们会从一个集合点到达另一个集合点(集合点按顺序A~Z~a~z..最多52个)...在每两个集合点间必须走最短的路径(当然可能有多条)..iSea可以最多在两个集合点间的最短路径上拿一个宝藏..每个宝藏只能被拿一次....问能否完成所有的路径..不行输出-1..否则输出iSea能拿到的最多的宝藏个数..

                   题解:

                             可见在最终的答案中..一个宝藏至多属于一条线段上...于是想到用宝藏和线段来做匹配...宝藏作为一侧..两两集合点的最短路径作为另一侧...来做最大匹配..那么关键就是构边了..什么时候某个宝藏可以和一条路劲相连呢?这个宝藏必定属于某两个集合点的最短路径上...所以枚举集合点,,每次用BFS找出图中点到当前起点的最短距离..图中点到当前终点的最短距离..若有一个宝藏的两者之和等于起点到终点的最短距离..那么这个宝藏就与这条路径的标号做边...

                   总结:

                             这类问题的特点是一个集合中的某个元素至多匹配另一个集合中的一个元素...问最多的数量...其实也很裸了...


Program:

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<cmath>
#include<stack>
#include<string.h>
#include<queue>
#define ll long long
#define esp 1e-5
#define MAXN 10005
#define MAXM 5000005
#define oo 1000007
using namespace std;     
struct node
{
       int y,next;
}line[MAXM];
int Lnum,_next[MAXN],match[MAXN],way[4][2]={{1,0},{-1,0},{0,1},{0,-1}},dis[2][MAXN];
bool had[105],used[MAXN];
char s[105][105];
void addline(int x,int y)
{
       line[++Lnum].next=_next[x],_next[x]=Lnum,line[Lnum].y=y;
}
int BFS(int t,int R,int C,int tp)
{
       int x,y,k,h;
       char c;
       queue<int> Q;
       memset(dis[tp],0x2f,sizeof(dis[tp]));
       while (!Q.empty()) Q.pop();   
       if (t<26) c='A'+t;
           else  c='a'+t-26;
       for (y=0;y<R;y++)
       {
             for (x=0;x<C;x++)
                 if (s[y][x]==c) break;
             if (x!=C) break;
       }
       dis[tp][y*C+x]=0;
       h=y*C+x;
       Q.push(y*C+x);
       while (!Q.empty())
       {
              x=Q.front();
              Q.pop();
              y=x/C,x%=C; 
              for (k=0;k<4;k++)
              {
                     int yy=y+way[k][0],xx=x+way[k][1];
                     if (yy<0 || xx<0 || yy>=R || xx>=C) continue;
                     if (s[yy][xx]=='#' || dis[tp][yy*C+xx]<oo) continue; 
                     Q.push(yy*C+xx),dis[tp][yy*C+xx]=dis[tp][y*C+x]+1;  
              }
       }
       return h;
}
bool build(int R,int C,int num)
{
       int i,j,x,y,t;
       for (i=0;i<num;i++)
          if (!had[i]) return false;
       Lnum=0,memset(_next,0,sizeof(_next)); 
       for (t=0;t<num-1;t++)
       {
               int p1=BFS(t,R,C,0),p2=BFS(t+1,R,C,1); 
               if (dis[0][p2]>oo) return false;
               for (y=0;y<R;y++) for (x=0;x<C;x++) 
                  if (s[y][x]=='*' && dis[0][y*C+x]+dis[1][y*C+x]==dis[0][p2])  
                     addline(t+1,y*C+x);
       }
       return true;
}
bool dfs(int x)
{
       for (int k=_next[x];k;k=line[k].next)
       {
               int y=line[k].y;
               if (used[y]) continue;
               used[y]=true;
               if (!match[y] || dfs(match[y]))
               {
                       match[y]=x;
                       return true;
               }
       }
       return false; 
}
int getmax(int n)
{
       int sum=0;
       memset(match,0,sizeof(match));
       for (int i=1;i<=n;i++)
       {
               memset(used,false,sizeof(used));
               sum+=dfs(i);
       }
       return sum;
}
int main()
{             
       int R,C,num,i,j,x;
       freopen("input.txt","r",stdin);
       freopen("output.txt","w",stdout); 
       while (~scanf("%d%d",&R,&C))
       {
                for (i=0;i<R;i++) scanf("%s",s[i]);
                num=0;
                memset(had,false,sizeof(had));
                for (i=0;i<R;i++)
                   for (j=0;j<C;j++)
                      if (s[i][j]>='A' && s[i][j]<='Z' || s[i][j]>='a' && s[i][j]<='z')
                      { 
                             if (s[i][j]>='A' && s[i][j]<='Z') x=s[i][j]-'A';
                                                         else  x=s[i][j]-'a'+26;
                             had[x]=true;
                             num=max(num,x+1);
                      }
                if (!build(R,C,num)) 
                {
                       printf("-1\n");
                       continue;
                }
                printf("%d\n",getmax(num-1));
       }
       return 0;  
}   


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值