m e n non possono essere entrambi pari. Se così fosse, $ \displaystyle~4\mid mn\mid 3^m+1 $, ma $ \displaystyle~3^m+1\equiv (-1)^m+1\equiv 2\pmod 4 $. Dunque posto $ \displaystyle~x=(m,n) $ abbiamo x dispari.
Dall'ipotesi $ \displaystyle~mn\mid 3^m+1\mid 9^m-1 $ e $ \displaystyle~mn\mid 3^n+1\mid 9^n-1 $. Toh, proprio 10 mesi fa era capitato sul forum
un lemma che ci permette di dire che $ \displaystyle~mn\mid 9^x-1 $, da cui otteniamo l'ipotesi più debole $ \displaystyle~x^2\mid 9^x-1 $.
Supponiamo $ \displaystyle~x>1 $ e sia $ \displaystyle~x=p_1^{e_1}\cdots p_u^{e_u} $ la fattorizzazione di x con $ \displaystyle~p_i<p_j~\forall i<j $.
Deve valere $ \displaystyle~p_1^{2e_1}\mid 9^{p_1^{e_1}\cdots p_u^{e_u}}-1~(*) $, in particolare $ \displaystyle~p_1\mid 9^{p_1^{e_1}\cdots p_u^{e_u}}-1 $ o applicando ripetutamente il piccolo teorema di Fermà $ \displaystyle~p_1\mid 9^{p_2^{e_2}\cdots p_u^{e_u}}-1 $
Dunque possiamo applicare
quest'altro lemma che insieme alla $ \displaystyle~(*) $ dà
$ \displaystyle~\upsilon_{p_1}((9^{p_2^{e_2}\cdots p_u^{e_u}})^{p_1^{e_1}}-1)=\upsilon_{p_1}(9^{p_2^{e_2}\cdots p_u^{e_u}}-1)+\upsilon_{p_1}(p_1^{e_1})\ge 2e_1 $ cioè $ \displaystyle~\upsilon_{p_1}(9^{p_2^{e_2}\cdots p_u^{e_u}}-1)\ge e_1 $.
Da $ \displaystyle~9^{p_2^{e_2}\cdots p_u^{e_u}}\equiv 1\pmod {p_1^{e_1}} $ abbiamo $ \displaystyle~ord_{p_1^{e_1}}(9)\mid (p_2^{e_2}\cdots p_u^{e_u},p_1^{e_1-1}(p_1-1))=(p_2^{e_2}\cdots p_u^{e_u},p_1-1)=1 $, infatti $ \displaystyle~p_1-1<p_i~\forall i $ quindi $ \displaystyle~p_1-1 $ non è divisibile per nessun $ \displaystyle~p_i $. Per cui $ \displaystyle~p_1^{e_1}\mid 9-1=8 $, assurdo.
Segue che $ \displaystyle~x=1 $, cioè $ \displaystyle~mn\mid 9-1=8 $, quindi (1,1), (1,2) e (2,1) sono le uniche soluzioni.