spread6.rb

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

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 seq(n)
	l = [n]
	while (n != 1)
		n = f2(n)
		l << n
	end
	return l
	
end

def norm(l)
	c = 100
	l2 = []
	c.times \
	{
		|x|
		i = (x.to_f / (c - 1) * (l.size - 1)).to_i
		l2 << l[i]
	}
	return l2
end

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

def coef(x, y)

        x1 = x2 = y1 = y2 = xy = 0.0

        n = x.size
        n.times \
        {
                |i|

                x1 += x[i]
                y1 += y[i]
                xy += x[i] * y[i]
                x2 += x[i] ** 2
                y2 += y[i] ** 2
        }

        d = n * x2 - x1 ** 2
        b = (y1 * x2 - x1 * xy) / d
        m = (n * xy - x1 * y1) / d

	b = 1 if (b.nan?)
	m = 1 if (m.nan?)
        return b, m
end

def spread(w)
	
	l2 = []
	200.times \
	{
		n = dense(w, 0.5).to_i(2)
		l = seq(n).map { |x| Math.log(x) }
		b, m = coef((0...l.size).to_a, l)
		l1 = []
		l.size.times \
		{
			|x|
			y1 = l[x]
			y2 = x * m + b
			l1 << (y1 - y2).abs
		}
		
		l2 << norm(l1)
	}

	l1 = l2.transpose.map { |x| avg(x) }

	return l1
	
end

(1..10).each \
{
	|x|
	
	w = x * 100
	l1 = spread(w)
	mx = l1.max
	l1.size.times { |x| puts([x, w, l1[x], l1[x] / mx].join("\t")) }
	puts
}