POJ1486-二分图必须边

对于这道题嘛，先把这道题转化成二分图匹配，左面是图，右面是节点，一图和一节点进行匹配，但这题还有点特别，看看样例二就知道了，（我比较菜，哈哈，有的我也解释不明白，凭感觉，奥利给）我们必须得去了解一个东西，必须边，就是如果没有这条边（就是删除这条边），匹配数就达不到没删除这条边时的匹配数，就说这条边不是必须边，就是可有可无，想要更好的理解这个概念，我们可以看下第二个样例

额。。不会作图啊，凑合看吧，左边代表的是图，右面代表的是节点，由于1节点和2节点在图1中所以1图连接1节点和2节点，其余的也是同理，
我们把1-1这个匹配删除，发现最大匹配数不变，对吧，这就是传说中的非必须边，最后发现四个都是非必须边，所以结果是none。

那么怎么实现呢，先进行一遍匈牙利，然后一条边一条边删除，看看匹配数能不能达到删除前的匹配数，能的话就不输出，不能的话输出（就是必须边）

当然这题有个细节，如果最后存二分图的时候，（dp[i][j])如果i代表节点，j代表图（就是四边形），match[j]代表的就是相应的节点，如果i代表图，j代表节点，match[j]最后代表的是图，所以最好用第一种表示方法，我就错在这，我真了。

下面上代码
下面展示一些 内联代码片。

//POJ-1486
#include<cstdio> 
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
const int maxn = 50;

int N;
struct Node//存四边形 
{
	int x1,x2,y1,y2;
}Node[maxn];
struct point//存节点 
{
	int X,Y;
}Point[maxn];

int match[maxn];
int vis[maxn],girl[maxn];
int dp[maxn][maxn];//存储二分图 

bool judge(int x,int y)//在图里头就连上 
{
	if(Point[y].X >= Node[x].x1&&Point[y].X <= Node[x].x2)
	{
		if(Point[y].Y >= Node[x].y1&&Point[y].Y <= Node[x].y2)
		{
			return true;
		}
		return false;
	}
	return false;
}

bool dfs(int n) //匈牙利 
{
	for(int r = 1; r <= N; r++)
	{
		if(vis[r]||!dp[n][r]) continue;
		vis[r] = 1;
		if(girl[r] == -1||dfs(girl[r]))
		{
			girl[r] = n;
			return true;
		} 
	}
	return false;
}

bool make() //判断这个边删除之后可不可以达到之前的匹配数 
{
	int s = 0;
	memset(girl,-1,sizeof(girl));
	for(int u = 1; u <= N; u++)
	{
		memset(vis,0,sizeof(vis));
		if(dfs(u)) s++;
		else break;
	}
	return s == N;
}

void solve() //慢慢删 
{
	bool flag = true;
	for(int i = 1; i <= N; i++)
	{
		int t = match[i];
		dp[t][i] = 0;
		if(!make())
		{
			if(!flag) printf(" ");
			printf("(%c,%d)",'A'+i-1,t);
			flag = false;
		}
		else continue;
		dp[t][i] = 1;
	}
	
	if(flag) printf("none\n");
	else printf("\n");
}

int main()
{
	int kase = 0;
	while(scanf("%d",&N) == 1&&N != 0)
	{
		memset(dp,0,sizeof(dp));
		memset(Node,0,sizeof(Node));
		
		int X1,X2,Y1,Y2;
		for(int i = 1; i <= N; i++)
		{
			scanf("%d%d%d%d",&X1,&X2,&Y1,&Y2);
			Node[i].x1 = X1,Node[i].x2 = X2;
			Node[i].y1 = Y1,Node[i].y2 = Y2;
		}
		
		int A,B;
		for(int t = 1; t <= N; t++)
		{
			scanf("%d%d",&A,&B);
			Point[t].X = A,Point[t].Y = B;
		}
		
		for(int i = 1; i <= N; i++)
		{
			for(int j = 1; j <= N; j++)
			{
				if(judge(j,i)) dp[i][j] = 1;
			}
		}
		
		memset(girl,-1,sizeof(girl));
		int tool = 0;
		for(int t = 1; t <= N; t++)
		{
			memset(vis,0,sizeof(vis));
			if(dfs(t)) tool++;
		} 
		
		for(int u = 1; u <= N; u++)
		{
			match[u] = girl[u];
		}
		printf("Heap %d\n",++kase);
		solve();
		puts("");
	}
	return 0;
} 

哦了，写的不好哈，理解，理解，我祝你幸福，奥利给。