Calkin-Wilf sequence
func calkin_wilf(n) is cached {
return 1 if (n == 1)
1/(2*floor(__FUNC__(n-1)) + 1 - __FUNC__(n-1))
}
func r2cw(r) {
var cfrac = r.as_cfrac
cfrac.len.is_odd || return nil
Num(cfrac.flip.map_kv {|k,v| (k.is_odd ? '0' : '1') * v }.join, 2)
}
with (20) {|n|
say "First #{n} terms of the Calkin-Wilf sequence:"
say calkin_wilf.map(1..n)
}
with (83116/51639) {|r|
say ("\n#{r.as_rat} is at index: ", r2cw(r))
}
Output:
First 20 terms of the Calkin-Wilf sequence:
[1, 1/2, 2, 1/3, 3/2, 2/3, 3, 1/4, 4/3, 3/5, 5/2, 2/5, 5/3, 3/4, 4, 1/5, 5/4, 4/7, 7/3, 3/8]
83116/51639 is at index: 123456789
Last updated