This sequence's twin numbers are given in A133135. - _Paul Curtz (bpcrtz(AT)free.fr), _, Aug 07 2008
This sequence's twin numbers are given in A133135. - _Paul Curtz (bpcrtz(AT)free.fr), _, Aug 07 2008
_Paul D. Hanna (pauldhanna(AT)juno.com), _, Dec 28 2004
proposed
approved
a(n) = 1 - Sum_{k=0, n-1, n} A001469(k) for n>0, with a(0)=1.
approved
proposed
This sequence's twin numbers are given in A133135. - Paul Curtz (bpcrtz(AT)free.fr), Aug 07 2008
sign,new
sign
sign,new
sign
Paul D . Hanna (pauldhanna(AT)juno.com), Dec 28 2004
a(n) = 1 - Sum_{k=0, n-1} A001469(k) for n>0, with a(0)=1.
sign,new
sign
(PARI) {a(n)=local(M=matrix(n+2, n+2)); M[1, 1]=1; if(n>0, M[2, 1]=1; M[2, 2]=1); for(r=3, n+2, for(c=1, r, M[r, c]=if(c==1, M[r-1, 1], if(c==r, 1, M[r, c]=M[r-1, c]-((matrix(r-1, r-1, i, j, M[i, j]))^-1)[r-1, c-1])))); return(if(n==0, 1, M[n+2, 2]))}
sign,new
sign
1, 2, 1, 4, -13, 142, -1931, 36296, -893273, 27927346, -1081725559, 50861556172, -2854289486309, 188475382997654, -14467150771771043, 1277417937676246672, -128570745743431055281, 14632875988040732946106, -1869882665740777942166543, 266593648798424693540514836
0,2
1-a(n+1) equals the n-th partial sum of the Genocchi numbers (A001469).
a(n) = 1 - Sum_{k=0,n-1} A001469(k) for n>0, with a(0)=1.
(PARI) {a(n)=local(M=matrix(n+2, n+2)); M[1, 1]=1; if(n>0, M[2, 1]=1; M[2, 2]=1); for(r=3, n+2, for(c=1, r, M[r, c]=if(c==1, M[r-1, 1], if(c==r, 1, M[r, c]=M[r-1, c]-((matrix(r-1, r-1, i, j, M[i, j]))^-1)[r-1, c-1])))); return(if(n==0, 1, M[n+2, 2]))}
sign
Paul D Hanna (pauldhanna(AT)juno.com), Dec 28 2004
approved