construct127.rb

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

end

def seqg(n)

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

def triangle1(c)

	l = []
	c.times \
	{
		|nw|
		ns = '1' * (nw + 2)
		l << ns.to_i(2)
	}
	return l
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.shift if (l[0] == '')
	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 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 key top right opaque; set title '#{t}'; ")
	f.puts("set ytics nomirror; set y2tics;")
	f.puts("plot \\")
	k.each \
	{ 
		|x| 
		next if (a.member?(x) && a[x].nil?)
		opt = a.fetch(x, '')
		opt += ' lw 2 ' if (!opt.include?('lw'))
		opt += ' with line ' if (!opt.include?('with'))
		f.puts("'$dat' using (column('#{x}')) #{opt},\\") 
	}
	f.puts
#	f.puts("reset; pause -1;")

end

def outafn(l, a = {}, fn = nil)
	outa(f = File.open(fn.nil? ? 'gnuplot.cmd' : fn, 'w'), l, a)
	f.close
end

def avg(l)
	l.inject { |t, x| t + x }.to_f / l.size
end


def data(l1)
	l2 = l1.map { |x| (d(x) - 0.5).abs }
	dmx = l2.max
	da = avg(l2)
	
	l2 = l1.map { |x| (e(x) - 0.5).abs }
	emx = l2.max
	ea = avg(l2)
	
	return {'dmx' => dmx, 'da' => da, 'emx' => emx, 'ea' => ea}
end

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 random(l1)
	return l1.map { |x| dense(x.length, 0.5) }
	
end

def seqg2(x)
	l = seqg(x)
	nw = l[0].to_s(2).length
	
	l1 = l[nw..-1].map { |x| x.to_s(2) }
	return l1
end


l = []
l2 = []
triangle1(200).each \
{
	|x|
	
	l << data(l1 = seqg2(x))
	l2 << data(random(l1))
}
outafn(l + [{}] + l2)