OFFSET
1,1
COMMENTS
A005091(a(n)) > 0. - Reinhard Zumkeller, Jan 07 2013
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
FORMULA
a(n) ~ n. - Charles R Greathouse IV, Sep 01 2015
MATHEMATICA
pfQ[n_]:=AnyTrue[Transpose[FactorInteger[n]][[1]], Mod[#, 4]==3&]; Select[ Range[100], pfQ] (* The program uses the AnyTrue function from Mathematica version 10 *) (* Harvey P. Dale, Jan 26 2016 *)
PROG
(Haskell)
import Data.List (findIndices)
a187811 n = a187811_list !! (n-1)
a187811_list = map (+ 1) $ findIndices (> 0) a005091_list
-- Reinhard Zumkeller, Jan 07 2013
(PARI) is(n)=if(n%4==3, return(1)); my(f=factor(n)[, 1]%4); for(i=1, #f, if(f[i]==3, return(1))); 0 \\ Charles R Greathouse IV, Sep 01 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jan 07 2013
STATUS
approved