johtek

-- {A, R, S, B, K, L, M, V, H}
P = {0.90, 0.90, 0.95, 0.95, 0.99, 0.99, 0.009, 0.0032}
T = {}
RBnK = {}
RSB = {}
nRSB = {}
H = {}
HM = {}
for i = 1, 100000 do
  local t = {}
  t.A = math.random() < P[1] and true or false
  t.R = t.A and math.random() < P[2] and true or false
  t.S = t.A and math.random() < P[3] and true or false
  t.B = math.random() < P[4] and true or false
  t.K = t.S and t.B and math.random() < P[5] and true or false
  t.L = t.K and math.random() < P[6] and true or false
  t.M = math.random() < P[7] and true or false
  t.V = math.random() < P[8] and true or false
  if t.M and not t.V then t.H = math.random() < 0.81 and true or false end
  if not t.M and t.V then t.H = math.random() < 0.92 and true or false end
  if t.M and t.V then t.H = math.random() < 0.97 and true or false end
  if not t.M and not t.V then t.H = math.random() < 0.0095 and true or false end
  T[#T+1] = t
  if t.R and t.B and not t.K then RBnK[#RBnK+1] = t end
  if t.R and t.S and t.B then RSB[#RSB+1] = t end
  if not t.R and t.S and t.B then nRSB[#nRSB+1] = t end
  if t.H then H[#H+1] = t end
  if t.H and t.M then HM[#HM+1] = t end
end

local a = 0
for k,v in ipairs(RBnK) do
  if v.A then a = a+1 end
end
print(a/#RBnK)
local a = 0
for k,v in ipairs(RSB) do
  if v.K then a = a+1 end
end
print(a/#RSB)
local a = 0
for k,v in ipairs(nRSB) do
  if v.K then a = a+1 end
end
print(a/#nRSB)
local a = 0
for k,v in ipairs(H) do
  if v.V then a = a+1 end
end
print(a/#H)
local a = 0
for k,v in ipairs(HM) do
  if v.V then a = a+1 end
end
print(a/#HM)

print("done")