A082223 - OEIS

A082223

a(n) = (n-th partial sum of A082220)/n = (n-th partial sum of the first column of array A082218)/n.

1, 3, 4, 4, 5, 6, 8, 13, 15, 19, 20, 26, 28, 34, 38, 45, 49, 57, 62, 71, 78, 85, 92, 101, 109, 117, 129, 137, 143, 152, 161, 168, 175, 185, 196, 207, 218, 230, 238, 246, 257, 269, 281, 294, 309, 322, 331, 347, 357, 371, 385, 397, 413, 432, 451, 468, 484, 504, 518

(list; graph; refs; listen; history; text; internal format)

OFFSET

1,2

LINKS

Table of n, a(n) for n=1..59.

PROG

(PARI) lista(nn) = { my(a=matrix(nn, nn)); my(b=vector(nn)); my(d=vector(nn)); S=Set(); for(s=2, nn+1, if(s%2, i0=1; i1=s-1; i2=1, i0=s-1; i1=1; i2=-1);

forstep(i=i0, i1, i2, j=s-i; ii=sum(k=1, j-1, a[i, k]); jj=sum(k=1, i-1, a[k, j]);

c=chinese(Mod(ii, j), Mod(jj, i));

t=component(c, 1)-lift(c); while(setsearch(S, t), t+=component(c, 1));

a[i, j]=t; S=setunion(S, [t]); b[i] = a[i, 1]; ); ); d[1]=b[1]; for(kk=1, nn, d[kk]=sum(ss=1, kk, b[ss])/kk); d; } \\ This is a modification of Max Alekseyev's PARI program from A082219. - Petros Hadjicostas, Feb 25 2021

CROSSREFS

Cf. A082218, A082219, A082220, A082221, A082222.

Sequence in context: A120172 A113454 A157726 * A292351 A098181 A322407

Adjacent sequences: A082220 A082221 A082222 * A082224 A082225 A082226

KEYWORD

nonn

AUTHOR

Amarnath Murthy and Meenakshi Srikanth (menakan_s(AT)yahoo.com), Apr 09 2003

EXTENSIONS

More terms from Max Alekseyev, Jun 08 2007

Edited by N. J. A. Sloane, Jun 11 2007

Name edited by Petros Hadjicostas, Feb 25 2021

STATUS

approved