%I #7 Aug 02 2021 12:14:48

%S 0,1,2,10,9,23,26,43,40,67,73,80,69,124,155,156,186,232,186,247,258,

%T 312,388,345,310,568,377,541,580,585,530,673,615,911,777,1048,830,948,

%U 702,1271,1024,1286,1403,1320,1159,1458,1388,1545,1590,1632,1714,1906,1883

%N a(n) = Sum_{k=1..n} (n^k mod (n+k)).

%p P:=proc(i) local a,n; for n from 0 by 1 to i do a:=sum('(n^k mod (n+k))','k'=1..n); print(a); od; end: P(100);

%o (PARI) a(n) = sum(k=1, n, n^k%(n+k)); \\ _Jinyuan Wang_, Aug 02 2021

%Y Cf. A156557, A156558, A156559

%K nonn,easy

%O 0,3

%A _Paolo P. Lava_ and _Giorgio Balzarotti_, Feb 10 2009