Как привязать Шаблон к main() для поиска k-ой статистики

#include <iostream>
#include <vector>
#include <cstring>
using namespace std;
 
template < class T>
 
int main()
{
    int n;
    cout << "Введите размерность: "; cin >> n;
    int *mas = new int[n];
    for (int i = 0; i < n; i++)
        cin >> mas[i];
 
    int k; cout << "Введите k-ый эл-т: "; cin >> k;
 
    system("pause");
    return 0;
}
 
template <class T>
T order_statistics(vector<T> a, unsigned n, unsigned k)
{
    unsigned l, r;
    using std::swap;
    for (l = 1, r = n;;)
    {
 
        if (r <= l + 1)
        {
            // текущая часть состоит из 1 или 2 элементов -
            //   легко можем найти ответ
            if (r == l + 1 && a[r] < a[l])
                swap(a[l], a[r]);
            return a[k];
        }
 
        // упорядочиваем a[l], a[l+1], a[r]
        unsigned mid = (l + r) >> 1;
        swap(a[mid], a[l + 1]);
        if (a[l] > a[r])
            swap(a[l], a[r]);
        if (a[l + 1] > a[r])
            swap(a[l + 1], a[r]);
        if (a[l] > a[l + 1])
            swap(a[l], a[l + 1]);
 
        // выполняем разделение
        // барьером является a[l+1], т.е. медиана среди a[l], a[l+1], a[r]
        unsigned
            i = l + 1,
            j = r;
        const T
            cur = a[l + 1];
        for (;;)
        {
            while (a[++i] < cur);
            while (a[--j] > cur);
            if (i > j)
                break;
            swap(a[i], a[j]);
        }
 
        // вставляем барьер
        a[l + 1] = a[j];
        a[j] = cur;
 
        // продолжаем работать в той части,
        //   которая должна содержать искомый элемент
        if (j >= k)
            r = j - 1;
        if (j <= k)
            l = i;
 
    }
}

template <class T>
T order_statistics(vector<T> a, unsigned n, unsigned k); 
 
int main()
{
    int n;
    cout << "Введите размерность: "; 
    cin >> n;
    vector<int> mas(n);
    for (int i = 0; i < n; i++)
        cin >> mas[i];
    int k; cout << "Введите k: "; cin >> k;
    int res=order_statistics(mas, n, k);
    cout<<res; 
    system("pause");
    return 0;
}