A354996 a(n) = 1 if n > A034699(n)^2, otherwise 0, where A034699(n) is the largest prime power factor of n.

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

Links

Antti Karttunen, Table of n, a(n) for n = 1..100000

Index entries for characteristic functions

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

Adjacent sequences: A354993 A354994 A354995 * A354997 A354998 A354999

Keyword

nonn

Author

Antti Karttunen, Jun 18 2022

Status

approved