construct180c.rb

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

end

def adj(x, m, n, s2, p2, p, x3)
	s2.replace(p[0...x] + (x3 % 2).to_s)
	m = (s2[x, 1] == p[x, 1])
	n = p2.reverse.to_i(2)
	p2[x, 1] = m ? '11' : '01'
	return m, n, x3
end

def terras121(p)
	
	p2 = ['01', '11'][p[0, 1].to_i]
	n = 1
	n3 = 0
	m = p2 == '11'
	x3 = (p2.reverse.to_i(2)) >> 1

	(1...p.length).each \
	{
		|x|

		if (!m)
			ns = n.to_s(2)
			ns[0, 1] = ''
			m1 = x3 - 3**n3
			n1 = f2(m1)
		end

		n3 = p[0...x].split('').select { |z| z == '1' }.size
		
		s2 = ''
		m, n, x3 = adj(x, m, n, s2, p2, p, m ? (f2(x3) + 3**n3) : (n1 + 3**n3))
	}

	n = p2.reverse.to_i(2)

	return n
end

def seq(n)
	l = [n]
	
	while (n != 1)
		n = f2(n)
		
		l << n
		
	end
	return l
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 len(ns, p)
	l = ns.split(p)
	l = [] if (l.nil?)
	l.shift if (l[0] == '')
	return l.map { |x| x.length }
end

def len1(ns)
	return len(ns, /0+/)
end

def len0(ns)
	return len(ns, /1+/)
end

def len01(ns)
	return len1(ns), len0(ns)
end

def e(ns)
	return len01(ns).flatten.size.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 sum(l)
	return l.inject { |t, x| t + x }
end

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

def add(a, x)
	a.replace({'x' => 0.0, 'x2' => 0.0, 'k' => x, 'n' => 0}) if (a.empty?)
	a['x'] += x - a['k']
	a['x2'] += (x - a['k']) ** 2
	a['n'] += 1
end

def meanvar(a, k = '')
	{"a#{k}" => a['k'] + a['x'] / a['n'], 
	"sd#{k}" => (a['x2'] - a['x'] ** 2 / a['n']) / (a['n'] - 1)}
end

def out(a)
	if ($x.nil?) then
		$x = true
		puts(a.keys.join("\t"))
	else
		puts(a.values.join("\t"))
	end
	$stdout.flush
end

e1 = {}
e2 = {}
t = 0
loop \
{

	p = dense(500, 0.64)
	n1 = terras121(p)


	ns = dense(500, 0.5)

	add(e1, e(n1.to_s(2)))
	add(e2, e(ns))

	t += 1
	out(meanvar(e1, 1).merge(meanvar(e2, 2)))
	$stdout.flush
}


# plot 'out.txt' using (column('a1')):(column('sd1')) with errorbar,'' using (column('a2')):(column('sd2')) with errorbar