在我进行matlab代码调整的时候，发现一个有趣的现象。有matlab的同学可以拿以下代码实验：
for b=[1:1:80];
a=b*2-1;r=a-b;
c=0;p=1;
for i=[1:1:b];
p=p*(b-i+1)/(a-i+1);
p1=p*r/(a-i);
q=i*p1;
c=c+q;
end
cb(b)=c;
end
cb
这串代码描述的牌堆剩余总牌数正好是黑牌的2倍-1，也就是剩余牌堆红牌仅仅比黑牌少1张的情况。
我们输出的cb是一个长度为80的向量，里面的每个分量都是1.0000，于是我们找到了洛神正好为1的情况：牌堆剩余黑牌数量恰好比红牌多1的时候，洛神期望为1.当牌堆黑牌比红牌多2或2以上的时候，期望大于1.当黑牌等于红牌或者少于红牌的时候，期望小于1.但是期望正好等于1的时候，需要严格数学公式推导，有兴趣的同学请帮忙推导一下，也就是a=2b-1的时候，该期望=1，楼主看这个公式非常蛋疼，就不想推了。
【完】