bitwise18b.rb

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

end


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

def len01(ns)
	nl = ns.length
	l1 = ns.split(/0+/).map { |x| x.length }
	
	l2 = ns.split(/1+/)[1..-1]
	l2 = [''] if (l2.nil?)
	
	l0 = l2.map { |x| x.length }
	return l0, l1
end

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

def slopediff(l, l2)
	return 0 if (l.nil? || l.empty?)
	
	l = l.map { |x| Math.log(x) }
	
	a, b = [l[0], l[-1]]
	
	return (0...l.size).map \
	{
		|x|
		
		y = a + x.to_f / (l.size - 1) * (b - a)
		
		d = (y - l[x]).abs / (l[x] == 0 ? 1 : l[x])
		l2 << {'x' => l[x], 'y' => y, 'd' => d}
		d
	}.max
	
end

def count(x, l2 = [])
	n = n1 = x['n']
	l = [n1]
	
	while (n != 1)
		n = f2(n)
		
		l << n
		
	end
	
	d = d(ns = n1.to_s(2))
	
	cg = l.index { |x| x < n1 }
	cg = 0 if (cg.nil?)
	cm = (0..cg).max_by { |x| l[x] }
	cg1 = l.rindex { |x| x >= n1 }
	cm1 = (0...l.size).max_by { |x| l[x] }

	nl = ns.length
	ms = slopediff(l[nl..-1], l2)
	
	return x.merge({'da' => (d - 0.5).abs,
		'ns' => ns,
		'cm' => cm,
		'cg' => cg, 
		'cm1' => cm, 
		'cg1' => cg,
		'c' => l.size,
		'nl' => nl,
		'd' => d,
		'e' => e(ns),
		'ms' => ms
		})
end

def out(a, a2 = {}, t = nil)

	fn = 'out.txt'
	k = a.keys
	
	if ($f1.nil?) then
		$f1 = File.open(fn, 'w')
		$f1.puts(k.join("\t")) 
		
		f2 = File.open('gnuplot.cmd', 'w')
		st = t.nil? ? 'unset title; ' : "set title '#{t}'; "
		f2.puts("set ytics nomirror; set y2tics; set colors classic; set key top left box opaque; #{st}; plot \\")
		(a.keys - ['t']).each \
		{
			|x|
			f2.puts("'#{fn}' using (column('t')):(column('#{x}')) with line title '#{x}' #{a2[x]}, \\")
		}
		f2.puts
		f2.close
	end
	$f1.puts(a.values_at(*k).join("\t"))
	$f1.flush
end

def next2(z, l, ks, t)
	
	fmt = {}
	ks.each { |x| fmt[x] = 'lw 2' }
	['da', 'd', 'e'].each { |x| fmt[x] = 'axes x1y2 dt 3' }
	
	out(z.select { |k, v| !['n', 'ns', 'p'].member?(k) }, 
		fmt, t.inspect) if (z['t'] % 10 == 0)
	
	p = z['p'] + 1
	l << [z, count({'n'=>z['n'] + 2 ** p, 'p'=>p, 't' => z['t']}), z.merge({'p'=>p})]
	return
end

def maxs(l)
	l.map { |x| x[-1, 1] == '-' ? x.dup.chop : x }
end

def hard(ks)

	$stderr.puts(ks.inspect)
	$f1 = nil
	l = []
	mx = maxs(ks)
	
	next2(count({'n' => 1, 'p' => 0, 't' => 0}), l, mx, ks)

	add = nil
	seen = {}
	f = File.open($fn, 'w')
	c = 20
	
	5e3.to_i.times \
	{
		|t|
		
		l2 = (0...l.size).map { {} }
		
		mx.each \
		{
			|k2|
			
			l1 = (0...l.size).map { |x| {k2 => l[x][1][k2], 'i' => x} }
			
			l1.sort_by! { |x| x[k2] }

			l1 = l1.reverse if (ks.member?(k2))
			
			(0...l1.size).each { |x| l2[l1[x]['i']][k2] = x }
		}
		l2.each { |x| x['z'] = x.values.inject(0) { |t, x| t + x } }
		
		l1 = (0...l2.size).sort_by { |x| l2[x]['z'] }
		
		a = {'t1' => t, 'l1' => l.size, 'c' => seen.size}.merge(l[l1[0]][1])

		if (!seen.member?(a['n'])) then
			seen[a['n']] = nil
			f.puts(a)
			f.flush
		end
		
		
		$stderr.puts(a.select { |k, | !['n'].member?(k)}) if ((t + 1) % 500 == 0)

		i = l1[rand([l1.size, c].min)]
		z = l.delete_at(i)

		[z[1], z[2]].each { |x| next2(x.merge({'t' => t}), l, mx, ks) }

		l.delete_at(l1[-1]) if (l.size > 2000) 
	}
	f.close
end

def out2(l)
	f = File.open('gnuplot2.cmd', 'w')
	f.puts('$dat << eof')
	f.puts(l[0].keys.join("\t"))
	l.each { |x| f.puts(x.values.join("\t")) }
	f.puts('eof')
	f.puts('plot \\')
	l[0].keys.each { |x| f.puts("$dat using (column('#{x}')) with line,\\") }
	f.close
	
end

def review()
	l = File.open($fn).readlines.map { |x| Kernel.eval(x) }
	l.sort_by! { |x| -x['ms'] }
	x= l[0]
	count({'n' => x['n']}, l2 = [])
	out2(l2)
	exit
	
	
end


review() if (File.exist?($fn = 'mixdb.txt'))

hard(['ms'])
#hard(['m', 'da-'])