Рекурсия: нахождение биномиальных коэффициентов

#include <vector>
#include <iostream>
 
using namespace std;
 
template <class T>
vector<vector<T> > perm(vector<T> source, int k)
{
    vector<vector<T> > result;
    if (k == 1) {
        result.reserve(source.size());
        for (int i = 0, len = source.size(); i < len; ++i) {
            vector<T> temp;
            temp.push_back(source[i]);
            result.push_back(temp);
        }
    }
    else if (k == source.size()) {
        result.push_back(source);
    }
    else {
        for (int i = 0, len = source.size() - k + 1; i < len; ++i) {
            T last = source.back();
            source.pop_back();
            vector<vector<T> > temp = perm(source, k - 1);
            for (int j = 0, len = temp.size(); j < len; ++j) {
                temp[j].push_back(last);
                result.push_back(temp[j]);
            }
        }
    }
    return result;
}
 
template <class T>
ostream& operator<<(ostream &stream, vector<T> vec)
{
    stream << "[";
    bool notFirst = false;
    for (typename vector<T>::iterator i = vec.begin();
         i != vec.end();
         ++i)
    {
        if (notFirst) {
            stream << ", ";
        }
        stream << *i;
        if (!notFirst) {
            notFirst = true;
        }
    }
    stream << "]";
    return stream;
}
 
int main()
{
    vector<int> array;
    array.reserve(5);
    for (int i = 1; i <= 5; ++i) {
        array.push_back(i);
    }
    vector<vector<int> > perms = perm(array, 2);
    cout << perms;
    return 0;
}

#include <cmath>
#include <iostream>
using namespace std;
 
//Ôóíêöèÿ ïîëó÷àåò è ñ÷èòàåò êîë-âî ðàçðÿäîâ ÷èñëà
//åñëè ïîä arr âûäåëåíà ïàìÿòü òî â íåãî çàíîñÿòüñÿ ðàçðÿäû
int * splitDigit(int val, int *arr, int &n);
void rekursBinom(int pos, int *arr, int n, int k);
 
int main()
{
    int val = 0;
    int k   = 0;
    int n   = 0;//Áóäåò ñîäåðæàòü ÷èñëî ðàçðÿäîâ ÷èñëà
    int *arr= 0;//Óêàçàòåëü íà ìàññèâ ñ ðàççðÿäàìè ÷èñëà
    cout<<"Enter value(val): ";cin>>val;
    cout<<"Enter number(k) : ";cin>>k;
    //Óçíàëè ÷èñëî ðàçðÿäîâ
    arr = splitDigit(val, arr, n);
    arr = new int[n];//Âûäåëèëè ïàìÿòü ïîä ìàññèâ
    //Çàïèñàëè ðàçðÿäû â ìàññèâ
    arr = splitDigit(val, arr, n);
    rekursBinom(0, arr, n, k);
    return 0;
}
 
int * splitDigit(int val, int *arr, int &n)
{
    if(val < 0)
        val *= -1;
    n = 1;// ÷èñëå äàæå 0 âñåãäà åñòü 1 öèôðà
    //Ïîäñ÷¸ò êîë-âà öèôð â ÷èñëå val
    //è åñëè arr âàëèäíûé òî çàïîëíåíèå åãî ðàçðàäàìè
    if(arr)
        arr[0] = val % 10;//Çàïèñûâàåì òåêóùèé ðàçðÿä
    while(0 < (val /= 10))
    {
        if(arr)
            arr[n] = val % 10;
        n = n + 1;//Óâåëè÷èâàåì ñ÷¸ò÷èê ðàçðÿäîâ
    }
    return arr;
}
 
void rekursBinom(int pos, int *arr, int n, int k)
{
    int i, j;
    int val = arr[pos];
    for(i = 0; i < n;     i++)
    {
        if(i == pos)
            continue;
        val = arr[pos];
        for(j = 0; j < k - 1; j++)
        {
            val *= 10;
            if(i + j < n)
                val += arr[i +     j];
            else
                val += arr[0 + n - j];  
        }
        cout<<val<<endl;
    }
    while((pos += 1) < n)
        rekursBinom(pos, arr, n, k);
}