Sequencia de farey

[Maria-luizah]

Encontrei a solução da sequência de farey para python, porem não entendi muito bem, alguém pode me explicar por favor?

def farey( n ):
    a, b, c, d = 0, 1, 1, n
    print "%d/%d" % (a,b)
    while (c < n):
        k = int((n + b)/d)
        a, b, c, d = c, d, k*c - a, k*d - b
        print "%d/%d" % (a,b)

[Jhonas]

explicação academica:

http://www.dma.uem.br/kit/textos/farey/farey.html

abraço