%I #19 Sep 08 2022 08:45:11
%S 1,5,9,13,17,25,29,37,41,49,53,61,73,81,89,97,101,109,113,121,125,137,
%T 149,157,169,173,181,193,197,229,233,241,257,269,277,281,289,293,313,
%U 317,337,349,353,361,373,389,397,401,409,421,433,449,457,461,509,521
%N Prime powers of the form 4n+1.
%H Michael De Vlieger, <a href="/A085759/b085759.txt">Table of n, a(n) for n = 1..10000</a>
%H Wyatt J. Desormeaux, Teresa W. Haynes, Michael A. Henning, <a href="https://doi.org/10.7151/dmgt.2222">Restrained domination in self-complementary graphs</a>, Preprint, Discussiones Mathematicae Graph Theory (2019), 1-13.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/PaleyGraph.html">Paley Graph</a>
%F a(n) ~ 2n log n. - _Charles R Greathouse IV_, Jul 12 2018
%t {1}~Join~Select[1 + 4 Range[130], PrimePowerQ] (* _Michael De Vlieger_, Aug 29 2019 *)
%o (PARI) list(lim)=my(v=List([1])); forprime(p=5,lim\=1, if(p%4==1, listput(v,p))); for(e=2,logint(lim,3), forprime(p=3,sqrtnint(lim,e), if(e%2==0 || p%4==1, listput(v,p^e)))); Set(v) \\ _Charles R Greathouse IV_, Jul 12 2018
%o (Magma) [1] cat [4*k+1:k in [1..140]|IsPrimePower(4*k+1)]; // _Marius A. Burtea_, Sep 07 2019
%Y Subsequence of A000961. A002144 is a subsequence.
%K nonn,easy
%O 1,2
%A _Lekraj Beedassy_, Jul 22 2003
%E Corrected and extended by _Ray Chandler_, Aug 10 2003
|