login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A271328 a(n) = A269347(3*n)/3. 5
1, 5, 10, 17, 28, 37, 50, 65, 82, 106, 122, 145, 170, 197, 228, 257, 294, 325, 362, 406, 442, 485, 530, 577, 628, 677, 730, 790, 842, 906, 962, 1025, 1090, 1161, 1228, 1297, 1376, 1445, 1522, 1606, 1682, 1765, 1850, 1937, 2028, 2117, 2210 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
a(n) is equal to n^2 + 1 with predictable regularity; in particular, the values of n for which a(n) does not equal n^2 + 1 are exactly those values n for which 3n is divisible by A269347(3*m) for some m with 1 < m < n. This is in part because 1 + sum{i=1...n}(2i - 1) = n^2 + 1; when computing a(n), each term has this form except when the named condition holds.
For example, a(5) does not equal 5^2 + 1 because 3(5) is divisible by A269347(6).
LINKS
FORMULA
a(n) = A269347(3*n)/3.
PROG
(PARI) lista(nn) = {nn *= 3; va = vector(nn); va[1] = 1; for (n=2, nn, va[n] = sum(k=1, n-1, k*((n % va[k])==0)); ); vector(#va\3, n, va[3*n]/3); } \\ Michel Marcus, Apr 04 2016
CROSSREFS
Sequence in context: A076598 A306011 A080341 * A086653 A215449 A174062
KEYWORD
easy,nonn
AUTHOR
Alec Jones, Apr 04 2016
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 11:59 EDT 2024. Contains 371254 sequences. (Running on oeis4.)