construct177.rb

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

end

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

def sum(l)
	return l.inject { |t, x| t + x }
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 midpt(ns)
	w2 = ns.length / 4
	l = ns.split('')
	i = j = 0
	while (i < w2 && j < l.size)
		i += l[j].to_i
		j += 1
	end
	return j.to_f / ns.length
end

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

def dl(ns)
	w2 = ns.length / 2
	nw1 = ns[-w2..-1]
	return d(nw1)
	
end

def dh(ns)
	w2 = ns.length / 2
	nw1 = ns[0..w2]
	return d(nw1)
	
end

def dp(x)
	return x
	return (x > 0.5) ? 1 : 0
end

def outd(f, l)
	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')
	return k
end

def outa(f, l, a = {}, t = '')
	
	k = outd(f, l)
	f.puts("set colors classic; set title '#{t}'; ")
	f.puts("set key bottom right;")
	f.puts("set ytics nomirror; set y2tics;")
	f.puts("plot \\")
	k, ct = [k - ['t'], "(column('t')):"] if (k.member?('t'))
	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') && !opt.include?('pt'))
		f.puts("'$dat' using #{ct}(column('#{x}')) #{opt} title '#{x}',\\") 
	}
	f.puts
	f.puts("set title; set ytics; set y2tics; set key;")
#	f.puts("reset; pause -1;")

end

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

def midpt(ns)
	w2 = ns.length / 4
	l = ns.split('')
	i = j = 0
	while (i < w2 && j < l.size)
		i += l[j].to_i
		j += 1
	end
	return j.to_f / ns.length
end

def avg(l)
	return l.empty? ? 0 : sum(l).to_f / l.size
end

def runavg(l, k, c)

	l1 = l.map { |x| x[k] }

	t = sum(l1[0...c])
	l2 = (['-'] * (c - 1)) + [t.to_f / c]
	while (l1.size > c)
		t -= l1.shift
		t += l1[c - 1]
		l2 << (t.to_f / c)
	end
	ka = "#{k}a"
	
	l.each_with_index { |x, i| x[ka] = l2[i] }
	
end

def tavg(l, k)
	t = 0.0
	l.size.times \
	{
		|j|
		t += l[j][k]
		
		l[j]["#{k}a"] = t / (j + 1)
	}
	
end

def sample(l, n)
	
	d = (l.size - 1).to_f / n
	l2 = []
	(n + 1).times \
	{
		|i|
		j = (i * d).to_i
		l2 << l[j]
	}
	return l2
end

def spaceno(l)

	l2 = []
	t = 0
	l.each \
	{
		|l1|
		l1.each { |x| x['t'] = (t += 1) }
		l2 += l1 + [{}]
		t += 500
	}
	return l2
end



srand(sr = ARGV[0].to_i)
ac = 1000


ns = dense(1000, 0.5)
l = seq(ns.to_i(2))
l2 = []
l.each \
{
	|x|
	ns = x.to_s(2)
	d = d(ns)
	dl = dl(ns)
	dh = dh(ns)
	dlh = dl - dh + 0.5
	mp = midpt(ns)
	a = {
		'd' => d,
		'dl' => dl,
		'dh' => dh,
		'dlh' => dlh,
		'mp' => mp,
		'y' => Math.log(x) / Math.log(2)}
	l2 << a
}

l3 = l2.map { |x| x.dup }
l4 = l2.map { |x| x.dup }.reverse
i = l4.index { |x| x['y'] >= 100 }
l4[0..i] = []

(l2[0].keys - ['y']).each \
{
	|k|
	runavg(l2, k, ac)
	tavg(l3, k)
	tavg(l4, k)
	(l2 + l3 + l4).each { |x| x.delete(k) }
}


outafn(spaceno([l2, l3, l4]), {'y' => 'axes x1y2 dt 3'}, {'ac' => ac, 'sr' => sr}.inspect)