|
%I #24 Oct 17 2016 11:41:45
%S 1,2,11,17,27,140,213,1934,13871,38405,171707,200938,389409,1633357,
%T 5460156,8405437,41983357,68780189,179376463,130292951546,
%U 393142151459,2100234982892,5942636062287
%N Start of record runs with lambda(k) = lambda(k+1) = ..., where lambda is Liouville's function A008836.
%e Lambda(1) = 1 is the first (record) run, so a(1) = 1.
%e Lambda(2) = lambda(3) = -1 is the second record run, so a(2) = 2.
%e Lambda(11) = lambda(12) = lambda(13) = -1 is the third record run, so a(3) = 11.
%t sz[n_] := Module[{t = LiouvilleLambda[n], k = n}, While[LiouvilleLambda[k++] == t]; k - n]; r = 0; Reap[For[n = 1, n <= 10^6, n++, t = sz[n]; If[t > r, r = t; Print[t, " ", n]; Sow[n]]]][[2, 1]] (* _Jean-François Alcover_, Oct 17 2016, adapted from PARI *)
%o (PARI) L(n)=(-1)^bigomega(n);
%o sz(n)=my(t=L(n),k=n);while(L(k++)==t,);k-n
%o r=0;for(n=1,1e9,t=sz(n);if(t>r,r=t;print(t" "n)))
%Y Cf. A066793, A066794, A175201, A233459, A008836.
%K nonn
%O 1,2
%A _Charles R Greathouse IV_, Dec 10 2013
|
