OFFSET
1,1
COMMENTS
Number of primes obtained from the sequence ‘prime plus its digit sum is perfect square’ is 150 for n = 1 to 3*10^5, while the sequence for ‘perfect cube’ yields only 11 primes for the same range of n. Hence, sequence for ‘square’ is framed.
Subsequence of primes of A066564. - Michel Marcus, Jun 02 2015
LINKS
K. D. Bajpai and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000 (terms 2..150 from Bajpai)
EXAMPLE
a(2) = 17 is prime. Digit sum of 17 = 8, 17 + 8 = 25 = 5^2.
a(5) = 467 is prime. Digit sum of 467 = 17, 467 + 17 = 484 = 22^2.
MAPLE
KD:= proc() local a, b, c, d; a:= ithprime(n); b:=add( i, i = convert((a), base, 10))(a); c:=a+b; d:=evalf(sqrt(c)); if d=floor(d) then return (a) :fi; end:seq(KD(), n=1..50000);
PROG
(PARI) for(n=2, 1e4, forprime(p=n^2-9*#digits(n^2), n^2, if(p+sumdigits(p) == n^2, print1(p", ")))) \\ Charles R Greathouse IV, Oct 08 2013
(Magma) [p: p in PrimesUpTo(6*10^5) | IsSquare(p+(&+Intseq(p)))]; // Vincenzo Librandi, Jun 02 2015
CROSSREFS
KEYWORD
nonn,base,less
AUTHOR
K. D. Bajpai, Oct 08 2013
EXTENSIONS
a(1) from Charles R Greathouse IV, Oct 08 2013
STATUS
approved