(1) x↦x2+1(mod5) では 1⟼2⟼0⟼1 と巡回する。a1≡1(mod5) だから、a3m+1,a3m+2,a3m+3 の剰余はそれぞれ 1,2,0 となる(m≧0)。従って 3∣n なら 5∣an である。
(2) 補助的に a0=0 とおくと、漸化式は a1=a02+1 でも成り立つ。任意の正整数 k について ak+r≡ar(modak)(r≧0)(1) が成り立つ。実際、r=0 では ak≡a0=0(modak) であり、ak+r≡ar(modak) なら漸化式から ak+r+1=ak+r2+1≡ar2+1=ar+1(modak) となるので、数学的帰納法で示される。
n=qk+r (0≦r<k) とする。n≧k なら (1) を繰り返して an≡ar(modak) となる。数列は aj+1−aj=aj2−aj+1>0 より増加するので、1≦r<k なら 0<ar<ak であり、この場合 ak∤an。一方、r=0 なら an≡a0=0(modak)。また n<k の場合も 0<an<ak で割り切れない。従って ak∣anすなわちr=0すなわちk∣n.
(3) (1) より、m≦n のとき gcd(am,an)=gcd(am,an−m) である。添字にユークリッドの互除法を繰り返し適用すると gcd(am,an)=agcd(m,n) を得る。ここで 8091=4⋅2022+3、2022=674⋅3 だから gcd(a2022,a8091)=a3=5. また a1,a2,a3≡1,2,5(mod25) であり、52+1≡1(mod25) だから、剰余は 1,2,5 と3項周期で繰り返す。3∣2022 より a2022≡5(mod25)、従って 25∤a2022。さらに 3∣8091 なので (1) より 5∣a8091、すなわち 25∣(a8091)2。
gcd(a2022,(a8091)2) の任意の素因数は a2022 と a8091 の双方を割るため、gcd(a2022,a8091)=5 から素因数は 5 に限られる。25∤a2022 であり、両数は5で割り切れるから gcd(a2022,(a8091)2)=5.