def METHOD_NAME(alpha: Tensor) -> Tensor:
    F = getF(alpha)
    samples_gamma = F.sample_gamma(
        alpha=alpha, beta=F.ones_like(alpha), dtype=dtype
    )
    sum_gamma = F.sum(samples_gamma, axis=-1, keepdims=True)
    samples_s = F.broadcast_div(samples_gamma, sum_gamma)
    return samples_s