POJ2536 二分图 匈牙利算法

还是二分图。

个人觉得如果把二分图的匈牙利算法学会了。

这些应用都是挺简单的。

#include<stdio.h>

#include<memory.h>
#include<math.h>
int n,m,s,v;
const int N=105;
int link[N];
bool visit[N];
struct
{
double x,y;
}h[N],g[N];
bool mat[N][N];
bool dfs(int u)
{
for(int i=1;i<=m;i++)
{
if(!visit[i]&&mat[u][i])
{
visit[i]=true;
if(link[i]==-1||dfs(link[i]))
{
link[i]=u;
return true;
}
}
}
return false;
}
int main()
{
while(scanf("%d%d%d%d",&n,&m,&s,&v)!=EOF)
{
memset(mat,0,sizeof(mat));
for(int i=1;i<=n;i++)
{
scanf("%lf%lf",&g[i].x,&g[i].y);
}
for(int i=1;i<=m;i++)
{
scanf("%lf%lf",&h[i].x,&h[i].y);
}
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
{
if(((h[j].y-g[i].y)*(h[j].y-g[i].y)+(h[j].x-g[i].x)*(h[j].x-g[i].x))/double(v*v)<double(s*s))
{
mat[i][j]=true;
}
}
memset(link,-1,sizeof(link));
int ans=0;
for(int i=1;i<=n;i++)
{
memset(visit,0,sizeof(visit));
if(dfs(i))
ans++;
}
printf("%d\n",n-ans);
}
return 0;
}