%I #15 Aug 03 2017 04:26:21
%S 4,5,6,11,19,25,34,36,75,82,87,90,94,237,604,609,614,1583,1592,10466,
%T 10467,10498,10504,10505,70501,70511,180227,180294,180358,180443,
%U 180447,466078,8103422,21058343,21058649,143052872,143052877,143053068
%N Numbers n such that A065863(n) = 1, i.e., prime(n) mod (n - Pi(n)) = 1.
%C A004648, A065134 and A065863 behave similarly; they grow relatively slowly and drop suddenly at unexpected values of n. Parity of A004648 behaves most regularly.
%C Each cluster of entries exceeds the previous cluster by a power of e.
%e For the cluster started at n = 10466 the remainders of A065863(n) are as follows: {9089, 9092, 9117, 9127, 9148, 9159, 1, 1, 9180, 9183, 9182, 9179, 9172, 9169, 9168, 9177, 9176, 9178, 9183, 9192, 43}. It behaves like A004648 or A065134.
%t Do[ If[ Mod[ Prime[n], n-PrimePi[n]] == 1, Print[n]], {n, 1, 150000000}]
%t (* Second program: *)
%t Position[Table[Mod[Prime[n], n - PrimePi[n]], {n, 10^6}], 1] // Flatten (* _Michael De Vlieger_, Jul 30 2017 *)
%Y Cf. A065860-A065863, A065133-A065135, A072608-A072610, A004648, A057809.
%K nonn
%O 1,1
%A _Labos Elemer_, Jun 26 2002
%E Edited by _Robert G. Wilson v_, Jun 27 2002
|