login
A354996
a(n) = 1 if n > A034699(n)^2, otherwise 0, where A034699(n) is the largest prime power factor of n.
3
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1
OFFSET
1
FORMULA
a(n) = [A034699(n)^2 < n], where [ ] is the Iverson bracket.
EXAMPLE
For n = 30, A034699(30) = 5 and as 5^2 < 30, a(30) = 1.
MATHEMATICA
a[n_] := If[n > Max[Power @@@ FactorInteger[n]]^2, 1, 0]; Array[a, 100] (* Amiram Eldar, Jun 18 2022 *)
PROG
(PARI)
A034699(n) = if(1==n, n, my(f=factor(n)); vecmax(vector(#f[, 1], i, f[i, 1]^f[i, 2])));
A354996(n) = (A034699(n)^2 < n);
CROSSREFS
Characteristic function of A051283.
Cf. A034699.
Sequence in context: A011745 A011744 A358754 * A351539 A354926 A011743
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 18 2022
STATUS
approved