#include <bits/stdc++.h>

using i64 = long long;

using Point = std::array<double, 3>;
std::vector<Point> a[11] {
    {},
    {{0.000000000, 0.000000000, 0.000000000},
    {1.000000000, 0.000000000, 0.000000000}},
    {{0.000000000, 0.000000000, 0.000000000},
    {1.000000000, 0.000000000, 0.000000000},
    {0.500000000, 0.866025404, 0.000000000}},
    {{0.500000000, 0.000000000, -0.353553391},
    {-0.500000000, 0.000000000, -0.353553391},
    {0.000000000, 0.500000000, 0.353553391},
    {0.000000000, -0.500000000, 0.353553391}},
    {{0.707106781, 0.000000000, 0.000000000},
    {-0.707106781, 0.000000000, 0.000000000},
    {0.000000000, 0.707106781, 0.000000000},
    {0.000000000, -0.707106781, 0.000000000},
    {0.000000000, 0.000000000, 0.707106781},
    {0.000000000, 0.000000000, -0.707106781}},
    {{0.000000000, 0.500000000, 0.809016994},
    {0.000000000, -0.500000000, 0.809016994},
    {0.000000000, 0.500000000, -0.809016994},
    {0.000000000, -0.500000000, -0.809016994},
    {0.500000000, 0.809016994, 0.000000000},
    {0.500000000, -0.809016994, 0.000000000},
    {-0.500000000, 0.809016994, 0.000000000},
    {-0.500000000, -0.809016994, 0.000000000},
    {0.809016994, 0.000000000, 0.500000000},
    {0.809016994, 0.000000000, -0.500000000},
    {-0.809016994, 0.000000000, 0.500000000},
    {-0.809016994, 0.000000000, -0.500000000}}
};

std::mt19937 rng;
std::vector<Point> conv(const std::vector<Point> &a, std::vector<Point> b) {
    const double Pi = std::acos(-1.0);
    double x = Pi * rng() / rng.max();
    double y = Pi * rng() / rng.max();
    double z = Pi * rng() / rng.max();
    std::vector<Point> c;
    for (auto &p : b) {
        double u = p[1], v = p[2];
        p[1] = u * std::cos(x) - v * std::sin(x);
        p[2] = u * std::sin(x) + v * std::cos(x);
        u = p[2], v = p[0];
        p[2] = u * std::cos(y) - v * std::sin(y);
        p[0] = u * std::sin(y) + v * std::cos(y);
        u = p[0], v = p[1];
        p[0] = u * std::cos(z) - v * std::sin(z);
        p[1] = u * std::sin(z) + v * std::cos(z);
    }
    for (auto p : a) {
        for (auto q : b) {
            c.push_back({p[0] + q[0], p[1] + q[1], p[2] + q[2]});
        }
    }
    return c;
}

int main() {
    std::ios::sync_with_stdio(false);
    std::cin.tie(nullptr);
    
    a[6] = conv(a[3], a[3]);
    a[7] = conv(a[3], a[4]);
    a[8] = conv(a[4], a[4]);
    a[9] = conv(a[4], a[5]);
    a[10] = conv(a[4], a[6]);
    
    int n;
    std::cin >> n;
    
    std::cout << a[n].size() << "\n";
    std::cout << std::fixed << std::setprecision(10);
    for (auto [x, y, z] : a[n]) {
        std::cout << x << " " << y << " " << z << "\n";
    }
    
    return 0;
}