Is there a known distribution for this permutation with replacement problem? Choose $t$ numbers from $n$ $(n>t)$ distinct numbers with replacement and the order of the $t$ numbers matters.
Say,
$P(X=1) = \dfrac{{numbers\ of\ unique\ t-set \ which\ has\ 1\ distinct\ number}}{n^t}$ (i.e. $n/n^t$),
$P(X=2)= \dfrac{{numbers\ of\ unique\ t-set \ which\ has\ 2\ distinct\ numbers}}{n^t}$,
$...$
$P(X=t) = \dfrac{{numbers\ of\ unique\ t-set \ which\ has\ t\ distinct\ numbers}}{n^t}$,
$P(X>t) = 0$.
For example, $n = 5$ and $t = 3$,
$P(X=1) = \dfrac{5}{125}$,
$P(X=2) = \dfrac{5 \times 4 \times 3}{125}= \dfrac{60}{125}$,
$P(X=3) = \dfrac{5 \times 4 \times 3}{125}= \dfrac{60}{125}$.
It seems to be very complicated for large $t$, is there a formula or even distribution function for that?
Unique means different taking orders into consideration. e.g. (1,1,3) = (1,1,3) but not equal to (1,3,1)
 A: Let $k$ be such that $1\leq k\leq t$. You want to count the number of $t$-tuples $a=(a_1,\dots,a_t)$, with, say, each $a_i\in\{1,2,\dots,n\}$, such that the unordered set $\{a_1,\dots,a_t\}$ has exactly $k$ elements. But your tuple $a$ can be thought as a function $a:\{1,\dots,t\}\to\{1,2,\dots,n\}$ (which sends each number $i$ to the number $a_i$), so what you want to count is the number of functions $f:\{1,\dots,t\}\to\{1,\dots,n\}$ such that the image of $f$ has exactly $k$ elements.
How to count these functions? first, choose the "target" $k$-subset of $\{1,\dots,n\}$, say $A$. This can be made in $\binom nk$ ways. Now for this set $A$ you want to count all the surjective functions from $\{1,\dots,t\}$ onto $A$, which in turn can be made in $k!\,\begin{Bmatrix}
t\\
k
\end{Bmatrix}\,$ ways, where $\begin{Bmatrix}
\vphantom{} \\
 \vphantom{}
\end{Bmatrix}$ denotes the Stirling numbers of the second kind; see here for an explanation. Therefore the total number of desired functions is 
$$k!\,\binom nk\begin{Bmatrix}
t\\
k
\end{Bmatrix}\,.$$
