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.

A279451 Number of nonequivalent ways to place 7 points on an n X n square grid so that no more than 2 points are on a vertical or horizontal straight line.

Original entry on oeis.org

0, 0, 0, 115, 11810, 326190, 4444935, 38675954, 246563232, 1248782460, 5296300670, 19499431941, 63958228738, 190528987506, 523151460045, 1339408935540, 3227223506896, 7372750196952, 16069268866908, 33586411339335, 67610793877650, 131569779776182, 248290280743571
Offset: 1

Views

Author

Heinrich Ludwig, Dec 22 2016

Keywords

Comments

Column 8 of A279453.
Rotations and reflections of placements are not counted. For numbers if they are to be counted see A279441.
For condition "no more than 2 points on straight lines at any angle", see A235458.

Crossrefs

Same problem but 2..6 points: A014409, A279447, A279448, A279449, A279450.

Programs

  • Mathematica
    Table[(n^14 - 91 n^12 + 420 n^11 + 693 n^10 - 10500 n^9 + 33647 n^8 - 45316 n^7 + 3682 n^6 + 62300 n^5 - 51996 n^4 - 28504 n^3 + 54384 n^2 - 18720 n)/40320 + Boole[OddQ@ n] (2 n^7 - 17 n^6 + 50 n^5 - 59 n^4 + 38 n^3 - 71 n^2 + 102 n - 45)/384, {n, 23}] (* or *)
    Rest@ CoefficientList[Series[x^4*(115 + 11005 x + 245015 x^2 + 2317550 x^3 + 12037814 x^4 + 39232894 x^5 + 85494738 x^6 + 129182670 x^7 + 135873108 x^8 + 97856368 x^9 + 44499480 x^10 + 9709722 x^11 - 1359254 x^12 - 1352974 x^13 - 257282 x^14 + 13866 x^15 + 7705 x^16 + 419 x^17 + x^18)/((1 - x)^15*(1 + x)^8), {x, 0, 23}], x] (* Michael De Vlieger, Dec 23 2016 *)
  • PARI
    concat(vector(3), Vec(x^4*(115 +11005*x +245015*x^2 +2317550*x^3 +12037814*x^4 +39232894*x^5 +85494738*x^6 +129182670*x^7 +135873108*x^8 +97856368*x^9 +44499480*x^10 +9709722*x^11 -1359254*x^12 -1352974*x^13 -257282*x^14 +13866*x^15 +7705*x^16 +419*x^17 +x^18) / ((1 -x)^15*(1 +x)^8) + O(x^30))) \\ Colin Barker, Dec 23 2016

Formula

a(n) = (n^14 -91*n^12 +420*n^11 +693*n^10 -10500*n^9 +33647*n^8 -45316*n^7 +3682*n^6 +62300*n^5 -51996*n^4 -28504*n^3 +54384*n^2 -18720*n)/40320 + IF(MOD(n, 2) = 1, 2*n^7 -17*n^6 +50*n^5 -59*n^4 +38*n^3 -71*n^2 +102*n -45)/384.
G.f.: x^4*(115 +11005*x +245015*x^2 +2317550*x^3 +12037814*x^4 +39232894*x^5 +85494738*x^6 +129182670*x^7 +135873108*x^8 +97856368*x^9 +44499480*x^10 +9709722*x^11 -1359254*x^12 -1352974*x^13 -257282*x^14 +13866*x^15 +7705*x^16 +419*x^17 +x^18) / ((1 -x)^15*(1 +x)^8). - Colin Barker, Dec 23 2016