Let c=100!. Suppose that n≥m+2. Then am+n=a(m+1)+(n−1) divides both c(am+am+1+⋯+an−1+an) and c(am+1+⋯+an−1), so it also divides the difference c(am+an). Notice that if n=m+1 then am+n divides c(am+am+1)=c(am+an), and if n=m then am+n divides both cam and 2cam=c(am+an). In either cases; am+n divides c(am+an).
Analogously, one can prove that if m>n,am−n=a(m−1)−(n−1) divides c(am−an), as it divides both c(an+1+⋯+am) and c(an+⋯+am−1).
From now on, drop the original divisibility statement and keep the statements " am+n divides c(am+an) " and " am−n divides c(am−an)." Now, all conditions are linear, and we can suppose without loss of generality that there is no integer D>1 that divides every term of the sequence; if there is such an integer D, divide all terms by D.
Having this in mind, notice that am=am+n−n divides c(am+n−an) and also c(am+n−am−an); analogously, an also divides c(am+n−am−an), and since am+n divides c(am+an), it also divides c(am+n−am−an). Therefore, c(am+n−am−an) is divisible by am,an, and am+n, and therefore also by lcm(am,an,am+n). In particular, cam+n≡c(am+an)(modlcm(am,an)).
Since am+n divides c(am+an),am+n≤c(am+an).
From now on, we divide the problem in two cases.
Case 1: there exist m,n such that lcm(am,an)>c2(am+an).
If lcm(am,an)>c2(am+an)>c(am+an) then both c(am+an) and cam+n are less than cam+n≤c2(am+an)<lcm(am,an). This implies cam+n=c(am+an)⟺am+n=am+an. Now we can extend this further: since gcd(am,am+n)=gcd(am,am+an)=gcd(am,an), it follows that
lcm(am,am+n)=gcd(am,an)amam+n=anam+nlcm(am,an)>anc2(am+an)2>anc2(2aman+an2)=c2(2am+an)=c2(am+am+n).
We can iterate this reasoning to obtain that akm+n=kam+an, for all k∈Z>0. In fact, if the condition lcm(am,an)>c2(am+an) holds for the pair (n,m), then it also holds for the pairs (m+n,m),(2m+n,m),…,((k−1)m+n,m), which implies akm+n=a(k−1)m+n+am=a(k−2)m+n+2am=⋯=an+kam.
Similarly, am+kn=am+kan.
Now, am+n+mn=an+(n+1)m=am+(m+1)n⟹an+(n+1)am=am+(m+1)an⟺man=nam. If d=gcd(m,n) then dnam=dman.
Therefore, since gcd(dm,dn)=1,dm divides am and dn divides an, which means that
an=dm⋅t=dtm and am=dn⋅t=dtn, for some t∈Z>0
which also implies
akm+n=dt(km+n) and am+kn=dt(m+kn), for all k∈Z>0.
Now let's prove that akd=tk=dt(kd) for all k∈Z>0. In fact, there exist arbitrarily large positive integers R,S such that kd=Rm−Sn=(n+(R+1)m)−(m+(S+1)n) (for instance, let u,v∈Z such that kd=mu−nv and take R=u+Qn and S=v+Qm for Q sufficiently large.)
Let x=n+(R+1)m and y=m+(S+1)n. Then kd=x−y⟺x=y+kd,ax=dtx, and ay=dty=dt(x−kd)=ax−tk. Thus ax divides c(ay+akd)=c(akd+ax−tk), and therefore also c(akd−tk). Since ax=dtx can be arbitrarily large, akd=tk=dt(kd). In particular, ad=t, so akd=kad.
Since kad=akd=a1+(kd−1) divides c(a1+akd−1),bk=akd−1 is unbounded. Pick p>akd−1 a large prime and consider apd=pad. Then
lcm(apd,akd−1)≥lcm(p,akd−1)=pakd−1
We can pick akd−1 and p large enough so that their product is larger than a particular linear combination of them, that is,
lcm(p,akd−1)=pakd−1>c2(pad+akd−1)=c2(apd+akd−1)
Then all the previous facts can be applied, and since gcd(pd,kd−1)=1,ak=akgcd(pd,kd−1)=kagcd(pd,kd−1)=ka1, that is, the sequence is linear. Also, since a1 divides all terms, a1=1.
Case 2: lcm(am,an)≤c2(am+an) for all m,n.
Suppose that am≤an; then lcm(am,an)=Man≤c2(am+an)≤2c2an⟹M≤2c2, that is, the factor in the smaller term that is not in the larger term is at most 2c2.
We prove that in this case the sequence must be bounded. Suppose on the contrary; then there is a term am that is divisible by a large prime power pd. Then every larger term an is divisible by a factor larger than 2c2pd. So we pick pd>(2c2)2, so that every large term an is divisible by the prime power pe>2c2. Finally, fix ak for any k. It follows from ak+n∣c(ak+an) that (an) is unbounded, so we can pick an and ak+n large enough such that both are divisible by pe. Hence any ak is divisible by p, which is a contradiction to the fact that there is no D>1 that divides every term in the sequence.