docs: def max_entropy_distribution(node_indices, number_of_nodes):
      """Return the maximum entropy distribution over a set of nodes.
  
      This is different from the network's uniform distribution because nodes
      outside ``node_indices`` are fixed and treated as if they have only 1
      state.
  
      Args:
          node_indices (tuple[int]): The set of node indices over which to take
              the distribution.
          number_of_nodes (int): The total number of nodes in the network.
  
      Returns:
          np.ndarray: The maximum entropy distribution over the set of nodes.
      """
      distribution = np.ones(repertoire_shape(node_indices, number_of_nodes))
  
      return distribution / distribution.size