Найти точку на плоскости, сумма расстояний от которой до остальных точек множества максимальна

#include <iostream>
#include <cmath>
using namespace std;
 
struct point
{ float x,y,z;
};
double dlina (point a)
{ return sqrt(a.x*a.x+a.y*a.y+a.z*a.z);}
int main()
{ point a[100];
int n;
cout<<"введите n";
cin>>n;
for (int i=0; i<n; i++)
{cout<<"введите координаты точки"<<i<<" от точки";
cin>>a[i].x>>a[i].y>>a[i].z;}
double max=0; int k=0;
for (int i=0; i<n; i++)
    if (dlina (a[i])>max)
{max=dlina(a[i]);
k=i;}
cout<<"k="<<k;
return 0; }

#include<iostream>
#include<windows.h>
#include<fstream>
#include<cmath>
using namespace std;
 
struct my_point
{
    int x;
    int y;
    int z;
};
 
int GetPoint( my_point *& , const char* );
double Distance( my_point* , my_point , int );
 
int main()
{
    int k;
    my_point* p;
    k = GetPoint(p,"Output.txt");
 
    double MaxSumLong=0, Dl;
    for(int i=0; i<k; i++)
    {
        Dl=Distance(p,p[i],k);
        if(Dl >  MaxSumLong)
             MaxSumLong=Dl;
    }
    cout<<"Max long "<<MaxSumLong<<endl;
    system("pause");
 
}
 
 
int GetPoint(my_point*& pArrayPpoint, const char* FileName)
{
    fstream File(FileName, ios::in);
    int k;
    File>>k;
    pArrayPpoint = new my_point [k];
    for(int i=0; i<k; i++)
    {
        File>>pArrayPpoint[i].x>>pArrayPpoint[i].y>>pArrayPpoint[i].z;
    }
    return k;
}
 
 double Distance( my_point* pArrayPpoint,my_point pl, int k)
{
    double SumDistance=0;
    for(int i=0; i<k; i++)
        SumDistance += sqrt(  pow( (double) pl.x -pArrayPpoint[i].x, 2) + pow( (double) pl.y-pArrayPpoint[i].y, 2) +  pow( (double) pl.z-pArrayPpoint[i].z, 2)  );
    return SumDistance;
}