Triângulo de Pascal (espécie de)

24

Quase todo mundo aqui está familiarizado com o Triângulo de Pascal. É formado por linhas sucessivas, onde cada elemento é a soma de seus dois vizinhos superior esquerdo e superior direito. Aqui estão as primeiras 5linhas (emprestadas do triângulo de Generate Pascal ):

    1
   1 1
  1 2 1
 1 3 3 1
1 4 6 4 1
  . . .

Recolher essas linhas à esquerda

1
1 1
1 2 1
1 3 3 1
1 4 6 4 1
. . .

Classifique-os em ordem crescente

1
1 1
1 1 2
1 1 3 3
1 1 4 4 6
. . .

Leia este triângulo por linhas

[1, 1, 1, 1, 1, 2, 1, 1, 3, 3, 1, 1, 4, 4, 6 ...]

Dada uma entrada n, imprima o nnúmero th nesta série. Este é o OEIS 107430 .

Regras

  • Você pode escolher a indexação com base em 0 ou 1. Indique qual em sua submissão.
  • Supõe-se que a entrada e a saída se encaixem no tipo inteiro nativo do seu idioma.
  • A entrada e saída podem ser fornecidas por qualquer método conveniente .
  • Um programa completo ou uma função são aceitáveis. Se uma função, você pode retornar a saída em vez de imprimi-la.
  • As brechas padrão são proibidas.
  • Isso é portanto todas as regras usuais de golfe se aplicam e o código mais curto (em bytes) vence.
AdmBorkBork
fonte
6
Título muito bom!
Luis Mendo
11
Conforme o link OEIS, a única alteração necessária para gerar essa sequência em vez de um coeficiente binomial é uma divisão inteira. Isso certamente cai em "trivial".
22418 Peter Peter
5
@ PeterTaylor Isso não parece um óbvio para mim. Existem muitas outras abordagens possíveis que podem levar a oportunidades interessantes de golfe, especialmente para idiomas que não possuem um binômio embutido.
Arnauld
4
@ PeterTaylor Também não estou convencido de que seja uma duplicata. Até agora, as respostas MATL, JavaScript e Pascal são bastante diferentes entre os dois desafios. No entanto, como meu voto está aberto, não vou votar ainda.
AdmBorkBork
4
Concordo totalmente com @AdmBorkBork. Então conte comigo como reabrir voto. Isso faz 3 agora. Quantos votos são necessários para reabrir?
Luis Mendo

Respostas:

9

JavaScript (ES6), 79 bytes

Indexado a 0.

f=(n,a=[L=1])=>a[n]||f(n-L,[...a.map((v,i)=>k=(x=v)+~~a[i-1-i%2]),L++&1?k:2*x])

Demo

Quão?

f = (                       // f = recursive function taking:
  n,                        //   n = target index
  a = [L = 1]               //   a[] = current row, L = length of current row
) =>                        //
  a[n] ||                   // if a[n] exists, stop recursion and return it
  f(                        // otherwise, do a recursive call to f() with:
    n - L,                  //   n minus the length of the current row
    [                       //   an array consisting of:
      ...a.map((v, i) =>    //     replace each entry v at position i in a[] with:
        k =                 //       a new entry k defined as:
        (x = v) +           //       v +
        ~~a[i - 1 - i % 2]  //       either the last or penultimate entry
      ),                    //     end of map()
      L++ & 1 ?             //     increment L; if L was odd:
        k                   //       append the last updated entry
      :                     //     else:
        2 * x               //       append twice the last original entry
    ]                       //   end of array update
  )                         // end of recursive call

Esse algoritmo gera diretamente as linhas classificadas do triângulo de Pascal. Ele atualiza n de acordo com o comprimento da linha anterior até que exista um [n] . Por exemplo, 6 iterações são necessárias para n = 19 :

 L | n  | a[]
---+----+------------------------
 1 | 19 | [ 1 ]
 2 | 18 | [ 1, 1 ]
 3 | 16 | [ 1, 1, 2 ]
 4 | 13 | [ 1, 1, 3, 3 ]
 5 |  9 | [ 1, 1, 4, 4, 6 ]
 6 |  4 | [ 1, 1, 5, 5, 10, 10 ]
                        ^^
Arnauld
fonte
Bom trabalho. Eu não tenho certeza se eu entendo exatamente como isso funciona. Minha tentativa acabou sendo muito mais longa que a sua.
Kamoroso94
@ kamoroso94 Adicionei uma explicação.
Arnauld 07/02
Eu amo isto! Realmente gostei de descobrir o que estava fazendo.
Shaggy
6

Oitava , 46 bytes

@(n)(M=sort(spdiags(flip(pascal(n)))))(~~M)(n)

Baseado em 1.

Experimente online!

Explicação

Considere n=4como um exemplo.

pascal(n) dá uma matriz Pascal:

 1     1     1     1
 1     2     3     4
 1     3     6    10
 1     4    10    20

As linhas do triângulo Pascal são as antidiagonais desta matriz. Por isso, é invertido verticalmente usandoflip(···)

 1     4    10    20
 1     3     6    10
 1     2     3     4
 1     1     1     1

que transforma antidiagonais em diagonais.

spdiags(···) extrai as diagonais (diferentes de zero), começando na parte inferior esquerda, e as organiza como colunas preenchidas com zero:

 1     1     1     1     0     0     0
 0     1     2     3     4     0     0
 0     0     1     3     6    10     0
 0     0     0     1     4    10    20

M=sort(···)classifica cada coluna desta matriz e atribui o resultado à variável M:

 0     0     0     1     0     0     0
 0     0     1     1     4     0     0
 0     1     1     3     4    10     0
 1     1     2     3     6    10    20

A indexação lógica (···)(~~M)agora é usada para extrair os não-zeros desta matriz na ordem das colunas principais (abaixo e depois). O resultado é um vetor de coluna:

 1
 1
 1
 1
···
10
10
20

Finalmente, a n-ésima entrada desse vetor é extraída usando (···)(n), o que neste caso fornece 1.

Luis Mendo
fonte
5

Python 2 , 86 78 72 bytes

-8 bytes graças a Rod

g=lambda n,r=[1]:r[n:]and r[n/2]or g(n-len(r),map(sum,zip([0]+r,r+[0])))

Experimente online!

Ungolfed

def g(n, row=[1]):
  if n < len(row):
    return row[n/2]
  else:
    next_row = map(sum, zip([0] + row, row + [0]))
    return g(n - len(row), next_row)

Experimente online!

A função calcula recursivamente a linha do triângulo de Pascal. Dada a linha atual como row, map(sum, zip([0] + row, row + [0])).
Em cada chamada né reduzido pelo comprimento da linha atual. Se a função chegar à linha direita, o nthnúmero mais baixo da linha deve ser retornado.
Como a primeira metade de uma linha está em ordem crescente e cada linha é simétrica, o número está no índice n/2(indexado a 0, divisão inteira).

ovs
fonte
4

Wolfram Language (Mathematica) , 55 bytes

A indexação é baseada em 1.

