#include <algorithm>
#include <array>
#include <assert.h>
#include <bitset>
#include <chrono>
#include <cmath>
#include <complex>
#include <cstring>
#include <functional>
#include <fstream>
#include <iomanip>
#include <iostream>
#include <istream>

#include <map>
#include <math.h>
#include <numeric>
#include <ostream>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <vector>

namespace asl
{
    
    template <typename numeric_type>
    class point
    {
    public:
        numeric_type x, y;

        point(numeric_type x = 0, numeric_type y = 0) : x(x), y(y) {}

        point<numeric_type> conj() const
        {
            return point(x, -y);
        }

        bool operator<(const point<numeric_type> &ri) const
        {
            if (x != ri.x)
                return x < ri.x;
            return y < ri.y;
        }

        point<numeric_type> &operator+=(const point<numeric_type> &ri)
        {
            x += ri.x;
            y += ri.y;
            return *this;
        }

        point<numeric_type> &operator-=(const point<numeric_type> &ri)
        {
            x -= ri.x;
            y -= ri.y;
            return *this;
        }

        point<numeric_type> &operator*=(const point<numeric_type> &ri)
        {
            numeric_type _x = x * ri.x - y * ri.y;
            numeric_type _y = x * ri.y + y * ri.x;
            x = _x;
            y = _y;
            return *this;
        }

        friend point<numeric_type> operator+(const point<numeric_type> &a, const point<numeric_type> &b)
        {
            return point<numeric_type>(a) += b;
        }

        friend point<numeric_type> operator-(const point<numeric_type> &a, const point<numeric_type> &b)
        {
            return point<numeric_type>(a) -= b;
        }

        friend point<numeric_type> operator*(const point<numeric_type> &a, const point<numeric_type> &b)
        {
            return point<numeric_type>(a) *= b;
        }

    };

    template <typename numeric_type>
    std::ostream &operator<<(std::ostream &out, const point<numeric_type> &p)
    {
        return out << p.x << " " << p.y;
    }

    template <typename numeric_type>
    std::istream &operator>>(std::istream &is, point<numeric_type> &p)
    {
        return is >> p.x >> p.y;
    }

    template <typename numeric_type>
    numeric_type norm(const point<numeric_type> &a)
    {
        return a.x * a.x + a.y * a.y;
    }

    template <typename numeric_type>
    point<numeric_type> conj(const point<numeric_type> &a)
    {
        return a.conj();
    }

    template <const int eps_exp = 9>
    int sign(double v)
    {
        const double eps = std::pow(10, -eps_exp);

        if (v < -eps)
            return -1;
        else if (v > eps)
            return +1;
        else
            return 0;
    }

#include <stdint.h>

    template <typename T>
    std::ostream &operator<<(std::ostream &out, const std::vector<T> &vec)
    {
        bool start = true;

        for (auto x : vec)
        {
            if (start)
            {
                start = false;
            }
            else
            {
                out << " ";
            }

            out << x;
        }
        return out;
    }

} 

#include <experimental/optional>

namespace asl
{
    template <typename numeric_type>
    class circle
    {
    public:
        point<numeric_type> p;
        numeric_type r;
    };

    template <typename numeric_type>
    std::istream &operator>>(std::istream &is, circle<numeric_type> &c)
    {
        return is >> c.p >> c.r;
    }

    template <typename numeric_type, typename float_type = double>
    std::vector<point<numeric_type>> tangent_circle_point(circle<numeric_type> c, point<numeric_type> p)
    {
        float_type sin2 = c.r * c.r / norm(p - c.p);
        if (sign(1 - sin2) < 0)
            return {};
        if (sign(1 - sin2) == 0)
            return {p};
        point<numeric_type> z(1 - sin2, sqrt(sin2) * sqrt(1 - sin2));
        return {p + (c.p - p) * conj(z), p + (c.p - p) * z};
    }

} 

#include <random>

#define endl '\n'

using namespace asl;
using namespace std;

int main()
{
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.precision(12);
    cout << fixed;

    point<double> p;
    circle<double> c;

    cin >> p >> c;

    auto res = tangent_circle_point(c, p);
    sort(res.begin(), res.end());

    cout << res << endl;

    return 0;
}





Measured memory_bytes: 3304