cp's OEIS Frontend

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.

Showing 1-2 of 2 results.

A090466 Regular figurative or polygonal numbers of order greater than 2.

Original entry on oeis.org

6, 9, 10, 12, 15, 16, 18, 21, 22, 24, 25, 27, 28, 30, 33, 34, 35, 36, 39, 40, 42, 45, 46, 48, 49, 51, 52, 54, 55, 57, 58, 60, 63, 64, 65, 66, 69, 70, 72, 75, 76, 78, 81, 82, 84, 85, 87, 88, 90, 91, 92, 93, 94, 95, 96, 99, 100, 102, 105, 106, 108, 111, 112, 114, 115, 117, 118
Offset: 1

Views

Author

Robert G. Wilson v, Dec 01 2003

Keywords

Comments

The sorted k-gonal numbers of order greater than 2. If one were to include either the rank 2 or the 2-gonal numbers, then every number would appear.
Number of terms less than or equal to 10^k for k = 1,2,3,...: 3, 57, 622, 6357, 63889, 639946, 6402325, 64032121, 640349979, 6403587409, 64036148166, 640362343980, ..., . - Robert G. Wilson v, May 29 2014
The n-th k-gonal number is 1 + k*n(n-1)/2 - (n-1)^2 = A057145(k,n).
For all squares (A001248) of primes p >= 5 at least one a(n) exists with p^2 = a(n) + 1. Thus the subset P_s(3) of rank 3 only is sufficient. Proof: For p >= 5, p^2 == 1 (mod {3,4,6,8,12,24}) and also P_s(3) + 1 = 3*s - 2 == 1 (mod 3). Thus the set {p^2} is a subset of {P_s(3) + 1}; Q.E.D. - Ralf Steiner, Jul 15 2018
For all primes p > 5, at least one polygonal number exists with P_s(k) + 1 = p when k = 3 or 4, dependent on p mod 6. - Ralf Steiner, Jul 16 2018
Numbers m such that r = (2*m/d - 2)/(d - 1) is an integer for some d, where 2 < d < m is a divisor of 2*m. If r is an integer, then m is the d-th (r+2)-gonal number. - Jianing Song, Mar 14 2021

References

  • Albert H. Beiler, Recreations In The Theory Of Numbers, The Queen Of Mathematics Entertains, Dover, NY, 1964, pp. 185-199.

Crossrefs

Cf. A057145, A001248, A177028 (A342772, A342805), A177201, A316676, A364693 (characteristic function).
Complement is A090467.
Sequence A090428 (excluding 1) is a subsequence of this sequence. - T. D. Noe, Jun 14 2012
Other subsequences: A324972 (squarefree terms), A324973, A342806, A364694.
Cf. also A275340.

Programs

  • Maple
    isA090466 := proc(n)
        local nsearch,ksearch;
        for nsearch from 3 do
            if A057145(nsearch,3) > n then
                return false;
            end if;
            for ksearch from 3 do
                if A057145(nsearch,ksearch) = n then
                    return true;
                elif A057145(nsearch,ksearch) > n then
                    break;
                end if;
            end do:
        end do:
    end proc:
    for n from 1 to 1000 do
        if isA090466(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Jul 28 2016
  • Mathematica
    Take[Union[Flatten[Table[1+k*n (n-1)/2-(n-1)^2,{n,3,100},{k,3,40}]]],67] (* corrected by Ant King, Sep 19 2011 *)
    mx = 150; n = k = 3; lst = {}; While[n < Floor[mx/3]+2, a = PolygonalNumber[n, k]; If[a < mx+1, AppendTo[ lst, a], (n++; k = 2)]; k++]; lst = Union@ lst (* Robert G. Wilson v, May 29 2014 and updated Jul 23 2018; PolygonalNumber requires version 10.4 or higher *)
  • PARI
    list(lim)=my(v=List()); lim\=1; for(n=3,sqrtint(8*lim+1)\2, for(k=3,2*(lim-2*n+n^2)\n\(n-1), listput(v, 1+k*n*(n-1)/2-(n-1)^2))); Set(v); \\ Charles R Greathouse IV, Jan 19 2017
    
  • PARI
    is(n)=for(s=3,n\3+1,ispolygonal(n,s)&&return(s)); \\ M. F. Hasler, Jan 19 2017
    
  • PARI
    isA090466(m) = my(v=divisors(2*m)); for(i=3, #v, my(d=v[i]); if(d==m, return(0)); if((2*m/d - 2)%(d - 1)==0, return(1))); 0 \\ Jianing Song, Mar 14 2021

Formula

Integer k is in this sequence iff A176774(k) < k. - Max Alekseyev, Apr 24 2018

Extensions

Verified by Don Reble, Mar 12 2006

A177202 Nontrivially polygonal-free numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 11, 13, 14, 17, 19, 23, 26, 29, 31, 37, 38, 41, 43, 47, 53, 59, 61, 62, 67, 71, 73, 74, 77, 79, 83, 86, 89, 97, 101, 103, 107, 109, 113, 119, 122, 127, 131, 134, 137, 139, 143, 146, 149, 151, 157, 158, 161, 163, 167, 173, 179
Offset: 1

Views

Author

Jonathan Vos Post, May 04 2010

Keywords

Comments

Positive integers which are not multiples of any nontrivially polygonal numbers A090466. Generalization of squarefree numbers: numbers that are not divisible by a square greater than 1 (A005117) where "square" is replaced by triangular, square, pentagonal, hexagonal, and so forth. Positive integers which are not positive integer multiples of numbers in the array of A057145 below the second row (which has every positive integer) and right of the 2nd column (which has every positive integer). That is, positive integers which are not positive integer multiples of any triangular number >3, any square >4, any pentagonal number >5, any hexagonal number >6, any heptagonal number >7, any octagonal number >8, any 9-gonal (nonagonal) number >9 and so forth. Properly includes all primes.

Crossrefs

Formula

Complement of A177201.

Extensions

Corrected and extended by Sean A. Irvine, Apr 09 2013
Showing 1-2 of 2 results.