login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A175886 Numbers that are congruent to {1, 12} mod 13. 13

%I #41 Sep 04 2022 22:06:16

%S 1,12,14,25,27,38,40,51,53,64,66,77,79,90,92,103,105,116,118,129,131,

%T 142,144,155,157,168,170,181,183,194,196,207,209,220,222,233,235,246,

%U 248,259,261,272,274,285,287,298,300,311,313,324,326,337,339,350

%N Numbers that are congruent to {1, 12} mod 13.

%C Cf. property described by Gary Detlefs in A113801: more generally, these numbers are of the form (2*h*n+(h-4)*(-1)^n-h)/4 (h, n natural numbers), therefore ((2*h*n+(h-4)*(-1)^n-h)/4)^2-1 == 0 (mod h); in this case, a(n)^2-1 == 0 (mod 13).

%H Bruno Berselli, <a href="/A175886/b175886.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 G.f.: x*(1+11*x+x^2)/((1+x)*(1-x)^2).

%F a(n) = (26*n+9*(-1)^n-13)/4.

%F a(n) = -a(-n+1) = a(n-1)+a(n-2)-a(n-3).

%F a(n) = a(n-2)+13.

%F a(n) = 13*A000217(n-1)+1 - 2*sum(a(i), i=1..n-1) for n>1.

%F Sum_{n>=1} (-1)^(n+1)/a(n) = (Pi/13)*cot(Pi/13). - _Amiram Eldar_, Dec 04 2021

%F E.g.f.: 1 + ((26*x - 13)*exp(x) + 9*exp(-x))/4. - _David Lovler_, Sep 04 2022

%t Select[Range[1, 350], MemberQ[{1, 12}, Mod[#, 13]]&] (* _Bruno Berselli_, Feb 29 2012 *)

%t CoefficientList[Series[(1 + 11 x + x^2) / ((1 + x) (1 - x)^2), {x, 0, 55}], x] (* _Vincenzo Librandi_, Aug 19 2013 *)

%t LinearRecurrence[{1,1,-1},{1,12,14},60] (* _Harvey P. Dale_, Oct 23 2015 *)

%o (Haskell)

%o a175886 n = a175886_list !! (n-1)

%o a175886_list = 1 : 12 : map (+ 13) a175886_list

%o -- _Reinhard Zumkeller_, Jan 07 2012

%o (Magma) [n: n in [1..350] | n mod 13 in [1, 12]]; // Bruno Berselli, Feb 29 2012

%o (Magma) [(26*n+9*(-1)^n-13)/4: n in [1..55]]; // _Vincenzo Librandi_, Aug 19 2013

%o (PARI) a(n)=(26*n+9*(-1)^n-13)/4 \\ _Charles R Greathouse IV_, Sep 24 2015

%Y Cf. A000217, A091998, A113801, A005408, A047209, A007310, A047336, A047522, A056020, A090771, A175885, A175887.

%Y Cf. A195045 (partial sums).

%K nonn,easy

%O 1,2

%A _Bruno Berselli_, Oct 08 2010 - Nov 17 2010

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 04:05 EDT 2024. Contains 371235 sequences. (Running on oeis4.)