{"id":"n-584d62f138731882","revision":0,"first_index":0,"fingerprint":"f6-c37f57020308301f","name":"theta_plus_one_squared of Triangle: row #n has n+1 terms. T(n,m) = 4^m (2n+1)! / ( (2n-2m)! (2m+1)! ).","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"],["3","16"],["5","160","144"],["7","560","3024","1024"],["9","1344","18144","36864","6400"],["11","2640","66528","337920","352000","36864"],["13","4576","185328","1757184","4576000","2875392","200704"],["15","7280","432432","6589440","32032000","50319360","21073920","1048576"],["17","10880","891072","19914752","155584000","456228864","477675520","142606336","5308416"],["19","15504","1674432","51597312","591219200","2786254848","5445500928","4064280576","907739136","26214400"],["21","21280","2930256","119070720","1881152000","13002522624","40841256960","56899928064","31770869760","5505024000","126877696"],["23","28336","4845456","251040768","5230016000","49843003392","229618622464","514131492864","535868669952","232128512000","32100057088","603979776"],["25","36800","7650720","492236800","13075040000","164317593600","1043721011200","3427543285760","5741450035200","4642570240000","1605002854400","181193932800","2835349504"],["27","46800","11625120","909342720","29995680000","480628961280","4025781043200","18228298383360","44783310274560","58197934080000","37557066792960","10599845068800","995207675904","13153337344"],["29","58464","17100720","1598238720","64096032000","1275394498560","13620559196160","81326254325760","275485211688960","525074694144000","544577468497920","286902473195520","67342386069504","5340254961664","60397977600"]],"degrees":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14],"row_sums":["1","19","309","4615","62761","795963","9599197","111503023","1258219857","13874496355","150149653381","1599780268119","16822209266809","174912142312075","1801070805985581"],"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"]],"recurrence":{"denominator":{"coeffs":[["-1","4","0"],["2","-44","0"],["0","112","-16"],["0","-80","96"],["0","0","-192"],["0","0","128"]]},"terms":[{"coeff":{"coeffs":[["1","-84","-592","-192"],["2","124","1424","896"],["0","32","784","-832"],["0","-80","-1760","-2432"],["0","0","-64","1280"],["0","0","128","1536"]]},"deriv_order":0,"offset":1,"sign":"none"},{"coeff":{"coeffs":[["0","4","40","-192","-128"],["0","0","-96","256","512"],["0","0","-160","704","-256"],["0","0","0","512","-2048"],["0","0","0","256","-1024"]]},"deriv_order":1,"offset":1,"sign":"none"}]},"schema":"polytool.recurrence.v1"}}],"properties":[{"id":184900,"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-d53f3b6943517760"]},{"id":184901,"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-d53f3b6943517760"]},{"id":184902,"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-d53f3b6943517760"]},{"id":184903,"position":3,"claim":"interlacing","from_index":1,"region_key":null,"provenance":"edges","claim_data":{"claim":"interlacing","from_index":1,"kind":"weak","orientation":"unordered"},"label":"Interlacing of consecutive rows","description":"Rows P_(m-1), P_m interlace (weak, unordered) for every m >= 1","region_exact":null,"edges":["e-d53f3b6943517760"]}],"region_exact":{},"incoming_edges":[{"id":"e-d53f3b6943517760","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-584d62f138731882","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":"d0a97e5e9fce78160f9b0222b3d79ce76bf2267aef33b7a5f9f470d32bd34387","verification_polytool":"replayed","verification_lean":"not_attempted","claims":[{"claim":"real_rooted","from_index":0,"region":{"upper":{"closed":true,"value":"0"}}},{"claim":"interlacing","from_index":1,"kind":"weak","orientation":"unordered"}],"premises":[{"position":0,"node_id":"n-9aec4822212db1dc","role":"sequence","name":"Triangle: row #n has n+1 terms. T(n,m) = 4^m (2n+1)! / ( (2n-2m)! (2m+1)! ).","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:  3,    16;\n  2:  5,   160,     144;\n  3:  7,   560,    3024,     1024;\n  4:  9,  1344,   18144,    36864,      6400;\n  5: 11,  2640,   66528,   337920,    352000,      36864;\n  6: 13,  4576,  185328,  1757184,   4576000,    2875392,     200704;\n  7: 15,  7280,  432432,  6589440,  32032000,   50319360,   21073920,    1048576;\n  8: 17, 10880,  891072, 19914752, 155584000,  456228864,  477675520,  142606336,   5308416;\n  9: 19, 15504, 1674432, 51597312, 591219200, 2786254848, 5445500928, 4064280576, 907739136, 26214400;\n","recurrence_sentence":"The row polynomials satisfy (2*n-1-(80*n^3-112*n^2+44*n-4)*x+(128*n^5-192*n^4+96*n^3-16*n^2)*x^2)*P(n,x) = (2*n+1-(80*n^3-32*n^2-124*n+84)*x+(128*n^5-64*n^4-1760*n^3+784*n^2+1424*n-592)*x^2+(1536*n^5+1280*n^4-2432*n^3-832*n^2+896*n-192)*x^3)*P(n-1,x) + (4*x-(160*n^2+96*n-40)*x^2+(256*n^4+512*n^3+704*n^2+256*n-192)*x^3-(1024*n^4+2048*n^3+256*n^2-512*n+128)*x^4)*d/dx P(n-1,x) for n >= 1, with P(0,x) = 1.","mathematica":"P[0, x_] := 1;\nP[n_, x_] := P[n, x] = Expand[Together[((2*n + 1 - (80*n^3 - 32*n^2 - 124*n + 84)*x + (128*n^5 - 64*n^4 - 1760*n^3 + 784*n^2 + 1424*n - 592)*x^2 + (1536*n^5 + 1280*n^4 - 2432*n^3 - 832*n^2 + 896*n - 192)*x^3)*P[n - 1, x] + (4*x - (160*n^2 + 96*n - 40)*x^2 + (256*n^4 + 512*n^3 + 704*n^2 + 256*n - 192)*x^3 - (1024*n^4 + 2048*n^3 + 256*n^2 - 512*n + 128)*x^4)*D[P[n - 1, x], x])/(2*n - 1 - (80*n^3 - 112*n^2 + 44*n - 4)*x + (128*n^5 - 192*n^4 + 96*n^3 - 16*n^2)*x^2)]];\nTable[CoefficientList[P[n, x], x], {n, 0, 9}] // Flatten","pari":"P(N) = {my(v = vector(N + 1)); v[1] = 1; for(n = 1, N, my(i = n + 1); v[i] = ((2*n + 1 - (80*n^3 - 32*n^2 - 124*n + 84)*x + (128*n^5 - 64*n^4 - 1760*n^3 + 784*n^2 + 1424*n - 592)*x^2 + (1536*n^5 + 1280*n^4 - 2432*n^3 - 832*n^2 + 896*n - 192)*x^3)*v[i - 1] + (4*x - (160*n^2 + 96*n - 40)*x^2 + (256*n^4 + 512*n^3 + 704*n^2 + 256*n - 192)*x^3 - (1024*n^4 + 2048*n^3 + 256*n^2 - 512*n + 128)*x^4)*deriv(v[i - 1], x))/(2*n - 1 - (80*n^3 - 112*n^2 + 44*n - 4)*x + (128*n^5 - 192*n^4 + 96*n^3 - 16*n^2)*x^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-9aec4822212db1dc","steps":["(\u03b8+1)\u00b2"],"edges":["e-d53f3b6943517760"],"text":"A085841 \u2192 (\u03b8+1)\u00b2"},"headline":"A085841 \u2192 (\u03b8+1)\u00b2","chain":null}