Lib.rs

This crate implements an iterator returning the tuples (k_n, d_n), where k_n is the nth term of the Kolakoski sequence and d_n = ∑_i in 1..=n (−1)^k_i is the “Kolakoski discrepancy function”.

use kolakoski_algorithms::Kolakoski;

println!("{:?}", Kolakoski::default().take(20).map(|(k, _)| k).collect::<Vec<_>>());
// [1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 1]
# assert_eq!(&[1, 2, 2, 1, 1, 2, 1, 2, 2, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 1][..], &Kolakoski::default().take(20).map(|(k, _)| k).collect::<Vec<_>>());

If you are interested in analysing the behaviour of the sequence for certain large values of the argument, you can avoid generating all preceding terms:

# use kolakoski_algorithms::Kolakoski;
println!("{:?}", Kolakoski::new(1_000_000_000).take(100).collect::<Vec<_>>());