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.

Previous Showing 21-30 of 33 results. Next

A307539 Heinz numbers of square integer partitions, where the Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Original entry on oeis.org

1, 2, 9, 125, 2401, 161051, 4826809, 410338673, 16983563041, 1801152661463, 420707233300201, 25408476896404831, 6582952005840035281, 925103102315013629321, 73885357344138503765449, 12063348350820368238715343, 3876269050118516845397872321
Offset: 0

Views

Author

Gus Wiseman, Apr 13 2019

Keywords

Examples

			The square partition (4,4,4,4) has Heinz number prime(4)^4 = 7^4 = 2401.
		

Crossrefs

Programs

  • Maple
    a:= n-> mul(ithprime(i), i=[n$n]):
    seq(a(n), n=0..20);  # Alois P. Heinz, Mar 03 2020
  • Mathematica
    Table[If[n==0,1,Prime[n]]^n,{n,0,10}]

Formula

a(n) = A330394(A088218(n)). - Alois P. Heinz, Mar 03 2020

A325230 Numbers of the form p^k * q, p and q prime, p > q, k > 0.

Original entry on oeis.org

6, 10, 14, 15, 18, 21, 22, 26, 33, 34, 35, 38, 39, 46, 50, 51, 54, 55, 57, 58, 62, 65, 69, 74, 75, 77, 82, 85, 86, 87, 91, 93, 94, 95, 98, 106, 111, 115, 118, 119, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 147, 155, 158, 159, 161, 162, 166, 177, 178
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2019

Keywords

