%I #22 Jan 03 2022 21:37:58
%S 1,2,3,7,13,14,25,26,37,38,49,50,61,62,73,74,85,86,97,98,109,110,121,
%T 122,133,134,145,146,157,158,169,170,181,182,193,194,205,206,217,218,
%U 229,230,241,242,253,254,265,266,277,278,289,290,301,302,313,314
%N {3, 7} together with the numbers congruent to {1, 2} mod 12.
%C The old definition of this sequence was "a(n) = least positive integer > a(n-1) and not equal to a(i)+a(j)+a(k) for 1<=i<=j<=k<=n". However, Ralf Stephan observes that this does not fit the terms shown. The present definition (due to Stephan) has been adopted as a temporary solution. - _N. J. A. Sloane_, Nov 24 2004
%C Regarding the old definition, see Comments at A047239. - _Clark Kimberling_, Oct 09 2019
%H Clark Kimberling, <a href="/A026472/b026472.txt">Table of n, a(n) for n = 1..10000</a>
%H <a href="/index/Rec#order_03">Index entries for linear recurrences with constant coefficients</a>, signature (1,1,-1).
%F From _Colin Barker_, Oct 10 2019: (Start)
%F G.f.: x*(1 + x + 3*x^3 + 5*x^4 - 3*x^5 + 5*x^6) / ((1 - x)^2*(1 + x)).
%F a(n) = a(n-1) + a(n-2) - a(n-3) for n>6.
%F a(n) = -(39/2) - (5*(-1)^n)/2 + 6*n for n>4.
%F (End)
%t p = {1, 2, 3, 7}; r = 12 Range[200]; Union[p, 1 + r, 2 + r] (* _Clark Kimberling_, Oct 10 2019 *)
%o (PARI) Vec(x*(1 + x + 3*x^3 + 5*x^4 - 3*x^5 + 5*x^6) / ((1 - x)^2*(1 + x)) + O(x^40)) \\ _Colin Barker_, Oct 20 2019
%Y Cf. A026476, A026474.
%K nonn
%O 1,2
%A _Clark Kimberling_
%E More terms from _Clark Kimberling_, Oct 10 2019
|