BASIC_1(p):
p1=p,
p0=1-p;

BASIC_2(p):
p2=p,
p1=0,
p0=1-p;

BASIC_1_1(p):
p11=p,
p10=0,
p01=0,
p00=1-p;

BASIC_1_2(p):
p12=p,
p11=0,
p10=0,
p02=0,
p01=0,
p00=1-p;

VIRTUALIZED_1_1(p,pVM):
p11=p*pVM*pVM,
p10=p*pVM*(1-pVM),
p01=p*(1-pVM)*pVM,
p00=(1-p)+p*(1-pVM)*(1-pVM);

VIRTUALIZED_2_1(p,pVM):
p21=p*pVM^3,
p20=p*pVM^2*(1-pVM),
p11=2*p*pVM*(1-pVM)*pVM,
p10=2*p*pVM*(1-pVM)^2,
p01=p*(1-pVM)^2*pVM,
p00=(1-p)+p*(1-pVM)^3;

VIRTUALIZED-M_1_1(p,pd,pm,r,pVM):
f_0_1=1-pVM,
f_1_1=pVM,
g_0_1=((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r),
g_1_1=1-g_0_1,
p11=p*f_1_1*f_1_1 + (1-p)*pd*g_1_1*g_1_1,
p10=p*f_1_1*f_0_1 + (1-p)*pd*g_1_1*g_0_1,
p01=p10,
p00=p*f_0_1*f_0_1 + (1-p)*pd*g_0_1*g_0_1+(1-p)*(1-pd);

VIRTUALIZED-M_2_1(p,pd,pm,r,pVM):
f_0_1=1-pVM,
f_1_1=pVM,
g_0_1=((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r),
g_1_1=1-g_0_1,
f_0_2=(1-pVM)^2,
f_1_2=2*pVM*(1-pVM),
f_2_2=1-f_0_2-f_1_2,
g_0_2=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^2,
g_1_2=2*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r)^2,
g_2_2=1-g_0_2-g_1_2,
p21=p*f_2_2*f_1_1 + (1-p)*pd*g_2_2*g_1_1,
p20=p*f_2_2*f_0_1 + (1-p)*pd*g_2_2*g_0_1,
p11=p*f_1_2*f_1_1 + (1-p)*pd*g_1_2*g_1_1,
p10=p*f_1_2*f_0_1 + (1-p)*pd*g_1_2*g_0_1,
p01=p*f_0_2*f_1_1 + (1-p)*pd*g_0_2*g_1_1,
p00=p*f_0_2*f_0_1 + (1-p)*pd*g_0_2*g_0_1+(1-p)*(1-pd);

VIRTUALIZED_2_2(p,pVM):
f_0_2=(1-pVM)^2,
f_1_2=2*pVM*(1-pVM),
f_2_2=1-f_0_2-f_1_2,
p22=p*f_2_2*f_2_2,
p21=p*f_2_2*f_1_2,
p20=p*f_2_2*f_0_2,
p12=p21,
p11=p*f_1_2*f_1_2,
p10=p*f_1_2*f_0_2,
p02=p20,
p01=p10,
p00=p*f_0_2*f_0_2+(1-p);

VIRTUALIZED_4_2(p,pVM):
f_0_2=(1-pVM)^2,
f_1_2=2*pVM*(1-pVM),
f_2_2=1-f_0_2-f_1_2,
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
p22=p*f_2_4*f_2_2,
p21=p*f_2_4*f_1_2,
p20=p*f_2_4*f_0_2,
p12=p*f_1_4*f_2_2,
p11=p*f_1_4*f_1_2,
p10=p*f_1_4*f_0_2,
p02=p*f_0_4*f_2_2,
p01=p*f_0_4*f_1_2,
p00=p*f_0_4*f_0_2+(1-p);

VIRTUALIZED-M_2_2(p,pd,pm,r,pVM):
f_0_2=(1-pVM)^2,
f_1_2=2*pVM*(1-pVM),
f_2_2=1-f_0_2-f_1_2,
g_0_2=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^2,
g_1_2=2*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r)^2,
g_2_2=1-g_0_2-g_1_2,
p22=p*f_2_2*f_2_2 + (1-p)*pd*g_2_2*g_2_2,
p21=p*f_2_2*f_1_2 + (1-p)*pd*g_2_2*g_1_2,
p20=p*f_2_2*f_0_2 + (1-p)*pd*g_2_2*g_0_2,
p12=p21,
p11=p*f_1_2*f_1_2 + (1-p)*pd*g_1_2*g_1_2,
p10=p*f_1_2*f_0_2 + (1-p)*pd*g_1_2*g_0_2,
p02=p20,
p01=p10,
p00=p*f_0_2*f_0_2 + (1-p)*pd*g_0_2*g_0_2+(1-p)*(1-pd);

