二章
书迷正在阅读:俺不是耕田的 , 某柯学的黑巫守则 , (柯南)今天大小姐又不准我离开GL , 怪物们的朋友[nph人外] , 【SVT】续约条件是集体住宿! , 我的100分男孩 , 转生成猪的我,突破只能靠双修 , 平行世界的管理者突然当起了腐女!? , 岁岁好 , 五月二十九 , 小坏蛋,再来一次好不好? , 奔向你的星眸
ndom importstatistics defsimute_spreadpop_size=1000,initial_seed=5,i_prob=0.03,resist_prob=0.1,steps=50: #状态:0=未接触,1=接触但未采纳,2=采纳「感染」 popution=[0]*pop_size #随机设置初始采纳者seed foriinrandom.samplerangepop_size,initial_seed: popution[=2 history=[] forstepinraeps: new_popution=popution.copy foriinrangepop_size: ifpopution[==2: #试着影响k随机接触者 tacts=random.samplerangepop_size,k=5 fortacts: ifpopution[c]==0: ifrandom.random #接触後有机会直接采纳 new_popution[c]=2 1 else: #有机会是只是「接触但未采纳」 new_popution[c]=1 elifpopution[c]==1: #已接触过但犹豫的人,再次接触有机