The answer is all pairs (1,1) and (P,P+i),(P,P−i), where P is a non-constant monic polynomial in C[x] and i is the imaginary unit. Notice that if P∣Q2+1 and Q∣P2+1 then P and Q are coprime and the condition is equivalent with PQ∣P2+Q2+1. Lemma. If P,Q∈C[x] are monic polynomials such that P2+Q2+1 is divisible by PQ, then degP=degQ. Proof. Assume for the sake of contradiction that there is a pair (P,Q) with degP=degQ. Among all these pairs, take the one with smallest sum degP+degQ and let (P,Q) be such pair. Without loss of generality, suppose that degP>degQ. Let S be the polynomial such that PQP2+Q2+1=S. Notice that P a solution of the polynomial equation X2−QSX+Q2+1=0, in variable X. By Vieta's formulas, the other solution is R=QS−P=PQ2+1. By R=QS−P, the R is indeed a polynomial, and because P,Q are monic, R=PQ2+1 is also monic. Therefore the pair (R,Q) satisfies the conditions of the Lemma. Notice that degR=2degQ−degP<degP, which contradicts the minimality of degP+degQ. This contradiction establishes the Lemma. By the Lemma, we have that deg(PQ)=deg(P2+Q2+1) and therefore PQP2+Q2+1 is a constant polynomial. If P and Q are constant polynomials, we have P=Q=1. Assuming that degP=degQ≥1, as P and Q are monic, the leading coefficient of P2+Q2+1 is 2 and the leading coefficient of PQ is 1 , which give us PQP2+Q2+1=2. Finally we have that P2+Q2+1=2PQ and therefore (P−Q)2=−1, i.e Q=P+i or Q=P−i. It's easy to check that these pairs are indeed solutions of the problem.