bounded products of some primes

[G. M.-S.]

What is the "best" way of def'ing a function for this in SageMath?

Input:  a list of positive integers [x1, … , xn] (actually distinct primes) and a bound M.

Output:  the list of all (e1, … , en) with ei ⩾ 0 such that x1^e1*…*xn^en ⩽ M.

TIA,

Guillermo