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)).
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
Keywords
Links
- Michel Marcus, Table of n, a(n) for n = 1..10000
- Mathematics Stack Exchange, Generating primes of the form floor(sqrt(3)*n)
Programs
-
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
Comments