The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified March 24 15:19 EDT 2023. Contains 361479 sequences. (Running on oeis4.)