1, 2, 4, 8, 16,… 33?


24

Виклик

Напишіть функцію / програму, яка виводить або n'-й елемент, або перші nелементи, у добре відому послідовність чисел:

         1, 2, 4, 8, 16 ...

О, зачекайте ... Я забув перші кілька номерів:

1, 1, 1, 1, 2, 4, 8, 16 ...

Чорт забираю, я додам ще кілька для хорошої міри:

1, 1, 1, 1, 2, 4, 8, 16, 33, 69, 146, 312, 673, 1463, 3202, 7050, 15605, 34705 ...

Цифри - це узагальнені каталонські числа, задані формулою (нульова індексація):

a(n+1)=a(n)+k=2n1a(k)a(n1k)

де

a(0)=a(1)=a(2)=a(3)=1

Це OEIS A004149 .

Ви можете вибрати, чи бажаєте ви мати послідовність нульової чи одноіндексованої. Зрозуміло, послідовність повинна бути однаковою, тому ви повинні переписати формулу, якщо у вас є одноіндексована.


Поправте мене , якщо я помиляюся тут, але модифікація для однієї-індексованої формули змінити a(n-1-k)на a(n-k), правильно?
Sumner18

Відповіді:


23

Пітон , 51 байт

f=lambda n,k=2:n<3or k<n and f(k)*f(n-k-2)+f(n,k+1)

Спробуйте в Інтернеті!

Трохи спрощує формулу:

a(n)=k=2n1a(k)a(n2k)

a(1)=a(0)=a(1)=a(2)=1


8
Поздравляю на 100к !!
Стюі Гріффін

Оскільки я також прийшов до цього рішення самостійно, я мушу сказати, що шлях до нього трохи
кумедний

10

Perl 6 , 44 байти

{1,1,1,1,{sum @_[2..*]Z*@_[@_-4...0,0]}...*}

Спробуйте в Інтернеті!

Блок анонімного коду, який повертає ледачу нескінченну послідовність значень. Це в значній мірі реалізує послідовність, як описано, за допомогою ярлика, в якому поштовий індекс перемножує всі елементи поки що після другого елемента з реверсом списку, починаючи з четвертого елемента, і додаючи додатковий 1в кінці.

Пояснення:

{                                          }  # Anonymous code block
                                       ...*   # Create an infinite sequence
 1,1,1,1,                                     # Starting with four 1s
         {                            }       # Where each new element is:
          sum                                   # The sum of
              @_[2..*]                          # The second element onwards
                      Z*                        # Zip multiplied with
                        @_[@_-4...0  ]          # The fourth last element backwards
                                   ,0           # And 1

10

05AB1E , 14 13 11 байт

$ƒˆ¯Âø¨¨¨PO

Спробуйте в Інтернеті!

Виводить n-й елемент, індексований 0.

$                # push 1 and the input
 ƒ               # repeat (input+1) times
  ˆ              #  add the top of the stack (initially 1) to the global array
   ¯             #  push the global array
    Â            #  and a reversed copy of it
     ø           #  zip the two together, giving a list of pairs
      ¨¨¨        #  drop the last 3 pairs
         P       #  take the product of each pair (or 1 if the list is empty)
          O      #  take the sum of those products
                 #  after the last iteration, this is implicitly output;
                 #  otherwise, it's added to the global array by the next iteration




6

05AB1E , 17 13 байт

4Å1λ£₁λ¨Â¦¦s¦¦*O+

Не коротший за відповідь 05AB1E , але я хотів спробувати рекурсивну функціональність нової версії 05AB1E як практику для себе. Можливо, може бути поле для гри в кілька байтів. EDIT: І справді можна побачити рекурсивну версію відповіді @Grimy 05AB1E нижче, що становить 13 байт .

Виводить перше нелементи: Спробуйте в Інтернеті .

Можна змінити на 0 н'-й предмет при заміні £на è: Спробуйте його в Інтернеті ;
або нескінченний список, видаливши £: Спробуйте його в Інтернеті .

Пояснення:

Це реалізує формулу, використану в описі виклику, як це:
а(н)=а(н-1)+к=2н-1(а(к)а(н-1-к))

а(0)=а(1)=а(2)=а(3)=1

   λ               # Create a recursive environment,
    £              # to output the first (implicit) input amount of results after we're done
