Elimination Odds Calculator for Algicosathlon, Etc.

import random
sims=1000000 # Number of simulations to run
num=21 # Number of contestants remaining
toadd=[0,1,2,3,4,6,8,10,12,15,20,25,30,35,40,45,50,60,70,80,100,125,150,175,200,250,300,400,500,600,750,1000] # Scores for next event
multFactor=2**12 # What to multiply the scores by
loseCount=list()
for i in range(num):
    loseCount.append(0)
for i in range(sims):
    swap=list()
    taken=list()
    scores=[362316,358504,337394,300163,288051,220352,177754,169334,148832,148667,135386,135052,128540,128060,123618,122407,105775,105671,97714,97235,87663] # Contestant scores
    for k in range(num):
        swap.append(-1)
        taken.append(False)
    for k in range(num):
        while swap[k]==-1 or taken[swap[k]]:
            swap[k]=random.randint(0,num-1)
        taken[swap[k]]=True
        scores[k]+=round(toadd[swap[k]]*multFactor)
    recordLow=scores[0]*1000000
    recordHolder=0
    for k in range(num):
        if scores[k]<recordLow:
            recordLow=scores[k]
            recordHolder=k
    loseCount[recordHolder]+=1
print(loseCount)