poj 1486 二分图最大匹配必须边

题意:就是求点和矩形的最大匹配。如果有这种唯一的对应关系则输出，没有的话输出none

思路：点和矩阵建立连接，然后求最大匹配。同时注意删边的技巧。

#include<iostream>
#include<cmath>
using namespace std;
struct point
{
  int x,y;
}P[100];
struct rec
{
  int minx;
  int miny;
  int maxx;
  int maxy;
}R[100];
bool g[100][100];
bool vis[100];
int link[100];
int N;
bool check(int i,int j)
{
  int xx=P[i].x;
  int yy=P[i].y;
  if(xx>R[j].minx&&xx<R[j].maxx&&yy>R[j].miny&&yy<R[j].maxy)
  return true;
  return false;
}
bool dfs(int u)
{
  for(int i=1;i<=N;i++)
  {
    if(!vis[i]&&g[u][i])
    {
      vis[i]=1;
      if(link[i]==-1||dfs(link[i]))
      {
        link[i]=u;
        return true;
      }
    }
  }
  return false;
}
int maxmatch()
{
  int num=0;
  memset(link,-1,sizeof(link));
  for(int i=1;i<=N;i++)
  {
    memset(vis,0,sizeof(vis));
    if(dfs(i))  num++;
  }
  return num;
}
int main()
{
  int h=1;
  while(scanf("%d",&N)!=EOF&&N)
  {
    for(int i=1;i<=N;i++)
    scanf("%d%d%d%d",&R[i].minx,&R[i].maxx,&R[i].miny,&R[i].maxy);
    for(int i=1;i<=N;i++)
    scanf("%d%d",&P[i].x,&P[i].y);
    memset(g,0,sizeof(g));
    for(int i=1;i<=N;i++)
      for(int j=1;j<=N;j++)
       {
         if(check(i,j))
         {
           g[i][j]=1;
           //printf("(%d,%d)\n",i,j);
         }
       }
       maxmatch();//先求出最大匹配
       //枚举边 i代表y图的节点 也就是矩形的标号
       bool flag=false;
       bool end=false;
       for(int i=1;i<=N;i++)
       {
         if(link[i]==-1)  continue;//不一定都是完美匹配
         int temp=link[i];
         link[i]=-1;//删边
         g[temp][i]=0;
         memset(vis,0,sizeof(vis));
         if(!dfs(temp))//从删掉这个边的x图节点temp开始找，看是否存在增广路，如果不存在则是必须边。 
         {
           if(!end)
           printf("Heap %d\n",h++);
           flag=true;
           if(!end)
           printf("(%c,%d)",'A'+i-1,temp);
           else
           printf(" (%c,%d)",'A'+i-1,temp);
           end=true;
           link[i]=temp;//如果是必须边的话要把link还原 
         }
         //link[i]=temp;
           g[temp][i]=1;//不是必须边的话，link已经改变不需要还原，只要把图g还原就可以 
       }
       if(!flag)
       {
         printf("Heap %d\n",h++);
         printf("none\n\n");
       }
       else
       printf("\n\n");
  }
  //system("pause");
}