4Å1                # Start this recursive list with [1,1,1,1], thus a(0)=a(1)=a(2)=a(3)=1
                   # Within the recursive environment, do the following:
      λ            #  Push the list of values in the range [a(0),a(n)]
       ¨           #  Remove the last one to make the range [a(0),a(n-1)]
        Â          #  Bifurcate this list (short for Duplicate & Reverse copy)
         ¦¦        #  Remove the first two items of the reversed list,
                   #  so we'll have a list with the values in the range [a(n-3),a(0)]
           s       #  Swap to get the [a(0),a(n-1)] list again
            ¦¦     #  Remove the first two items of this list as well,
                   #  so we'll have a list with the values in the range [a(2),a(n-1)]
              *    #  Multiply the values at the same indices in both lists,
                   #  so we'll have a list with the values [a(n-3)*a(2),...,a(0)*a(n-1)]
               O   #  Take the sum of this list
               +  #  And add it to the a(n-1)'th value
                   # (afterwards the resulting list is output implicitly)

13 байт версії @Grimy (переконайтеся , що upvote своєї відповіді , якщо у вас ще немає!):

1λ£λ1šÂ¨¨¨øPO

Виводить перше нелементи: Спробуйте в Інтернеті.

Можна знову змінити на 0-індексування або нескінченний список замість цього:
- (на основі 0) індексація 1λèλ1šÂ¨¨¨øPO: спробуйте в Інтернеті ;
- Нескінченний список λλ1šÂ¨¨¨øPO: Спробуйте в Інтернеті . (Зверніть увагу, що тут збережено 2 байти замість 1, оскільки починається рекурсивне середовищеа(0)=1 за замовчуванням.)

Пояснення:

Він замість цього реалізує формулу, знайдену @xnor для своєї відповіді Python, як це:
a(n)=k=2n1(a(k)a(n2k))

a(1)=a(0)=a(1)=a(2)=1

 λ             # Create a recursive environment,
  £            # to output the first (implicit) input amount of results after we're done
1              # Start this recursive list with 1, thus a(0)=1
               # Within the recursive environment, do the following:
   λ           #  Push the list of values in the range [a(0),a(n)]
    1š         #  Prepend 1 in front of this list
      Â        #  Bifurcate the list (short for Duplicate & Reverse copy)
       ¨¨¨     #  Remove (up to) the last three value in this reversed list
          ø    #  Create pairs with the list we bifurcated earlier
               #  (which will automatically remove any trailing items of the longer list)
           P   #  Get the product of each pair (which will result in 1 for an empty list)
            O  #  And sum the entire list
               # (afterwards the resulting list is output implicitly)

1
Interesting that this can solve a(1200) in 40 second on tio, while other recursive approaches time out for numbers n than 100...
Stewie Griffin

1
I also made (but didn't publish) a recursive version. It's 13 bytes for the first n terms, or 11 bytes for an infinite list. Special-casing a(n-1) costs a lot of bytes and isn't needed (see for example xnor's formula).
Grimmy

@Grimy Do you mind if I add your recursive solutions to my answer (crediting you of course)? I will leave my original answer as well. But it's nice to see the differences between the original formula and xnor's byte-saving formula. :)
Kevin Cruijssen

1
Sure, that's fine!
Grimmy

@StewieGriffin Yeah, I was also impressed by the speed of these recursive infinite functions. Maybe one of Elixir's strengths, and definitely due to the builtin lazy-loading. It calculates n=100 in 0.65 seconds, but when I disable lazy-loading, it will time out after 60 seconds instead, even for n=25.
Kevin Cruijssen





2

Japt, 19 17 16 bytes

Outputs the nth term, 1-indexed.

@Zí*Zz2)Ťx}g4Æ1

Try it

@Zí*Zz2)Ťx}g4Æ1     :Implicit input of integer U
@                    :Function taking an array as an argument via parameter Z
 Zí                  :  Interleave Z with
    Zz2              :  Z rotated clockwise by 180 degrees (simply reversing would be a bye shorter but would modify the original array)
   *                 :  Reduce each pair by multiplcation
       )             :  End interleave
        Å            :  Slice off the first element
         ¤           :  Slice off the first 2 elements
          x          :  Reduce by addition
           }         :End function
            g        :Pass the following as Z, push the result back to it and repeat until it has length U
             4Æ1     :Map the range [0,4) to 1s
                     :Implicit output of the last element

