A277851 - OEIS

%I #8 Nov 25 2016 12:34:01

%S 3,6,7,8,11,13,14,18,19,20,21,22,23,25,26,29,30,32,34,35,38,39,40,43,

%T 44,45,47,51,53,54,56,58,62,67,68,69,70,71,72,74,75,77,78,80,82,83,87,

%U 89,90,91,92,94,97,98,99,102,103,104,105,106,107,108,110,111,112,115,117,118,119,120,122,123,125,126

%N Numbers not occurring in A075771 (= q(n) + r(n), with n^2 = prime(n)*q(n) + r(n), 0 <= r(n) < prime(n)).

%C The lower bound A075771(n) >= n^2/prime(n) ensures that a given number can't occur beyond a certain index in that sequence.

%t mx = 6300; Take[ Complement[ Range@ mx, Array[Plus @@ IntegerDigits[#^2, Prime[#]] &, mx]], mx/90] (* _Robert G. Wilson v_, Nov 25 2016 *)

%o (PARI) A277851_vec(N,L=N^2\prime(N),A=A075771)=setminus([1..L],Set(vector(N,n,A(n))))

%Y Cf. A075771, A277852, A277853, A277854.

%K nonn

%O 1,1

%A _M. F. Hasler_, Nov 25 2016