Solution. Denote bn=rad(an)an. Since rad(an) divides rad(an+1) we have bn+1∣bn+1. If there are indices i<j with bi<2022<bi+1, we will be done by “continuity”. If, to the contrary, this does not happen, there are two possible cases.
* bn<2022 for all n big enough. Since an increases indefinitely, then so does rad(an), so at some moment rad(an) receives a new prime p>2022. This means that p∤an and p∣an+1=bnbn+1an, so p∣bn+1 and hence bn≥2022, a contradiction.
* bn>2022 for all n. We can assume WLOG that b0 is the smallest term of the sequence (bn). Suppose that bi+1=bi+1 for all 0≤i<n. Then
rad(a0)=⋯=rad(an−1)=R.
But for every prime p≤n there is a multiple of p among us sus b0,…,bn−1, so p∣ak for some k and consequently p∣R. Since not every prime divides R, there must be an index n such that bn<bn−1+1, i.e. bn−1+1=dbn for some d>1 and rad(an+1)=dR, so gcd(d,R)=1.
Recall that b0≤bn=db0+n, which reduces to n≥(d−1)b0. By above, this means that all primes up to (d−1)b0 divide R, but d does not divide R, so d>(d−1)b0, which is impossible. □