2021-04-05 14:48:01 +00:00
|
|
|
package connections
|
|
|
|
|
|
|
|
import "math/rand"
|
|
|
|
|
|
|
|
// https://www.keithschwarz.com/darts-dice-coins/
|
2021-04-05 17:10:03 +00:00
|
|
|
type Sampler struct {
|
2021-04-05 14:48:01 +00:00
|
|
|
randomGenerator *rand.Rand
|
|
|
|
probabilities []float64
|
|
|
|
alias []int
|
|
|
|
}
|
|
|
|
|
|
|
|
type workList []int
|
|
|
|
|
|
|
|
func (wl *workList) push(e int) {
|
|
|
|
*wl = append(*wl, e)
|
|
|
|
}
|
|
|
|
|
|
|
|
func (wl *workList) pop() int {
|
|
|
|
l := len(*wl) - 1
|
|
|
|
n := (*wl)[l]
|
|
|
|
*wl = (*wl)[:l]
|
|
|
|
return n
|
|
|
|
}
|
|
|
|
|
2021-04-05 17:10:03 +00:00
|
|
|
func NewSampler(probabilities []float64, source rand.Source) *Sampler {
|
2021-04-05 17:13:49 +00:00
|
|
|
sampler := &Sampler{}
|
2021-04-05 14:48:01 +00:00
|
|
|
var (
|
|
|
|
small workList
|
|
|
|
large workList
|
|
|
|
)
|
|
|
|
n := len(probabilities)
|
2021-04-05 17:13:49 +00:00
|
|
|
sampler.randomGenerator = rand.New(source)
|
|
|
|
sampler.probabilities = make([]float64, n)
|
|
|
|
sampler.alias = make([]int, n)
|
2021-04-05 14:48:01 +00:00
|
|
|
// Compute scaled probabilities.
|
|
|
|
p := make([]float64, n)
|
|
|
|
for i := 0; i < n; i++ {
|
|
|
|
p[i] = probabilities[i] * float64(n)
|
|
|
|
}
|
|
|
|
for i, pi := range p {
|
|
|
|
if pi < 1 {
|
|
|
|
small = append(small, i)
|
|
|
|
} else {
|
|
|
|
large = append(large, i)
|
|
|
|
}
|
|
|
|
}
|
|
|
|
for len(large) > 0 && len(small) > 0 {
|
|
|
|
l, g := small.pop(), large.pop()
|
2021-04-05 17:13:49 +00:00
|
|
|
sampler.probabilities[l] = p[l]
|
|
|
|
sampler.alias[l] = g
|
2021-04-05 14:48:01 +00:00
|
|
|
p[g] = (p[g] + p[l]) - 1
|
|
|
|
if p[g] < 1 {
|
|
|
|
small.push(g)
|
|
|
|
} else {
|
|
|
|
large.push(g)
|
|
|
|
}
|
|
|
|
}
|
|
|
|
for len(large) > 0 {
|
|
|
|
g := large.pop()
|
2021-04-05 17:13:49 +00:00
|
|
|
sampler.probabilities[g] = 1
|
2021-04-05 14:48:01 +00:00
|
|
|
}
|
|
|
|
for len(small) > 0 {
|
|
|
|
l := small.pop()
|
2021-04-05 17:13:49 +00:00
|
|
|
sampler.probabilities[l] = 1
|
2021-04-05 14:48:01 +00:00
|
|
|
}
|
2021-04-05 17:13:49 +00:00
|
|
|
return sampler
|
2021-04-05 14:48:01 +00:00
|
|
|
}
|
|
|
|
|
2021-04-05 17:10:03 +00:00
|
|
|
func (g *Sampler) Next() int {
|
2021-04-05 14:48:01 +00:00
|
|
|
n := len(g.alias)
|
|
|
|
i := g.randomGenerator.Intn(n)
|
|
|
|
if g.randomGenerator.Float64() < g.probabilities[i] {
|
|
|
|
return i
|
|
|
|
}
|
|
|
|
return g.alias[i]
|
|
|
|
}
|