|
%I #15 Aug 11 2014 22:23:53
%S 1,2,3,4,5,7,11,17,24,31,59,89,127,184,277,461,669,709,1097,1787,1995,
%T 3259,4999,5381,8807,15299,17351,30133,48593,52711,60810,91081,167449,
%U 192263
%N Fixed points of A245821 and A245822.
%C If n is a member, then the n-th prime, p_n (= A000040(n)) is also a member. If p_n is a member, then its index n is also a member. Thus the sequence is completely determined by its nonprime terms: 1, 4, 24, 184, 669, 1995, 60810, ..., and is obtained as a union of prime recurrences that start with those values: A007097 U A057450 U ..., etc.
%o (PARI)
%o default(primelimit, 2^30);for(n=1, 2^18, if(A245821(n) == n, print1(n, ", "))); \\ Other code as in A245821.
%o (Scheme, with _Antti Karttunen_'s IntSeq-library)
%o (define A245823 (FIXED-POINTS 1 1 A245821))
%Y A007097 and A057450 are subsequences.
%Y Cf. A000040, A245821, A245822.
%K nonn
%O 1,2
%A _Antti Karttunen_, Aug 02 2014
|
