  For a $\psi$-mixing process $\xi_0,\xi_1,\xi_2,...$ we consider the number
$\mathcal{N}_N$ of multiple returns $\{\xi_{q_{i,N}(n)}\in\Gamma_N,\,
i=1,...,\ell\}$ to a set $\Gamma_N$ for $n$ until either a fixed number $N$ or
until the moment $\tau_N$ when another multiple return
$\{\xi_{q_{i,N}(n)}\in\Delta_N,\, i=1,...,\ell\}$ takes place for the first
time where $\Gamma_N\cap\Delta_N=\emptyset$ and $q_{i,N},\, i=1,...,\ell$ are
certain functions of $n$ taking on nonnegative integer values when $n$ runs
from 0 to $N$. The dependence of $q_{i,N}(n)$'s on both $n$ and $N$ is the main
novelty of the paper. Under some restrictions on the functions $q_{i,N}$ we
obtain Poisson distributions limits of $\mathcal{N}_N$ when counting is until
$N$ as $N\to\infty$ and geometric distributions limits when counting is until
$\tau_N$ as $N\to\infty$. We obtain also similar results in the dynamical
systems setup considering a $\psi$-mixing shift $T$ on a sequence space
$\Omega$ and studying the number of multiple returns $\{
T^{q_{i,N}(n)}\omega\in A^a_n,\, i=1,...,\ell\}$ until the first occurrence of
another multiple return $\{ T^{q_{i,N}(n)}\omega\in A^b_m,\, i=1,...,\ell\}$
where $A^a_n,\, A_m^b$ are cylinder sets of length $n$ and $m$ constructed by
sequences $a,b\in\Omega$, respectively, and chosen so that their probabilities
have the same order.
