login
A395053
a(n) = 1 if the greatest prime dividing A276085(n) is larger than the greatest prime dividing n, otherwise 0, where A276085 is fully additive with a(p) = p#/p. a(1) = 0 by convention.
4
0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1
OFFSET
1
FORMULA
a(n) = [n>1 and A395052(n) > A006530(n)], where [ ] is the Iverson bracket.
PROG
(PARI)
A276085(n) = { my(f=factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1, primepi(f[k, 1]-1), prime(i))); };
A395053(n) = if(n<=1 || isprime(n), 0, my(f = factor(n), gpf = vecmax(f[, 1]), x = A276085(n)); forprime(p=2, gpf, x /= (p^valuation(x, p))); (x>1));
CROSSREFS
Characteristic function of A395054.
Cf. also A395043.
Sequence in context: A011693 A144198 A144200 * A011695 A011706 A011697
KEYWORD
nonn,easy
AUTHOR
Antti Karttunen, Apr 15 2026
STATUS
approved