SG18.py

import scipy.special as special
from math import comb
K = 18   # number of sections
def survival_prob(r , K):
    probability = 2 ** (-K-1) * comb(K, r-1) * special.hyp2f1(1, 1+K, 2 + K - r, 1/2)
    return probability

def main():
    for r in range(16):
        if r <= K:
            p = survival_prob(r + 1, K)
        else:
            p = 1.0
        print(f'Probability of survival for player {r+1} is {p}!')        
    
if __name__ == "__main__":
    main()