|
| |
|
|
A071821
|
|
Numbers whose largest prime factor is of the form 4k+1.
|
|
2
|
|
|
|
5, 10, 13, 15, 17, 20, 25, 26, 29, 30, 34, 37, 39, 40, 41, 45, 50, 51, 52, 53, 58, 60, 61, 65, 68, 73, 74, 75, 78, 80, 82, 85, 87, 89, 90, 91, 97, 100, 101, 102, 104, 106, 109, 111, 113, 116, 117, 119, 120, 122, 123, 125, 130, 135, 136, 137, 143, 145, 146, 148, 149
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
|
OFFSET
|
1,1
|
|
|
COMMENTS
|
Subsequence of A009003. - M. F. Hasler, Feb 06 2009
|
|
|
LINKS
|
Robert Israel, Table of n, a(n) for n = 1..10000
|
|
|
FORMULA
|
Numbers k such that A006530(k) == 1 (mod 4).
|
|
|
MAPLE
|
filter:= proc(n)
max(numtheory:-factorset(n)) mod 4 = 1
end proc:
select(filter, [$1..200]); # Robert Israel, Sep 11 2020
|
|
|
MATHEMATICA
|
Select[Range[2, 150], Mod[FactorInteger[#][[-1, 1]], 4] == 1 &] (* Amiram Eldar, May 04 2022 *)
|
|
|
PROG
|
(PARI) for(n=2, 200, if((component(component(factor(n), 1), omega(n))-1)%4==0, print1(n, ", ")))
(PARI) for( n=2, 99, vecmax(factor(n)[, 1])%4==1 && print1(n", ")) \\ M. F. Hasler, Feb 06 2009
|
|
|
CROSSREFS
|
Cf. A004431, A006530, A009000, A009003, A073503, A083025.
Sequence in context: A057100 A304436 A009003 * A334678 A201012 A360020
Adjacent sequences: A071818 A071819 A071820 * A071822 A071823 A071824
|
|
|
KEYWORD
|
easy,nonn
|
|
|
AUTHOR
|
Benoit Cloitre, Jun 07 2002
|
|
|
STATUS
|
approved
|
| |
|
|
