Although computationally slow, this might do the trick:

import numpy as np

def draw_random_samples(len_deck, n_simulations, n_cards):
    
    values = np.random.random((n_simulations, len_deck))
    indices = np.argsort(values, axis=-1)
    indices = indices[...,:n_cards]
    
    return indices