login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A157068
Number of integer sequences of length n+1 with sum zero and sum of absolute values 38.
2
2, 114, 3612, 80180, 1374690, 19234194, 227605448, 2335932504, 21186110970, 172295622730, 1271112537684, 8588601364668, 53573492643034, 310601807143530, 1683493452034320, 8573748834211984, 41210997268585158, 187693442844729174, 812839595630249540
OFFSET
1,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (39,-741,9139,-82251,575757,-3262623,15380937, -61523748,211915132,-635745396,1676056044,-3910797436,8122425444,-15084504396, 25140840660,-37711260990,51021117810,-62359143990,68923264410,-68923264410, 62359143990,-51021117810,37711260990,-25140840660,15084504396,-8122425444, 3910797436,-1676056044,635745396,-211915132,61523748,-15380937,3262623,-575757, 82251,-9139,741,-39,1).
FORMULA
a(n) = T(n,19); T(n,k) = Sum_{i=1..n} binomial(n+1, i)*binomial(k-1, i-1)*binomial(n-i+k, k).
From G. C. Greubel, Jan 25 2022: (Start)
a(n) = (n+1)*binomial(n+18, 19)*Hypergeometric3F2([-18, -n, 1-n], [2, -n-18], 1).
a(n) = (35345263800/38!)*n*(n+1)*(778817392288148379660189696000000 + 1984223956005743569581323059200000*n + 3392214823876583668626122342400000*n^2 + 3227079634641025484578928197632000*n^3 + 2701114821085872776574503662387200*n^4 + 1477486663626167257723210367631360*n^5 + 794678697494482855499280703586304*n^6 + 289485264886342590944226501328896*n^7 + 112195641614805001937808853208064*n^8 + 29309532027252333838411983247872*n^9 + 8732100429652853130168723017472*n^10 + 1708566742801697011435174735872*n^11 + 408081704870580048838437092992*n^12 + 61460345467484307832839519168*n^13 + 12123027157132710911533327584*n^14 + 14298582910205269163512480328n^15 + 238150505845545646647030204*n^16 + 222226805381963345901159308n^17 + 3179819458407554816818235*n^18 + 235823049245552968253250*n^19 + 29394217444775030780985*n^20 + 17315150592375085755608n^21 + 190160234133314656140*n^22 + 8844512620448927880*n^23 + 864030358357843740*n^24 + 31339517913669420*n^25 + 2745580274521866*n^26 + 76036376515644*n^27 + 6015727425006*n^28 + 122857968168*n^29 + 8831668028*n^30 + 125358408*n^31 + 8231808*n^32 + 72522*n^33 + 4371*n^34 + 18*n^35 + n^36).
G.f.: 2*x*(1 + 18*x + 324*x^2 + 2754*x^3 + 23409*x^4 + 124848*x^5 + 665856*x^6 + 2496960*x^7 + 9363600*x^8 + 26218080*x^9 + 73410624*x^10 + 159056352*x^11 + 344622096*x^12 + 590780736*x^13 + 1012766976*x^14 + 1392554592*x^15 + 1914762564*x^16 + 2127513960*x^17 + 2363904400*x^18 + 2127513960*x^19 + 1914762564*x^20 + 1392554592*x^21 + 1012766976*x^22 + 590780736*x^23 + 344622096*x^24 + 159056352*x^25 + 73410624*x^26 + 26218080*x^27 + 9363600*x^28 + 2496960*x^29 + 665856*x^30 + 124848*x^31 + 23409*x^32 + 2754*x^33 + 324*x^34 + 18*x^35 + x^36)/(1-x)^39. (End)
MATHEMATICA
A103881[n_, k_]:= (n+1)*Binomial[n+k-1, k]*HypergeometricPFQ[{1-n, -n, 1-k}, {2, 1-n - k}, 1];
A157068[n_]:= A103881[n, 19];
Table[A157068[n], {n, 50}] (* G. C. Greubel, Jan 25 2022 *)
PROG
(Sage)
def A103881(n, k): return sum( binomial(n+1, i)*binomial(k-1, i-1)*binomial(n-i+k, k) for i in (0..n) )
def A157068(n): return A103881(n, 19)
[A157068(n) for n in (1..50)] # G. C. Greubel, Jan 25 2022
CROSSREFS
Sequence in context: A224871 A230471 A140986 * A008271 A362575 A209184
KEYWORD
nonn
AUTHOR
R. H. Hardin, Feb 22 2009
STATUS
approved