——一个业余民科之谈
好久没有写些东西了，似乎是从研究生毕业之后，工作的这两年一直没有写。我觉得工作了让人更加实际，以前我会写一些诗，姑且认为是诗吧。工作之后，成了比996更惨的打工资深社畜，每天7点醒，12点睡，中间除了吃喝拉撒，就是工作了，实在没有闲心写些东西。后来我换了部门来到市场部，工作清闲了许多，但是当我一天心血来潮，想再写点什么，却发现憋了三个小时，只写了几句现代诗。
海洋中，一只鱼儿的游动带起了涟漪，
天空中，一只鸟儿的飞行搅动着气旋，
嘘，打破了宁静。森林里，干燥的木柴被意外地点燃，它想努力地燃烧，温暖它的朋友，却又担心热情，灼伤了友谊。
哎，笨拙的木柴。
我一直也没什么文采，做的只是用诗的形式传递我的想法，重在于思想。
幸运的是我获得了出国留学的机会，重新从社畜回到学生，有很多感慨，虽然将来还会再面临工作压力，但也能够在自由的国度中畅游一小段人生旅途。
这两年，我过的并不顺利，经历过失恋，经历过追求而不得，经历过疾病，所幸有亲人、朋友的关心和陪伴，我走过来了。所有的这些都让当时的我险些被击垮，但迈过去之后，苦难带给了我新的收获。首先是，我知道了自己不是超人，以前总觉得自己很幸运，所有的病痛、苦难都不会降临在我身上，然而我现在明白了，在不幸方面一切都是平等，这也让我开始给父母买了保险。其次是，我相信缘分了，我想起那些离我们而去的熟悉的面孔，虽然不多，但也让我觉得生命十分脆弱，人生充满不确定性。
以前我总觉得工作好像停不下来，但是我真的因为住院停了四个月，一切都照旧，井然有序，除了我自己在医院里，每天早晚在母亲的陪伴下散步，迎接太阳，又送走太阳。
当我发现离开我，即使是我自己从头策划的工作都可以正常继续下去，有一种虚无感上来了。好在这种念头随着我身体恢复，也逐渐消失，取而代之的是，我更加珍惜眼前的缘分，有时候我觉得自己好像一部移动的摄像机在不断地记录生活。
从概率上来讲每个人诞生的几率是非常渺茫，人生也许可以排出数千次的生殖液体，每次排出的体液里含有数以亿计的拥有一半染色体的细胞，其中只有一个能遇到另一半，组成受精卵，孕育出新生命。我们被不同的父母带到了世界各地，在有限的生命的前提下，能够遇见实属不易。这种不易，我从中学会了珍惜相遇、坦诚相待。
上面是我想说的概率论的第一层含义，从名字上来看概率论确实可能是计算概率，但我还想讲第二层含义，这一层和宿命更相关。
我曾经看过一个视频，带给我很大震撼，刚才我尝试去找，却没有找到了。我大概描述一下内容：科学家让一群人选择他自己喜欢的数字，选择好之后，统计前面n-1个人的数字，能够推导出最后一个人的数字满足一定的数学规律。很有意思的视频，如果哪一天我找到了再分享出来。
给我的震撼在于，我们以为自己在做绝对的随机选择，而且是由自由意志所支配的选择，却仍然能够在统计意义上被确定，就好似有一种神秘力量在支配我们的选择。其实从身边大家的感知也能知道，我们的随机没有那么的随机。前几天和朋友玩狼人杀，一个小伙伴替自己做好人辩解的台词是：“我上一把是狼人，这一把还是狼人的概率很低。”一听似乎是对的，但一想两次抽牌是独立事件，不会互相干扰，这是概率给独立事件的定义。但为什么我们都会有这样的潜意识认为两件独立的事件是有联系呢？
我的概率学的不好，也并非统计学出身，既然是闲谈，我说说自己的想法：概率和统计是两回事，概率上是独立的，统计上这些独立事件却仍需要满足一定的规律。就好像仍硬币，每次都是独立的1/2的概率正面朝上，但是仍了大量的硬币之后，做统计就是呈现了一个趋势。随着不断投掷更多的硬币，正面朝上的硬币数占总数的比例逐渐趋近于1/2，这是大数定律讲的内容。
因此我想讲的概率与命运，其实是概率里面的统计意义与命运。
这里我想自己做一个实验，这是我从统计热力学课堂上学到的。假如我们有400个人，每个人手里有5万块钱，可以通过石头剪刀布和其他人猜拳，赢的人拿走输的人一万块。在此过程，可以随机找任何人进行猜拳，唯一注意的是钱不能为负数，最少为0元。大家猜猜最终的结果是什么样？是1、所有人的钱仍保持不变？还是2、所有人钱是不均匀的？我猜测大家应该都会觉得是不均匀的，那么这种不均匀的是随机的还是满足一定的分布规律呢？这里做一个模拟实验，在电脑里给400个对象，他们初始有5万块钱，互相比赛，最后的结果满足一种统计规律——玻尔兹曼分布（Boltzmann Distribution）。如下图。
看看这个分布，以及我们做的假设，其实和真实社会很类似，社会的财富分布也呈现该分布特征——少部分拥有大量的财富，贫困人占大多数。
代码非常简单，提供给大家，有兴趣可以在python里面实现。
所以我想说的是，虽然我们看似在做自由选择，却仍然被笼罩在一个神秘力量的支配中，不自觉的参与到了大的社会实验中。
我想把这种神秘力量叫做命运。
import numpy as np
from matplotlib import pyplot as plt
# For the convenience of coding, I arrange the people in a 2-dimensional square.
# So the population in this numerical experiment is equal to row * col.
row = 20
col = 20
# Initializing the money of each people
bolt = np.ones((row,col))*5
# Repeating the process of exchanging money between two random people.
for i in np.arange(1000000):
x1 = np.random.randint(row)
y1 = np.random.randint(col)
x2 = np.random.randint(row)
y2 = np.random.randint(col)
# bolt[x1, y1] -= 1
# bolt[x2, y2] += 1
rand = np.random.rand()
# rand = 1
if bolt[x1,y1] - rand >= 0:
bolt[x1,y1] -= rand
bolt[x2,y2] += rand
else:
pass
# Making a plot to show the result.
n = 30
plt.hist(bolt.flatten(), bins = np.arange(0,n))
plt.title("Boltzmann distribution")
plt.xlim(0, n)
plt.xticks(np.linspace(0, n, num=n+1))
plt.show()
