Консультация # 189359: Уважаемые эксперты! Пожалуйста,помогите с реализацией на Си, желательно с комментариями строк: Задано множество m точек в трёхмерном пространстве. Найти такую из них, что шар заданного радиуса с центром в этой точке содержит максимальное число точек множества. ...
Уважаемые эксперты! Пожалуйста,помогите с реализацией на Си, желательно с комментариями строк:
Задано множество m точек в трёхмерном пространстве. Найти такую из них, что шар заданного радиуса с центром в этой точке содержит максимальное число точек множества.
/*
Задано множество m точек в трёхмерном пространстве.
Найти такую из них, что шар заданного радиуса с центром в этой точке
содержит максимальное число точек множества
*/
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <malloc.h>
#include <math.h>
#include <windows.h>
typedef struct _point
{
double x;
double y;
double z;
}point;
// Вычисляет расстояние между точками
double distance(point a, point b)
{
return sqrt((a.x - b.x)*(a.x - b.x) + (a.y - b.y)*(a.y - b.y) + (a.z - b.z)*(a.z - b.z));
}
int main()
{
SetConsoleCP(1251);
SetConsoleOutputCP(1251);
int i, j, m, r, n, max, imax;
printf("Введите число точек m: ");
scanf("%d", &m);
printf("Введите радиус r: ");
scanf("%d", &r);
point* p = (point*)malloc(m*sizeof(point)); //массив точек
srand(time(NULL)); //инициируем генератор псевдослучайных чисел
for(i=0; i<m; i++) //заполним массив координат случайными дробными числами
{ // в интервале [-100,99]
// чтобы было число дробное, делим на случайное число [1,9]
// +1, чтобы не было деления на 0
p[i].x = (double)((rand()%200) - 100) / ((rand()%9)+1);
p[i].y = (double)((rand()%200) - 100) / ((rand()%9)+1);
p[i].z = (double)((rand()%200) - 100) / ((rand()%9)+1);
}
for(i=0; i<m; i++) //по всем точкам
{
for(n=j=0; j<m; j++) //считаем расстояние от i-й до всех остальных
{ // и считаем те, которые ближе, чем радиус r
if (j!=i) //исключаем из расчета i-ю точку
n += distance(p[i], p[j]) <= r; //если расстояние меньше r, то прибавится 1
}
if (i==0)
{
max = n; //количество для первой точки сохраняем, как максимальное
imax = i; //сохраним индекс центра
}
else if (max < n) //для всех остальных сравниваем с максимумом
{
max = n; //новый максимум
imax = i; //индекс нового центра
}
}
if (max == 0) //если окажется, что не найдено ни одной точки, ближе, чем на r
printf("Точки не найдены\n");
else //выводим результат
printf("Центр шара: x = %.2f, y = %.2f, z = %.2f\nЧисло точек внутри = %d\n",
p[imax].x, p[imax].y, p[imax].z, max);
free(p); //освободим память под массив точек
return 0;
}
Консультировал: Лысков Игорь Витальевич (Старший модератор)
Дата отправки: 14.05.2016, 12:24
5
нет комментария ----- Дата оценки: 15.05.2016, 20:29
Команда портала RFPRO.RU благодарит Вас за то, что Вы пользуетесь нашими услугами. Вы только что прочли очередной выпуск рассылки. Мы старались.
Пожалуйста, оцените его. Если совет помог Вам, если Вам понравился ответ, Вы можете поблагодарить автора -
для этого в каждом ответе есть специальные ссылки. Вы можете оставить отзыв о работе портале. Нам очень важно знать Ваше мнение.
Вы можете поближе познакомиться с жизнью портала, посетив наш форум, почитав журнал,
который издают наши эксперты. Если у Вас есть желание помочь людям, поделиться своими знаниями, Вы можете зарегистрироваться экспертом.
Заходите - у нас интересно!