#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) {}

        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;
        }

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

        point<numeric_type> &operator/=(const numeric_type &ri)
        {
            x /= ri;
            y /= ri;
            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;
        }

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

        friend point<numeric_type> operator/(const point<numeric_type> &a, const 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, typename float_type = double>
    float_type abs(const point<numeric_type> &a)
    {
        return sqrtl(norm(a));
    }

    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>> intersect_circle_circle(const circle<numeric_type> &c, const circle<numeric_type> &d)
    {
        float_type dist = abs(c.p - d.p);
        if (sign(dist - c.r - d.r) > 0)
            return {}; 
        if (sign(dist - std::abs(c.r - d.r)) < 0)
            return {}; 
        float_type a = (c.r * c.r - d.r * d.r + dist * dist) / (2 * dist);
        float_type h = sqrt(c.r * c.r - a * a);
        auto v = (d.p - c.p) / dist;
        if (sign(h) == 0)
            return {c.p + v * a};                                
        return {c.p + v * a + point<numeric_type>(0, 1) * v * h, 
                c.p + v * a - point<numeric_type>(0, 1) * v * h};
    }

} 

#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;

    circle<double> c1, c2;
    cin >> c1 >> c2;

    auto res = intersect_circle_circle(c1, c2);

    sort(res.begin(), res.end());
    if (res.size() == 1)
        res.push_back(res[0]);

    cout << res << endl;

    return 0;
}





Measured memory_bytes: 3304