adapt.rb

def dense(w, d)
	w2 = w - 1
	
	a = (0...w2).to_a
	s = '0' * w2
	(1..(d * w - 1)).map { a.delete_at(rand(a.size)) }.each { |x| s[x, 1] = '1' }
	return ('1' + s)
	
end

def len(ns, p)
	l = ns.split(p)
	l = [] if (l.nil?)
	l.shift if (l[0] == '')
	return l.map { |x| x.length }
end

def len1(ns)
	return len(ns, /0+/)
end

def len0(ns)
	return len(ns, /1+/)
end

def len01(ns)
	return len1(ns), len0(ns)
end

def hist(l)
	a = {}
	l.each { |x| a[x] = a.fetch(x, 0) + 1 }
	return a
end

def hdiff(h1, h2)
	d = 0
	(h1.keys | h2.keys).each { |x| d += (h1.fetch(x, 0) - h2.fetch(x, 0)).abs }
	return d
end

def h2(ns)
	l1, l0 = len01(ns)
	return hist(l0), hist(l1)
end

def hdiff2(h2a, h2b)
	return hdiff(h2a[0], h2b[0]) + hdiff(h2a[1], h2b[1])
end

def hdiff1(ns, i)
	
	h2a = h2(ns)
	t = c = 0
	h1 = nil
	loop \
	{
		ns2 = dense(ns.length, 0.5)
		h2b = h2(ns2)
		d = hdiff2(h2a, h2b)
		t += d
		c += 1
		h = t.to_f / c
		e = h1.nil? ? nil : (h - h1) / h
		puts([c, h, i].join("\t"))
		h1 = h
		break if (c == 20)
	}
	puts

end

50.times \
{
	|i|
	hdiff1(dense(100 + 10 * i, 0.5), i)
}