1

Haskell, 65 bytes

f a|a<4=1|z<-g[2..a]=sum$zipWith(*)z$reverse(1:g[0..a-4])
g=map f

Try it online!

You can use either f to get a single element of a sequence, or pass a list of values to g and get all the indexes for that list.


1

Forth (gforth), 99 81 bytes

: f recursive dup 4 > if 0 over 3 do over 1- i - f i f * + loop else 1 then nip ;

Try it online!

Output is nth term and input is 1-indexed

Edit: Saved 17 bytes by switching to xnor's formula. Saved another 1 byte by using 1-indexed

Code Explanation

: f                     \ start a new word definition
  recursive             \ mark that this word will be recursive
  dup 4 >               \ duplicate the input and check if it is greater than 4
  if                    \ if it is:
    0 over              \ create an accumulator and copy n to top of stack
    3 do                \ start counted loop from 3 to n-1
      over 1- i - f     \ recursively calculate f(n-1-i)
      i f               \ recursively calculate f(i)
      * +               \ multiply results and add to accumulator
    loop                \ end the counted loop        
  else                  \ otherwise, if n < 5
    1                   \ put 1 on the stack
  then                  \ end the if block
  nip                   \ drop n from the stack
;                       \ end the word definition

1

Charcoal, 26 bytes

F⁵⊞υ¹FN⊞υΣ✂E⮌υ×κ§υλ³→I§υ±⁴

Try it online! Link is to verbose version of code. Prints the 0-indexed nth number, although it calculates using 1-indexing internally. Explanation:

F⁵⊞υ¹

Start with a[0] = a[1] = a[2] = a[3] = a[4] = 1. Yes, this is 1-indexed, but then with an extra zeroth value. That's code golf for you.

FN

Calculate an additional n terms. This is overkill, but it makes finding the desired term easier when n<5.

⊞υΣ✂E⮌υ×κ§υλ³

For each term, compute the next term as the sum of the terms so far termwise multiplied by the reverse of the terms so far, excluding three terms.

This is a no-op used to trick Charcoal into parsing the 2-argument form of Slice, otherwise I would have to use a less golfy way of removing three terms.

I§υ±⁴

Output the 4th last term.


1

Pyth, 30 bytes

J*4]1VQ=+J+eJsPP*M.t,PJ_PJ0;<J

Try it online!

Returns the first n elements of the sequence.

J*4]1VQ=+J+eJsPP*M.t,PJ_PJ0;<JQ # Full program, last Q = input (implicitly added)
J*4]1                  # J = 4 * [1] (=[1,1,1,1])
VQ                     # for N in range(Q):
  =+J                  #  J +=
     +eJ               #   J[-1] + 
        s              #    sum(                           )
           *M          #     map(__operator_mul,          )
             .t      0 #      transpose(          , pad=0)
               ,       #       [       ,         ]
                PJ     #         J[:-1] 
                  _PJ  #                 J[1::-1]
<JQ                    # J[::Q]

Alternative: Replace < with @ to return the n-th element of the sequence, 0-indexed.



1

Octave, 73 bytes

g=(1:4).^0;for(i=3:(n=input('')))g(i+2)=g(4:i+1)*g(i-(2:i-1))';end;g(end)

Try it online!

-2 bytes thanks to Stewie Griffin. Once more, the imperative approach wins over the functional recursive approach. That one is shown below.

Octave, 75 bytes

f(f=@(a)@(n){@()sum(arrayfun(@(k)a(a)(k)*a(a)(n-2-k),2:n-1)),1}{2-(n>3)}())

Try it online!

Captcha wanted to verify I was a human when posting this. To be honest, I'm not so sure.


I can't see any obvious ways to shorten the loop approach... It looks pretty well golfed! Also, it's not often I see zero-based indexing in Octave :)
Stewie Griffin

@StewieGriffin Since the recursion has some offsets it doesn't really matter if you pick zero- or one indexing. I think maybe I could shave some bytes of if I did 2-indexing, but that seemed like cheating. Anyway, your intuition was right - somehow, this was indeed shorter in an anonymous recursive way. I think the main advantage is that it handles creating the four initial values very well because it just returns 1 for n<4.
Sanchises

1
@StewieGriffin Of course, good old matrix multiplication. Well done!
Sanchises


Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.