|
%I #17 Feb 26 2021 02:10:14
%S 1,5,6,4,9,11,20,48,31,55,30,92,52,112,94,150,113,193,152,242,218,232,
%T 246,308,301,317,441,353,311,413,431,385,399,515,570,592,614,674,542,
%U 558,697,761,785,853,969,907,745,1099,837,1057,1085,1009,1245,1439,1477
%N First column of array A082218.
%o (PARI) lista(nn) = { my(a=matrix(nn, nn)); S=Set();
%o for(s=2, nn+1, if(s%2, i0=1; i1=s-1; i2=1, i0=s-1; i1=1; i2=-1);
%o forstep(i=i0, i1, i2, j=s-i;
%o ii=sum(k=1, j-1, a[i, k]); jj=sum(k=1, i-1, a[k, j]);
%o c=chinese(Mod(ii, j), Mod(jj, i));
%o t=component(c, 1)-lift(c); while(setsearch(S, t), t+=component(c, 1));
%o a[i, j]=t; S=setunion(S, [t]);
%o if(j==1, print1(", ", sum(k=1, j, a[i, j])/j) ); ))} \\ This is a modification of _Max Alekseyev_'s PARI program from A082219. - _Petros Hadjicostas_, Feb 25 2021
%Y Cf. A082218, A082219, A082221, A082222, A082223.
%K nonn
%O 1,2
%A _Amarnath Murthy_ and Meenakshi Srikanth (menakan_s(AT)yahoo.com), Apr 09 2003
%E More terms from _Max Alekseyev_, Jun 08 2007
%E Edited by _N. J. A. Sloane_, Jun 11 2007
%E Name edited by _Petros Hadjicostas_, Feb 25 2021
|
