|
%I #6 Nov 04 2020 05:57:52
%S 1,2,3,7,9,14,19,21,24,26,30,34,40,42,48,51,122,128,174,207,211,247,
%T 264,296,300,310,354,368,408,438,448,454,464,490,494,520,538,728,746,
%U 798,828,866,1008,1010,1026,1076,1188,1198,1262,1370,1376,1403,1474,1542
%N Numbers whose position is not changed by the rearrangement defined by A096070; or k such that A096070(k) = k.
%t Module[{n, m = 1550, s = 0, v = {}}, Reap[For[n = 1, n <= m, n++, k = 1; While[MemberQ[v, k] || PrimeQ[s + k] != EvenQ[n], k++]; If[n == k, Sow[k]]; s = s + k; v = Union[v, {k}]]][[2, 1]]] (* _Jean-François Alcover_, Nov 03 2020, after PARI *)
%o (PARI) {m=1550;s=0;v=Set([]);for(n=1,m,k=1;while(setsearch(v,k)>0||isprime(s+k)!=(n%2==0),k++);if(n==k,print1(k,","));s=s+k;v=setunion(v,Set(k)))}
%Y Cf. A096070.
%K nonn
%O 1,2
%A _Amarnath Murthy_, Jun 20 2004
%E Edited, corrected and extended by _Klaus Brockhaus_, Jun 21 2004
|
