Dice roll probability sim

roll.py

import random

print(
    """
This script simulates this scenario:
you have a die of N faces
you have Y quantity of that die
what is the average of X rolls,
taking the highest value of all
dies for that roll?
"""
)

sides = int(input("Enter the number of die faces: "))
advantage = int(input("Enter the number of dice to roll"))
rolls = int(input("Enter the number of simulations: "))

sum_ = 0

for i in range(rolls):
    sum_ += max([random.randrange(0, sides) for x in range(advantage)])
res = sum_ / rolls

print(f"The expected value is: {res}")

Inspired by Matt Parker.

made in 5 minutes.