cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A236069 Primes p such that f(f(p)) is prime where f(x) = x^4 + 1.

This page as a plain text file.
%I A236069 #10 Jan 26 2014 03:22:40
%S A236069 3,79,83,107,211,401,491,881,1013,1061,1367,1637,1669,1811,2029,2309,
%T A236069 2399,2459,2671,2713,2963,3109,3203,3407,3593,3709,3733,3929,4219,
%U A236069 4457,4513,4639,4703,4729,5417,5641,6047,6113
%N A236069 Primes p such that f(f(p)) is prime where f(x) = x^4 + 1.
%F A236069 a(n) = (A235982(n)-1)^(1/4).
%e A236069 881 is prime and (881^4+1)^4+1 is also prime. So, 881 is a member of this sequence.
%o A236069 (Python)
%o A236069 import sympy
%o A236069 from sympy import isprime
%o A236069 {print(p) for p in range(10**4) if isprime(p) and isprime((p**4+1)**4+1)}
%o A236069 (PARI) isok(p) = isprime(p) && (q = p^4+1) && isprime(q^4+1); \\ _Michel Marcus_, Jan 19 2014
%Y A236069 Cf. A235982.
%K A236069 nonn
%O A236069 1,1
%A A236069 _Michel Marcus_ and _Derek Orr_, Jan 19 2014