login
a(n) is the index of A210487 when the n-th prime appears for the first time.
1

%I #17 Apr 17 2013 04:18:31

%S 3,4,2,20,60,92,246,752,1289,2084,3383

%N a(n) is the index of A210487 when the n-th prime appears for the first time.

%e In A210487, 2 appears for the first time as A210487(3), so a(1)=3;

%e 3 appears for the first time as A210487(4), so a(2)=4;

%e 5 appears for the first time as A210487(2), so a(3)=2.

%p N := 60 ;

%p A210488 := [seq(0,i=1..N)] ;

%p for n from 2 do

%p a := A210487(n) ;

%p if isprime(a) then

%p idx := numtheory[pi](a) ;

%p if idx <= N then

%p if op(idx,A210488) = 0 then

%p A210488 := subsop(idx=n,A210488) ;

%p print(A210488) ;

%p end if;

%p end if;

%p end if;

%p end do: # _R. J. Mathar_, Apr 17 2013

%t nn = 11; a = Table[0, {nn}]; k = 1; While[Times @@ a == 0, k++; p1 = Prime[k]; id = PrimePi[Min[Table[Last[FactorInteger[Prime[k]^2 - Prime[j]^2]][[1]], {j, k - 1}]]]; If[id <= nn && a[[id]] == 0, a[[id]] = k]]; a

%Y Cf. A000040, A210487, A147971.

%K nonn,hard

%O 1,1

%A _Lei Zhou_, Jan 23 2013