longclaps

1. Ищется функция G, обратная функции распределения.
   
2. Берутся N равномерно распределённых на [0, 1] случайных чисел x (их даёт обычный random()).
   
3. G(x) будут иметь нужное распределение.
   

from random import randint

def shift_rand(a, b, x, f):
    t = randint(a, (b - x) * f + x)
    return t if t <= x else (t - x) // f + x

нужный_процент = 0.33 // например
for (i = 0; i < n ; i++) {
   x = random_от_нуля_до_единицы();
   if (x > нужный_процент) x = 0;
   else x /= нужный_процент; // ненулевые равномерно распределены на [0..1]
}