Solution:
Observe that n=1 is peculiar and that every prime is peculiar. Consider dn, where d>1 is a divisor of n. Then dn(dn+1) divides n(n+1), equivalent to n+d dividing d2(n+1). Since n≡−d(modn+d), we obtain that n+d divides d3−d2, hence n+d≤d3−d2 and in particular n<d3. Now taking d to be the smallest prime divisor of n, we deduce that n cannot have more than two prime factors, counting multiplicities. Hence from now on we work with n=pq, where p and q are primes. In fact, p=q, as otherwise p+1∣p2−p, i.e. p+1∣2, contradiction.
Let us first prove that n is square-free. Assume that n=pkh, where gcd(p,h)=1. Now pk−1h(pk−1h+1)∣pkh(pkh+1), which is equivalent to (pk−1h+1)∣p(pkh+1). Since gcd(p,pk−1h+1)=1, we must have
(pk−1h+1)∣(pkh+1)=pkh+p−p+1
so (pk−1h+1)∣(p−1). If k≥2, this is a contradiction. Hence, n is square-free.
Let us now show that n has at most 2 different prime factors. Write n=p1p2⋯pk. Now we must have
p2⋯pk(p2⋯pk+1)∣p1⋯pk(p1⋯pk+1)
which is equivalent to
(p2⋯pk+1)∣p1(p1⋯pk+1)
Now
p1(p1⋯pk+1)−p12(p2⋯pk+1)=p1−p12
Hence,
(p2⋯pk+1)∣p12−p1
Pick p1 to be the smallest prime factor of n. Then p2⋯pk>p12 if k>2, which is a contradiction.
The only relations of concern are p(p+1)∣n(n+1) and q(q+1)∣n(n+1). The former yields p(p+1)∣pq(pq+1), i.e. p+1∣q(pq+1), that is, p+1∣q(q−1); analogously q+1∣p(p−1). Without loss of generality assume p>q. If q does not divide p+1, then since p+1 and q are relatively prime, we obtain from p+1∣q(q−1) that p+1 divides q−1 and hence p<p+1≤q−1<q, contradiction. Hence q divides p+1. Also, p does not divide q+1 due to p>q (unless p=q+1, i.e. p=3,q=2, but then p+1∤q(q−1) ) and now from q+1∣p(p−1) we deduce that q+1 divides p−1.
Write p+1=mq for some positive integer m. Then q+1 divides p−1=mq−2=m(q+1)−m−2, so q+1 divides m+2. On the other hand, p+1∣q(q−1) implies that m divides q−1. If the ratio of the latter is greater than 1, then m≤2q−1, but q+1 dividing m+2 from above insists on q≤m+1, so m≤2q−1≤2m, contradiction. Hence m=q−1 and we have proven the following: if n>1 is composite and peculiar, then n=pq, where p>q are primes and p=q2−q−1.
Hence from now on we work with n=ab, with a<b primes. The only relations of concern are a(a+1)∣n(n+1) and b(b+1)∣n(n+1). The former yields a(a+1)∣ab(ab+1), i.e. a+1∣b(ab+1), that is, a+1∣b(b−1); analogously b+1∣a(a−1). If a does not divide b+1, then since b+1 and a are relatively prime, we obtain from b+1∣a(a−1) that b+1 divides a−1 and hence b<b+1≤a−1<a, contradiction. Hence a divides b+1. Writing b+1=ac, from b+1∣a(a−1) we have that c divides a−1.
Now from b+1=ac with c∣a−1 and a+1∣b(b−1) we derive that b=a2−a−1 in a different (and longer) manner, without using the primality of b until the final step. Note that since b=ac−1≡−c−2(moda+1), we have a+1∣(c+1)(c+2). Since c∣a−1, we may write a−1=cd, and note that cd≡−2(moda+1). We then have
0≡−d(c+1)(c+2)=−c2d−3cd−2d≡2(c−d)+6(moda+1)
Since 1≤c,d≤a−1, we have ∣2(c−d)+6∣≤2(a−2)+6=2a+2. Furthermore, we have c≥2 (otherwise b=ac−1<a, contradiction), and so d≤2a−1, thus 2(c−d)+6≥11−a>−(a+1). From these inequalities and equivalence (1), we obtain 2(c−d)+6∈{0,a+1,2(a+1)}. It is immediate to verify that these three possible values correspond to the following three cases, respectively:
(i) d=c+3 and a=c2+3c+1;
(ii) d=2,c=2a−1,a=2c+1; and
(iii) d=1,c=a−1.
Substituting these values in b=ac−1, we obtain in case (i), b=c3+3c2+c−1=(c+1)(c2+2c−1), and in case (ii), b=2c2+c−1=(c+1)(2c−1). In both cases these decompositions contradict the primality of b, as c≥2. Thus we are left only with case (iii), where b=a2−a+1, as claimed.
In conclusion, a prime p can only divide the peculiar numbers p, p(p2−p−1) (if p2−p−1 is prime) and pq where p=q2−q−1 (if q is prime). In particular, the greatest common divisor of more than three peculiar numbers must be 1. (Note that three peculiar numbers may indeed have a common divisor, e.g. 5,15=3⋅5 and 95=5⋅19 are all peculiar.)