POJ 2536 二分匹配 匈牙利算法 || 网络流 (一)

2014-11-24 02:42:07 · 作者: · 浏览: 4

题意很简单,不叙述了,我分别用二分匹配和最大流实现了一下这道题,算是对前面所学知识的一次综合运用吧。

这是二分匹配:


[cpp]
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define PI acos(-1.0)
#define Max 2005
#define inf 1<<28
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define REP(i,s,t) for(int i=(s);i<=(t);++i)
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
#define mp(a,b) make_pair(a,b)
using namespace std;

int Map[200][200] ;
int link[200] ;
bool vis[200] ;
int n , m ,s, v;
int dfs(int now)
{
for (int i = 1 ;i <= m ;i ++)
{
if(Map[now][i])
if(!vis[i])
{
vis[i] = 1 ;
if(link[i] == -1 || dfs(link[i]))
{
link[i] = now ;
return 1 ;
}
}
}
return 0 ;
}
struct pp
{
double x, y;
}a[200],b[200] ;
bool Can(pp x ,pp y)
{
double dis = sqrt((x.x - y.x) * (x.x - y.x) + (x.y - y.y) * (x.y - y.y)) ;
double t = dis / v ;
if(t <= s)
return 1 ;
return 0;
}
void build()
{
REP(i,1,n)
{
REP(j,1,m)
{
Map[i][j] = Can(a[i],b[j]) ;
}
}
}
int main()
{
while(cin >> n >> m >> s >> v )
{
REP(i,1,n)cin >> a[i].x >> a[i].y ;
REP(i,1,m)cin >> b[i].x >> b[i].y ;
build() ;
mem(link ,-1) ;
int ans = 0 ;
REP(i,1,n)
{
mem(vis,0) ;
ans += dfs(i) ;
}
cout << n - ans < }
return 0;
}

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define PI acos(-1.0)
#define Max 2005
#define inf 1<<28
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define REP(i,s,t) for(int i=(s);i<=(t);++i)
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
#define mp(a,b) make_pair(a,b)
using namespace std;

int Map[200][200] ;
int link[200] ;
bool vis[200] ;
int n , m ,s, v;
int dfs(int now)
{
for (int i = 1 ;i <= m ;i ++)
{
if(Map[now][i])
if(!vis[i])
{
vis[i] = 1 ;
if(link[i] == -1 || dfs(link[i]))
{
link[i] = now ;
return 1 ;
}
}
}
return 0 ;
}
struct pp
{
double x, y;
}a[200],b[200] ;
bool Can(pp x ,pp y)
{
double dis = sqrt((x.x - y.x) * (x.x - y.x) + (x.y - y.y) * (x.y - y.y)) ;
double t = dis / v ;
if(t <= s)
return 1 ;
return 0;
}
void build()
{
REP(i,1,n)
{
REP(j,1,m)
{
Map[i][j] = Can(a[i],b[j]) ;
}
}
}
int main()
{
while(cin >> n >> m >> s >> v )
{
REP(i,1,n)cin >> a[i].x >> a[i].y ;
REP(i,1,m)cin >> b[i].x >> b[i].y ;
build() ;
mem(link ,-1) ;
int ans = 0 ;
REP(i,1,n)
{
mem(vis,0) ;
ans += dfs(i) ;
}
cout << n - ans < }
return 0;
}

这是网络流:


[cpp]
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define PI acos(-1.0)
#define Max 2005
#define inf 1<<28
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define REP(i,s,t) for(int i=(s);i<=(t);++i)
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
#define mp(a,b) make_pair(a,b)
using namespace std;

int n , m ,s, v;
struct kdq
{
int s , e, l ,next ;
}ed[200000] ;
int head[200], num ;
void add(int s ,int e ,int l)
{
ed[num].s = s ;
ed[num].e = e ;
ed[num].l = l ;
ed[num].next = head[s] ;
head[s] = num ++ ;

ed[num].s = e ;
ed[num].e = s ;
ed[num].l = 0 ;