VIRTUALIZED-M_4_2(p,pd,pm,r,pVM):
f_0_2=(1-pVM)^2,
f_1_2=2*pVM*(1-pVM),
f_2_2=1-f_0_2-f_1_2,
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
g_0_2=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^2,
g_1_2=2*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r)^2,
g_2_2=1-g_0_2-g_1_2,
g_0_4=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^4,
g_1_4=4*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))^3/(1-(1-pm)*r)^4,
g_2_4=1-g_0_4-g_1_4,
p22=p*f_2_4*f_2_2 + (1-p)*pd*g_2_4*g_2_2,
p21=p*f_2_4*f_1_2 + (1-p)*pd*g_2_4*g_1_2,
p20=p*f_2_4*f_0_2 + (1-p)*pd*g_2_4*g_0_2,
p12=p*f_1_4*f_2_2 + (1-p)*pd*g_1_4*g_2_2,
p11=p*f_1_4*f_1_2 + (1-p)*pd*g_1_4*g_1_2,
p10=p*f_1_4*f_0_2 + (1-p)*pd*g_1_4*g_0_2,
p02=p*f_0_4*f_2_2 + (1-p)*pd*g_0_4*g_2_2,
p01=p*f_0_4*f_1_2 + (1-p)*pd*g_0_4*g_1_2,
p00=p*f_0_4*f_0_2 + (1-p)*pd*g_0_4*g_0_2+(1-p)*(1-pd);

VIRTUALIZED_4_4(p,pVM):
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
p22=p*f_2_4*f_2_4,
p21=p*f_2_4*f_1_4,
p20=p*f_2_4*f_0_4,
p12=p21,
p11=p*f_1_4*f_1_4,
p10=p*f_1_4*f_0_4,
p02=p20,
p01=p10,
p00=p*f_0_4*f_0_4+(1-p);

VIRTUALIZED_8_4(p,pVM):
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
f_0_8=(1-pVM)^8,
f_1_8=8*pVM*(1-pVM)^7,
f_2_8=1-f_0_8-f_1_8,
p22=p*f_2_8*f_2_4,
p21=p*f_2_8*f_1_4,
p20=p*f_2_8*f_0_4,
p12=p*f_1_8*f_2_4,
p11=p*f_1_8*f_1_4,
p10=p*f_1_8*f_0_4,
p02=p*f_0_8*f_2_4,
p01=p*f_0_8*f_1_4,
p00=p*f_0_8*f_0_4+(1-p);

VIRTUALIZED-M_4_4(p,pd,pm,r,pVM):
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
g_0_4=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^4,
g_1_4=4*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))^3/(1-(1-pm)*r)^4,
g_2_4=1-g_0_4-g_1_4,
p22=p*f_2_4*f_2_4 + (1-p)*pd*g_2_4*g_2_4,
p21=p*f_2_4*f_1_4 + (1-p)*pd*g_2_4*g_1_4,
p20=p*f_2_4*f_0_4 + (1-p)*pd*g_2_4*g_0_4,
p12=p21,
p11=p*f_1_4*f_1_4 + (1-p)*pd*g_1_4*g_1_4,
p10=p*f_1_4*f_0_4 + (1-p)*pd*g_1_4*g_0_4,
p02=p20,
p01=p10,
p00=p*f_0_4*f_0_4 + (1-p)*pd*g_0_4*g_0_4 + (1-p)*(1-pd);

VIRTUALIZED-M_8_4(p,pd,pm,r,pVM):
f_0_4=(1-pVM)^4,
f_1_4=4*pVM*(1-pVM)^3,
f_2_4=1-f_0_4-f_1_4,
f_0_8=(1-pVM)^8,
f_1_8=8*pVM*(1-pVM)^7,
f_2_8=1-f_0_8-f_1_8,
g_0_4=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^4,
g_1_4=4*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))^3/(1-(1-pm)*r)^4,
g_2_4=1-g_0_4-g_1_4,
g_0_8=(((1-pm)*(1-r)+pm*(1-pVM))/(1-(1-pm)*r))^8,
g_1_8=4*pm*pVM*((1-pm)*(1-r)+pm*(1-pVM))^7/(1-(1-pm)*r)^8,
g_2_8=1-g_0_8-g_1_8,
p22=p*f_2_8*f_2_4 + (1-p)*pd*g_2_8*g_2_4,
p21=p*f_2_8*f_1_4 + (1-p)*pd*g_2_8*g_1_4,
p20=p*f_2_8*f_0_4 + (1-p)*pd*g_2_8*g_0_4,
p12=p*f_1_8*f_2_4 + (1-p)*pd*g_1_8*g_2_4,
p11=p*f_1_8*f_1_4 + (1-p)*pd*g_1_8*g_1_4,
p10=p*f_1_8*f_0_4 + (1-p)*pd*g_1_8*g_0_4,
p02=p*f_0_8*f_2_4 + (1-p)*pd*g_0_8*g_2_4,
p01=p*f_0_8*f_1_4 + (1-p)*pd*g_0_8*g_1_4,
p00=p*f_0_8*f_0_4 + (1-p)*pd*g_0_8*g_0_4 + (1-p)*(1-pd);