Mathematica有三个函数可以干类似的事，根据同一个元素能否重复取，以及列表是否有序分类来看的话：
Tuples 可重复取 有顺序
Permutations 不可重复取 有顺序
Subsets 不可重复取 无顺序
那这样来看的话，好像缺少一个函数，可重复取，无顺序
例如对于 f[{1, 2, 3}, {3}] 应该给出
{{1, 1, 1}, {1, 1, 2}, {1, 1, 3}, {1, 2, 2}, {1, 2, 3}, {1, 3, 3}, {2, 2, 2}, {2, 2, 3}, {2, 3, 3}, {3, 3, 3}}
我能想到的办法是
Union@Subsets[{1, 1, 1, 2, 2, 2, 3, 3, 3}, {3}]
或者
Union[Sort /@ Tuples[{1, 2, 3}, {3}]]
但是这两个方法对于比较大的列表，中间结果就会很占内存，有什么更好的办法吗？

SE帖子《How to obtain all subsets of a list with repetition?》（编号 258495 ）似乎说的就是这个问题。注意 thorimur 的评论。

除了2楼的方法又想到一个，由于计算顺序的问题会有一个警告，暂时没想好怎么规避
f[l_List,n_Integer]:=
Module[{i},
i[0]=1;
Flatten[Table@@Prepend[{i[
#],i[#
-1],Length@l}&/@Range@n,l[[i/@Range@n]]],n-1]]

关于4楼的计算次序控制，我能想到最简单的是，Block：
f[l_List, n_Integer] := Module[{i}, i[0] = 1;
Flatten[
Block[{Part, Table},
Table @@Prepend[{i[
#], i[#
- 1], Length@l} & /@ Range@n, l[[i /@ Range@n]]]],
n - 1]]

另一种调法：
f[l_List, n_Integer] := Module[{i}, i[0] = 1;
Flatten[
Table[l[[#]], ##2] & @@
Prepend[Function[index, {i[index], i[index - 1], Length@l}] /@
Range@n, i /@ Range@n], n - 1]]

GroupTheory`Tools`Multisets[Range[3],3]
MSE：fastest-way-to-generate-a-a-a-b-a-c-b-b-b-c-c-c-from-a-b-c

Clear[f]
f[n_,k_]:=SymmetrizedIndependentComponents[Table[n,{k}],Symmetric@Range@k];
f[3,3]