|
|
A296084
|
|
a(1) = 0 and for n > 1, a(n) = 1 if tau(n)-1 divides sigma(n)-1, 0 otherwise. Here tau = A000005, sigma = A000203.
|
|
5
|
|
|
0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1
|
|
COMMENTS
|
Also number of numbers of the form i*n with 1 <= i <= n and tau(i*n) = 3 (equivalently, i*n is the square of a prime). - N. J. A. Sloane, Nov 11 2020
|
|
LINKS
|
|
|
FORMULA
|
|
|
MATHEMATICA
|
Join[{0}, Table[If[Divisible[DivisorSigma[1, n]-1, DivisorSigma[0, n]-1], 1, 0], {n, 2, 120}]] (* Harvey P. Dale, Dec 15 2018 *)
|
|
PROG
|
(PARI) A296084(n) = if(1==n, 0, !((sigma(n)-1)%(numdiv(n)-1)));
(Python)
from math import prod
from sympy import factorint
f = factorint(n).items()
return int(not (prod((p**(e+1)-1)//(p-1) for p, e in f)-1)%(prod(e+1 for p, e in f)-1)) if n>1 else 0 # Chai Wah Wu, Oct 14 2023
|
|
CROSSREFS
|
Characteristic function of A284288.
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|