login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A304822 a(n) = A304821(n+3)/A304821(n) - 1. 1

%I

%S 1,1,2,5,1,7,1,3,5,11,1,13,7,5,1,17,1,19,5,7,11,23,1,5,13,3,1,29,5,31,

%T 1,11,17,1,1,37,19,13,1,41,1,43,11,5,23,47,1,7,1,17,13,53,1,1,1,19,29,

%U 59,1,61,31,1,1,13,11,67,17,23,1,71,1,73

%N a(n) = A304821(n+3)/A304821(n) - 1.

%C Conjecture: This sequence consists of 1's and primes only.

%C a(n) divides n+1. - _A.H.M. Smeets_, Jul 02 2018

%H A.H.M. Smeets, <a href="/A304822/b304822.txt">Table of n, a(n) for n = 1..20000</a>

%t Drop[#, 3] &@ Fold[Append[#1, {#2, #2/#1 - 1} & @@ {#1[[-3, 1]], #1[[-3, 1]] + LCM[#1[[-3, 1]], #2 + 1]}] &, {{6, 0}, {6, 0}, {6, 0}}, Range@ 72] [[All, -1]] (* _Michael De Vlieger_, May 20 2018 *)

%o (PARI) Generator(n)={b1=6;b2=6;b3=6;list=[];for(k=4,n,b4=b1+lcm(k-2,b1);a=b4/b1-1;list=concat(list,a);b1=b2;b2=b3;b3=b4);print(list)}

%o (Python)

%o from math import gcd

%o n,an0,an1,an2 = 3,6,6,6

%o while n-3 < 200:

%o n += 1

%o an0,an1,an2,an3 = an2+an2*(n-2)//gcd(an2,n-2),an0,an1,an2

%o # an0 = A304821(n), an3 = A304821(n-3)

%o print(n-3,an0//an3-1)

%o # _A.H.M. Smeets_, Jul 02 2018

%Y Cf. A135504, A217662, A217663, A304821.

%K nonn

%O 1,3

%A _Pedja Terzic_, May 19 2018

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified June 17 05:30 EDT 2021. Contains 345080 sequences. (Running on oeis4.)