Uva Oj 1388 - Graveyard

最新推荐文章于 2024-12-10 05:00:00 发布

原创最新推荐文章于 2024-12-10 05:00:00 发布 · 649 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#1Y #UvaOJ #ACM #大一 #计算几何

本文介绍了一种通过构建特殊结构体来实现点匹配优化的方法。该算法首先定义了原始点和新点，并通过排序确定新点与原始点之间的最优距离。接着记录下这些距离，并在前n个距离中寻找最小总成本，最终输出优化后的结果。

依旧是训练指南上的例题,不过这个事我自己写的

和lrj中的没有关系

建了个结构体来标记这个点是原来的点还是新的点

然后排序

接下来找到序列中的新出现的点

如果这个点左边(右边)是原本点,则记录距离,并更新到z中

如果两边都不是原本点,那么原本点一定不会舍近求远,无视掉就好了

接下来去z的前n项求和,表示n个原本点的最优选择分别是某n个新点

然后转化输出就好了

#include<stdio.h>
#include<vector>
#include<stdlib.h>
struct tmd{int x,y;};
tmd q[3010],z[3010];
int cmp(const void *a,const void *b){
	tmd *p=(tmd *)a;
	tmd *q=(tmd *)b;
	tmd x=*p,y=*q;
	return x.x-y.x;
	}
int min(int a,int b){
	if(a>b)return b;
	return a;
	}
int main(){
	int n,m;
	while(scanf("%d%d",&n,&m)!=EOF){
		int i;
		for(i=0;i<n;i++)q[i]=(tmd){(n+m)*i,1};
		for(i=0;i<m+n;i++)q[i+n]=(tmd){n*i,-1};
		qsort(q,m+n+n,sizeof(tmd),cmp);
		int rs=0;
		int pz=0;
		for(i=0;i<m+n+n;i++){
			if(q[i].y==-1){
				int temp=n*(m+n);
				if(i>0&&q[i-1].y==1)temp=min(temp,abs(q[i-1].x-q[i].x));
				if(i+1<m+n+n&&q[i+1].y==1)temp=min(temp,abs(q[i+1].x-q[i].x));
				if(temp!=n*(m+n))z[pz++]=(tmd){temp,0};
			//	printf("z[%d].x=%d,y=%d\n",pz-1,z[pz-1].x,z[pz-1].y);
				}
		//	printf("q[%d].x=%d,y=%d\n",i,q[i].x,q[i].y);
			}
		qsort(z,pz,sizeof(tmd),cmp);
	//	for(i=0;i<pz;i++)printf("z[%d].x=%d,y=%d\n",i,z[i].x,z[i].y);
		for(i=0;i<min(n,pz);i++)rs+=z[i].x;
		double res=(double)rs/(n*(n+m))*10000;
		printf("%.4lf\n",res);
		}
	return 0;
	}