{"id":"n-64b0c861287f610b","revision":0,"first_index":0,"fingerprint":"f6-4b7f5e0d58f0c946","name":"theta_plus_one_squared of basis_map_A013609 of Pascal-(1,3,1) array.","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","8"],["7","56","36"],["20","240","432","128"],["61","976","2808","2176","400"],["182","3640","14544","19328","8800","1152"],["547","13128","67068","126848","99600","31104","3136"],["1640","45920","285984","705280","803200","428544","100352","8192"],["4921","157472","1157040","3523328","5298400","4124160","1630720","303104","20736"],["14762","531432","4500576","16326912","30628800","31732992","18289152","5677056","870912","51200"],["44287","1771480","16989588","71567616","161205600","209827584","161535360","72818688","18475776","2406400","123904"],["132860","5845840","62626896","300529536","791011200","1243759104","1208288256","727965696","266996736","57036800","6443008","294912"],["398581","19131888","226443816","1219657600","3675740400","6787929600","8000563200","6096642048","2985050880","917401600","168757248","16809984","692224"],["1195742","62178584","805780656","4814057600","16353408800","34727605632","48261609984","44745646080","27731206656","11354163200","2991042560","482181120","42917888","1605632"],["3587227","200884712","2828976876","18567486592","70212258800","168675214464","270495927744","296550653952","223976358144","116044134400","40625767424","9338978304","1338068992","107577344","3686400"]],"degrees":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14],"row_sums":["1","10","99","820","6421","47646","341431","2379112","16219881","108623794","716766283","4670930844","30115219069","192374600134","1218969561375"],"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","8"]],"recurrence":{"denominator":{"coeffs":[["4","-4","-6","-5","2","0"],["0","24","-4","-16","-4","0"],["0","0","44","20","-6","-2"],["0","0","0","24","8","-2"],["0","0","0","0","4","0"]]},"terms":[{"coeff":{"coeffs":[["12","84","222","273","150","12","-24"],["0","72","396","792","648","114","-60"],["0","0","132","540","678","228","-48"],["0","0","0","72","216","138","-12"],["0","0","0","0","12","24","0"]]},"deriv_order":0,"offset":2,"sign":"none"},{"coeff":{"coeffs":[["8","24","12","-50","-48","-6","-4"],["0","48","128","120","12","-16","-12"],["0","0","88","200","140","10","-12"],["0","0","0","48","80","28","-4"],["0","0","0","0","8","8","0"]]},"deriv_order":0,"offset":1,"sign":"none"}]},"schema":"polytool.recurrence.v1"}}],"properties":[{"id":179092,"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-1489b6d0566d06f7"]},{"id":179093,"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-1489b6d0566d06f7"]},{"id":179094,"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-1489b6d0566d06f7"]}],"region_exact":{},"incoming_edges":[{"id":"e-1489b6d0566d06f7","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-64b0c861287f610b","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":"0b798266c9696fcdc43c674aa1a6cd7d48501ce333e39e4c591c5f0ee0e9a7d8","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-c1d69b4cd7e1b71e","role":"sequence","name":"basis_map_A013609 of Pascal-(1,3,1) array.","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":[],"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,      8;\n  2:     7,     56,      36;\n  3:    20,    240,     432,      128;\n  4:    61,    976,    2808,     2176,      400;\n  5:   182,   3640,   14544,    19328,     8800,     1152;\n  6:   547,  13128,   67068,   126848,    99600,    31104,     3136;\n  7:  1640,  45920,  285984,   705280,   803200,   428544,   100352,    8192;\n  8:  4921, 157472, 1157040,  3523328,  5298400,  4124160,  1630720,  303104,  20736;\n  9: 14762, 531432, 4500576, 16326912, 30628800, 31732992, 18289152, 5677056, 870912, 51200;\n","recurrence_sentence":"The row polynomials satisfy (4+(24*n-4)*x+(44*n^2-4*n-6)*x^2+(24*n^3+20*n^2-16*n-5)*x^3+(4*n^4+8*n^3-6*n^2-4*n+2)*x^4-(2*n^3+2*n^2)*x^5)*P(n,x) = (12+(72*n+84)*x+(132*n^2+396*n+222)*x^2+(72*n^3+540*n^2+792*n+273)*x^3+(12*n^4+216*n^3+678*n^2+648*n+150)*x^4+(24*n^4+138*n^3+228*n^2+114*n+12)*x^5-(12*n^3+48*n^2+60*n+24)*x^6)*P(n-2,x) + (8+(48*n+24)*x+(88*n^2+128*n+12)*x^2+(48*n^3+200*n^2+120*n-50)*x^3+(8*n^4+80*n^3+140*n^2+12*n-48)*x^4+(8*n^4+28*n^3+10*n^2-16*n-6)*x^5-(4*n^3+12*n^2+12*n+4)*x^6)*P(n-1,x) for n >= 2, with P(0,x) = 1, P(1,x) = 2 + 8*x.","mathematica":"P[0, x_] := 1;\nP[1, x_] := 2 + 8*x;\nP[n_, x_] := P[n, x] = Expand[Together[((12 + (72*n + 84)*x + (132*n^2 + 396*n + 222)*x^2 + (72*n^3 + 540*n^2 + 792*n + 273)*x^3 + (12*n^4 + 216*n^3 + 678*n^2 + 648*n + 150)*x^4 + (24*n^4 + 138*n^3 + 228*n^2 + 114*n + 12)*x^5 - (12*n^3 + 48*n^2 + 60*n + 24)*x^6)*P[n - 2, x] + (8 + (48*n + 24)*x + (88*n^2 + 128*n + 12)*x^2 + (48*n^3 + 200*n^2 + 120*n - 50)*x^3 + (8*n^4 + 80*n^3 + 140*n^2 + 12*n - 48)*x^4 + (8*n^4 + 28*n^3 + 10*n^2 - 16*n - 6)*x^5 - (4*n^3 + 12*n^2 + 12*n + 4)*x^6)*P[n - 1, x])/(4 + (24*n - 4)*x + (44*n^2 - 4*n - 6)*x^2 + (24*n^3 + 20*n^2 - 16*n - 5)*x^3 + (4*n^4 + 8*n^3 - 6*n^2 - 4*n + 2)*x^4 - (2*n^3 + 2*n^2)*x^5)]];\nTable[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten","pari":"P(N) = {my(v = vector(N + 1)); v[1] = 1; v[2] = 2 + 8*x; for(n = 2, N, my(i = n + 1); v[i] = ((12 + (72*n + 84)*x + (132*n^2 + 396*n + 222)*x^2 + (72*n^3 + 540*n^2 + 792*n + 273)*x^3 + (12*n^4 + 216*n^3 + 678*n^2 + 648*n + 150)*x^4 + (24*n^4 + 138*n^3 + 228*n^2 + 114*n + 12)*x^5 - (12*n^3 + 48*n^2 + 60*n + 24)*x^6)*v[i - 2] + (8 + (48*n + 24)*x + (88*n^2 + 128*n + 12)*x^2 + (48*n^3 + 200*n^2 + 120*n - 50)*x^3 + (8*n^4 + 80*n^3 + 140*n^2 + 12*n - 48)*x^4 + (8*n^4 + 28*n^3 + 10*n^2 - 16*n - 6)*x^5 - (4*n^3 + 12*n^2 + 12*n + 4)*x^6)*v[i - 1])/(4 + (24*n - 4)*x + (44*n^2 - 4*n - 6)*x^2 + (24*n^3 + 20*n^2 - 16*n - 5)*x^3 + (4*n^4 + 8*n^3 - 6*n^2 - 4*n + 2)*x^4 - (2*n^3 + 2*n^2)*x^5)); 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-d8f9ad9ce0e039d0","steps":["basis A013609","(\u03b8+1)\u00b2"],"edges":["e-471a582f60d8a81d","e-1489b6d0566d06f7"],"text":"A081578 \u2192 basis A013609 \u2192 (\u03b8+1)\u00b2"},"headline":"A081578 \u2192 basis A013609 \u2192 (\u03b8+1)\u00b2","chain":null}
