VIJOS 1212 Way Selection

題目大意：

非常顯然的二分圖最大匹配。

解題思路：

上匈牙利。

關於匈牙利算法。

就是關於尋找增廣路的算法。

僞代碼如下：

bool 尋找從k出發的對應項出的可增廣路
{
    while (從鄰接表中列舉k能關聯到頂點j)
    {
        if (j不在增廣路上)
        {
            把j加入增廣路;
            if (j是未蓋點 或者 從j的對應項出發有可增廣路)
            {
                修改j的對應項爲k;
                則從k的對應項出有可增廣路,返回true;
            }
        }
    }
    則從k的對應項出沒有可增廣路,返回false;
}

void 匈牙利hungary()
{
    for i->1 to n
    {
        if (則從i的對應項出有可增廣路)
            匹配數++;
    }
    輸出 匹配數;
}

本題代碼：

#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int r,a,t;
double px[1001],py[1001];
int lk[1001],y[1001],mp[1001][1001];
int ans=0;
bool pd(double x,double y,double v,int k)
{
	if(v*t>=sqrt((x-px[k])*(x-px[k])+(y-py[k])*(y-py[k])))return 1;
	return 0;
}
bool find(int v)
{
	for(int i=1;i<=a;i++)
	   if(mp[v][i]&&!y[i])
	   {
	   	   y[i]=1;
	   	   if(!lk[i]||find(lk[i]))
	   	     {
			  lk[i]=v;
	   	      return 1;
	   	  }
	   }
	   return 0;
}
int main()
{
    double xx,yy,v;
	cin>>r>>a>>t;
	for(int i=1;i<=a;i++)
	   cin>>px[i]>>py[i];
	for(int i=1;i<=r;i++)
	   {
	   	  cin>>xx>>yy>>v;
	   	  for(int j=1;j<=a;j++)
	   	  {
	   	  	 if(pd(xx,yy,v,j))mp[i][j]=1;
	   	  }
	   }
	for(int i=1;i<=r;i++)
	{
		memset(y,0,sizeof(y));
		if(find(i))ans++;
	}
	cout<<ans;
	return 0;
}