After showing around the formula that I was very proud of, I was surprised by how people though it was complex. So for you enjoyment and clarity’s sake:
f(i, j) = z i^(a.log j)
where z = 1/Z(a.log j) for the distribution to sum up to 1, and a a case-dependent constant. Doesn’t it look better this way?