## Infrastructural
from decimal import Decimal as D
import pandas as pa
DF = pa.DataFrame

u_n = D("1")  # utility of nlife  # arbitrary constant here
u_x = u_n * D("1e5")  # utility of xlife  # very humble low guess
p_x = D("1e-6")  # base probability of xlife  # also somewhat low

p_n = 1 - p_x  # probability of nlife  # see assumptions

columns = 'u_n u_x p_n p_x min_p_d'.split()
c_a_x = lambda: p_d > p_n * (u_n / u_x)
min_p_d = p_n * (u_n / u_x)

DF([(u_n, u_x, p_n, p_x, min_p_d)], columns=columns)

values = [
    # u_x, p_x
    [D("1e2"), D("1e-10")],
    [D("1e10"), D("1e-1")],
]
DF([[u_n, u_x, 1 - p_x, p_x, (1 - p_x) * (u_n / u_x)]
    for u_x, p_x in values], columns=columns)