This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A341276 #20 Apr 01 2021 01:43:19 %S A341276 1,6,16,32,53,81,113,153,197,248,304,368,434,510,590,676,767,867,969, %T A341276 1081,1195,1317,1445,1581,1717,1864,2016,2174,2336,2508,2680,2864, %U A341276 3050,3244,3444,3650,3857,4077,4301,4531,4763,5007,5251,5507,5765 %N A341276 a(n) = 1 + 3*n*(n+1) - Sum_{k=1..n} d(k), where d(k) is the number of divisors of k. %C A341276 Number of regions of the plane delimited by the 3n curves: x=k, y=k, x*y=k for k=1..n with x>0 and y>0. %F A341276 a(n) = 1 + A028896(n) - A006218(n). %e A341276 The curves with equations x=1, y=1 and x*y=1 are delimiting a(1)=6 regions in the (+; +) quadrant of the plane. With the addition of the curves with equations x=2, y=2 and x*y=2, the number of delimited regions reaches a(2)=16. %t A341276 Table[1+3n(n+1)-Sum[DivisorSigma[0,k],{k,n}],{n,0,44}] (* _Stefano Spezia_, Feb 08 2021 *) %o A341276 (PARI) a(n)=1+3*n*(n+1)-sum(k=1,n,n\k) %o A341276 (Python) %o A341276 from sympy import integer_nthroot %o A341276 def A341276(n): return 1+3*n*(n+1)-2*sum(n//k for k in range(1,integer_nthroot(n,2)[0]+1))+integer_nthroot(n,2)[0]**2 # _Chai Wah Wu_, Mar 31 2021 %Y A341276 Cf. A000005, A006218, A028896. %K A341276 nonn %O A341276 0,2 %A A341276 _Luc Rousseau_, Feb 08 2021