|
%I #19 Jun 26 2022 05:53:02
%S 23,37,53,59,61,83,103,107,113,127,137,139,149,151,197,211,223,227,
%T 229,331,347,349,353,359,383,421,439,461,479,491,509,523,529,541,557,
%U 563,569,607,631,739,751,757,761,769,797,809,811,821,823,827,829,839,851
%N Numbers k such that the continued fraction for m/k contains a term >= 3 for every 1 <= m < k.
%C Composite terms begin a(33) = 529, a(53) = 851, a(77) = 1357, .... - _Charles R Greathouse IV_, Mar 11 2014
%H Amiram Eldar, <a href="/A055114/b055114.txt">Table of n, a(n) for n = 1..10000</a>
%e 23 is in sequence because continued fractions for 1/23, 2/23, ..., 22/23 each contain a term >= 3.
%t q[n_] := AllTrue[Range[1, n - 1], Max[ContinuedFraction[#/n]] > 2 &]; Select[Range[2, 1000], q] (* _Amiram Eldar_, Jun 25 2022 *)
%o (PARI) drop(v)=v[2..#v]
%o is(n)=for(k=1,n-1,if(vecmax(drop(contfrac(k/n)))<3, return(0))); n>1 \\ _Charles R Greathouse IV_, Mar 11 2014
%Y Cf. A082409.
%K nonn
%O 1,1
%A _David W. Wilson_, Jun 16 2000
|
