{"id":"n-5ddd616bf530de21","revision":0,"first_index":0,"fingerprint":"f6-856ad640314d7fe1","name":"theta_plus_one_squared of T(n, k) = Sum_{j=0..n-k} (-1)^j*binomial(n - k + 1, j)*(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.","polytool_family":null,"interpretation":null,"rows_observed":15,"degree_pattern":"affine","degree_slope":1,"degree_intercept":0,"palindromic":0,"gamma_positive":0,"nonnegative":1,"rr_by_factors":false,"rows":[["1"],["2","4"],["6","24","9"],["24","144","126","16"],["120","960","1350","480","25"],["720","7200","14040","8640","1550","36"],["5040","60480","151200","134400","45150","4536","49"],["40320","564480","1723680","2016000","1020600","208656","12446","64"],["362880","5806080","20956320","30481920","20853000","6713280","889350","32640","81"],["3628800","65318400","272160000","474163200","410886000","183708000","40107480","3582720","82782","100"],["39916800","798336000","3772137600","7664025600","8062362000","4652303040","1450772400","223872000","13851486","204600","121"],["479001600","10538035200","55683936000","129330432000","160299216000","113834730240","46698442560","10593792000","1188757944","51915600","495374","144"],["6227020800","149448499200","873339667200","2283240960000","3261402144000","2757324810240","1411481111040","429526056960","72982629720","6078072000","189939750","1179360","169"],["87178291200","2266635571200","14515185484800","42194292940800","68289661440000","67083695078400","41316052377600","15911652556800","3699609273360","480900420000","30177627480","681670080","2768558","196"],["1307674368000","36614882304000","254996501760000","815988805632000","1476400685760000","1654208075520000","1193619736108800","559003808563200","166864561783920","30289915656000","3060335678400","146382163200","2408352414","6422136","225"]],"degrees":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14],"row_sums":["1","6","39","310","2935","32186","400855","5586246","86095551","1453637482","26677781647","528698754662","11251242090439","255875725500474","6192503780072295"],"nonn":true,"definitions":[{"position":0,"source":"derived:theta_plus_one_squared","source_kind":"derived","equality_edge_id":null,"recurrence":{"first_index":0,"initial_polynomials":[["1"],["2","4"]],"recurrence":{"denominator":{"coeffs":[["0"],["0"],["1"]]},"terms":[{"coeff":{"coeffs":[["0","4","0"],["0","0","0"],["-1","0","-7"],["-2","0","-2"],["-1","-4","9"]]},"deriv_order":0,"offset":2,"sign":"none"},{"coeff":{"coeffs":[["0","1","5","0"],["0","0","0","0"],["0","1","-4","-13"],["0","2","4","-14"],["0","0","4","4"]]},"deriv_order":1,"offset":2,"sign":"none"},{"coeff":{"coeffs":[["0","0","1","1","0"],["0","0","0","0","0"],["0","0","-1","-4","-3"],["0","0","0","-4","-4"]]},"deriv_order":2,"offset":2,"sign":"none"},{"coeff":{"coeffs":[["0","0"],["1","0"],["3","0"],["2","4"]]},"deriv_order":0,"offset":1,"sign":"none"},{"coeff":{"coeffs":[["0","1","0"],["0","0","0"],["0","-2","-4"],["0","0","-4"]]},"deriv_order":1,"offset":1,"sign":"none"}]},"schema":"polytool.recurrence.v1"}}],"properties":[{"id":182650,"position":0,"claim":"real_rooted","from_index":0,"region_key":"(-inf, 0]","provenance":"edges","claim_data":{"claim":"real_rooted","from_index":0,"region":{"upper":{"closed":true,"value":"0"}}},"label":"Real-rooted","description":"Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]","region_exact":null,"edges":["e-cbf3127edaff7909"]},{"id":182651,"position":1,"claim":"one_signed_coefficients","from_index":0,"region_key":null,"provenance":"derived_from_edges","claim_data":{"alternating":false,"claim":"one_signed_coefficients","from_index":0},"label":"Coefficients of one sign","description":"Every nonzero row, n >= 0: the coefficients of P_n have one sign","region_exact":null,"edges":["e-cbf3127edaff7909"]},{"id":182652,"position":2,"claim":"hadamard_kernel","from_index":0,"region_key":null,"provenance":"derived_from_edges","claim_data":{"claim":"hadamard_kernel","from_index":0},"label":"Hadamard kernel (Malo)","description":"Rows P_n, n >= 0, are Hadamard kernels (preserve real-rootedness)","region_exact":null,"edges":["e-cbf3127edaff7909"]}],"region_exact":{},"incoming_edges":[{"id":"e-cbf3127edaff7909","operation":{"index_shift":0,"operation":"preserver_image","operator":"{\"shape\":\"euler_polynomial\",\"coefficients\":[{\"numerator\":[\"1\"],\"denominator\":[\"1\"]},{\"numerator\":[\"2\"],\"denominator\":[\"1\"]},{\"numerator\":[\"1\"],\"denominator\":[\"1\"]}]}"},"operation_detail":"euler_polynomial, shift 0","conclusion_id":"n-5ddd616bf530de21","premise_count":1,"theorem":"Laguerre's theorem on multiplier sequences: if f is a real polynomial whose roots are all real and <= 0, then (f(k))_(k >= 0) is a multiplier sequence, so x^k -> f(k) x^k maps every real-rooted polynomial to a real-rooted polynomial or zero (Laguerre; see T. Craven and G. Csordas, Multiplier sequences for fields, Illinois J. Math. 21 (1977); Polya--Schur, Theorem 40 of P. Branden, arXiv:1410.6601: sum_k f(k) x^k / k! = f(theta) e^x = C e^x x^m prod (1 + x/a_j), a_j > 0); f(k) has one sign for k >= 0, so coefficients of one sign stay of one sign; applied to f(n, .) for each n","lean_route":null,"lean_declarations":["RealRooted.isMultiplierSequence_eval_of_roots_nonpos","RealRooted.isPFMultiplierSequence_eval_of_roots_nonpos","RealRooted.isMultiplierSequence_iff_preserves"],"certificate_hash":"b3ab8590b401a6aae5f4794b9542fba1303d52446e66f00944c67bb463d84d34","verification_polytool":"replayed","verification_lean":"not_attempted","claims":[{"claim":"real_rooted","from_index":0,"region":{"upper":{"closed":true,"value":"0"}}}],"premises":[{"position":0,"node_id":"n-3b2b9e31bd278a34","role":"sequence","name":"T(n, k) = Sum_{j=0..n-k} (-1)^j*binomial(n - k + 1, j)*(n - k + 1 - j)^n. Triangle read by rows, T(n, k) for 1 <= k <= n.","uses":{"claim":"real_rooted","from_index":0,"region":{"upper":{"closed":true,"value":"0"}}},"uses_text":"Every row P_n, n >= 0, is zero or real-rooted with all roots in (-inf, 0]"}]}],"used_by":[],"used_by_total":0,"links":[{"id":"l-045b9c29a3b5d839","kind":"fingerprint_related","node_a_id":"n-5ddd616bf530de21","node_b_id":"n-7016994b0a4cefab","fingerprint":"f6-856ad640314d7fe1","oeis_id":null,"in_catalog":null,"target":null,"row_alignment":null,"rows_matched":null,"oeis_nonn":null,"details":{"fingerprint":"f6-856ad640314d7fe1"},"other_id":"n-7016994b0a4cefab","other":{"id":"n-7016994b0a4cefab","name":"theta_plus_one_squared of shift_x_plus_one of eulerian","first_index":1,"fingerprint":"f6-856ad640314d7fe1","polytool_family":null,"degree_pattern":"affine","palindromic":0,"gamma_positive":0,"nonnegative":1,"rr_by_factors":0,"preview_rows":[["1"],["2","4"],["6","24","9"]],"card_rows":[["1"],["2","4"],["6","24","9"],["24","144","126","16"]],"row_sums_preview":["1","6","39","310","2935","32186","400855","5586246"],"regions":["(-inf, 0]"],"region_exact":{},"oeis_ids":[],"derivation":{"root":"n-ff270bcb802a85af","steps":["x \u21a6 x+1","(\u03b8+1)\u00b2"],"edges":["e-16a5ff84df4996fb","e-578b675dca3d0df1"],"text":"A008292 \u2192 x \u21a6 x+1 \u2192 (\u03b8+1)\u00b2"},"headline":"A008292 \u2192 x \u21a6 x+1 \u2192 (\u03b8+1)\u00b2","chain":null}},{"id":"l-462b5da64379a971","kind":"fingerprint_related","node_a_id":"n-2ee019a0ed5b4470","node_b_id":"n-5ddd616bf530de21","fingerprint":"f6-856ad640314d7fe1","oeis_id":null,"in_catalog":null,"target":null,"row_alignment":null,"rows_matched":null,"oeis_nonn":null,"details":{"fingerprint":"f6-856ad640314d7fe1"},"other_id":"n-2ee019a0ed5b4470","other":{"id":"n-2ee019a0ed5b4470","name":"theta_plus_one_squared of shift_x_plus_one of Euler's triangle: triangle of Eulerian numbers T(n,k) (n>=0, 0 <= k <= n) read by rows.","first_index":1,"fingerprint":"f6-856ad640314d7fe1","polytool_family":null,"degree_pattern":"irregular","palindromic":0,"gamma_positive":0,"nonnegative":1,"rr_by_factors":0,"preview_rows":[["1"],["1"],["2","4"]],"card_rows":[["1"],["1"],["2","4"],["6","24","9"]],"row_sums_preview":["1","1","6","39","310","2935","32186","400855"],"regions":["(-inf, 0]"],"region_exact":{},"oeis_ids":[],"derivation":{"root":"n-efd1fdf07e5468d8","steps":["x \u21a6 x+1","(\u03b8+1)\u00b2"],"edges":["e-aed75c9f11c92985","e-f2325e238de872b6"],"text":"A173018 \u2192 x \u21a6 x+1 \u2192 (\u03b8+1)\u00b2"},"headline":"A173018 \u2192 x \u21a6 x+1 \u2192 (\u03b8+1)\u00b2","chain":null}}],"oeis":[],"bfiles":[],"oeis_texts":[{"oeis_id":"","target":"node","first_row":0,"oeis_offset":0,"reversed":false,"trim_zeros":"none","regular_triangle":false,"rows_shown":10,"polytool_commit":"c48629374597092cc41ef4d6273f37991516f984","triangle_example":"Triangle begins:\n  0:       1;\n  1:       2,        4;\n  2:       6,       24,         9;\n  3:      24,      144,       126,        16;\n  4:     120,      960,      1350,       480,        25;\n  5:     720,     7200,     14040,      8640,      1550,        36;\n  6:    5040,    60480,    151200,    134400,     45150,      4536,       49;\n  7:   40320,   564480,   1723680,   2016000,   1020600,    208656,    12446,      64;\n  8:  362880,  5806080,  20956320,  30481920,  20853000,   6713280,   889350,   32640,    81;\n  9: 3628800, 65318400, 272160000, 474163200, 410886000, 183708000, 40107480, 3582720, 82782, 100;\n","recurrence_sentence":"The row polynomials satisfy n^2*P(n,x) = -(n^4+2*n^3+n^2+(4*n^4-4)*x-(9*n^4-2*n^3-7*n^2)*x^2)*P(n-2,x) + ((2*n^3+n^2+1)*x+(4*n^4+4*n^3-4*n^2+5)*x^2+(4*n^4-14*n^3-13*n^2)*x^3)*d/dx P(n-2,x) - ((n^2-1)*x^2+(4*n^3+4*n^2-1)*x^3+(4*n^3+3*n^2)*x^4)*d^2/dx^2 P(n-2,x) + (2*n^3+3*n^2+n+4*n^3*x)*P(n-1,x) - ((2*n^2-1)*x+(4*n^3+4*n^2)*x^2)*d/dx P(n-1,x) for n >= 2, with P(0,x) = 1, P(1,x) = 2 + 4*x.","mathematica":"P[0, x_] := 1;\nP[1, x_] := 2 + 4*x;\nP[n_, x_] := P[n, x] = Expand[Together[(-(n^4 + 2*n^3 + n^2 + (4*n^4 - 4)*x - (9*n^4 - 2*n^3 - 7*n^2)*x^2)*P[n - 2, x] + ((2*n^3 + n^2 + 1)*x + (4*n^4 + 4*n^3 - 4*n^2 + 5)*x^2 + (4*n^4 - 14*n^3 - 13*n^2)*x^3)*D[P[n - 2, x], x] - ((n^2 - 1)*x^2 + (4*n^3 + 4*n^2 - 1)*x^3 + (4*n^3 + 3*n^2)*x^4)*D[P[n - 2, x], {x, 2}] + (2*n^3 + 3*n^2 + n + 4*n^3*x)*P[n - 1, x] - ((2*n^2 - 1)*x + (4*n^3 + 4*n^2)*x^2)*D[P[n - 1, x], x])/(n^2)]];\nTable[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten","pari":"P(N) = {my(v = vector(N + 1)); v[1] = 1; v[2] = 2 + 4*x; for(n = 2, N, my(i = n + 1); v[i] = (-(n^4 + 2*n^3 + n^2 + (4*n^4 - 4)*x - (9*n^4 - 2*n^3 - 7*n^2)*x^2)*v[i - 2] + ((2*n^3 + n^2 + 1)*x + (4*n^4 + 4*n^3 - 4*n^2 + 5)*x^2 + (4*n^4 - 14*n^3 - 13*n^2)*x^3)*deriv(v[i - 2], x) - ((n^2 - 1)*x^2 + (4*n^3 + 4*n^2 - 1)*x^3 + (4*n^3 + 3*n^2)*x^4)*deriv(deriv(v[i - 2], x), x) + (2*n^3 + 3*n^2 + n + 4*n^3*x)*v[i - 1] - ((2*n^2 - 1)*x + (4*n^3 + 4*n^2)*x^2)*deriv(v[i - 1], x))/(n^2)); v};\nrow(n, p) = {my(c = Vecrev(p)); c};\nT(N) = {my(v = P(N)); concat(vector(#v, i, row(i - 1, v[i])))};\nT(9)","reading":{"first_row":0,"offset":0,"reversed":false,"trim_zeros":"none","regular_triangle":false}}],"derivation":{"root":"n-3b2b9e31bd278a34","steps":["(\u03b8+1)\u00b2"],"edges":["e-cbf3127edaff7909"],"text":"A090582 \u2192 (\u03b8+1)\u00b2"},"headline":"A090582 \u2192 (\u03b8+1)\u00b2","chain":null}