

A306302


Number of regions into which a figure made up of a row of n adjacent congruent rectangles is divided upon drawing diagonals of all possible rectangles (a(0)=0 by convention).


47



0, 4, 16, 46, 104, 214, 380, 648, 1028, 1562, 2256, 3208, 4384, 5924, 7792, 10052, 12744, 16060, 19880, 24486, 29748, 35798, 42648, 50648, 59544, 69700, 80992, 93654, 107596, 123374, 140488, 159704, 180696, 203684, 228624, 255892, 285152, 317400, 352096, 389576
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

0,2


COMMENTS

Assuming that the rectangles have vertices at (k,0) and (k,1), k=0..n, the projective map (x,y) > ((1y)/(x+1),y/(x+1)) maps their partition to the partition of the right isosceles triangle described by Alekseyev et al. (2015), for which Theorem 13 gives the number of regions, line segments, and intersection points.  Max Alekseyev, Apr 10 2019


LINKS



FORMULA

a(n) = Sum_{i=1..n, j=1..n, gcd(i,j)=1} (n+1i)*(n+1j) + n^2 + 2*n.  N. J. A. Sloane, Apr 11 2020
a(n) = 2n(n+1) + Sum_{i=2..n} (n+1i)*(2n+2i)*phi(i).  Chai Wah Wu, Aug 16 2021


MAPLE

z := proc(n)
local a, b, r ;
r := 0 ;
for a from 1 to n do
for b from 1 to n do
if igcd(a, b) = 1 then
r := r+(n+1a)*(n+1b);
end if;
end do:
end do:
r ;
end proc:
a := n> z(n)+n^2+2*n;
[seq(a(n), n=1..50)];


MATHEMATICA

z[n_] := Sum[(n  i + 1)(n  j + 1) Boole[GCD[i, j] == 1], {i, n}, {j, n}];
a[0] = 0;
a[n_] := z[n] + n^2 + 2n;


PROG

(Python)
from sympy import totient
def A306302(n): return 2*n*(n+1) + sum(totient(i)*(n+1i)*(2*n+2i) for i in range(2, n+1)) # Chai Wah Wu, Aug 16 2021


CROSSREFS

See A331755 for the number of vertices, A331757 for the number of edges.


KEYWORD

nonn


AUTHOR



EXTENSIONS



STATUS

approved