Examples

			The sequence of terms together with their prime indices begins:
    6: {1,2}
   10: {1,3}
   14: {1,4}
   15: {2,3}
   18: {1,2,2}
   21: {2,4}
   22: {1,5}
   26: {1,6}
   33: {2,5}
   34: {1,7}
   35: {3,4}
   38: {1,8}
   39: {2,6}
   46: {1,9}
   50: {1,3,3}
   51: {2,7}
   54: {1,2,2,2}
   55: {3,5}
   57: {2,8}
   58: {1,10}
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local F;
       F:= sort(ifactors(n)[2],(a,b)-> a[1]Robert Israel, Apr 14 2019
  • Mathematica
    Select[Range[100],PrimeOmega[#/Power@@FactorInteger[#][[-1]]]==1&]
  • Python
    from sympy import factorint
    A325230_list = [n for n, m in ((n, factorint(n)) for n in range(2,10**6)) if len(m) == 2 and m[min(m)] == 1] # Chai Wah Wu, Apr 16 2019

A325231 Numbers of the form 2 * p or 3 * 2^k, p prime, k > 1.

Original entry on oeis.org

6, 10, 12, 14, 22, 24, 26, 34, 38, 46, 48, 58, 62, 74, 82, 86, 94, 96, 106, 118, 122, 134, 142, 146, 158, 166, 178, 192, 194, 202, 206, 214, 218, 226, 254, 262, 274, 278, 298, 302, 314, 326, 334, 346, 358, 362, 382, 384, 386, 394, 398, 422, 446, 454, 458, 466
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2019

Keywords

Comments

Also numbers n such that the sum of prime indices of n minus the greater of the number of prime factors of n counted with multiplicity and the largest prime index of n is 1. A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798, and their sum is A056239.

Examples

			The sequence of terms together with their prime indices begins:
    6: {1,2}
   10: {1,3}
   12: {1,1,2}
   14: {1,4}
   22: {1,5}
   24: {1,1,1,2}
   26: {1,6}
   34: {1,7}
   38: {1,8}
   46: {1,9}
   48: {1,1,1,1,2}
   58: {1,10}
   62: {1,11}
   74: {1,12}
   82: {1,13}
   86: {1,14}
   94: {1,15}
   96: {1,1,1,1,1,2}
  106: {1,16}
  118: {1,17}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Total[primeMS[#]]-Max[Length[primeMS[#]],Max[primeMS[#]]]==1&]
  • Python
    from sympy import isprime
    A325231_list = [n for n in range(6,10**6) if ((not n % 2) and isprime(n//2)) or (bin(n)[2:4] == '11' and bin(n).count('1') == 2)] # Chai Wah Wu, Apr 16 2019

A343938 Twice the number of prime factors of n minus the sum of prime indices of n, both counted with multiplicity.

Original entry on oeis.org

0, 1, 0, 2, -1, 1, -2, 3, 0, 0, -3, 2, -4, -1, -1, 4, -5, 1, -6, 1, -2, -2, -7, 3, -2, -3, 0, 0, -8, 0, -9, 5, -3, -4, -3, 2, -10, -5, -4, 2, -11, -1, -12, -1, -1, -6, -13, 4, -4, -1, -5, -2, -14, 1, -4, 1, -6, -7, -15, 1, -16, -8, -2, 6, -5, -2, -17, -3, -7
Offset: 1

Views

Author

Gus Wiseman, Aug 04 2021

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.

Examples

			For n = 1050 we have 5 prime indices {1,2,3,3,4}, so a(1050) = 10 - 13 = -3.
		

Crossrefs

First appearances are the elements of A174090 except for 3.
Positions of zeros are A340387.
Positions of nonpositive terms are A344291.
Positions of nonnegative terms are A344296.
A001222 counts prime indices.
A056239 adds up prime indices.
A112798 lists prime indices.

Programs

  • Mathematica
    Table[2*PrimeOmega[n]-Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimePi[p]]],{n,100}]

Formula

Totally additive with a(prime(k)) = 2 - k.

A352747 Array read by ascending antidiagonals. A(n, k) = F(k, n) mod n for n >= 1 and k >= 0, where F(n, k) = A352744(n, k) are the Fibonacci numbers, A(0, k) = 1 for k >= 0.

Original entry on oeis.org

1, 0, 1, 1, 0, 1, 1, 0, 0, 1, 2, 0, 1, 0, 1, 3, 1, 2, 0, 0, 1, 5, 3, 0, 1, 1, 0, 1, 1, 1, 3, 3, 0, 0, 0, 1, 5, 0, 3, 3, 2, 2, 1, 0, 1, 3, 2, 6, 5, 3, 1, 1, 0, 0, 1, 4, 1, 7, 5, 1, 3, 0, 0, 1, 0, 1, 0, 9, 8, 4, 4, 3, 3, 3, 2, 0, 0, 1, 5, 1, 4, 6, 1, 3, 5, 3, 2, 1, 1, 0, 1
Offset: 0

Views

Author

Peter Luschny, Apr 08 2022

Keywords

Comments

This array aims the study of the divisibility properties of the Fibonacci numbers A352744. The identity F(n, k) = (-1)^k*F(1 - n, -k) from A352744 shows that negative indices do not add to the divisibility properties of F(n, k).
All rows A(n, .) are pure periodic sequences. The length of the periods is given by (1, A270313). For n > 0 the length of the period of row A(n, .) is <= n.
The period length is 1 for n in (1, A023172) and n for n in (1, A074215), as observed by Robert Israel in A270313. In particular, if n is a power of 2 or a prime (A174090), then the period length is n.
The indices of the zero-free rows are in A353280. A zero-free row A(n, .) means that n will not divide F(k, n) whatever value k takes. For that it is sufficient to check that period(A(n, .)) is zero-free.
If period(A(n, .)) = [k | 0 <= k < n] we call n a 'Fibonacci friend'. In other words, in this case F(k, n) mod n = k for 0 <= k < n. A Fibonacci friend does not have to be prime (since 1 is a Fibonacci friend), but if it is prime then it is congruent to {1, 4} mod 5 (A045468), and all such primes are Fibonacci friends.
To say that n is a Fibonacci friend is equivalent to saying that A(n, n) = 0 and that n divides F(n, n). Fibonacci friends are the indices of the zeros in A002752.
Integers n > 0 that divide Sum{k=0..n-1} (F(k, n) mod n) are congruent to {0, 1, 3, 5} mod 6 (A301729).

Examples

			Array starts (periods are indicated with () ):
[n\k] 0   1   2   3   4  5  6   7   8   9  10  11  12
----------------------------------------------------------
[ 0] (1), 1,  1,  1,  1, 1, 1,  1,  1,  1,  1,  1,  1, ...
[ 1] (0), 0,  0,  0,  0, 0, 0,  0,  0,  0,  0,  0,  0, ...
[ 2] (1,  0), 1,  0,  1, 0, 1,  0,  1,  0,  1,  0,  1, ...
[ 3] (1,  0,  2), 1,  0, 2, 1,  0,  2,  1,  0,  2,  1, ...
[ 4] (2,  1,  0,  3), 2, 1, 0,  3,  2,  1,  0,  3,  2, ...
[ 5] (3), 3,  3,  3,  3, 3, 3,  3,  3,  3,  3,  3,  3, ...
[ 6] (5,  1,  3), 5,  1, 3, 5,  1,  3,  5,  1,  3,  5, ...
[ 7] (1,  0,  6,  5,  4, 3, 2), 1,  0,  6,  5,  4,  3, ...
[ 8] (5,  2,  7,  4,  1, 6, 3,  0), 5,  2,  7,  4,  1, ...
[ 9] (3,  1,  8,  6,  4, 2, 0,  7,  5), 3,  1,  8,  6, ...
[10] (4,  9), 4,  9,  4, 9, 4,  9,  4,  9,  4,  9,  4, ...
[11] (0,  1,  2,  3,  4, 5, 6,  7,  8,  9, 10), 0,  1, ...
[12] (5), 5,  5,  5,  5, 5, 5,  5,  5,  5,  5,  5,  5, ...
		

Crossrefs

Programs

  • Maple
    f := n -> combinat:-fibonacci(n + 1):
    F := proc(n, k) option remember; (n-1)*f(k-1) + f(k) end:
    A := (n, k) -> ifelse(n = 0, 1, modp(F(k, n), n)):
    for n from 0 to 12 do seq(A(n, k), k = 0..10) od;
  • Mathematica
    F[n_, k_] := (n - 1)*Fibonacci[k] + Fibonacci[k + 1];
    A[n_, k_] := If[n == 0, 1, Mod[F[k, n], n]];
    Table[A[n, k], {n, 0, 12}, {k, 0, 10}] // TableForm
  • SageMath
    def F(n, k): return (n - 1)*fibonacci(k) + fibonacci(k + 1)
    def A(n,k): return mod(F(k, n), n)
    for n in range(13): print([A(n,k) for k in range(13)])

Formula

A(n, 0) = A(n, n) = A002752(n).
Clearly 0 <= A(n, k) < n for all k and n > 0.

A325223 Sum of the prime indices of n minus the greater of the number of prime factors of n counted with multiplicity and the largest prime index of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 2, 1, 0, 1, 0, 1, 2, 0, 0, 2, 0, 2, 2, 1, 0, 1, 3, 1, 3, 2, 0, 3, 0, 0, 2, 1, 3, 2, 0, 1, 2, 2, 0, 3, 0, 2, 4, 1, 0, 1, 4, 4, 2, 2, 0, 3, 3, 3, 2, 1, 0, 3, 0, 1, 4, 0, 3, 3, 0, 2, 2, 4, 0, 2, 0, 1, 5, 2, 4, 3, 0, 2, 4, 1, 0, 4, 3, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Apr 12 2019

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n.
Also the number of squares in the Young diagram of the integer partition with Heinz number n after the first row or the first column, whichever is larger, is removed. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			88 has 4 prime indices {1,1,1,5} with sum 8 and maximum 5, so a(88) = 8 - max(4,5) = 3.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[primeMS[n]]-Max[Length[primeMS[n]],Max[primeMS[n]]],{n,100}]

Formula

a(n) = A056239(n) - max(A001222(n), A061395(n)) = A056239(n) - A263297(n).

A383961 Square array read by upward antidiagonals: T(n,k) is the n-th number whose largest odd divisor is its k-th divisor, n >= 1, k >= 1.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 8, 7, 9, 15, 16, 11, 10, 20, 18, 32, 13, 12, 21, 50, 36, 64, 17, 14, 27, 81, 45, 30, 128, 19, 22, 28, 88, 63, 42, 105, 256, 23, 24, 33, 98, 75, 54, 135, 60, 512, 29, 25, 35, 104, 99, 66, 165, 84, 120, 1024, 31, 26, 39, 136, 117, 70, 189, 108, 140, 90
Offset: 1

Views

Author

Omar E. Pol, May 16 2025

Keywords

Comments

This is a permutation of the positive integers.
From Peter Munn, May 18 2025: (Start)
Numbers with the same factorization pattern of their sequence of divisors (see A290110) and the same parity appear here in the same column.
For example, each column k > 2 includes the subsequence 2^(k-2) * p for all prime p > 2^(k-2).
(End)

Examples

			The corner 15 X 15 of the square array is as follows:
      1,  3,  6,  15,  18,  36,  30, 105,  60, 120,  90, 315,  816, 1360, 180, ...
      2,  5,  9,  20,  50,  45,  42, 135,  84, 140, 126, 324,  880, 1520, 210, ...
      4,  7, 10,  21,  81,  63,  54, 165, 108, 168, 150, 432,  912, 1632, 252, ...
      8, 11, 12,  27,  88,  75,  66, 189, 132, 220, 198, 440, 1040, 1760, 270, ...
     16, 13, 14,  28,  98,  99,  70, 195, 156, 240, 216, 495, 1056, 1824, 300, ...
     32, 17, 22,  33, 104, 117,  72, 200, 162, 260, 234, 520, 1104, 1840, 330, ...
     64, 19, 24,  35, 136, 147,  78, 231, 204, 308, 264, 525, 1120, 1904, 378, ...
    128, 23, 25,  39, 152, 153, 100, 255, 225, 340, 280, 528, 1144, 2000, 390, ...
    256, 29, 26,  40, 176, 171, 102, 273, 228, 364, 294, 560, 1232, 2080, 396, ...
    512, 31, 34,  44, 184, 175, 110, 285, 276, 380, 306, 585, 1248, 2128, 462, ...
   1024, 37, 38,  51, 208, 207, 114, 297, 348, 405, 312, 616, 1392, 2208, 468, ...
   2048, 41, 46,  52, 232, 243, 130, 345, 372, 460, 336, 624, 1456, 2288, 510, ...
   4096, 43, 48,  55, 242, 245, 138, 351, 400, 476, 342, 675, 1458, 2320, 546, ...
   8192, 47, 49,  56, 248, 261, 144, 357, 441, 480, 350, 680, 1488, 2464, 570, ...
  16384, 53, 58,  57, 296, 272, 154, 375, 444, 500, 408, 693, 1496, 2480, 588, ...
  ...
		

Crossrefs

Column 1 gives A000079.
Column 2 gives A065091.
Column 3 consists of (A001248 U A091629 U A100484)\{4}.
Column 4 consists of numbers >= 15 in (A001749 U A030078 U A046388 U A070875).
Row 1 gives A383402.

Programs

  • Mathematica
    f[n_] := If[OddQ[n], DivisorSigma[0, n], FirstPosition[Divisors[n], n/2^IntegerExponent[n, 2]][[1]]]; seq[m_] := Module[{t = Table[0, {m}, {m}], v = Table[0, {m}], c = 0, k = 1, i, j}, While[c < m*(m + 1)/2, i = f[k]; If[i <= m, j = v[[i]] + 1; If[j <= m - i + 1, t[[i]][[j]] = k; v[[i]]++; c++]]; k++]; Table[t[[j]][[i - j + 1]], {i, 1, m}, {j, 1, i}] // Flatten]; seq[11] (* Amiram Eldar, May 16 2025 *)

A067133 n is a term if the phi(n) numbers in [0,n-1] and coprime to n form an arithmetic progression.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 11, 13, 16, 17, 19, 23, 29, 31, 32, 37, 41, 43, 47, 53, 59, 61, 64, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241
Offset: 1

Views

Author

Amarnath Murthy, Jan 09 2002

Keywords

Comments

The sequence consists of primes, powers of 2 and 6. Sketch of proof: Let k be the common difference of the arithmetic progression. If n is odd, then 1 and 2 are coprime to n, so k=1 and n is prime. If n==0 (mod 4), then n/2-1 and n/2+1 are coprime to n, so k=2 and n is a power of 2. If n==2 (mod 4), then n/2-2 and n/2+2 are coprime to n, so k divides 4 and n is either 2 or 6.
From Bernard Schott, Jan 08 2021: (Start)
This sequence is the answer to the 2nd problem, proposed by Romania, during the 32nd International Mathematical Olympiad in 1991 at Sigtuna (Sweden) (see the link IMO Compendium and reference Kuczma).
These phi(m) numbers coprimes to m form an arithmetic progression with at least 3 terms iff m = 5 or m >= 7. (End)

Examples

			8 is a term as phi(8) = 4 and the coprime numbers 1,3,5,7 form an arithmetic progression. 17 is a member as phi(17) = 16 and the numbers 1 to 16 form an arithmetic progression.
		

References

  • Marcin E. Kuczma, International Mathematical Olympiads, 1986-1999, The Mathematical Association of America, 2003, pages 6 and 61-62.

Crossrefs

Equals A000040 U A000079 U {6}.
Equals A174090 U {6}.

Programs

  • Mathematica
    rps[ n_ ] := Select[ Range[ 0, n-1 ], GCD[ #, n ]==1& ]; difs[ n_ ] := Drop[ n, 1 ]-Drop[ n, -1 ]; Select[ Range[ 1, 250 ], Length[ Union[ difs[ rps[ # ] ] ] ]<=1& ]
  • PARI
    isok(n) = {my(v = select(x->gcd(x, n)==1, [1..n]), dv = vector(#v-1, k, v[k+1] - v[k])); if (#dv, if (vecmin(dv) != vecmax(dv), return(0))); return(1)} \\ Michel Marcus, Jan 08 2021
    
  • Python
    from sympy import primepi
    def A067133(n):
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x): return int(n+x-(x>5)+(0 if x<=1 else 1-primepi(x))-x.bit_length())
        return bisection(f,n,n) # Chai Wah Wu, Sep 19 2024

Extensions

Edited by Dean Hickerson, Jan 15 2002

A270306 Smallest number which is representable as a sum of n but no fewer consecutive nonnegative integers, or 0 for n > 1 if no such number exists.

Original entry on oeis.org

0, 1, 6, 10, 20, 0, 28, 44, 0, 0, 88, 0, 104, 0, 0, 136, 272, 0, 304, 0, 0, 0, 368, 0, 0, 0, 0, 0, 464, 0, 496, 592, 0, 0, 0, 0, 1184, 0, 0, 0, 1312, 0, 1376, 0, 0, 0, 1504, 0, 0, 0, 0, 0, 1696, 0, 0, 0, 0, 0, 1888, 0, 1952, 0, 0, 2144, 0, 0, 4288, 0, 0, 0, 4544
Offset: 1

Views

Author

Martin Renner, Mar 14 2016

Keywords

Examples

			a(2) = 1 = 0 + 1 (two but no fewer terms)
a(3) = 6 = 1 + 2 + 3 (three but no fewer terms)
a(4) = 10 = 1 + 2 + 3 + 4 (four but no fewer terms)
a(5) = 20 = 2 + 3 + 4 + 5 + 6 (five but no fewer terms)
a(6) = 0, since no number needs six but no fewer terms
a(7) = 28 = 1 + 2 + 3 + 4 + 5 + 6 + 7 (seven but no fewer terms)
		

Crossrefs

Formula

a(A111774(n)) = 0.
a(A174090(n)) = A163172(n).

A337939 Irregular triangle T(n, m) read by rows: row n gives the distinct length ratios diagonal/side of regular n-gons, DSR(n, k), for n >= 2, k = 1, 2, ..., floor(n/2), expressed by the coefficients in the power basis of the Galois group Gal(Q(rho(n))/Q), where rho(n) = 2*cos(Pi/n), for n >= 2. T(1, 1) is set to 1.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 2, 1, 0, 1, -1, 0, 1, 1, 0, 1, -1, 0, 1, 0, -2, 0, 1, 1, 0, 1, -1, 0, 1, 1, 1, 1, 0, 1, -1, 0, 1, 0, -2, 0, 1, -4, 0, 2, 1, 0, 1, -1, 0, 1, 0, -2, 0, 1, 1, 0, -3, 0, 1, 1, 0, 1, -1, 0, 1, 0, -2, 0, 1, 0, 0, 1, 0, 2
Offset: 1

Views

Author

Wolfdieter Lang, Jan 15 2021

Keywords

Comments

The length of row n is given in A338431(n), for n >= 1.
The length of the sublists t(n, k) of the power basis coefficients of DSR(n, k), for k = 1, 2, ..., floor(n/2), is 1 if n = 1, for n >= 2 it is k except for n = n(j) = A111774(j) for which the final A219839(n) sublists have fewer than k members.
Trailing vanishing coefficients of the delta(n) = A055034(n) power base elements <1 = rho(n)^0, rho(n)^1, ..., rho(n)^{delta(n)-1}> are not recorded. The coefficients of the minimal polynomial C(n, x) of rho(n) = 2*cos(Pi/n) of degree delta(n) are given in A187360. C(n, rho(n)) = 0 is used to eliminate all powers of rho with exponent >= delta(n).
The length ratios DSR(n, k) := diagonal(n, k)/side(n) of regular n-gons, for n >= 2, and k = 1, 2, ..., floor(n/2) (distinct diagonals, starting with the side for k = 1, in increasing order) are given by DSR(n, k) = S(k-1, rho(n)), with the Chebyshev S polynomials (A049310). See the W. Lang link.
For n = 2, the degenerate case, diagonal/side = side/side = 1 for k = 1. For n = 1 (a point) diagonal/side is undetermined, and T(1, 1) is set to 1.
For the power basis sublists t(n, k), for k = 1, 2, ..., delta(n), only the k coefficients of S(k-1, x) are present (trailing vanishing coefficients are not recorded). For k = delta(n)+1, ..., floor(n/2) less than k coefficients appear due to elimination via C(n, rho(n)) = 0. E.g., for n = 6 with delta(6) = 2 the only coefficient for k = 3 is 2 (coefficient of rho^0). This appears for n = n(j) = A111774(j), because then floor(n/2) - delta(n) = A219839(n) > 0.
Because A219839(n) = 0 means that n is from A174090, i.e., a prime or a power of 2 (complement of A111774), these rows n have all the sublists t(n, k) with the k coefficients of S(k-1, x), hence they are identical (but the basis differs). See especially the table for the pairs of consecutive numbers n with identical coefficients, like (2, 3), (4, 5), (16, 17), (256, 267), (65536, 65537), ?... (cf. Fermat primes A019434).

Examples

			The irregular triangle T(n, m) begins: (For n >= 4 the bar divides the DSR(n, k) power basis coefficients, the sublists t(n, k), for k = 1, 2, ..., floor(n/2))
n \ m  1   2 3    4  5  6   7  8 9 10   11 12 12  13 14  15 16 17 18 19 20 ...
1:     1
2:     1
3:     1
4:     1 | 0 1
5:     1 | 0 1
6:     1 | 0 1 |  2
7:     1 | 0 1 | -1  0  1
8:     1 | 0 1 | -1  0  1 | 0 -2 0  1
9:     1 | 0 1 | -1  0  1 | 1  1
10:    1 | 0 1 | -1  0  1 | 0 -2 0  1 | -4  0  2
11:    1 | 0 1 | -1  0  1 | 0 -2 0  1 |  1  0 -3   0  1
12:    1 | 0 1 | -1  0  1 | 0 -2 0  1 |  0  0  1 | 0  2
13:    1 | 0 1 | -1  0  1 | 0 -2 0  1 |  1  0 -3   0  1 | 0  3  0 -4  0  1
...
n = 14: 1 | 0 1 | -1 0 1 | 0 -2 0 1 | 1 0 -3 0 1 | 0 3 0 -4 0 1 | 6 0 -8 0 2,
n = 15: 1 | 0 1 | -1 0 1 | 0 -2 0 1 | 0 4 1 -1 | 1 -2 0 1 | -1 1 1,
n = 16 and n = 17: 1 | 0 1 | -1  0 1 | 0 -2 0 1 | 1 0 -3 0 1 | 0 3 0 -4 0 1 | -1 0 6 0 -5 0 1 | 0 -4 0 10 0 -6 0 1,
--------------------------------------------------------------------------------
n = 5: DSR(5, 1) = 1 = side(5)/side(5), DSR(5, 2) = 1*rho(5) = A001622 (golden section).
n = 8: DSR(8, 1) = 1 = side(8)/side(8), DSR(8, 2) = 1*rho(8) = sqrt(2+sqrt(2)) = A179260, DSR(8, 3) = -1 + rho(8)^2 = 1 + sqrt(2) = A014176, DSR(8, 4) = -2*rho(8) + 1*rho(8)^3 = sqrt(2)*rho(8) = A121601.
		

Crossrefs

Formula

T(1, 1) = 1, and in row n, for n >= 2, the power base coefficients of Gal(Q(2*cos(Pi/n))/Q) for DSR(n, k) := diagonal(n, k)/side(n) of regular n-gons, for k = 1, 2, ..., floor(n/2), are listed as t(n, k) in this order, with trailing vanishing coefficients omitted.
Previous Showing 21-30 of 33 results. Next