OFFSET
1,4
COMMENTS
If A067131(n) = 2, then a(n) = d(n)*(d(n)-1)/2, where d(n) is the number of divisors of n.
a(p) = 1 for all primes p.
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..65537
FORMULA
EXAMPLE
The divisors of 18 are 1,2,3,6,9,18. There are 2 sets of these divisors, (1,2,3) and (3,6,9), that have their terms in arithmetic progression and that each have the maximal number (3) of such divisors of 18. So a(18) = 2.
PROG
(PARI) A160752(n) = if(1==n, n, my(d=divisors(n), m=1, counts=vector(#d)); for(i=1, (#d-1), for(j=(i+1), #d, my(c=1, k=d[j], s=(d[j]-d[i])); while(!(n%k), k+=s; c++); counts[c]++; m = max(m, c))); (counts[m])); \\ Antti Karttunen, Sep 21 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, May 25 2009
EXTENSIONS
Extended by Ray Chandler, Jun 15 2009
STATUS
approved