%数学建模案例选集
%游戏网站奖项设置
clc
clear
format long
n = 4;p = 1/n;q = 1-p;
a = 40;b = 100;
so = 5000;
for i = 1:23
    s(i) = so + 100*i;
    pmax = 0;
    for M = 1:1000
        ko = ceil((s(i)-so+M*a)/(a+b));
        PM = 1-binocdf(ko-1,M,p);
        if PM > pmax
            pmax = PM;
            Mmax = M;
        end
    end
    pstar(i) = pmax;
    mstar(i) = Mmax;
end
ans = [s',mstar',pstar',ceil(10000*pstar)']