Find all integer coefficient polynomials $Q$ such that 
[list]
[*] $Q(n)\ge 1$   $\forall n\in \mathbb{Z}_+$.
[*] $Q(mn)$ and $Q(m)Q(n)$ have the same number of prime divisors $\forall m,n\in\mathbb{Z}_+$.
[/list]
1. **Restate the problem and notation:**
   We need to find all integer coefficient polynomials \( Q \) such that:
   - \( Q(n) \geq 1 \) for all \( n \in \mathbb{Z}_+ \).
   - \( Q(mn) \) and \( Q(m)Q(n) \) have the same number of prime divisors for all \( m, n \in \mathbb{Z}_+ \).

   Let \( f(x) = Q(x) \). We aim to show that the only solutions are of the form \( f(x) = cx^d \) for \( c \in \mathbb{Z}^+ \) and \( d \in \mathbb{Z}^+ \cup \{0\} \).

2. **Introduce the function \(\omega(n)\):**
   Define \(\omega(n)\) as the number of distinct prime divisors of \( n \). We need to show that \( \omega(f(mn)) = \omega(f(m)f(n)) \).

3. **Degree of the polynomial \( f \):**
   Let \( d = \deg(f) \). We will use the average order of \(\omega(f(n))\) to derive a contradiction if \( f \) is not of the form \( cx^d \).

4. **Main Claim:**
   Suppose \( f \) is nonconstant and satisfies \( f(n) \geq 1 \) for all positive integers \( n \) and \( \omega(f(ab)) \leq \omega(f(a)) + \omega(f(b)) + C \) for some constant \( C \). Then \( f(0) = 0 \).

5. **Proof of the Main Claim:**
   - Assume \( f(x) = x^r g(x) \) where \( g(0) \neq 0 \).
   - Consider the number of prime divisors:
     \[
     \omega(f(ab)) = \omega(ab \cdot g(ab)) \geq \omega(ab) + \omega(g(ab)) - \omega(\gcd(g(ab), ab))
     \]
     \[
     \omega(f(a)f(b)) = \omega(ab \cdot g(a)g(b)) \leq \omega(ab) + \omega(g(a)g(b)) \leq \omega(ab) + \omega(g(a)) + \omega(g(b))
     \]
   - If \( g \) is nonconstant, then by the claim (taking \( C = \omega(f(0)) \)), \( g(0) = 0 \), which is a contradiction. Hence, \( g \) must be constant.

6. **Lemma:**
   \( p_k < e^{O(k)} \) for all \( k \).

   - Mimic Erdős' proof of the infinitude of primes.
   - For each \( i \), write \( f(i) = a \cdot b^{d+1} \) where \( 0 \leq \nu_p(a) \leq d \) for any prime \( p \).
   - Considering all \( i \) such that \( f(i) < p_k \):
     - There are at most \( (p_k)^{\frac{1}{d+1}} \) choices of \( b \).
     - There are at most \( (d+1)^k \) choices of \( a \).
   - Thus,
     \[
     (p_k)^{\frac{1}{d+1}} \cdot (d+1)^k \gg (p_k)^{\frac{1}{d}}
     \]
     which gives the desired result after rearranging.

7. **Bad and Good Primes:**
   - Call a prime \( p \) *bad* if it is less than \( 2020d \), divides \( f(0) \), or the leading coefficient \( c \) of \( f \); otherwise, call it *good*.
   - Take \( N \) very large so that \( p_N \) exceeds all bad primes.
   - For each \( i \) where \( p_i \) is good, pick residues \( a_i, b_i \) such that:
     \[
     p \mid f(a_i b_i) \text{ and } p \nmid f(a_i) f(b_i)
     \]
   - Use the Chinese Remainder Theorem (CRT) to pick \( a, b < p_1 p_2 \ldots p_N \) such that:
     \[
     a \equiv a_i \pmod{p_i},\ b \equiv b_i \pmod{p_i} \text{ for all } i \leq N
     \]

8. **Summing the Inequality:**
   - For any constant \( k \), we have:
     \[
     N \leq \omega(f(a + kP)) + \omega(f(b + kP)) + C \quad \text{where } P = p_1 p_2 \ldots p_N
     \]
   - Summing from \( k = 1 \) to \( M - 1 \) (pick \( M = e^{N^3} \)):
     \[
     M(N - C) \leq \sum_{k = 1}^M \omega(f(a + kP)) + \omega(f(b + kP))
     \]

9. **Analyzing the Sum:**
   - If \( p \leq p_N \) and \( p \) is bad, then \( p \) may divide \( f(a + kP) \), but there are finitely many such primes (at most \( t \)).
   - If \( p \leq p_N \) and \( p \) is good, then \( p \) cannot divide \( f(a + kP) \) due to our construction.
   - If \( p_N < p < \sqrt{MP} \), then the leading coefficient of \( g(x) = f(a + Px) \) is not zero modulo \( p \). Thus, at most:
     \[
     d \left\lceil \frac{M}{p} \right\rceil < \frac{dM}{p} + d
     \]
     numbers are divisible by \( p \).
   - If \( p > \sqrt{MP} \), then most terms can have at most \( 2d \) such divisors due to size reasons. Thus, these primes contribute at most \( 2dM + O_{M, N}(1) \) to the sum.

10. **Final Summation:**
    \[
    M(N - C) \leq 2 \left( t + \sum_{p < \sqrt{MP}} \left( \frac{dM}{p} + d \right) + dM + O_{M, N}(1) \right)
    \]
    \[
    \leq 2 \left( (M \log \log (MP) + A + o_{M, N}(1) + d \sqrt{MP}) + 2dM + O_{M, N}(1) \right)
    \]
    \[
    \leq 2 \left( M \log \log (e^{N^5}) + O(e^{N^3/2} + O(N^2)) \right)
    \]
    \[
    \leq 2 \left( 5M \log N + O(e^{N^3/2 + O(N^2)}) \right)
    \]
    which is a contradiction if \( N \) is large enough.

Thus, the only possible form for \( f(x) \) is \( cx^d \).

The final answer is \( \boxed{ f(x) = cx^d } \) for \( c \in \mathbb{Z}^+ \) and \( d \in \mathbb{Z}^+ \cup \{0\} \).