#include <iostream>

#include <cassert>
#include <array>
#include <initializer_list>
#include <algorithm>

template<typename T, std::size_t d>
class basic_vec : public std::array<T, d> {
public:
    static constexpr std::size_t dimensions = d;

    basic_vec();
    basic_vec(std::initializer_list<T> l);
    template<std::size_t n>
    basic_vec(const basic_vec<T, n>& other);

    T length() const;
};

template<typename T, std::size_t d>
basic_vec<T, d>::basic_vec() : std::array<T, d>() {
    this->fill(static_cast<T>(0));
}

template<typename T, std::size_t d>
basic_vec<T, d>::basic_vec(std::initializer_list<T> l) {
    assert(l.size() <= d);
    this->fill(static_cast<T>(0));
    std::copy(l.begin(), l.end(), this->begin());
}

template<typename T, std::size_t d>
template<std::size_t n>
basic_vec<T, d>::basic_vec(const basic_vec<T, n>& other) {
    static_assert(d >= n);
    this->fill(static_cast<T>(0));
    for (std::size_t i = 0;  i < n; i++) {
        this->at(i) = other.at(i);
    }
}

template<typename T, std::size_t d>
T basic_vec<T, d>::length() const {
    T res = static_cast<T>(0);
    for (const T& val : *this) {
        res += val * val;
    }

    return res;
}

template<typename T, std::size_t d>
basic_vec<T, d> operator+(const basic_vec<T, d>& rhs, const basic_vec<T, d>& lhs) {
    basic_vec<T, d> out;
    auto rit = rhs.cbegin(), lit = lhs.cbegin();
    std::generate(out.begin(), out.end(), [rit, lit] () mutable {
        return *(rit++) + *(lit++);
    });

    return out;
}

template<typename T, std::size_t d>
basic_vec<T, d> operator*(const T& rhs, const basic_vec<T, d>& lhs) {
    basic_vec<T, d> out;
    auto lit = lhs.cbegin();

    std::generate(out.begin(), out.end(), [rhs, lit] () mutable {
        return rhs * *(lit++);
    });

    return out;
}

template<typename T, std::size_t d>
basic_vec<T, d> operator-(const basic_vec<T, d>& rhs, const basic_vec<T, d>& lhs) {
    return rhs + static_cast<T>(-1) * lhs;
}

template<typename T, std::size_t d>
T operator*(const basic_vec<T, d>& rhs, const basic_vec<T, d>& lhs) {
    T res = static_cast<T>(0);
    for (std::size_t i = 0; i < d; i++) {
        res += rhs[i] * lhs[i];
    }

    return res;
}

template<typename T, std::size_t d>
std::ostream& operator<<(std::ostream& os, const basic_vec<T, d>& v) {
    os << "<";
    for (std::size_t i = 0; i < d -1; i++) {
        os << v[i] << ", ";
    }
    os << v[d-1] << ">";
    return os;
}


template<typename T, std::size_t d>
class vec: public basic_vec<T, d> {
public:
    vec(std::initializer_list<T> l) : basic_vec<T, d>(l) {}

    template<std::size_t n>
    vec(const basic_vec<T, n>& other) : basic_vec<T, d>(other) {}
};

template<typename T>
class vec3 : public basic_vec<T, 3> {
public:
    vec3(std::initializer_list<T> l) : basic_vec<T, 3>(l) {}

    template<std::size_t n>
    vec3(const basic_vec<T, n>& other) : basic_vec<T, 3>(other) {}

    static vec3<T> cross(const vec3<T>& rhs, const vec3<T>& lhs);
};

template<typename T>
vec3<T> vec3<T>::cross(const vec3<T>& rhs, const vec3<T>& lhs) {
    vec3<T> res;
    // TODO: cross product

    return res;
}


template<typename T>
class vec2: public basic_vec<T, 2> {
public:
    vec2(std::initializer_list<T> l) : basic_vec<T, 2>(l) {}

    template<std::size_t n>
    vec2(const basic_vec<T, n>& other) : basic_vec<T, 2>(other) {}

    T angle();
    static vec3<T> cross(const vec2<T>& rhs, const vec2<T>& lhs);
};

template<typename T>
vec3<T> vec2<T>::cross(const vec2<T>& rhs, const vec2<T>& lhs) {
    vec3<T> res;
    // TODO: cross product
}


int main(int argc, char *argv[]) {
    vec3<double> v{1, 2, 3};
    vec3<double> u{3, 4, 5};

    std::cout << v << std::endl;
    std::cout << v.length() << std::endl;
    std::cout << v + u << std::endl;
    std::cout << v - u << std::endl;
    std::cout << 2.0 * u << std::endl;
    std::cout << v * u << std::endl;

    return 0;
}