/*
匈牙利算法
二分匹配 最小点覆盖=最大匹配。
即踢掉最小点覆盖
*/
#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<algorithm>
#include<iostream>
#include<queue>
#include<map>
#include<math.h>
using namespace std;
typedef long long ll;
//typedef __int64 int64;
const int maxn = 1015;
const int inf = 0x7fffffff;
const double pi=acos(-1.0);
const double eps = 1e-8;
struct Edge{
int v,next;
}edge[ maxn<<4 ];
struct Node{
int x,y;
}sum[ maxn<<1 ];
int cnt,head[ maxn<<1 ];
int vis[ maxn<<1 ],mylink[ maxn<<1 ];
void init(){
cnt = 0;
memset( head,-1,sizeof( head ) );
}
void addedge( int a,int b ){
edge[ cnt ].v = b;
edge[ cnt ].next = head[ a ];
head[ a ] = cnt++;
}
bool km( int u ){
for( int i=head[ u ];i!=-1;i=edge[ i ].next ){
int v = edge[ i ].v;
if( vis[ v ]==0 ){
vis[ v ] = 1;
if( mylink[ v ]==-1||km( mylink[v] ) ){
mylink[ v ] = u;
return true;
}
}
}
return false;
}
int solve( int n,int m ){
memset( mylink,-1,sizeof( mylink ) );
int ans = 0;
for( int i=1;i<=n;i++ ){
memset( vis,0,sizeof( vis ) );
if( km(i) )
ans++;
}
return ans;
}
bool Judge( int i,int j ){
int x1 = sum[i].x;
int y1 = sum[i].y;
int x2 = sum[j].x;
int y2 = sum[j].y;
if( x1==x2 ){
if( y1==y2 ) return true;
if( y1==y2+1 ) return true;
}
if( x1+1==x2 ){
if( y1==y2 ) return true;
if( y1==y2+1 ) return true;
}
return false;
}
int main(){
int n,m;
while( scanf("%d%d",&n,&m)==2,n+m ){
init();
int cc = 1;
for( int i=1;i<=n;i++ )
scanf("%d%d",&sum[cc].x,&sum[cc].y),cc++;
for( int j=1;j<=m;j++ )
scanf("%d%d",&sum[cc].x,&sum[cc].y),cc++;
for( int i=1;i<=n;i++ ){
for( int j=n+1;j<cc;j++ ){
if( Judge( i,j ) ){
addedge( i,j );
}
}
}
printf("%d\n",cc-1-solve( n,m ));
}
return 0;
}
HDU4619+匈牙利
最新推荐文章于 2021-09-12 13:26:16 发布