Ускорить программу сортировки списка в порядке возрастания

#include <iostream>
#include <cmath>
using namespace std;
 
int main ()
{
    long long n;
    cin >> n;
    long long arr[n][3];
    float arr1[n];
    for(int i=0; i<n; i++)
    {
        cin>>arr[i][0]>>arr[i][1]>>arr[i][2];
        arr1[i]=sqrt(arr[i][0]*arr[i][0]+arr[i][1]*arr[i][1]+arr[i][2]*arr[i][2]);
    }
    for(int i=0; i<n; i++)
    {
        for(int xi=0; xi<n; xi++)
        {
            if(arr1[i]<arr1[xi])
            {
                swap(arr1[i],arr1[xi]);
                swap(arr[i][3],arr[xi][3]);
                swap(arr[i][2],arr[xi][2]);
                swap(arr[i][1],arr[xi][1]);
                swap(arr[i][0],arr[xi][0]);
            }
            else if(arr1[i]==arr1[xi])
            {
                if(arr[i][3]<arr[xi][3])
                {
                    swap(arr1[i],arr1[xi]);
                    swap(arr[i][3],arr[xi][3]);
                    swap(arr[i][2],arr[xi][2]);
                    swap(arr[i][1],arr[xi][1]);
                    swap(arr[i][0],arr[xi][0]);
                }
                else  if(arr[i][3]==arr[xi][3])
                {
                    if(arr[i][2]<arr[xi][2])
                    {
                        swap(arr1[i],arr1[xi]);
                        swap(arr[i][3],arr[xi][3]);
                        swap(arr[i][2],arr[xi][2]);
                        swap(arr[i][1],arr[xi][1]);
                        swap(arr[i][0],arr[xi][0]);
                    }
                    else  if(arr[i][2]==arr[xi][2])
                    {
                        if(arr[i][1]<arr[xi][1])
                        {
                            swap(arr1[i],arr1[xi]);
                            swap(arr[i][3],arr[xi][3]);
                            swap(arr[i][2],arr[xi][2]);
                            swap(arr[i][1],arr[xi][1]);
                            swap(arr[i][0],arr[xi][0]);
                        }
                        else if(arr[i][1]==arr[xi][1])
                        {
                            if(arr[i][0]<arr[xi][0])
                            {
                                swap(arr1[i],arr1[xi]);
                                swap(arr[i][3],arr[xi][3]);
                                swap(arr[i][2],arr[xi][2]);
                                swap(arr[i][1],arr[xi][1]);
                                swap(arr[i][0],arr[xi][0]);
                            }
                        }
                    }
                }
            }
        }
    }
    for(int i=0; i<n; i++)
    {
        cout<<arr[i][0]<<" "<<arr[i][1]<<" "<<arr[i][2]<<endl;
 
    }
    return 0;
}

#include <cmath>
#include <iostream>
#include <map>
#include <tuple>
 
double Foo(const int x, const int y, const int z)
{
    return sqrt(x*x + y*y + z*z);
}
 
int main()
{
    size_t N;
    std::cin >> N;
    int x, y, z;
    std::multimap<double, std::tuple<int, int, int>> m;
    for (size_t i = 0; i < N && std::cin >> x >> y >> z; i++)
        m.emplace(Foo(x, y, z), std::make_tuple(x, y, z));
    for (const auto &p : m)
        std::cout << std::get<0>(p.second) << " " << std::get<1>(p.second) << " " << std::get<2>(p.second) << std::endl;
}

#include <algorithm>
#include <cmath>
#include <iostream>
#include <vector>
 
struct A
{
    A(const int x_, const int y_, const int z_) : x(x_), y(y_), z(z_), p(sqrt(x*x + y*y + z*z)) {}
    int x, y, z;
    double p;
};
 
int main()
{
    size_t N;
    std::cin >> N;
    int x, y, z;
    std::vector<A> v;
    v.reserve(N);
    for (size_t i = 0; i < N && std::cin >> x >> y >> z; i++)
        v.emplace_back(x, y, z);
    std::sort(v.begin(), v.end(), [](const A &lhs, const A &rhs) { return lhs.p < rhs.p; });
    for (const auto &a : v)
        std::cout << a.x << " " << a.y << " " << a.z << std::endl;
}

#include <algorithm>
#include <cmath>
#include <iostream>
#include <vector>
 
struct A
{
    A(const int x_, const int y_, const int z_) : x(x_), y(y_), z(z_), p(sqrt(x*x + y*y + z*z)) {}
    int x, y, z;
    double p;
};
 
int main()
{
    size_t N;
    std::cin >> N;
    int x, y, z;
    std::vector<A> v;
    v.reserve(N);
    for (size_t i = 0; i < N && std::cin >> x >> y >> z; i++)
        v.emplace_back(x, y, z);
    std::sort(v.begin(), v.end(), [](const A &lhs, const A &rhs) { return lhs.p < rhs.p; });
    for (const auto &a : v)
        std::cout << a.x << " " << a.y << " " << a.z << std::endl;
}

#include <algorithm>
#include <cmath>
#include <iostream>
#include <vector>
 
struct A
{
    A(const int x_, const int y_, const int z_) : x(x_), y(y_), z(z_), p(sqrt(x*x + y*y + z*z)) {}
    int x, y, z;
    double p;
};
 
bool operator<(const A &lhs, const A &rhs)
{
    if (fabs(lhs.p - rhs.p) > 1e-08)
        return lhs.p < rhs.p;
    if (lhs.x != rhs.x)
        return lhs.x < rhs.x;
    if (lhs.y != rhs.y)
        return lhs.y < rhs.y;
    return lhs.z < rhs.z;
}
 
int main()
{
    size_t N;
    std::cin >> N;
    int x, y, z;
    std::vector<A> v;
    v.reserve(N);
    for (size_t i = 0; i < N && std::cin >> x >> y >> z; i++)
        v.emplace_back(x, y, z);
    std::sort(v.begin(), v.end());
    for (const auto &a : v)
        std::cout << a.x << " " << a.y << " " << a.z << std::endl;
}