|
|
A116199
|
|
a(n) = the number of positive divisors of n which are coprime to sigma(n) = A000203(n).
|
|
2
|
|
|
1, 2, 2, 3, 2, 1, 2, 4, 3, 2, 2, 2, 2, 2, 2, 5, 2, 2, 2, 2, 4, 2, 2, 1, 3, 2, 4, 1, 2, 2, 2, 6, 2, 2, 4, 9, 2, 2, 4, 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 6, 2, 2, 2, 1, 4, 2, 4, 2, 2, 2, 2, 2, 6, 7, 4, 2, 2, 2, 2, 4, 2, 4, 2, 2, 6, 2, 4, 2, 2, 2, 5, 2, 2, 2, 4, 2, 2, 2, 2, 2, 2, 2, 4, 2, 2, 1, 2, 6, 2, 9, 2, 2, 2, 2, 4
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
In the first 1000 terms, only 69 are odd. - Harvey P. Dale, Jul 16 2016
|
|
LINKS
|
|
|
EXAMPLE
|
The sum of the positive divisors of 12 is 1+2+3+4+6+12 = 28. There are 2 positive divisors (1 and 3) of 12 which are coprime to 28. So a(12) = 2.
|
|
MAPLE
|
with(numtheory): a:=proc(n) local div, ct, j: div:=divisors(n): ct:=0: for j from 1 to tau(n) do if igcd(div[j], sigma(n))=1 then ct:=ct+1 else ct:=ct fi od: ct: end: seq(a(n), n=1..140); # Emeric Deutsch, May 05 2007
|
|
MATHEMATICA
|
pdc[n_]:=Module[{s=DivisorSigma[1, n]}, Count[Divisors[n], _?(CoprimeQ[ #, s]&)]]; Array[pdc, 110] (* Harvey P. Dale, Jul 16 2016 *)
|
|
PROG
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|