poj:4091:The Closest M Points

poj:4091:The Closest M Points


题目

描写叙述

每到饭点,就又到了一日几度的小L纠结去哪吃饭的时候了。由于有太多太多好吃的地方能够去吃,而小L又比較懒不想走太远,所以小L会先找到距离他近期的M家餐馆然后再做筛选。

小L如今所在的位置和每家餐馆的位置用同一笛卡尔坐标系中的点表示,而点与点之间的距离为欧几里得距离,对于点p = (p1, p2,..., pn)和点q = (q1,q2,..., qn),两者的距离定义例如以下

现给出在K维空间中小L所处的位置的坐标以及n个餐馆的位置,请帮助小L完毕他的需求。


输入
第1行包括两个整数n和K。1≤n≤5000,1≤K≤5。


接下来n行,每行包括K个数,表示每一个餐馆的坐标。
接下来1行。包括一个数t,1≤t≤10000。表示小L询问的数目。
每次询问包括两行。

第1行包括K个数,表示小L所在的坐标。第2行包括一个数M,1≤M≤10。


全部坐标值不会超过10000。
输入数据包括多组数据。请逐个处理直到文件结束。

输出
对于每个询问,输出m+1行:
第1行输出:”the closest M points are:”,当中M在输入中给出。
接下来M行输出距离近期的M家餐馆的坐标,依照由近及远的顺序输出。


输出数据保证答案唯一。保证从小L的位置到近期的M+1家餐馆位置各不同样。这说明例如以下输入数据:
2 2
1 1
3 3
1
2 2
1
不会存在。

例子输入
3 2
1 1
1 3
3 4
2
2 3
2
2 3
1
例子输出
the closest 2 points are:
1 3
3 4
the closest 1 points are:
1 3
解题方案,这个题目比較简单
我们转换问题这是求出k最大值或者最小值问题
暴力法是先排序再求出前m个最小值
快捷法是建立m规模的小顶堆。然后去全部数据去更新一边,最后全部的堆的元素就是结果

代码(暴力法代码)
#include <iostream>
#include <fstream>
#include <list>
#include <math.h>
#include <algorithm>
using namespace std;


typedef pair<double,int> distance_restaurant;



void read_data();
// data 为餐馆地址表
list<distance_restaurant> main_solution(int ** data,int n,int k,int *dizhi);


int main()
{
	read_data();
	system("pause");
	return 0;
}

void read_data()
{
	ifstream reader;
	reader.open("data.txt");
	int n,k;
	reader>>n;
	reader>>k;
	int ** data = new int *[n];
	for( int i=0; i<n; i++ )
	{
		data[i] = new int[k];
		for(int j=0;j<k;j++)
			reader>>data[i][j];
	}

	int t;
	reader>>t;
	while( t>0 )
	{
		t-- ;

		int * dizhi = new int[k] ;
		for( int i=0;i<k;i++ )
			reader>>dizhi[i] ;
		int m;
		reader>>m;
		cout<<"the closest "<<m<<" points are:"<<endl;
		list<distance_restaurant> result = main_solution(data,n,k,dizhi);
		for( list<distance_restaurant>::iterator it = result.begin();m>0 && it != result.end() ;m-- ,it++)
		{
			for(int j=0;j<k;j++)
			{
				cout<<data[it->second][j]<<" ";
			}
			cout<<endl;
		}

	}
}


// data 为餐馆地址表
list<distance_restaurant> main_solution(int ** data,int n,int k,int *dizhi)
{
	list<distance_restaurant> result;
	distance_restaurant elem;
	for( int i=0;i<n;i++)
	{
		elem.second = i;
		elem.first = 0;
		for(int j=0;j<k;j++)
		{
			elem.first += (data[i][j]-dizhi[j]) * (data[i][j]-dizhi[j]);
		}
		elem.first = sqrt(elem.first);
		result.push_back(elem);
	}

	stable_sort(result.begin(),result.end());
	return result;
}



原文地址:https://www.cnblogs.com/zsychanpin/p/7027205.html