OFFSET
1,2
COMMENTS
Values of n for which all numbers in row A238747(n) are even. Also, numbers n such that A000005(n^m) is a perfect square for all nonnegative integers m; numbers n such that A181819(n) is a perfect square; numbers n such that A182860(n) is odd.
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 3, 33, 314, 3119, 31436, 315888, 3162042, 31626518, 316284320, 3162915907, ... . Apparently, the asymptotic density of this sequence exists and equals 0.3162... . - Amiram Eldar, Nov 28 2023
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..10000
EXAMPLE
The prime signature of 36 = 2^2 * 3^2 is {2,2}. One distinct integer (namely, 2) appears in the prime signature, and it appears an even number of times (2 times). Hence, 36 appears in the sequence.
The prime factorization of 1260 = 2^2 * 3^2 * 5^1 * 7^1. Exponent 2 occurs twice (an even number of times), as well as exponent 1, thus 1260 is included. It is also the first term k > 1 in this sequence for which A182850(k) = 4, not 3. - Antti Karttunen, Feb 06 2016
MATHEMATICA
q[n_] := n == 1 || AllTrue[Tally[FactorInteger[n][[;; , 2]]][[;; , 2]], EvenQ]; Select[Range[200], q] (* Amiram Eldar, Nov 28 2023 *)
PROG
(Scheme)
(define (square? n) (not (zero? (A010052 n))))
;; Requires also MATCHING-POS macro from my IntSeq-library - Antti Karttunen, Feb 06 2016
(PARI) is(n) = {my(e = factor(n)[, 2], m = #e); if(m%2, return(0)); e = vecsort(e); forstep(i = 1, m, 2, if(e[i] != e[i+1], return(0))); 1; } \\ Amiram Eldar, Nov 28 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
Matthew Vandermast, May 08 2014
STATUS
approved