Another update:
As explained in the paper below, you can use rook polynomials to solve such problems.
Playing with a full deck of 52 cards we will call "one" 18 times,
we will call "two" 17 times, and we will call "three" 17 times.
The forbidden positions in the 52 by 52 board consist of three "independent"
complete rectangles; one $18\times 4$ and the other two $17\times 4$.
The rook polynomial for a full $m\times n$ rectangle with $m\geq n$ is
$$\sum_{k=0}^n{m\choose k}\, {n!\over (n-k)!}\, x^k. $$
Multiply the polynomials for these three rectangles to give us the rook polynomial for our problem
$$R(x)=(1+73440x^4+19584x^3+1836x^2+72x)(1+57120x^4+16320x^3+1632x^2+68x)^2.$$
The number of winning deck orders is
$$\int_0^\infty x^N R(-1/x) \exp(-x)\,dx $$
so the probability is this divided by $N!$, i.e.
$$\mathbb{P}(\text{win})= 24532967512/3004641364725= 0.008165023553.$$
Reference: F. F. Knudsen and I. Skau,
On the Asymptotic Solution of a Card-Matching Problem,
Mathematics Magazine 69 (1996), 190-197.
Update: The solution below is for a simplified version of the problem where you work
with a deck of size 12: four each of ace, deuce, and trey.
This is a problem in generalized derangements and joriki's answer here
tells you what to do. In general, the number of deck orders that lead to a win is
$$\int_0^\infty L_{n_1}(x)\cdots L_{n_r}(x)\,\mathrm e^{-x}\mathrm dx.$$
In this problem, we have $r=3$ and $n_1=n_2=n_3=4$. The fourth Laguerre
polynomial is $L_4(x)=(x^4-16x^3+72x^2-96x+24)/24$. Raising this to the third
power and integrating against $\exp(-x)$ gives $346$.
That is, there are $346$ ways to order the deck that give a win.
Divide this by the total number of orders $12!/(4!)^3$, to give
$$\mathbb{P}(\text{win})=173/17325=0.00998.$$