/**
 *  author:  fractal
 *  timus:   288481RF
 *  created: 03/25/21 01:16
**/

#include <bits/stdc++.h>
using namespace std;

#define F first
#define S second
#define mp make_pair
#define pb push_back
#define pf push_front
#define ppb pop_back
#define ppf pop_front
#define speed ios_base::sync_with_stdio(0),cin.tie(0),cout.tie(0)
#define sz(x) (int)x.size()
#define len(x) (int)strlen(x)
#define all(x) x.begin(), x.end()
#define debug cerr << "OK\n";
#define ub upper_bound
#define lb lower_bound
#define nl printf("\n");
#define clbuff fflush(stdin);

mt19937 bruh(chrono::steady_clock::now().time_since_epoch().count());
mt19937_64 rofl(chrono::steady_clock::now().time_since_epoch().count());

typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef vector<int> vi;
typedef vector<ll> vll;
typedef vector<pii> vpii;
typedef vector<pll> vpll;
typedef set<int> si;
typedef set<ll> sll;
typedef set<pii> spii;
typedef set<pll> spll;
typedef multiset <int> msi;
typedef multiset <ll> msll;
typedef map <int, int> mi;
typedef map <ll, ll> mll;

const int N = 1e5 + 2;
const int M = 1e5;
const int mod = 0;
const int inf = 2e9 + 3;
const ll INF = 1e18;
const ld pi2 = 2.0 * 3.14159265359;
const ld pi = 3.14159265359;

const int dx[4] = {1, -1, 0, 0};
const int dy[4] = {0, 0, -1, 1};

void files(string s = "main") {
    #ifndef PC
        freopen((s + ".in").c_str(), "r", stdin);
        freopen((s + ".out").c_str(), "w", stdout);
    #endif
}

int add(int a, int b) {
    if (a + b < 0) return a + b + mod;
    if (a + b >= mod) return a + b - mod;
    return a + b;
}

int mul(int a, int b) {
    return a * 1LL * b % mod;
}

int binpow(int a, int n) {
    int ret = 1;
    while (n) {
        if (n & 1) ret = mul(ret, a);
        a = mul(a, a);
        n >>= 1;
    }
    return ret;
}

struct point {
    ld x, y;

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

    point norm() {
        return point(-y, x);
    }

    ld l() {
        return sqrt(x * x + y * y);
    }

    point unit() {
        return point(x / this->l(), y / this->l());
    }

    point rot(ld alph) {
        alph = alph * pi / 180.0;
        ld c = cos(alph);
        ld s = sin(alph);
        return point(x * c - y * s, x * s + y * c);
    }

    point operator*(const ld c) {
        return point(x * c, y * c);
    }

    point operator/(const ld c) {
        return point(x / c, y / c);
    }

    point operator+(const point a) {
        return point(x + a.x, y + a.y);
    }


    point operator-(const point a) {
        return point(x - a.x, y - a.y);
    }

    ll operator*(const point a) {
        return x * a.y - y * a.x;
    }

    ll operator%(const point a) {
        return x * a.x + y * a.y;
    }

    bool operator<(const point a) {
        if (x == a.x)
            return y < a.y;
        return x < a.x;
    }

    bool operator==(const point a) {
        return x == a.x && y == a.y;
    }

    ld polar() {
        if (x >= 0 && y >= 0)
            return asin(y / this->l());
        if (x >= 0 && y < 0)
            return pi2 - asin(-y / this->l());
        if (x < 0 && y >= 0)
            return pi - asin(y / this->l());
        if (x < 0 && y < 0)
            return pi + asin(-y / this->l());
    }
};

ld dist(point a, point b) {
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

int n;
point p[N];
vector <point> q;

void add(point a) {
    while (sz(q) && q.back() == a)
        q.ppb();
    while (sz(q) > 1) {
        int l = sz(q) - 1;
        if ((a - q[l]) * (q[l] - q[l - 1]) > 0)
            break;
        if ((a - q[l]) * (q[l] - q[l - 1]) == 0 && (a - q[l]) % (q[l] - q[l - 1]) < 0)
            break;
        q.ppb();
    }
    q.pb(a);
}

int main() {
    speed;
    files("tangent");
    point p, q;
    ld r;
    cin >> p.x >> p.y >> r;
    cin >> q.x >> q.y;
    ld d = dist(q, p);
    if (d == r) {
        cout << 1 << '\n';
        cout << fixed << setprecision(10) << q.x << " " << q.y << '\n';
        return 0;
    }
    if (d < r) {
        cout << 0 << '\n';
        return 0;
    }
    ld l = sqrt(d * d - r * r);
    ld c = (l * l - r * r + d * d) / (2.0 * d);
    ld h = sqrt(l * l - c * c);
    point go = (p - q).unit();
    point v = q + go * c;
    point a1 = v + (go.norm() * h);
    point a2 = v - (go.norm() * h);
    cout << 2 << '\n';
    cout << fixed << setprecision(10) << v.x << " " << v.y << '\n';
    cout << fixed << setprecision(10) << c << " " << h << '\n';
    cout << fixed << setprecision(10) << a1.x << " " << a1.y << '\n';
    cout << fixed << setprecision(10) << a2.x << " " << a2.y << '\n';
}