%I #21 Oct 27 2023 21:17:12
%S 1,3,8,12,21,29,40,52,67,83,100,116,140,160,185,210,237,264,298,327,
%T 363,397,435,472,514,557,602,644,690,741,791,837,897,950,1009,1063,
%U 1126,1185,1253,1313,1381,1450,1521,1593,1667,1739,1820,1894,1973,2054,2140
%N a(n) = Sum_{k=1..n} (-1)^(k+1) * floor((n/k)^2).
%F a(n) = A309081(n^2).
%t a[n_] := Sum[(-1)^(k + 1)*Floor[(n/k)^2], {k, 1, n}]; Array[a, 50] (* _Amiram Eldar_, Dec 20 2021 *)
%o (PARI) a(n) = sum(k=1, n, (-1)^(k+1)*(n^2\k^2));
%o (Python)
%o def A350221(n): return (m:=n**2)+sum(m//k**2 if k&1 else -(m//k**2) for k in range(2,n+1)) # _Chai Wah Wu_, Oct 27 2023
%Y Cf. A000290, A153818, A309081, A344720, A350222, A350224.
%K nonn
%O 1,2
%A _Seiichi Manyama_, Dec 20 2021