如何使用C++ ranges实现numpy.ndindex？

#include <iostream>
#include <vector>
#include <range/v3/all.hpp>

namespace views = ranges::views;

template <typename T>
auto ndindex(const std::vector<T>& shape)
{
    // 生成N维数组的索引序列
    auto ranges = shape | views::transform([](T size) {
        return views::iota(static_cast<T>(0), size);
    });
    
    // 使用ranges::cartesian_product生成笛卡尔积
    return views::cartesian_product(ranges);
}

int main()
{
    std::vector<int> shape = {3, 4, 2}; // 3维数组，每个维度大小分别为3、4和2
    
    // 生成索引序列
    auto indices = ndindex(shape);
    
    // 输出索引序列
    for (const auto& index : indices)
    {
        // index是一个std::tuple，表示一个索引
        std::cout << "(";
        std::apply([](auto... args) {((std::cout << args << ", "), ...);}, index);
        std::cout << ")" << std::endl;
    }
    
    return 0;
}