有关于毕导THU尴尬定理的代码实现

此篇文章灵感来自于毕导29号更新的视频，感兴趣的童鞋请戳这里：https://www.bilibili.com/video/BV1RT4y1j7pP

在视频当中，毕导向我们阐释了尴尬定理的含义，以及解决尴尬定理通式的传统计算机解决问题的分治方法，在这里就直接摆出递推式了

在这里，我用MATLAB加以实现，代码如下所示

m=50;%手动输入k的取值范围，假定是1至50个人上厕所
f=zeros(1,m);%初始化占位函数，全赋值为1
a=ones(1,m);%初始化坑位利用率函数，全赋值为1
f(1)=1;%1人入厕所，自然是一个坑
f(2)=1;
f(3)=2;
for i=[4:1:m]%输入递归函数
    if mod(i,2)~=0
        f(i)=2*f((i+1)/2)-1;
    else
        f(i)=f(i/2)+f((i/2)+1)-1;
    end
end

for i=[1:1:m]%计算坑位利用率函数
    a(i)=f(i)/i;
end

figure(1)

x=[1:1:m];
y=f(1:1:m);
plot(x,y);
xlabel('坑位数量')
ylabel('最多容纳人数')

figure(2)
x=[1:1:m];
y=a;
plot(x,y);
xlabel('坑位数量')
ylabel('坑位利用率')

 

我们来看一下效果如何：

结果是特别amazing