if sel<0:
                    k=3
                    T=0
                    M=0
                    P=0
                    Q=0
                    dY=0
                    R=0
                    z=0
                else:
                    k=sel
                    tau=min(par[3],par[8])
                    v0=sum((bm,2*bp))//4
                    v1=sum((2*bm,bp))//4
                    v2=sum((bm,bp))
                    if g==0:
                        vg=v0
                    elif g==1:
                        vg=v1
                    else:
                        vg=v2
                    ell=sum((e,bh,-gamma[g]*u,vg,par[13]//12))
                    qp=sum((par[13],ell))
                    T=sum((tau,de,d))
                    th=T
                    to=par[8-5*h]
                    tau2=min((th,to))
                    Q=max(0,qp-5*(tau2-tau))
                    dY=sum((rm[g]*bm,rp[g]*bp,sigma[g]*(2-u),qp//12,rho*max(0,T-D)))
                    if qp>101:
                        R=0
                        z=1
                    else:
                        n=n 1
                        if g==0:
                            na=sum((a,1))
                            nb=b
                            nc=c
                        elif g==1:
                            na=a
                            nb=sum((b,1))
                            nc=c
                        else:
                            na=a
                            nb=b
                            nc=sum((c,1))
                        if h==0:
                            ch=(T,g,M,P,min(2,sum((u,1))),par[8],par[9],par[10],par[11],par[12])
                        else:
                            ch=(par[3],par[4],par[5],par[6],par[7],T,g,M,P,min(2,sum((u,1))))
                        child=(na,nb,nc) ch (Q,sum((par[14],dY)),sum((par[15],1 if k==2 else 0)),n)
                        ck=child[0:13]
                        dead=False
                        for old in NL:
                            if old[0:13]==ck and dom(old,child):
                                dead=True
                                break
                        if dead:
                            z=2
                            R=0
                        else:
                            rmv=[]
                            for old in NL:
                                if old[0:13]==ck and dom(child,old):
                                    rmv.append(old)
                            R=len(rmv)
                            if R>0:
                                z=4
                                NL=[o for o in NL if not (o[0:13]==ck and dom(child,o))]
                            else:
                                z=3
                            NL.append(child)
                H=mod(sum((37*H,19*j,17*par[16],13*i,11*h,7*k,5*T,3*M,2*P,Q,23*dY,29*R,31*z)),9973)