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

inline int rd() {
    int x = 0;
    bool f = 0;
    char c = getchar();
    for (; !isdigit(c); c = getchar()) f |= (c == '-');
    for (; isdigit(c); c = getchar()) x = x * 10 + (c ^ 48);
    return f ? -x : x;
}

#define pb push_back

const double eps = 1e-4;
#define z(x)   (abs((x)) <= eps)
#define sqr(x) ((x) * (x)) 

#define letp const P

struct P {
    double x, y, z;
    P operator + (letp &p) const {return {x + p.x, y + p.y, z + p.z};} 
    P operator - (letp &p) const {return {x - p.x, y - p.y, z - p.z};} 
    bool operator == (letp &p) const {return x == p.x && y == p.y && z == p.z;}
    P rotx(double ang) const { 
        double cosa = cos(ang), sina = sin(ang);
        return {x, y * cosa - z * sina, y * sina + z * cosa};
    }
    P roty(double ang) const { 
        double cosa = cos(ang), sina = sin(ang);
        return {x * cosa - z * sina, y, x * sina + z * cosa};
    } 
    P rotz(double ang) const { 
        double cosa = cos(ang), sina = sin(ang);
        return {x * cosa - y * sina, x * sina + y * cosa, z};
    } 
    inline double dis(letp &p) {
        return sqrt(sqr(x - p.x) + sqr(y - p.y) + sqr(z - p.z));
    }
};

inline bool check(vector<P> &s, int n) {
    for (auto x : s) {
        int cnt = 0;
        for (auto y : s) { 
            if (x == y) continue;
            double d = x.dis(y);
            if (z(d)) return false;
            if (z(d - 1.0)) ++cnt;
        }
        if (cnt != n) return false;
    }
    return true;
}

inline vector<P> random_rotate(vector<P> &s) {
    vector<P> res = s;
    double ang = 360.0 * rand() / RAND_MAX;
    for (auto &x : res) x = x.rotx(ang); 
    ang = 360.0 * rand() / RAND_MAX;
    for (auto &x : res) x = x.roty(ang); 
    ang = 360.0 * rand() / RAND_MAX;
    for (auto &x : res) x = x.rotz(ang); 
    return res;
}

vector<P> ans[11];

int main() {
    ans[1].pb({0, 0, 0}); ans[1].pb({1, 0, 0});
    ans[2] = ans[1]; ans[2].pb({0.5, sqrt(3) / 2, 0});
    ans[3] = ans[2]; ans[3].pb({0.5, sqrt(3) / 6, sqrt(6) / 3});
    ans[4].pb({0, 0, 0}); ans[4].pb({1, 0, 0}); 
    ans[4].pb({0, 1, 0}); ans[4].pb({1, 1, 0});
    ans[4].pb({0.5, 0.5, sqrt(2) / 2}); ans[4].pb({0.5, 0.5, -sqrt(2) / 2});
    for (int i = 5; i <= 10; ++i) {
        int x = 1, nw = 100;
        for (int j = 1; j < i; ++j)
            if (ans[j].size() * ans[i - j].size() < nw) {
                nw = ans[j].size() * ans[i - j].size(); x = j;
            }
        vector<P> p;
        while (true) {
            ans[i].clear();
            p = random_rotate(ans[x]);
            for (auto a : p)
                for (auto b : ans[i - x]) ans[i].pb(a + b);
            if (check(ans[i], i)) break;
        }
    }

    int n = rd();
    printf("%d\n", (int)ans[n].size());
    for (auto x : ans[n]) printf("%.10lf %.10lf %.10lf\n", x.x, x.y, x.z);
    return 0;
}