OFFSET
1,2
COMMENTS
Numbers n such that A033570(2n) is semiprime.
LINKS
Harvey P. Dale, Table of n, a(n) for n = 1..1000
FORMULA
a(n) = A130800(n)/2.
EXAMPLE
10 is in this sequence because 4*10+1=41 and 6*10+1=61 are primes.
MAPLE
A255607:=n->`if`(isprime(4*n+1) and isprime(6*n+1), n, NULL): seq(A255607(n), n=1..600); # Wesley Ivan Hurt, Feb 28 2015
MATHEMATICA
Select[Range[600], PrimeQ[4 # + 1] && PrimeQ[6 # + 1] &]
Select[Range[600], AllTrue[{4#, 6#}+1, PrimeQ]&] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Sep 22 2020 *)
PROG
(Magma) [n: n in [1..600] | IsPrime(6*n+1) and IsPrime(4*n+1)];
(PARI) for(n=1, 10^3, if(isprime(4*n+1)&&isprime(6*n+1), print1(n, ", "))) \\ Derek Orr, Mar 01 2015
(PARI) select( is_A255607(n)=isprime(4*n+1)&&isprime(6*n+1), [1..555]) \\ M. F. Hasler, Dec 13 2019
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Feb 28 2015
STATUS
approved