import itertools

# Capacitores
E6_C = [1.0, 1.5, 2.2, 3.3, 4.7, 6.8]

# resistencias de 5% de error
E24_R = [
    1.0, 1.1, 1.2, 1.3, 1.5, 1.6, 1.8, 2.0, 2.2, 2.4, 2.7, 3.0,
    3.3, 3.6, 3.9, 4.3, 4.7, 5.1, 5.6, 6.2, 6.8, 7.5, 8.2, 9.1
]

E24_R_custom = [
    1.0, 1.1, 1.2, 1.3, 1.5, 1.6, 1.8, 2.0, 2.2, 2.4, 2.7, 3.0,
    3.3, 3.6, 3.9, 4.3, 4.7, 5.1, 5.6, 6.2, 6.8, 7.5, 8.2
]

# resistencias de 2% de error
E48_R = [
    1.00, 1.05, 1.10, 1.15, 1.21, 1.27, 1.33, 1.40, 1.47, 1.54, 1.62, 1.69,
    1.78, 1.87, 1.96, 2.05, 2.15, 2.26, 2.37, 2.49, 2.61, 2.74, 2.87, 3.01,
    3.16, 3.32, 3.48, 3.65, 3.83, 4.02, 4.22, 4.42, 4.64, 4.87, 5.11, 5.36,
    5.62, 5.90, 6.19, 6.49, 6.81, 7.15, 7.50, 7.87, 8.25, 8.66, 9.09, 9.53
]

# resistencias de 1% de error
E96_R = [
    1.00, 1.02, 1.05, 1.07, 1.10, 1.13, 1.15, 1.18, 1.21, 1.24, 1.27, 1.30,
    1.33, 1.37, 1.40, 1.43, 1.47, 1.50, 1.54, 1.58, 1.62, 1.65, 1.69, 1.74,
    1.78, 1.82, 1.87, 1.91, 1.96, 2.00, 2.05, 2.10, 2.15, 2.21, 2.26, 2.32,
    2.37, 2.43, 2.49, 2.55, 2.61, 2.67, 2.74, 2.80, 2.87, 2.94, 3.01, 3.09,
    3.16, 3.24, 3.32, 3.40, 3.48, 3.57, 3.65, 3.74, 3.83, 3.92, 4.02, 4.12,
    4.22, 4.32, 4.42, 4.53, 4.64, 4.75, 4.87, 4.99, 5.11, 5.23, 5.36, 5.49,
    5.62, 5.76, 5.90, 6.04, 6.19, 6.34, 6.49, 6.65, 6.81, 6.98, 7.15, 7.32,
    7.50, 7.68, 7.87, 8.06, 8.25, 8.45, 8.66, 8.87, 9.09, 9.31, 9.53, 9.76
]


SERIE_R = E24_R_custom
SERIE_C = E6_C

# décadas
DECADAS_R = [1, 10, 100, 1_000, 10_000, 100_000]
DECADAS_C = [1e-9, 10e-9, 100e-9, 1e-6]  # pF–uF


# lista de resistencias y capacitores disponibles
resistencias = [v * d for v, d in itertools.product(SERIE_R, DECADAS_R)]
capacitores  = [v * d for v, d in itertools.product(SERIE_C, DECADAS_C)]

# a = 1/(RC)
def find_lower_error(a, res_list=resistencias, cap_list=capacitores):
    resultados = []

    for R, C in itertools.product(res_list, cap_list):
        R_teorica = 1 / (C * a)

        if R_teorica == 0:
            continue

        error_porcentual = abs((R - R_teorica) / R_teorica) * 100
        resultados.append((error_porcentual, R, C, R_teorica))

    # ordenar por error porcentual
    resultados.sort(key=lambda x: x[0])
    return resultados

if __name__ == "__main__":
    top = find_lower_error(250)

    print(f"    R [Ω], R_ideal [Ω],    C [F], Err [%]")
    for i, (err, R, C, R_ideal) in enumerate(top[:10], 1):
        print(f"{R:9.2f}, {R_ideal:11.2f}, {C:8.2e}, {err:7.2f}")

    top_cap_list=[c for (err, R, c, R_ideal) in top[:10]]
    top_2 = find_lower_error(500, cap_list=top_cap_list)

    print("")
    print(f"    R [Ω], R_ideal [Ω],    C [F], Err [%]")
    for i, (err, R, C, R_ideal) in enumerate(top_2[:10], 1):
        print(f"{R:9.2f}, {R_ideal:11.2f}, {C:8.2e}, {err:7.2f}")