对于这道题嘛,先把这道题转化成二分图匹配,左面是图,右面是节点,一图和一节点进行匹配,但这题还有点特别,看看样例二就知道了,(我比较菜,哈哈,有的我也解释不明白,凭感觉,奥利给)我们必须得去了解一个东西,必须边,就是如果没有这条边(就是删除这条边),匹配数就达不到没删除这条边时的匹配数,就说这条边不是必须边,就是可有可无,想要更好的理解这个概念,我们可以看下第二个样例
额。。不会作图啊,凑合看吧,左边代表的是图,右面代表的是节点,由于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;
}
哦了,写的不好哈,理解,理解,我祝你幸福,奥利给。