fzu 2194 星系碰撞 二分图最大匹配

fzu 2194

福州大学第十二届程序竞赛的题，曾经有一个ak的机会在我面前，我却没有好好珍惜，如果上天给我一次重来的机会，我就要写二分图，这个题一看n最大为5e4，立马就被吓跑了，赛完看这题时间给30s.....

这题说同一个星系两个星球距离必定大于5，如果两个星球距离<=5，那么肯定就是不同的星系，把所有距离小于等于5的星球用无向图连接，首先用二分图染色判断是否是二分图，不是二分图说明这些星球不能被两个星系所划分，如果是二分图，求最大匹配即所有不能在一个星系的星球有多少个，然后总星球减去最大匹配的一半就是答案。

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
const int maxn=5e4+10;
vector<int>G[maxn];
int n,left[maxn],vis[maxn];
struct node
{
	int x,y;
}a[maxn];
bool cmp(node t1,node t2)
{
	if(t1.x==t2.x)
	return t1.y<t2.y;
	return t1.x<t2.x;
}
void init()
{
	memset(left,0,sizeof(left));
	memset(vis,0,sizeof(vis));
	for(int i=0;i<maxn;i++)
	G[i].clear();
}
int check(int i,int j)
{
	ll tmp1=ll(1)*(a[i].x-a[j].x)*(a[i].x-a[j].x);
	ll tmp2=ll(1)*(a[i].y-a[j].y)*(a[i].y-a[j].y); 
	return tmp1+tmp2<=25;
}
int dfs1(int u,int color)
{
	for(int i=0;i<G[u].size();i++)
	{
		int v=G[u][i];
		if(!vis[v])
		{
			vis[v]=3-color;
			if(!dfs1(v,3-color))
			return 0;
		}
		else if(vis[v]==color)
		return 0;
	}
	return 1;
}
int dfs(int u)
{
	for(int i=0;i<G[u].size();i++)
	{
		int v=G[u][i];
		if(!vis[v])
		{
			vis[v]=1;
			if(!left[v]||dfs(left[v]))
			{
				left[v]=u;
				return 1;
			}
		}
	}
	return 0;
}
int main()
{
	while(~scanf("%d",&n))
	{
		init();
		for(int i=1;i<=n;i++)
		scanf("%d%d",&a[i].x,&a[i].y);
		sort(a+1,a+1+n,cmp);
		for(int i=2;i<=n;i++)
		{
			for(int j=i-1;j;j--)
			if(a[i].x-a[j].x>5)break;
			else
			{
				if(check(i,j))
				G[i].push_back(j),G[j].push_back(i);
			}
		}
		vis[1]=1;
		if(!dfs1(1,1))//判断二分图 
		puts("-1");
		else
		{
			int ans=0;
			for(int i=1;i<=n;i++)
			{
				memset(vis,0,sizeof(vis));
				ans+=dfs(i);//二分图最大匹配 
			}
			printf("%d\n",n-ans/2);
		}
	}
}