[go: up one dir, main page]

login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
Revision History for A036569 (Underlined text is an addition; strikethrough text is a deletion.)

Showing entries 1-10 | older changes
A036569 Increments used in Sedgewick-Incerpi upper bound for shell sort.
(history; published version)
#20 by N. J. A. Sloane at Mon Oct 19 12:27:44 EDT 2020
STATUS

proposed

approved

#19 by Joerg Arndt at Thu Oct 15 06:38:06 EDT 2020
STATUS

editing

proposed

#18 by Joerg Arndt at Thu Oct 15 06:38:02 EDT 2020
MATHEMATICA

With[{r = Floor[Sqrt[2 k + Sqrt[2 k]]]},

Select[Range[r], # != b &])]]; (* &])]]; (* _Morgan Owens_, Oct 08 2020 *)

STATUS

proposed

editing

#17 by Morgan Owens at Thu Oct 15 06:09:04 EDT 2020
STATUS

editing

proposed

#16 by Morgan Owens at Thu Oct 15 06:07:25 EDT 2020
MATHEMATICA

With[{r = Floor[Sqrt[2 k + Sqrt[2 k]]]},

With[{b = r (r + 1)/2 - k + 1},

Times @@ (A036567 /@

Times @@ (A036567 /@ Select[Range[0, r - 1], # != b &])] &])]]; (* Morgan Owens_, Oct 08 2020 *)

] (* Morgan Owens, Oct 08 2020 *)

Array[A036569, 10]

STATUS

proposed

editing

Discussion
Thu Oct 15 06:09
Morgan Owens: Tweaked the function to conform to the changes to A036567. Added a command to explicitly output some terms.
#15 by Morgan Owens at Thu Oct 08 05:49:33 EDT 2020
STATUS

editing

proposed

Discussion
Thu Oct 08 11:46
Wesley Ivan Hurt: Mmca doesn't produce any terms for me.
Fri Oct 09 05:40
Morgan Owens: I don't have any problem:

In[1]:= A036567[0] = 3;
A036567[q_] := 
 A036567[q] = 
  With[{prev = A036567 /@ Range[0, q - 1]}, 
   Block[{n = Ceiling[(5/2)^(q + 1)]}, 
    While[Nand @@ ((# == 1 &) /@ GCD[prev, n]), n++]; 
    n]];(*_Morgan Owens_,Oct 08 2020*)

In[3]:= A036569[k_] := 
 With[{r = Floor[Sqrt[2 k + Sqrt[2 k]]]}, 
  With[{b = r (r + 1)/2 - k}, 
   Times @@ (A036567 /@ 
      Select[Range[0, 
        r - 1], # != b &])]] (*_Morgan Owens_,Oct 08 2020*)

In[4]:= A036569[8]

Out[4]= 861
#14 by Morgan Owens at Thu Oct 08 05:49:25 EDT 2020
MATHEMATICA

With[{r = Floor[Sqrt[2 k + Sqrt[2 k]]]},

]

] (* Morgan Owens, Oct 08 2020 *)

#13 by Michel Marcus at Wed Oct 07 07:35:34 EDT 2020
STATUS

proposed

editing

#12 by Morgan Owens at Wed Oct 07 05:41:02 EDT 2020
STATUS

editing

proposed

Discussion
Wed Oct 07 07:35
Michel Marcus: please sign new program : see 1st bullet of https://oeis.org/wiki/Style_Sheet#Mathematica
#11 by Morgan Owens at Fri Oct 02 21:33:44 EDT 2020
MATHEMATICA

a[0] = 3;

a[q_] := a[q] =

With[{prev = a /@ Range[0, q - 1]},

Block[{n = Ceiling[(5/2)^(q + 1)]},

While[Nand @@ ((# == 1 &) /@ GCD[prev, n]), n++]; n]];

Times @@ (aA036567 /@ Select[Range[0, r - 1], # != b &])]

Discussion
Fri Oct 02 21:36
Morgan Owens: Added a Mathematica implementation of the sequence; calls on an implementation of A036567 (q.v) to avoid clutter.

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 30 07:09 EDT 2024. Contains 375532 sequences. (Running on oeis4.)