Let $U_1$ and $U_2$ be $n$-qubit unitaries, and denote by $P_{U_1U_2}(y \mid x) = |\langle y | U_1U_2 | x \rangle|^2$ the probability of measuring $y \in \{0,1\}^n$ on input $x \in \{0,1\}^n$. Suppose $P_{U_1U_2}(y \mid x)$ is efficiently sampleable in the sense that there exists a classical probabilistic polynomial time algorithm $A$ such that on input $x$, $A$ outputs $y$ with probability $P_{U_1U_2}(y \mid x)$. Is it the case that $P_{(U_1U_2)^\dagger}(y \mid x) = |\langle y | (U_1U_2)^\dagger | x \rangle|^2 = |\langle y | U_2^\dagger U_1^\dagger | x \rangle|^2$ is also efficiently sampleable?
I feel like this should be the case but I cannot prove it. Perhaps there is a simple reason why or a paper on the topic.