MATLAB模拟一百万次(十万次方差有点大)
j=0;
for i=1:1000000
A=randperm(54);
a=A(1:13);
b=A(14:26);
c=A(27:40);
d=A(41:54);
e=[1,2];
a1=prod(ismember(e,a));
b1=prod(ismember(e,b));
c1=prod(ismember(e,c));
d1=prod(ismember(e,d));
if a1||b1||c1||d1
j=j+1;
end
end
fprintf('%d',j/i)
概率约为0.236,实际值是338/1431≈ 0.236198462613557