OFFSET
1,1
COMMENTS
A156207 read without the 1's and without the composites. - R. J. Mathar, Sep 07 2016
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
Given a number n with digits d1d2d3...dm, a(n) = d1*1+d2*2+d3*3+...+dm*m.
If a(n) is prime, list it.
EXAMPLE
For n=19 we have 1*1 + 2*9 = 19 prime and the sequence.
MAPLE
f:= proc(n) local L, i, a;
L:= convert(n, base, 10);
a:= add(L[-i]*i, i=1..nops(L));
if isprime(a) then a else NULL fi
end proc:
map(f, [$1..1000]); # Robert Israel, Sep 07 2016
PROG
(PARI) g1(n) = for(j=1, n, if(isprime(g(j)), print1(g(j)", ")))
g(n) = v=Vec((Str(n))); sum(x=1, length(v), x*eval(v[x]))
CROSSREFS
KEYWORD
base,nonn,less
AUTHOR
Cino Hilliard, Feb 05 2009, Feb 08 2009
EXTENSIONS
Definition clarified. - R. J. Mathar, Sep 07 2016
STATUS
approved