OFFSET
3,1
COMMENTS
See A335057 for illustrations.
LINKS
Lars Blomberg, Table of n, a(n) for n = 3..270
FORMULA
Empirically for n <= 270.
For n > 3 select the row in the table below for which d = n mod m. Then a(n) = (a*n^2+b*n+c)/denom.
+=============================================+
| d | m | a | b | c | denom |
+---------------------------------------------+
| 1, 5 | 6 | 5 | 0 | -17 | 12 |
| 3 | 6 | 5 | -16 | 27 | 12 |
| 2, 10 | 12 | 5 | -15 | 22 | 12 |
| 4, 8 | 12 | 5 | -15 | 40 | 12 |
| 0 | 60 | 5 | -31 | 0 | 12 |
| 6, 18, 42, 54 | 60 | 5 | -31 | 126 | 12 |
| 12, 24, 36, 48 | 60 | 5 | -31 | 144 | 12 |
| 30 | 60 | 5 | -31 | -18 | 12 |
+=============================================+
PROG
(PARI) bc=[[5, 0, -17, 12], [5, -16, 27, 12], [5, -15, 22, 12], [5, -15, 40, 12], [5, -31, 0, 12], [5, -31, 126, 12], [5, -31, 144, 12], [5, -31, -18, 12]];
m=[[1, 6, 1], [5, 6, 1], [3, 6, 2], [2, 12, 3], [10, 12, 3], [4, 12, 4], [8, 12, 4], [0, 60, 5], [6, 60, 6], [18, 60, 6], [42, 60, 6], [54, 60, 6], [12, 60, 7], [24, 60, 7], [36, 60, 7], [48, 60, 7], [30, 60, 8]];
ix(n)=for(i=1, length(m), x=m[i]; if(n%x[2]==x[1], return(x[3]))); -1
a(n)=if(n==3, return(3)); x=bc[ix(n)]; (x[1]*n^2+x[2]*n+x[3])/x[4]
vector(200, x, a(x+2))
CROSSREFS
KEYWORD
nonn
AUTHOR
Lars Blomberg, May 24 2020
STATUS
approved