login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Complement of A005563.
3

%I #28 Feb 04 2022 14:42:56

%S 1,2,4,5,6,7,9,10,11,12,13,14,16,17,18,19,20,21,22,23,25,26,27,28,29,

%T 30,31,32,33,34,36,37,38,39,40,41,42,43,44,45,46,47,49,50,51,52,53,54,

%U 55,56,57,58,59,60,61,62,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,81,82,83,84,85,86,87,88

%N Complement of A005563.

%C From _Andres Cicuttin_, Apr 18 2016: (Start)

%C Defining the probability for integer value x in a Poisson distribution of integer mean = q > 0 as P(q,x) = e^(-q)*(q^x)/x! then it is conjectured that P(n, a(n)+1)- P(n, a(n)) <= P(n, k+1)- P(n, k), for every k > 0.

%C That is, a(n) is the position where the Poisson distribution with mean n has its minimum discrete difference (not proved, but tested up to n = 20*10^3).

%C (Very qualitative) Plot of a Poisson Distribution with mean q = n. The vertical line above a(n) indicates the place where the distribution has its minimum (negative) discrete difference.

%C P

%C ^

%C | *

%C | * *

%C | * *

%C | * *

%C | * *

%C | * | *

%C | * | *

%C | * | *

%C | * | *

%C *-------------+----+----------------------> x

%C n a(n)

%C For instance, if n = 8 then

%C P(8, a(8)+1) - P(8, a(8)) = P(8,11) - P(8,10) = -0.027071

%C If we now calculate the discrete difference in a(n)+1 we then obtain

%C P(8,a(8)+2) - P(8,a(8)+1) = P(8,12)- P (8,11) = -0.0240634

%C and in a(n)-1

%C P(8,a(8)) - P(8,a(8)-1) = P(8,10) - P(8,9) = -0.0248154

%C Both previous values are larger than the minimum obtained at a(n). (End)

%C Numbers k such that sqrt(k+1) is not an integer. - _Wesley Ivan Hurt_, Feb 03 2022

%F (See the Mathematica code.)

%t a=1; b=2;

%t F[n_]:=a*n^2+b*n;

%t R[n_]:=(n/a+((b-1)/(2a))^2)^(1/2);

%t G[n_]:=n-1+Ceiling[R[n]-(b-1)/(2a)];

%t Table[F[n], {n,60}]

%t Table[G[n], {n,100}]

%Y Cf. A005563.

%K nonn

%O 1,2

%A _Clark Kimberling_, Jan 03 2011