Skip to content

Instantly share code, notes, and snippets.

@vznvzn
Created January 16, 2020 03:14
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save vznvzn/cc01501d7abb234d42a48c3193c0842a to your computer and use it in GitHub Desktop.
Save vznvzn/cc01501d7abb234d42a48c3193c0842a to your computer and use it in GitHub Desktop.
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 f2(n)
return n.odd? ? (n * 3 + 1) / 2 : n / 2
end
def seq(n)
n1 = n
l = []
begin
l << n
n = f2(n)
end while (n != 1)
return l
end
def d(s)
c = s.split('').select { |x| x == '1' }.size
d = c.to_f / s.length
return d
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 len01(ns)
return len(ns, /0+/), len(ns, /1+/)
end
def e(ns)
return len01(ns).flatten.size.to_f / ns.length
end
def outa(f, l, a = {}, c = '')
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 right opaque; ")
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 #{c}(column('#{x}')):(column('c')) with line #{a.fetch(x, '')} linecolor variable lw 2 title '#{x}',\\")
}
f.puts
# f.puts("reset; pause -1;")
end
def stat(k, l)
t = l.inject { |a, x| a + x }
t2 = l.inject(0) { |a, x| a + (x ** 2) }
c = l.size
a = t.to_f / c
z = t2.to_f / c - a ** 2
sd = Math.sqrt(z < 0 ? 0 : z)
raise if (sd.nan?)
return {"#{k}_a" => a, "#{k}_s" => sd}
end
def convert(l)
x2 = nil
l2 = []
l.each_with_index \
{
|x, j|
next if (x == x2)
l2 << [j, x]
x2 = x
}
return l2
end
def scale(l, m)
l1 = l.map { |x| x[m] }
mn = l1.min
mx = l1.max
l.each { |x| x[m] = (x[m].to_f - mn) / (mx - mn) }
end
def review(l, c)
l1, l0 = len01(l.map { |x| x.to_s(2) }.join)
convert(l1.sort).each \
{
|x, y|
puts([x, y, c].join("\t"))
}
puts
end
ns = dense(200, 0.5)
l = seq(ns.to_i(2))
l1 = []
l2 = []
(1...l.size).each \
{
|x|
next if (!(l[x - 1].odd? && l[x].even?))
l1 << l[x - 1]
l2 << l[x]
}
review(l1, 1)
review(l2, 2)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment