Sumas absolutas de coeficientes polinomiales de Sidi

28

Fondo

El polinomio Sidi de grado n - o el (n + 1) th polinomio Sidi - se define de la siguiente manera.

definición de polinomios de Sidi

Los polinomios de Sidi tienen varias propiedades interesantes, pero también sus coeficientes. Este último forma la secuencia OEIS A075513 .

Tarea

Escriba un programa completo o una función que, dado un número entero no negativo n , imprima o devuelva la suma absoluta de los coeficientes del polinomio Sidi de grado n , es decir

definición de salida prevista

Estas sumas forman la secuencia OEIS A074932 .

Si prefiere la indexación basada en 1, se puede tomar un número entero positivo n lugar y calcular la suma absoluta de los coeficientes de la n º Sidi polinomio.

Como se trata de , debe hacer su código lo más corto posible. Se aplican todas las reglas estándar.

Casos de prueba (basados ​​en 0)

 n           Σ

 0           1
 1           3
 2          18
 3         170
 4        2200
 5       36232
 6      725200
 7    17095248
 8   463936896
 9 14246942336

Casos de prueba (basados ​​en 1)

 n           Σ

 1           1
 2           3
 3          18
 4         170
 5        2200
 6       36232
 7      725200
 8    17095248
 9   463936896
10 14246942336
Dennis
fuente

Respuestas:

46

Python 2 , 43 bytes

f=lambda n,k=1:k/n or n*f(n,k+1)+k*f(n-1,k)

Pruébalo en línea!

Un enfoque diferente

Desde que publiqué este desafío, traté de encontrar una solución recursiva para este problema. Si bien no pude usar nada más que lápiz y papel, logré convertir la fórmula del golf en un problema práctico, al menos para ciertas definiciones prácticas, lo que hizo que fuera más fácil de analizar.

Imagine un programa de juegos con k + m candidatos que funciona de la siguiente manera.

En la ronda 1, todos los candidatos deben realizar una tarea determinada lo más rápido posible. Los k candidatos que realizan la tarea más rápido ganan 1 k $ (un kilodólar) cada uno y avanzan a la ronda 3.

En la ronda 2, los m candidatos restantes tienen una segunda oportunidad de unirse a los otros k . A cada candidato se le hace una pregunta. Si responden la pregunta correctamente, ganan 1 k $ y avanzan a la ronda 3. Sin embargo, si no responden la pregunta, son eliminados del juego. Esto significa que la ronda 3 tendrá entre k y k + m candidatos, dependiendo de cuántos puedan responder sus preguntas.

La ronda 3 consiste en m más concursos que son similares a la ronda 1. En cada competencia, los participantes tienen que cumplir una determinada tarea. A diferencia de la ronda 1, solo un candidato recibe un premio, pero todos los candidatos pueden participar en el próximo concurso. Cada concurso paga el doble que el concurso anterior; el primero paga 2 k $ y el último 2 m k $ .

Tenga en cuenta que, dado que todos los premios son poderes de dos, saber cuánto dinero en premios ganó un candidato significa que sabemos si avanzaron a la ronda 3 y cuáles de los concursos de la ronda 3 ganaron.

Suponga que está viendo el programa de juegos y que la ronda 1 ya terminó, por lo que sabe qué k candidatos ya llegaron a la ronda 3 y cuáles m candidatos aún están atrapados en la ronda 2. ¿De cuántas maneras se puede distribuir el dinero restante del premio?

Una vez que sabemos cuáles de los candidatos m de la segunda ronda han avanzado a la ronda 3, es fácil calcular los posibles resultados para este escenario específico. Si j candidatos avanzan, hay k + j total de candidatos en la ronda 3 y, por lo tanto, k + j posibles resultados para cada concurso. Con m concursos individuales en la ronda 3, esto produce (k + j) m resultados para todos los m concursos.

Ahora, j puede tomar cualquier valor entre 0 y m , dependiendo de qué candidatos responden correctamente en la ronda 2. Para cada valor fijo de j , hay m C j diferentes combinaciones de j candidatos que podrían haber avanzado a la ronda 3. Si llamamos el número total de resultados posibles para k candidatos de la ronda 3 ym candidatos de la ronda 2 g (m, k) , obtenemos la siguiente fórmula.

fórmula para g

Si arreglamos k = 1 , obtenemos el siguiente caso especial, que constituye nuestro nuevo enfoque para resolver el problema original.

relación entre Sigma yg

Una fórmula recursiva

Ahora, suponga que se quedó dormido durante los comerciales después de la ronda 1, y se despertó justo a tiempo para ver quién ganó el último concurso de la ronda 3 y, por lo tanto, el gran premio de 2 m k $ . No tiene ninguna otra información, ni siquiera cuánto dinero en premios ganó ese candidato en total. ¿De cuántas maneras se puede distribuir el dinero restante del premio?

Si el ganador fue uno de los m candidatos de la ronda 2, ya sabemos que debe haber avanzado a la ronda 3 . Por lo tanto, efectivamente tenemos k + 1 candidatos en la ronda 3, pero solo m - 1 candidatos en la ronda 2. Dado que conocemos al ganador del último concurso, solo hay m - 1 con resultados inciertos, por lo que hay g (m - 1, k + 1) posibles resultados.

Si el ganador es uno de los k candidatos que se saltó la ronda 2, el cálculo se vuelve un poco más complicado. Como antes, solo quedan m - 1 rondas, pero ahora todavía tenemos k candidatos en la ronda 3 ym candidatos en la ronda 2. Dado que el número de candidatos de la ronda 2 y el número de concursos de la ronda 3 son diferentes, los posibles resultados no pueden se calculará con una simple invocación de g . Sin embargo, después de que el candidato de la primera ronda 2 ha respondido, correcta o incorrectamente, el número de candidatos de la segunda ronda coincide una vez más con los concursos m - 1 de la ronda 3. Si el candidato avanza, hay k + 1 ronda 3 candidatos y por lo tanto g (m - 1, k + 1)posibles resultados; Si el candidato es eliminado, el número de candidatos de la ronda 3 permanece en k y hay g (m - 1, k) posibles resultados. Dado que el candidato avanza o no, hay g (m - 1, k + 1) + g (m - 1, k) posibles resultados que combinan estos dos casos.

Ahora, si sumamos los resultados potenciales para todos los candidatos k + m que podrían haber ganado el gran premio, el resultado debe coincidir con g (m, k) . Hay m concursantes de la ronda 2 que conducen a g (m - 1, k + 1) resultados potenciales cada uno, yk concursantes de la ronda 3 que conducen a g (m - 1, k + 1) + g (m - 1, k) unos. En resumen, obtenemos la siguiente identidad.

fórmula recursiva para g

Junto con el caso base

caso base para g

Estas dos fórmulas caracterizan completamente la función g .

Una implementación de golf

Mientras

g=lambda m,k=1:0**m or(m+k)*g(m-1,k+1)+k*g(m-1,k)

(49 bytes, 0**mproduce 1 una vez que m cae a 0 ) o incluso

g=lambda m,k=1:m<1 or(m+k)*g(m-1,k+1)+k*g(m-1,k)

(48 bytes, devuelve True en lugar de 1 ) serían soluciones válidas, todavía hay bytes para guardar.

Si definimos una función f que toma el número n de candidatos de la ronda 1 en lugar del número m de candidatos de la ronda 2 como primer argumento, es decir,

definición de f en términos de g

obtenemos la fórmula recursiva

fórmula recursiva para f

con estuche base

caso base para f

Finalmente tenemos

relación entre Sigma y f

entonces la implementación de Python

f=lambda n,k=1:k/n or n*f(n,k+1)+k*f(n-1,k)

( k/nproduce 1 una vez que n = k ) resuelve la tarea en cuestión con indexación basada en 1.

Dennis
fuente
8

Mathematica, 33 32 bytes

Guardado un byte gracias a JungHwan Min .

Binomial[#,r=0~Range~#].(r+1)^#&
alephalpha
fuente
4

Pyth - 13 12 bytes

Solo implementa la fórmula sin la (-1)^kparte.

sm*.cQd^hdQh

Test Suite .

Maltysen
fuente
3

MATL , 12 bytes

t:XnG:QG^*sQ

La entrada está basada en 0.

Pruébalo en línea!

Explicación

Considere la entrada 5como un ejemplo.

t      % Take n implicitly. Duplicate
       % STACK: 5, 5
:      % Range [1 2 ...n]
       % STACK: 5, [1 2 3 4 5]
Xn     % N-choose-k, vectorized
       % STACK: [5 10 10 5 1]
G:Q    % Push [2 3 ... n+1]
       % STACK: [5 10 10 5 1], [2 3 4 5 6]
G^     % Raise to n
       % STACK: [5 10 10 5 1], [32 243 1024 3125 7776]
*      % Multiply, element-wise
       % STACK: [160 2430 10240 15625 7776]
s      % Sum of array
       % STACK: 36231
Q      % Add 1. Display implicitly
       % STACK: 36232
Luis Mendo
fuente
2

R, 36 bytes

sum(choose(n<-scan(),0:n)*(0:n+1)^n)

La vectorización de R es útil aquí cuando se aplica la fórmula.

Billywob
fuente
2

J , 19 bytes

+/@((!{:)*>:^{:)@i.

Utiliza indexación basada en uno.

Pruébalo en línea!

Explicación

+/@((!{:)*>:^{:)@i.  Input: integer n
                 i.  Range [0, 1, ..., n-1]
   (           )@    Operate on that range
             {:        Get the last value, n-1
          >:           Increment, range becomes [1, 2, ..., n]
            ^          Exponentiate. [1^(n-1), 2^(n-1), ..., n^(n-1)]
    ( {:)              Get the last value, n-1
     !                 Binomial coefficient. [C(n-1, 0), C(n-1, 1), ..., C(n-1, n-1)]
         *             Multiply
+/@                  Reduce by addition
millas
fuente
1

Máxima, 39 bytes

f(n):=sum(binomial(n,k)*(k+1)^n,k,0,n);
rahnema1
fuente
1

PARI / GP, 35 bytes

n->sum(k=0,n,binomial(n,k)*(k+1)^n)
alephalpha
fuente
0

Axioma, 39 bytes

f(n)==sum(binomial(n,i)*(i+1)^n,i=0..n)

código de prueba y resultados

(35) -> [[i,f(i)] for i in 0..9]
   (35)
   [[0,1], [1,3], [2,18], [3,170], [4,2200], [5,36232], [6,725200],
    [7,17095248], [8,463936896], [9,14246942336]]
RosLuP
fuente
0

Jalea , 9 bytes

cR×R‘*ƊS‘

Pruébalo en línea!

Cómo funciona

cR×R‘*ƊS‘ - Main link. Argument: n (integer)        e.g.   5
 R        - Range from 1 to n                              [1, 2, 3, 4, 5]
c         - Binomial coefficient                           [5, 10, 10, 5, 1]
      Ɗ   - Last three links as a monad:
   R      -   Link 1: Range from 1 to n                    [1, 2, 3, 4, 5]
    ‘     -   Link 2: Increment                            [2, 3, 4, 5, 6]
     *    -   Link 3: To the power of n                    [32, 243, 1024, 3125, 7776]
  ×       - Multiply, pairwise                             [160, 2430, 10240, 15625, 7776]
       S  - Sum                                            36231
        ‘ - Increment                                      36232

fuente