login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A323388
a(n) = b(n+1)/b(n) - 1 where b(1)=3 and b(k) = b(k-1) + lcm(floor(sqrt(3)*k), b(k-1)).
3
1, 5, 1, 1, 5, 1, 13, 5, 17, 19, 1, 11, 1, 5, 1, 29, 31, 1, 17, 1, 19, 13, 41, 43, 1, 23, 1, 1, 17, 53, 1, 19, 29, 1, 31, 1, 13, 67, 23, 71, 1, 37, 1, 1, 79, 1, 83, 1, 43, 1, 1, 13, 31, 1, 1, 1, 1, 1, 103, 1, 107, 109, 1, 1, 1, 29, 1, 1, 1, 61, 1, 1
OFFSET
1,2
COMMENTS
Conjectures:
1. This sequence consists only of 1's and primes.
2. Every odd prime of the form floor(sqrt(3)*m) greater than 3 is a term of this sequence.
3. At the first appearance of each prime of the form floor(sqrt(3)*m), it is larger than any prime that has already appeared.
The 2nd and 3rd conjectures are proved at the Mathematics Stack Exchange link. - Sungjin Kim, Jul 17 2019
PROG
(PARI) Generator(n)={b1=3; list=[]; for(k=2, n, b2=b1+lcm(floor(sqrt(3)*k), b1); a=b2/b1-1; list=concat(list, a); b1=b2); print(list)}
(PARI) lista(nn)={my(b1=3, b2, va=vector(nn)); for(k=2, nn+1, b2=b1+lcm(sqrtint(3*k^2), b1); va[k-1]=b2/b1-1; b1=b2); va}; \\ Michel Marcus, Aug 20 2022
CROSSREFS
Cf. A184796 (primes of the form floor(sqrt(3)*m)).
Sequence in context: A010333 A378898 A131777 * A260877 A353305 A237888
KEYWORD
nonn
AUTHOR
Pedja Terzic, Jan 13 2019
STATUS
approved