cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A065156 Numbers n such that some Lucas number (A000204) is divisible by n.

This page as a plain text file.
%I A065156 #36 Oct 05 2020 01:04:32
%S A065156 1,2,3,4,6,7,9,11,14,18,19,22,23,27,29,31,38,41,43,44,46,47,49,54,58,
%T A065156 59,62,67,71,76,79,81,82,83,86,94,98,101,103,107,116,118,121,123,124,
%U A065156 127,129,131,134,139,142,151,158,161,162,163,166,167,179,181,191,199
%N A065156 Numbers n such that some Lucas number (A000204) is divisible by n.
%C A065156 From _A.H.M. Smeets_, Sep 20 2020 (Start)
%C A065156 For the Fibonacci numbers, each natural number divides some Fibonacci number (see A001177).
%C A065156 If, for some number m, m divides some Lucas number L_i (=A000204(i)), then, the smallest i satisfies i <= m. (End)
%H A065156 A.H.M. Smeets, <a href="/A065156/b065156.txt">Table of n, a(n) for n = 1..20000</a> (terms 1..1000 from T. D. Noe)
%H A065156 B. Avila and T. Khovanova, <a href="http://arxiv.org/abs/1403.4614">Free Fibonacci Sequences</a>, arXiv preprint arXiv:1403.4614 [math.NT], 2014 and <a href="https://cs.uwaterloo.ca/journals/JIS/VOL17/Avila/avila4.html">J. Int. Seq. 17 (2014) # 14.8.5</a>
%F A065156 Equals {1,2,4} union {p^e | p in A140409 and e > 0} union {2*p^e | p in A140409 and e > 0} union {4*p | p in A053032} union {4*p*q | p, q in A053032}. - _A.H.M. Smeets_, Sep 20 2020
%t A065156 test[ n_ ] := For[ a=1; b=3, True, t=b; b=Mod[ a+b, n ]; a=t, If[ b==0, Return[ True ] ]; If[ a==2&&b==1, Return[ False ] ] ]; Select[ Range[ 200 ], test ]
%t A065156 Take[Flatten[Divisors/@LucasL[Range[200]]]//Union,70] (* _Harvey P. Dale_, Jun 07 2020 *)
%o A065156 (Python)
%o A065156 a, n = 0, 0
%o A065156 while n < 1000:
%o A065156     a, f0, f1, i = a+1, 1, 2, 1
%o A065156     if f1%a == 0:
%o A065156         n = n+1
%o A065156         print(n,a)
%o A065156     else:
%o A065156         while f0%a != 0 and i <= a:
%o A065156             f0, f1, i = f0+f1, f0, i+1
%o A065156         if i <= a:
%o A065156             n = n+1
%o A065156             print(n,a) # _A.H.M. Smeets_, Sep 20 2020
%Y A065156 Complement of A064362. Cf. A000204.
%Y A065156 Cf. A001177, A053032, A140409.
%K A065156 nonn,easy
%O A065156 1,2
%A A065156 _Dean Hickerson_, Oct 18 2001