import math; from decimal import Decimal m=200; s=0.2; t=0.5; m1s=int(m*(1-s)); mt=int(m*t) attackrisk1_inv = (1 - Decimal(math.comb(m1s, mt))/Decimal(math.comb(m, mt))) for c in [2,3,4,5,6,7,8,9,10,20,30,40,50,60,70,80,90,100,200,300,400,500,600,700]: result=[] for c_ in range(1,200): fc_ = float(c_)/100 if (1 - Decimal(math.comb(int(m1s*fc_), int(mt*fc_)))/Decimal(math.comb(int(m*fc_), int(mt*fc_))))**c - attackrisk1_inv<0: result=[c,fc_] continue print("c:",result[0]," c':", result[1])