review139b.rb

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

end

def outin()
        l = File.open('out.txt').readlines
        k = l.shift.split
        l2 = l.map { |x| Hash[[k, x.split.map { |x| Kernel.eval(x) }].transpose] }
	return l2
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 left 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?)
		f.puts("'$dat' using (column('#{x}')) with line #{a.fetch(x, '')} lw 2,\\") 
	}
	f.puts
#	f.puts("reset; pause -1;")

end


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

def len01x(ns)
	return lenx(ns, /0+/), lenx(ns, /1+/)
end

def seq(n)

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

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

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

def dp(l, b)
	c = c1 = c2 = 0
	l.each \
	{
		|x|
		d = d(x)
		d1 = d < 0.5 ? 0 : 1
		x1 = x[-1, 1].to_i
		c += 1 if (d1 == x1)
		c1 += d1
		c2 += x1
	}
	return {"mr#{b}" => c.to_f / l.size, "dp#{b}" => c1.to_f / l.size} #, "pl#{b}" => c2.to_f / l.size}
end

l = outin()
l.sort_by! { |x| x['cg'] }

l3 = []
(1..100).each \
{
	|i|

	x = l[-i]
	l1 = seq(x['n'])
	l2 = l1.map { |x| x.to_s(2) }
	
	cm = x['cm']
	cg = x['cg']

	sl1 = l2[0...cm]
	sr1 = l2[cm..-1]
	v1 = dp(sl1.select { |x| x[-1,1] = '1'}, '1l')
	v2 = dp(sl1.select { |x| x[-1,1] = '0'}, '0l')
	v3 = dp(sr1.select { |x| x[-1,1] = '1'}, '1r')
	v4 = dp(sr1.select { |x| x[-1,1] = '0'}, '0r')
	
	sl = sl1.map { |x| x[-50..-1] }.join
	sr = sr1.map { |x| x[-50..-1] }.join
	
	dl = d(sl)
	dr = d(sr)
	drl = dr - dl
	
	el = e(sl)
	er = e(sr)
	erl = er - el
	
	a = v1.merge(v2).merge(v3).merge(v4)
	a = a.merge({'dl' => dl, 'dr' => dr, 'drl' => drl})
	a = a.merge({'el' => el, 'er' => er, 'erl' => erl})
	
	l3 << a
	$stderr.puts(i)
}

outa(f = File.open('gnuplot.cmd', 'w'), l3, {'drl' => 'axes x1y2', 'erl' => 'axes x1y2'})
f.close