(##&@@@Sort/@Table[n~Binomial~k,{n,0,#},{k,0,n}])[[#]]&

Experimente online!

Explicação

Provavelmente é possível jogar golfe, não sou um usuário muito experiente do Mathematica.

Table[n~Binomial~k,{n,0,#},{k,0,n}]

Para cada n ∈ [0, Entrada] ℤ ℤ , gere a tabela de binômios com cada k ∈ [0, n] ∩ ℤ .

Sort/@

Classifique cada um. Usa uma abreviação para Map[function,object]- function/@object.

(##&@@@...)[[#]]

Achate a lista resultante e recupere o elemento cujo índice na lista é a entrada.

Mr. Xcoder
fonte
3

APL (Dyalog) , 26 25 bytes

1 byte salvo graças a @ngn

{⍵⊃0~⍨∊(⍋⊃¨⊂)¨↓⍉∘.!⍨⍳1+⍵}

Experimente online!

Uriel
fonte
{⍵[⍋⍵]}->(⍋⊃¨⊂)
ngn 07/02
3

R , 58 bytes

function(n)(m=apply(outer(0:n,0:n,choose),1,sort))[m>0][n]

Experimente online!

Calcula n choose kpara cada n,kem [0,1,...,n]como uma matriz, classifica as linhas ascendente (*), e remove os zeros, em seguida, selecciona o nelemento th.

(*) Isso também os transforma em colunas, mas é melhor, pois R armazena uma matriz como um vetor coluna a coluna , o que nos permite indexar diretamente nela, preservando a ordem.

Giuseppe
fonte
3

Haskell , 143 132 125 123 bytes

((p>>=s.h)!!)
p=[1]:map(\r->zipWith(+)(0:r)(r++[0]))p
h r=splitAt(div(length r)2)r
s(a,b)=reverse b!a
(h:t)!b=h:(b!t)
x!_=x

A primeira linha é uma função sem ponto que pega um índice (com base em 0) e retorna o número apropriado na sequência. Experimente online!

Este é o meu primeiro programa Haskell! Tenho certeza que pode ficar muito mais curto. Dicas são apreciadas.

Economizou 2 bytes graças a nimi

Ungolfed

pascalRows = [1] : map (\row -> zipWith (+) (0:row) (row++[0])) pascalRows
halves row = splitAt (div (length row) 2) row
joinSorted (first, second) = interleave (reverse second) first
interleave [] _ = []
interleave longer shorter = (head longer) : (interleave shorter (tail longer))
f n = (concatMap (joinSorted.halves) pascalRows) !! n
DLosc
fonte
Você ainda tem uma ifunção s, que foi renomeada para !, eu acho. Se você usar uma função infix você pode soltar a ()volta reverse b: s(a,b)=reverse b!a.
nimi
@ nimi Ah, obrigado - eu mudei no TIO, mas perdi um lugar no código aqui. E obrigado pela dica de parênteses.
7608 DLosc
3

JavaScript, 57 bytes

f=(i,r=1)=>i<r?i>1?f(i-2,--r)+f(i<r?i:r-1,r):1:f(i-r,r+1)

Indexado a 0.

Como isso acontece:

Etapa 0:

c=(i,r)=>i?r&&c(i-1,r-1)+c(i,r-1):1
f=(i,r=1)=>i<r?c(i>>1,r-1):f(i-r,r+1)

Este código é fácil de entender:

  • A função ccalcula a fórmula de combinação de uso: C (n, k) = C (n-1, k) + C (n-1, k-1); ou 1 se k == 0 ou k == n
  • função ftente descobrir o número e o índice da linha e chame a função c para obter o resultado.

Passo 1:

c=(i,r)=>i>1?--r&&c(i-2,r)+c(i,r):1
f=(i,r=1)=>i<r?c(i,r):f(i-r,r+1)

Nesta etapa, tentamos modificar a chamada da função cpara c(i,r)que ela seja igual ao parâmetro de f.

Passo 2:

c=(i,r)=>i>1?--r&&c(i-2,r)+c(i<r?i:r-1,r):1
f=(i,r=1)=>i<r?c(i,r):f(i-r,r+1)

Testamos i<rse estamos usando função fou função c. É por isso que nós pensamos em manter o i<rcontrole durante a recursão da função c.

Etapa 3:

f=(i,r=1)=>i<r?i>1?--r&&f(i-2,r)+f(i<r?i:r-1,r):1:f(i-r,r+1)

Nesta etapa, mesclamos essas duas funções em uma.

Depois de mais um pouco de golfe, finalmente conseguimos a resposta descrita acima.

tsh
fonte
2

Gelatina , 13 bytes

0rcþ`ZṢ€Ẏḟ0⁸ị

Experimente online!

Usando o algoritmo Dyalog de Uriel.

1 indexado.

Explicação:

0rcþ`ZṢ€Ẏḟ0⁸ị
0r            Return inclusive range from 0 to n
    `         Call this dyad with this argument on both sides
   þ           Outer product with this dyad
  c             Binomial coefficient
     Z        Zip
       €      Call this link on each element
      Ṣ        Sort
        Ẏ     Concatenate elements
         ḟ0   Remove 0s
           ⁸ị Take the nth element
Erik, o Outgolfer
fonte
Você poderia adicionar uma explicação? Não consigo descobrir o que þestá fazendo aqui.
Shaggy
11
@ Shagy É produto externo, vou adicionar uma explicação.
Erik the Outgolfer
2

JavaScript (Node.js) , 65 bytes

Nem mesmo uma matriz é usada. Indexado a 0.

f=(n,i=0,g=x=>x?x*g(x-1):1)=>n>i?f(n-++i,i):g(i)/g(c=n>>1)/g(i-c)

Experimente online!

Explicação:

f=(n,i=0,                 )=>                                     // Main Function
         g=x=>x?x*g(x-1):1                                        // Helper (Factorial)
                             n>i?                                 // Is n > i?
                                 f(n-++i,i):                      // If so, call function
                                                                  // f(n-i-1, i+1) to skip
                                                                  // . i+1 terms
                                            g(i)/g(c=n>>1)/g(i-c) // If not, since sorting 
                                                                  // . the binomial coeffs
                                                                  // . equals to writing
                                                                  // . the first floor(i/2)
                                                                  // . coefficients twice
                                                                  // . each, so a shortcut
Shieru Asakoto
fonte
1

Pascal , 373 bytes

function t(n,k,r:integer):integer;begin if(n<k)then t:=r-1 else t:=t(n,k+r,r+1)end;
function s(n,k:integer):integer;begin if(k=0)then s:=n else s:=s(n+k,k-1)end;
function f(n,k:integer):integer;begin if((k<1)or(k>n))then f:=0 else if n=1 then f:=1 else f:=f(n-1,k-1)+f(n-1,k)end;
function g(n:integer):integer;var k:integer;begin k:=t(n,0,1);g:=f(k,(n-s(0,k-1)+2)div 2)end;

g é a função

Experimente online!

Uriel
fonte
n=1 thenpode ser n=1then.
Jonathan Frech
Da mesma forma, parece que if(k=0)thenpode se tornar if k=0then.
Salsicha
se algum número sempre for maior que 0, você deve usar em wordvez de integer.
tsh
1

Java 8, 187 bytes

n->{int r=~-(int)Math.sqrt(8*n+1)/2+1,a[]=new int[r],k=r,x=0;for(;k-->0;a[k]=p(r,k))x+=k;java.util.Arrays.sort(a);return a[n-x];}int p(int r,int k){return--r<1|k<2|k>r?1:p(r,k-1)+p(r,k);}

Explicação:

Experimente online.

n->{                   // Method with integer as both parameter and return-type
  int r=~-(int)Math.sqrt(8*n+1)/2+1,
                       //  Calculate the 1-indexed row based on the input
      a[]=new int[r],  //  Create an array with items equal to the current row
      k=r,             //  Index integer
      x=0;             //  Correction integer
  for(;k-->0;          //  Loop down to 0
    a[k]=p(r,k))       //   Fill the array with the Pascal's Triangle numbers of the row
    x+=k;              //   Create the correction integer
  java.util.Arrays.sort(a);
                       //  Sort the array
  return a[n-x];}      //  Return the `n-x`'th (0-indexed) item in this sorted array

// Separated recursive method to get the k'th value of the r'th row in the Pascal Triangle
int p(int r,int k){return--r<1|k<2|k>r?1:p(r,k-1)+p(r,k);}
Kevin Cruijssen
fonte
1

MATL , 11 bytes

:qt!XnSXzG)

Baseado em 1.

Experimente online! Ou verifique todos os casos de teste .

Explicação

Considere a entrada 4como um exemplo. ;é o separador de linhas para matrizes ou vetores de coluna.

:     % Implicit input: n. Push the row vector [1 2 ... n]          
      S STACK: [1 2 3 4]
q     % Subtract 1, emlement-wise: gives [0 1 ... n-1]
      % STACK: [0 1 2 3]
t!    % Duplicate and transpose into a column vector
      % STACK: [0 1 2 3], [0; 1; 2; 3]
Xn    % Binomial coefficient, element-wise with broadcast. Gives an
      % n×n matrix where entry (i,j) is binomial(i,j), or 0 for i<j
      % STACK: [1 1 1 1;
                0 1 2 3;
                0 0 1 3;
                0 0 0 1]
S     % Sort each column
      % STACK: [0 0 0 1;
      %         0 0 1 1;
      %         0 1 1 3;
      %         1 1 2 3]
Xz    % Keep only nonzeros. Gives a column vector
      % STACK: [1; 1; 1; 1; 1; 2; 1; 1; 3; 3]
G)    % Get the n-th element. Implicitly display
      % STACK: 1
Luis Mendo
fonte
1

Lote, 128 bytes

@set/as=2,t=r=m=i=1
:l
@if %1 geq %t% set/as+=r,t+=r+=1&goto l
@for /l %%i in (%s%,2,%1)do @set/ar-=1,m=m*r/i,i+=1
@echo %m%

Indexado a 0.

Neil
fonte
Você pode adicionar uma explicação, por favor? Não consigo seguir a lógica aqui.
AdmBorkBork
@AdmBorkBork As três primeiras linhas calculam a linha re a coluna %1-(s-2)do dia %1da série. A quarta linha usa isso para calcular o coeficiente binomial (n k) = n!/(n-k)!k!= n(n-1)...(n+1-k)/(1)(2)...k= (n/1)((n-1)/2)...((n+1-k)/k). Onde está o MathJax quando eu precisar?
Neil
1

APL (Dyalog Classic) , 17 bytes

⎕⊃∊i!⍨,\⌊.5×i←⍳99

Experimente online!

Indexação baseada em 0

observe que (49!98) > 2*53, ou seja, o coeficiente binomial 98 acima de 49 é maior que 2 53 , portanto, nesse momento, o Dyalog já começou a perder precisão devido ao ponto flutuante IEEE

ngn
fonte
@Abigail veja aqui e aqui
ngn
1

05AB1E , 10 bytes

Indexado a 0

ÝεDÝc{}˜sè

Experimente online!

Explicação

Ý             # push range [0 ... input]
 ε    }       # apply to each element
  DÝc         # N choose [0 ... N]
     {        # sort
       ˜      # flatten result to a list
        sè    # get the element at index <input>
Emigna
fonte
1

Gelatina , 11 bytes

Ḷc€`Ṣ€Fḟ0ị@

Experimente online!

Um link monádico que pega o índice e retorna um número inteiro - usa indexação baseada em 1.

Quão?

Realiza o desafio da mesma maneira que está escrito, com mais da direita do triângulo de Pascal (zeros), que é então jogado fora ...

Ḷc€`Ṣ€Fḟ0ị@ - Link: integer, i    e.g. 1   or    9
Ḷ           - lowered range            [0]       [0,1,2,3,4,5,6,7,8]
   `        - repeat left as right arg [0]       [0,1,2,3,4,5,6,7,8]
 c€         - binomial choice for €ach [[1]]     [[1,0,0,0,0,0,0,0,0],[1,1,0,0,0,0,0,0,0],[1,2,1,0,0,0,0,0,0],[1,3,3,1,0,0,0,0,0],[1,4,6,4,1,0,0,0,0],[1,5,10,10,5,1,0,0,0],[1,6,15,20,15,6,1,0,0],[1,7,21,35,35,21,7,1,0],[1,8,28,56,70,56,28,8,1]]
    Ṣ€      - sort €ach                [[1]]     [[0,0,0,0,0,0,0,0,1],[0,0,0,0,0,0,0,1,1],[0,0,0,0,0,0,1,1,2],[0,0,0,0,0,1,1,3,3],[0,0,0,0,1,1,4,4,6],[0,0,0,1,1,5,5,10,10],[0,0,1,1,6,6,15,15,20],[0,1,1,7,7,21,21,35,35],[1,1,8,8,28,28,56,56,70]]
      F     - flatten                  [1]       [0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,1,1,0,0,0,0,0,0,1,1,2,0,0,0,0,0,1,1,3,3,0,0,0,0,1,1,4,4,6,0,0,0,1,1,5,5,10,10,0,0,1,1,6,6,15,15,20,0,1,1,7,7,21,21,35,35,1,1,8,8,28,28,56,56,70]
       ḟ0   - filter discard zeros     [1]       [1,1,1,1,1,2,1,1,3,3,1,1,4,4,6,1,1,5,5,111,1,6,6,15,15,21,1,7,7,21,21,35,35,1,1,8,8,28,28,56,56,70]
         ị@ - index into (sw@p args)    1         3 --------------^
Jonathan Allan
fonte
1

Vermelho , 206 bytes

f: func[n][t: copy[[1]]l: 0
while[l < n][a: copy last t insert append a 0 0 b: copy[]repeat i k:(length? a)- 1[append b a/(i) + a/(i + 1)]append t reduce[b]l: l + k]foreach p t[sort p]pick split form t{ }n]

Baseado em 1

Experimente online!

Explicação:

f: func [n] [
    t: copy [[1]]                       ; start with a list with one sublist [1]
    l: 0                                ; there are 1 items overall
    while [l < n] [                     ; while the number of items is less than the argument
        a: copy last t                  ; take the last sublist 
        insert append a 0 0             ; prepend and append 0 to it  
        b: copy []                      ; prepare a list for the sums  
        repeat i k: (length? a) - 1 [   ; loop throught the elements of the list
            append b a/(i) + a/(i + 1)] ; and find the sum of the adjacent items
        append t reduce [b]             ; append the resulting list to the total list
        l: l + k                        ; update the number of the items
    ]
    foreach p t [sort p]                ; sort each sublist
    v: pick split form t { } n          ; flatten the list and take the n-th element
]
Galen Ivanov
fonte
1

Perl, 48 bytes

Inclui +1parap

perl -pe '$_-=$%until$_<++$%;$./=$_/--$%for 1..$_/2;$_=$.' <<< 19

Usa a indexação base 0.

Ton Hospel
fonte
1

J, 46 bytes

f=:](([-2!]){/:~@(i.!<:)@])[:<.2&!@,:^:_1

Indexado a 0

Experimente online!

Notas:

  • <.2&!@,:^:_1fornece o número da linha relevante do triângulo de Pascal arredondando o inverso de y choose 2.
  • /:~@(i.!<:)@] calcula a linha e a classifica.
  • [-2!] fornece o índice para a linha.
Kyle Miller
fonte
Olá. Bem vindo ao site! Esta é uma primeira resposta nice :)
DJMcMayhem
1

Julia , 70 bytes

f(x)=map(n->binomial(n-1,ceil(Int,x/2-(n^2-n)/4-1)),round(Int,√(x*2)))

Baseado em 1

Explicação:

primeiro encontre o número da linha, depois o número da coluna e calcule o binômio

Jimmy Chen
fonte
Bem-vindo ao PPCG!
Martin Ender
sim rosto thx feliz
Jimmy Chen
0

Pitão, 15 bytes

@u+GSm.cHdhHhQY

Indexado a 0

Tente

Explicação

@u+GSm.cHdhHhQY
 u          hQY   Reduce on [0, ..., input], starting with the empty list...
  +G              ... append to the accumulator...
    Sm.cHdhH      ... the sorted binomial coefficients.
@              Q  Take the 0-indexed element.

fonte
0

Limpo , 80 bytes

import StdEnv

\n=flatten[sort[prod[j+1..i]/prod[1..i-j]\\j<-[0..i]]\\i<-[0..]]!!n

Experimente online!

Como uma função lambda.

Furioso
fonte
0

Ruby , 56 bytes

->n{a=0;n-=a until n<a+=1;[*2..a].combination(n/2).size}

Baseado em 0

Primeiro, obtenha a linha e a coluna no triângulo e depois calcule o coeficiente binomial correspondente a essa posição.

Experimente online!

GB
fonte
0

Na realidade , 8 bytes

Em grande parte baseado na resposta de Jonathan Allan's Jelly . Usa indexação 0.

;r♂╣♂SΣE

Experimente online!

Ungolfing

          Implicit input n.
;         Duplicate n.
 r        Lowered range. [0..n-1].
  ♂╣      Pascal's triangle row of every number.
    ♂S    Sort every row.
      Σ   Sum each row into one array.
       E  Get the n-th element of the array (0-indexed).
          Implicit return.
Sherlock9
fonte
É suposto produzir um único número; o dia nda série. Isso produz uma matriz.
recursivo em
Ops. Fixo. Obrigado @recursive
Sherlock9
0

C (gcc) , 140 123 bytes

F(n){n=n?n*F(~-n):1;}f(n,j,k,c){for(c=j=0;j<n;j++)for(k=0;k<=j/2;k++)if(++c==n||j&1|k-j/2&&++c==n)return F(j)/F(k)/F(j-k);}

Experimente online!

Jonathan Frech
fonte
@ceilingcat Obrigado.
Jonathan Frech 14/09