La verdad que http://www.vitutor.com/pro/1/a_6.html es todo lo que necesite para entenderlo..
Pero al caso, << es para para correr el numero binario tantos como son necesarios hacia el lado que indica el simbolo, o sea, >> es su operacion contraria.
Si estuvieramos en C++ esto significaria el desplazamiento de flujo que funciona distinto dependiendo del tipo/clase

En cuanto al algoritmo de generacion, la parte mas complicada es la de seleccionar la posicion a la cual incrementar al siguiente elemento del grupo..

Fijate si con la lectura que te pase lo entendes sino te lo comento..

*El operador >> sirve para lo que se denomina 'Desplazamiento de Bits' Esto sirve para lo que su nombre dice. Me explico con un ejemplo: tenemos una variable de tipo entero (int) cuyo valor es 12. Dicho numero en binario es 1100. Si hacemos la siguiente operación...

varialbe = variable >> 1

entonces estamos desplazando todos los bits un lugar, dando como resultado el numero en binario 0110. El 1 que estaba en cuarto lugar ha pasado al tercer lugar. El 1 que estaba en tercer lugar ha pasado al segundo lugar... y así sucesivamente. Si hubiéramos hecho esto...

variable = variable >> 2

entonces daría como resultado 0011 en binario

*El operador & es parecido al |. Lo que hace es multiplicar bit a bit dos valores. Por ejemplo, tenemos los valores 011010 y 110011 (en binario). Si hacemos 011010 & 110011 obtendremos 010010.
Empezamos de derecha a izquierda: El primer bit está determinado por la operación 1 * 0, por lo que vale 0. El segundo bit está determinado por la operación 1 * 1, por lo que vale 1 y así sucesivamente

Bien, he regresado xD
Vengo a explicarte por que se usa asi el tema de la posicion:
Examinemos el codigo si?
Lo primero que hacemos es guardar el valor de i, que es el indicador del numero de combinacion general en la que estamos.
Luego entrmos en el for que hace las combinaciones.
El orden se determina por

            gp[gp_len -j -1] = list[n % list_len];
            n /= list_len;
 

gp -> grupo, gp_len -> largo del grupo, j-> posicion actual en el grupo, -1 porque el indice esta basado en ceros
list es la lista de elementos permutables en el grupo.
Entonces:
para cada elemento de gp le en la posicion donde estamos en el bucle (-1) se da el valor del elemento del grupo en relacion al largo del grupo de elementos intercambiables en el round actual (round = valor de i), luego se ajusta n con el largo del grupo deseado.

Esto es asi dado que cada permutacion tiene asociado un numero, si "0000" es 0, "0001" es 1, y "AAAA" es X

No se explicarme muy bien, tendria que agarrar latex y escribirte la formula algebraica de generacion de permutaciones.

Basicamente es que tenes un grupo vacio de universo hasta X elementos.
De tal manera que vas a seleccionar dependiendo del numero un elemento Y para meterlo y puede ser biyectiva(es decir tener dos valores iguales dentro del grupo X partiendo del mismo elemento en Y)