HDU 5046 Airport(dlx)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5046

题意：n个城市修建m个机场，使得每个城市到最近进场的最大值最小。

思路：二分+dlx搜索判定。

 

const int INF=1000000005;

const int N=4444;





int m;



struct node

{

	int L[N],R[N],D[N],U[N],e;

	int col[N];

	int H[N],num[N];





	int visit[N],KK;



	void init(int n)

	{

		KK=0;

		int i;

		for(i=0;i<=n;i++)

		{

			if(i==n) L[i]=0;

			else L[i]=i+1;

			if(i==0) R[i]=n;

			else R[i]=i-1;



			num[i]=0;



			H[i]=-1;

			D[i]=U[i]=i;

		}

		e=n+1;

	}



	void add(int r,int c)

	{

		U[e]=c;

		D[e]=D[c];

		U[D[c]]=e;

		D[c]=e;

		if(H[r]<0) H[r]=L[e]=R[e]=e;

		else

		{

			L[e]=L[H[r]];

			R[e]=H[r];

			R[L[H[r]]]=e;

			L[H[r]]=e;

		}

		num[c]++;

		col[e]=c;

		e++;

	}



	void remove(int c)

	{

		int i;

		for(i=D[c];i!=c;i=D[i])

		{

			R[L[i]]=R[i];

			L[R[i]]=L[i];

		}

	}



	void resume(int c)

	{

		int i;

		for(i=U[c];i!=c;i=U[i])

		{

			L[R[i]]=R[L[i]]=i;

		}

	}



	int astar()

	{

		KK++;

		int ans=0,i,j,k;

		for(i=L[0];i!=0;i=L[i]) if(KK!=visit[i])

		{

			visit[i]=KK;

			ans++;

			for(j=D[i];j!=i;j=D[j]) for(k=L[j];k!=j;k=L[k])

			{

				visit[col[k]]=KK;

			}

		}

		return ans;

	}



	int DFS(int cnt)

	{

		if(L[0]==0) return 1;

		if(cnt+astar()>m) return 0;



		int tmp=INF,c,i;

		for(i=L[0];i!=0;i=L[i]) if(num[col[i]]<tmp)

		{

			tmp=num[col[i]];

			c=i;

		}

		for(i=D[c];i!=c;i=D[i])

		{

			remove(i);

			int j;

			for(j=L[i];j!=i;j=L[j]) remove(j);

			if(DFS(cnt+1)) return 1;



			for(j=R[i];j!=i;j=R[j]) resume(j);

			resume(i);

		}

		return 0;

	}

}A;



int a[66][2];

i64 dis[66][66];

int n;

i64 d[4444];





inline i64 dist(int i,int j)

{

	return ((i64)abs(a[i][0]-a[j][0]))+abs(a[i][1]-a[j][1]);

}



int ok(i64 M)

{

	A.init(n);

	int i,j;

	for(i=1;i<=n;i++) for(j=1;j<=n;j++) if(dis[i][j]<=M)

	{

		A.add(i,j);

	}

	return A.DFS(0);

}



i64 cal()

{

	int num=0;

	int i,j;

	for(i=1;i<=n;i++) for(j=1;j<=n;j++) dis[i][j]=dist(i,j),d[num++]=dis[i][j];

	sort(d,d+num);

	int M=unique(d,d+num)-d;

	int low=0,high=M-1;

	i64 ans=0;



	while(low<=high)

	{

		int mid=(low+high)>>1;

		if(ok(d[mid])) ans=d[mid],high=mid-1;

		else low=mid+1;

	}

	return ans;

}



int main()

{

	int num=0;

	int T=getInt();

	while(T--)

	{

		n=getInt();

		m=getInt();

		int i;

		for(i=1;i<=n;i++) scanf("%d%d",&a[i][0],&a[i][1]);

		printf("Case #%d: ",++num);

		cout<<cal()<<endl;

	}

}