Pure Python Gumbel Dice

d12.py

import dice
print(dice.roll(12))

d20.py

import dice
print(dice.roll(20))

d6.py

import dice
print(dice.roll(6))

dice.py

import math
from random import random

def gumbel(eps=1e-6):
    U = random()
    return -math.log(-math.log(U+eps)+eps)

def roll(sides):
    G = [gumbel() for i in range(sides)]
    return [i for g,i in zip(G, range(sides)) if g == max(G)][0]+1