snoi多校模拟赛 1.18 t2 clique

B最大团(clique.pas/c/cpp)

TL:3S  ML:512MB

【Description】

给定二维坐标上的N个点，如果两个点之间的距离大于K，则他们不能同时被选取。

求最大团的大小。也即，选出最多点，使得这些点两两之间的距离不大于K。

【Input】

多组测试，第一行一个整数T，表示数据的组数。

之后T组数据，每组数据第一行两个非负整数N, K。

之后N行，每行两个整数，表示点的坐标。

【Output】

共T行，每行一个整数，表示最大团的大小。

【Sample Input】

       2

5 20

0 0

0 2

100 100

100 110

100 120

4

1

0 0

0 1

1 0

1 1

【Sample Output】

3

2

【Hint】

20%：N <= 20

40%：N<=50

100%：N<=100，0<=K<=10000，所有输入坐标的绝对值不超过10000。T<=10。

这道题需要的是几何建模而不是图论建模，如果你想的是完全图，那么恭喜GG。

每次如果枚举两个点，以K为半径画圆，两圆的交部分中的所有点与两圆的距离都在K以内，所以我们只需要知道两圆交部分内的点的距离是否满足。

我们连接两圆圆心，分成两部分，两部分内的点都不会距离大于K，但是跨连线就有可能，这时候就需要图论建模了。

二分图，显而易见。

我们枚举找出超距离的点加入二分图内，跑一边最大匹配，用原来的点数+2-最大匹配数即为这次的答案，取max即可

特判：仅有1个点

ans初值赋为-1即可避免特判。

真的要%想出这个方法的julao，太玄了。

代码：

#include <bits/stdc++.h>
using namespace std;
#define ll long long
int read()
{
	char c;int sum=0,f=1;
	c=getchar();  
    while(c<'0' || c>'9'){if(c=='-')f=-1;c=getchar();}  
    while(c>='0' && c<='9'){sum=sum*10+c-'0';c=getchar();}  
    return sum*f;  
}
struct P
{
	int x,y;
}p[305];
P operator-(P a,P b){return (P){a.x-b.x,a.y-b.y};}
int dis(P a){return abs(a.x*a.x+a.y*a.y);}
int delta(P a,P b){return a.x*b.y-a.y*b.x;}
int T,N,K,ans,match[305],e[305][305],vis[305],L,R,l[305],r[305];
bool dfs(int now)
{
	for(int i=1;i<=R;i++)
	    if(e[now][i]&&!vis[i])
	    {
	    	vis[i]=1;
	    	if(!match[i]||dfs(match[i]))
	    	{
	    		match[i]=now;
	    		return 1;
			}
		}
	return 0;
}
int hungry()
{
	int res=0;
	memset(match,0,sizeof(match));
	for(int i=1;i<=R;i++) match[i]=0;
	for(int i=1;i<=L;i++)
	{
		for(int j=1;j<=R;j++) vis[j]=0;
		if(dfs(i)) res++;
	}
	return res;
}
int main()
{
	freopen("clique.in","r",stdin);freopen("clique.out","w",stdout);
	T=read();
	while(T--)
	{
		ans=-1;
		N=read();K=read();
		for(int i=1;i<=N;i++) p[i].x=read(),p[i].y=read();
		for(int i=1;i<=N;i++)
		    for(int j=i+1;j<=N;j++)
		        if(dis(p[i]-p[j])<=K*K)
		        {
		        	L=R=0;
		        	for(int k=1;k<=N;k++)
		        		if(k!=i&&k!=j&&dis(p[i]-p[k])<=dis(p[i]-p[j])&&dis(p[j]-p[k])<=dis(p[i]-p[j]))
		        		{
		        			if(delta(p[j]-p[i],p[k]-p[i])<=0) l[++L]=k;
		        			else r[++R]=k;
						}
					for(int k=1;k<=L;k++)
					    for(int o=1;o<=R;o++)
						    e[k][o]=0;
		            for(int k=1;k<=L;k++)
		                for(int o=1;o<=R;o++)
		                    if(dis(p[l[k]]-p[r[o]])>K*K)
		                        e[k][o]=1;
		            ans=max(ans,L+R-hungry());
				}
		printf("%d\n",ans+2);
	}
}