Código C++ – Ordenamiento Radix

#include <iostream>
#include <cmath>
using namespace std;
#define NUMELTS 20

void radixsort(int x[], int n)
{
    int front[10], rear[10];

    struct {
        int info;
        int next;
    } node[NUMELTS];

    int exp, first, i, j, k, p, q, y;

    /* Inicializar una lista vinculada */
    for (i = 0; i < n - 1; i++)
    {
        node[i].info = x[i];
        node[i].next = i + 1;
    } /* fin del for */

    node[n - 1].info = x[n - 1];
    node[n - 1].next = -1;
    first = 0; /* first es la cabeza de la lista vinculada */

    for (k = 1; k < 5; k++)
    {
        /* Suponer que tenemos números de cuatro dígitos */
        for (i = 0; i < 10; i++)
        {
            /* Inicializar colas */
            rear[i] = -1;
            front[i] = -1;
        } /* fin del for */

        /* Procesar cada elemento en la lista */
        while (first != -1)
        {
            p = first;
            first = node[first].next;
            y = node[p].info;
            /* Extraer el k-ésimo dígito */
            exp = pow(10, k - 1); /* elevar 10 a la (k-1)-ésima potencia */
            j = (y / exp) % 10;
            /* Insertar y en queue[j] */
            q = rear[j];
            if (q == -1)
                front[j] = p;
            else
                node[q].next = p;
            rear[j] = p;
        } /* fin del while */

        /* En este punto, cada registro está en su cola basándose en el dígito k
           Ahora formar una lista única de todos los elementos de la cola.
           Encontrar el primer elemento. */
        for (j = 0; j < 10 && front[j] == -1; j++)
            ;
        first = front[j];

        /* Vincular las colas restantes */
        while (j <= 9)
        { /* Verificar si se ha terminado */
            /* Encontrar el elemento siguiente */
            for (i = j + 1; i < 10 && front[i] == -1; i++)
                ;
            if (i <= 9)
            {
                p = i;
                node[rear[j]].next = front[i];
            } /* fin del if */
            j = i;
        } /* fin del while */
        node[rear[p]].next = -1;
    } /* fin del for */

    /* Copiar de regreso al archivo original */
    for (i = 0; i < n; i++)
    {
        x[i] = node[first].info;
        first = node[first].next;
    } /* fin del for */
} /* fin de radixsort */

int main(void)
{
    int x[50] = {0}, i;
    static int n;

    cout << "Cadena de números enteros:\n";
    for (n = 0;; n++)
    {
        cin >> x[n];
        if (x[n] == -1)
            break;
    }

    if (n)
        radixsort(x, n);
    for (i = 0; i < n; i++)
        cout << x[i] << endl;

    return 0;
}