construct51.rb

def f2(n)
	return n.odd? ? (n * 3 + 1) / 2 : n / 2

end

def adj(ns)

	ns1 = ns.dup
	ns1[0, 1] = '10'
	
	ns2 = ns.dup
	ns2[0, 1] = '11'
	
	return ns1, ns2
end

def seq(n)

	n1 = n
	l = []
	begin
		l << n
		n = f2(n)
	end while (n != 1)
	
	return l
end

def adv(l)
	
	l2 = []

	l.each \
	{
		|x|
		l2 += adj(x)
	}

	return l2
end

def glide(ns)
	n1 = n = ns.to_i(2)
	ns.length.times \
	{
		n = f2(n)
		return if (n < n1)
	}
	return ns
	
end

def filter(l)

	l2 = []
	l.each \
	{
		|x|
		l2 << x if (glide(x))
	}
	return l2
end

def d(s)
        c = s.split('').select { |x| x == '1' }.size
	d = c.to_f / s.length
        return d
        
end

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

def len01(ns)
	return len(ns, /0+/), len(ns, /1+/)
end

def e(ns)
	return len01(ns).flatten.size.to_f / ns.length
end

def stat(l)
	return 0, 0, 0 if (l.empty?)
	
	t = l.inject { |a, x| a + x }
	t2 = l.inject { |a, x| a + x ** 2 }
	
        c = l.size
        a = t.to_f / c
        z = t2.to_f / c - a ** 2
        sd = Math.sqrt(z < 0 ? 0 : z)
        raise if (sd.nan?)
        return {'a' => a, 'sd' => sd, 'mn' => l.min, 'mx' => l.max}
end

def stats(x, l)
	l2 = []
	l.each \
	{
		|x|
		n = x.to_i(2)
		l2 << {'e' => e(x), 'd' => d(x) }
	}
	h2 = {'x' => x}
	l2[0].keys.each \
	{
		|k|
		h = stat(l2.map { |x| x[k] })
		h2 = h2.merge(Hash[[h.keys.map { |x| "#{k}_{#{x}}" }, h.values].transpose])
		
	}
	$stderr.puts([x, l.size].inspect)
	return h2
end

def outa(f, l, a = {}, t = '')
	f.puts('$dat << eof')
	k = l[0].keys
	
	f.puts(k.join("\t"))
	l.each { |x| f.puts(x.values.join("\t")) }

	f.puts('eof')
	f.puts("set colors classic; set ytics nomirror; set y2tics; set title '#{t}'; plot \\")
	k.each \
	{ 
		|x| 
		next if (a.member?(x) && a[x].nil?)
		f.puts("'$dat' using (column('#{x}')) with line #{a.fetch(x, '')} lw 2,\\") 
	}
	f.puts
#	f.puts("reset; pause -1;")

end


l = ['1']
l2 = []
23.times { |x| l = filter(adv(l)); l2 << stats(x, l); }

outa(f = File.open('gnuplot.cmd', 'w'), l2, {'x' => nil})
f.close