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 61-70 of 144 results. Next

A342529 Number of compositions of n with distinct first quotients.

Original entry on oeis.org

1, 1, 2, 3, 7, 13, 19, 36, 67, 114, 197, 322, 564, 976, 1614, 2729, 4444, 7364, 12357, 20231, 33147, 53973, 87254, 140861, 227535, 368050, 589706, 940999, 1497912, 2378260, 3774297, 5964712, 9416411, 14822087, 23244440, 36420756
Offset: 0

Views

Author

Gus Wiseman, Mar 17 2021

Keywords

Comments

The first quotients of a sequence are defined as if the sequence were an increasing divisor chain, so for example the first quotients of (6,3,1) are (1/2,1/3).

Examples

			The composition (2,1,2,3) has first quotients (1/2,2,3/2) so is counted under a(8).
The a(1) = 1 through a(5) = 13 compositions:
  (1)  (2)    (3)    (4)      (5)
       (1,1)  (1,2)  (1,3)    (1,4)
              (2,1)  (2,2)    (2,3)
                     (3,1)    (3,2)
                     (1,1,2)  (4,1)
                     (1,2,1)  (1,1,3)
                     (2,1,1)  (1,2,2)
                              (1,3,1)
                              (2,1,2)
                              (2,2,1)
                              (3,1,1)
                              (1,1,2,1)
                              (1,2,1,1)
		

Crossrefs

The version for differences instead of quotients is A325545.
The version for equal first quotients is A342495.
The unordered version is A342514, ranked by A342521.
The strict unordered version is A342520.
A000005 counts constant compositions.
A000009 counts strictly increasing (or strictly decreasing) compositions.
A000041 counts weakly increasing (or weakly decreasing) compositions.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A003238 counts chains of divisors summing to n - 1 (strict: A122651).
A167865 counts strict chains of divisors > 1 summing to n.

Programs

  • Mathematica
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],UnsameQ@@Divide@@@Partition[#,2,1]&]],{n,0,15}]

Extensions

a(21)-a(35) from Alois P. Heinz, Jan 16 2025

A056503 Number of periodic palindromic structures of length n using a maximum of two different symbols.

Original entry on oeis.org

1, 2, 2, 4, 4, 7, 8, 14, 16, 26, 32, 51, 64, 100, 128, 198, 256, 392, 512, 778, 1024, 1552, 2048, 3091, 4096, 6176, 8192, 12324, 16384, 24640, 32768, 49222, 65536, 98432, 131072, 196744, 262144, 393472, 524288, 786698, 1048576, 1573376, 2097152, 3146256, 4194304
Offset: 1

Views

Author

Keywords

Comments

For example, aaabbb is not a (finite) palindrome but it is a periodic palindrome. Permuting the symbols will not change the structure.
A periodic palindrome is just a necklace that is equivalent to its reverse. The number of binary periodic palindromes of length n is given by A164090(n). A binary periodic palindrome can only be equivalent to its complement when there are an equal number of 0's and 1's. - Andrew Howroyd, Sep 29 2017
Number of cyclic compositions (necklaces of positive integers) summing to n that can be rotated to form a palindrome. - Gus Wiseman, Sep 16 2018

Examples

			From _Gus Wiseman_, Sep 16 2018: (Start)
The sequence of palindromic cyclic compositions begins:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (111)  (22)    (113)    (33)      (115)
                    (112)   (122)    (114)     (133)
                    (1111)  (11111)  (222)     (223)
                                     (1122)    (11113)
                                     (11112)   (11212)
                                     (111111)  (11122)
                                               (1111111)
(End)
		

References

  • M. R. Nester (1999). Mathematical investigations of some plant interaction designs. PhD Thesis. University of Queensland, Brisbane, Australia. [See A056391 for pdf file of Chap. 2]

Crossrefs

Programs

  • Mathematica
    (* b = A164090, c = A045674 *)
    b[n_] := (1/4)*(7 - (-1)^n)*2^((1/4)*(2*n + (-1)^n - 1));
    c[0] = 1; c[n_] := c[n] = If[EvenQ[n], 2^(n/2-1) + c[n/2], 2^((n-1)/2)];
    a[n_?OddQ] := b[n]/2; a[n_?EvenQ] := (1/2)*(b[n] + c[n/2]);
    Array[a, 45] (* Jean-François Alcover, Oct 08 2017, after Andrew Howroyd *)
    Table[Length[Select[Join@@Permutations/@IntegerPartitions[n],Function[q,And[Array[OrderedQ[{q,RotateRight[q,#]}]&,Length[q]-1,1,And],Array[SameQ[RotateRight[q,#],Reverse[RotateRight[q,#]]]&,Length[q],1,Or]]]]],{n,15}] (* Gus Wiseman, Sep 16 2018 *)

Formula

a(2n+1) = A164090(2n+1)/2 = 2^n, a(2n) = (A164090(2n) + A045674(n))/2. - Andrew Howroyd, Sep 29 2017

Extensions

a(17)-a(45) from Andrew Howroyd, Apr 07 2017

A275527 Number of distinct classes of permutations of length n under reversal and complement to n+1.

Original entry on oeis.org

1, 1, 1, 4, 12, 64, 360, 2544, 20160, 181632
Offset: 1

Views

Author

Olivier Gérard, Jul 31 2016

Keywords

Comments

Let us consider two permutations to be equivalent if they can be obtained from each other by cyclic rotation (12345->(23451,34512,45123,51234) or n+1-complement (31254->35412), or a combination of those two transformations (they commute with each other). a(n) is the number of classes.
We obtain the same number of classes if the transformations are (addition of a constant modulo n and reversal (12345->54321)) but not the same set of representatives.
It seems probable that a(2n+1) = (2n)!/2
This sequence may be related to A113247 (and A113248) as they share a common dissection 1, 4, 64, 2544, 181632. The fact that they count permutation classes for the major index is a further indication.
Number of path necklaces, defined as equivalence classes of (labeled, undirected) Hamiltonian paths under rotation of the vertices. The cycle version is A000939. - Gus Wiseman, Mar 02 2019

Examples

			Examples of permutation representatives. The representative is chosen to be the first of the class in lexicographic order.
n=4 case addition mod n and reversal
1234, 1243, 1324, 1423.
n=4 case rotation and complement
1234, 1243, 1324, 1342.
.
n=5 case addition mod n and reversal
12345, 12354, 12435, 12453, 12534, 13245, 13425, 13452, 13524, 14235, 14523, 15234.
n=5 case rotation and complement
12345, 12354, 12435, 12453, 12534, 13245, 13425, 13452, 13524, 14235, 14325, 14352.
		

Crossrefs

Cf. A000939, A000940, A002619, A089066, A262480 (other symmetry classes of permutations).
Cf. A193651 (inspiration for a(2n)).

Programs

  • Mathematica
    rotgra[g_,m_]:=Sort[Sort/@(g/.k_Integer:>If[k==m,1,k+1])];
    Table[Length[Select[Union[Sort[Sort/@Partition[#,2,1]]&/@Permutations[Range[n]]],#==First[Sort[Table[Nest[rotgra[#,n]&,#,j],{j,n}]]]&]],{n,8}] (* Gus Wiseman, Mar 02 2019 *)

Formula

(Conjecture). If n odd a(n)=((n - 1))!/2. If n even a(n)= 1/2 (n - 2)!! (1 + ( n - 1)!!).

A320778 Inverse Euler transform of the Euler totient function phi = A000010.

Original entry on oeis.org

1, 1, 0, 1, 0, 2, -3, 4, -4, 4, -9, 14, -19, 30, -42, 50, -76, 128, -194, 286, -412, 598, -909, 1386, -2100, 3178, -4763, 7122, -10758, 16414, -25061, 38056, -57643, 87568, -133436, 203618, -311128, 475536, -726355, 1109718, -1697766, 2601166, -3987903, 6114666
Offset: 0

Views

Author

Gus Wiseman, Oct 22 2018

Keywords

Comments

The Euler transform of a sequence q is the sequence of coefficients of x^n, n > 0, in the expansion of Product_{n > 0} 1/(1 - x^n)^q(n). The constant term 1 is sometimes taken to be the zeroth part of the Euler transform.

Crossrefs

Number theoretical functions: A000005, A000010, A000203, A001055, A001221, A001222, A008683, A010054.
Inverse Euler transforms: A059966, A320767, A320776, A320777, A320779, A320780, A320781, A320782.

Programs

  • Maple
    # The function EulerInvTransform is defined in A358451.
    a := EulerInvTransform(n -> ifelse(n=0, 1, NumberTheory:-Totient(n))):
    seq(a(n), n = 0..43); # Peter Luschny, Nov 21 2022
  • Mathematica
    EulerInvTransform[{}]={};EulerInvTransform[seq_]:=Module[{final={}},For[i=1,i<=Length[seq],i++,AppendTo[final,i*seq[[i]]-Sum[final[[d]]*seq[[i-d]],{d,i-1}]]];
    Table[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
    EulerInvTransform[Array[EulerPhi,30]]

A323866 Number of aperiodic toroidal necklaces of positive integers summing to n.

Original entry on oeis.org

1, 1, 1, 3, 5, 12, 18, 42, 72, 145, 262, 522, 960, 1879, 3531, 6831, 13013, 25148, 48177, 93186, 179507, 347509, 671955, 1303257, 2527162, 4910681, 9545176, 18579471, 36183505, 70540861, 137603801, 268655547, 524842088, 1026067205, 2007118657, 3928564113
Offset: 0

Views

Author

Gus Wiseman, Feb 04 2019

Keywords

Comments

The 1-dimensional (Lyndon word) case is A059966.
We define a toroidal necklace to be an equivalence class of matrices under all possible rotations of the sequence of rows and the sequence of columns. An n X k matrix is aperiodic if all n * k rotations of its sequence of rows and its sequence of columns are distinct.

Examples

			Inequivalent representatives of the a(6) = 18 toroidal necklaces:
  [6] [1 5] [2 4] [1 1 4] [1 2 3] [1 3 2] [1 1 1 3] [1 1 2 2] [1 1 1 1 2]
.
  [1] [2] [1 1]
  [5] [4] [1 3]
.
  [1] [1] [1]
  [1] [2] [3]
  [4] [3] [2]
.
  [1] [1]
  [1] [1]
  [1] [2]
  [3] [2]
.
  [1]
  [1]
  [1]
  [1]
  [2]
		

Crossrefs

Programs

  • GAP
    List([0..30], A323866); # See A323861 for code; Andrew Howroyd, Aug 21 2019
  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    ptnmats[n_]:=Union@@Permutations/@Select[Union@@(Tuples[Permutations/@#]&/@Map[primeMS,facs[n],{2}]),SameQ@@Length/@#&];
    apermatQ[m_]:=UnsameQ@@Join@@Table[RotateLeft[m,{i,j}],{i,Length[m]},{j,Length[First[m]]}];
    neckmatQ[m_]:=m==First[Union@@Table[RotateLeft[m,{i,j}],{i,Length[m]},{j,Length[First[m]]}]];
    Table[If[n==0,1,Length[Union@@Table[Select[ptnmats[k],And[apermatQ[#],neckmatQ[#]]&],{k,Times@@Prime/@#&/@IntegerPartitions[n]}]]],{n,0,10}]

Extensions

Terms a(16) and beyond from Andrew Howroyd, Aug 21 2019

A329316 Irregular triangle read by rows where row n gives the sequence of lengths of components of the Lyndon factorization of the reversed first n terms of A000002.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 3, 1, 4, 1, 1, 4, 1, 2, 4, 1, 1, 2, 4, 1, 1, 1, 2, 4, 1, 3, 2, 4, 1, 1, 3, 2, 4, 1, 1, 1, 3, 2, 4, 1, 3, 3, 2, 4, 1, 9, 4, 1, 1, 9, 4, 1, 2, 9, 4, 1, 16, 1, 1, 16, 1, 1, 1, 16, 1, 3, 16, 1, 1, 3, 16, 1, 5, 16, 1, 6, 16, 1, 1, 6, 16, 1, 2, 6
Offset: 0

Views

Author

Gus Wiseman, Nov 11 2019

Keywords

Comments

There are no repeated rows, as row n has sum n.
We define the Lyndon product of two or more finite sequences to be the lexicographically maximal sequence obtainable by shuffling the sequences together. For example, the Lyndon product of (231) with (213) is (232131), the product of (221) with (213) is (222131), and the product of (122) with (2121) is (2122121). A Lyndon word is a finite sequence that is prime with respect to the Lyndon product. Equivalently, a Lyndon word is a finite sequence that is lexicographically strictly less than all of its cyclic rotations. Every finite sequence has a unique (orderless) factorization into Lyndon words, and if these factors are arranged in lexicographically decreasing order, their concatenation is equal to their Lyndon product. For example, (1001) has sorted Lyndon factorization (001)(1).
It appears that some numbers (such as 10) never appear in the sequence.

Examples

			Triangle begins:
   1: (1)
   2: (1,1)
   3: (1,1,1)
   4: (3,1)
   5: (4,1)
   6: (1,4,1)
   7: (2,4,1)
   8: (1,2,4,1)
   9: (1,1,2,4,1)
  10: (3,2,4,1)
  11: (1,3,2,4,1)
  12: (1,1,3,2,4,1)
  13: (3,3,2,4,1)
  14: (9,4,1)
  15: (1,9,4,1)
  16: (2,9,4,1)
  17: (16,1)
  18: (1,16,1)
  19: (1,1,16,1)
  20: (3,16,1)
For example, the reversed first 13 terms of A000002 are (1221221211221), with Lyndon factorization (122)(122)(12)(1122)(1), so row 13 is (3,3,2,4,1).
		

Crossrefs

Row lengths are A329317.
The non-reversed version is A329315.

Programs

  • Mathematica
    lynQ[q_]:=Array[Union[{q,RotateRight[q,#]}]=={q,RotateRight[q,#]}&,Length[q]-1,1,And];
    lynfac[q_]:=If[Length[q]==0,{},Function[i,Prepend[lynfac[Drop[q,i]],Take[q,i]]][Last[Select[Range[Length[q]],lynQ[Take[q,#]]&]]]];
    kolagrow[q_]:=If[Length[q]<2,Take[{1,2},Length[q]+1],Append[q,Switch[{q[[Length[Split[q]]]],q[[-2]],Last[q]},{1,1,1},0,{1,1,2},1,{1,2,1},2,{1,2,2},0,{2,1,1},2,{2,1,2},2,{2,2,1},1,{2,2,2},1]]]
    kol[n_Integer]:=Nest[kolagrow,{1},n-1];
    Table[Length/@lynfac[Reverse[kol[n]]],{n,100}]

A334265 Numbers k such that the k-th composition in standard order is a reversed Lyndon word.

Original entry on oeis.org

0, 1, 2, 4, 5, 8, 9, 11, 16, 17, 18, 19, 21, 23, 32, 33, 34, 35, 37, 39, 41, 43, 47, 64, 65, 66, 67, 68, 69, 71, 73, 74, 75, 77, 79, 81, 83, 85, 87, 91, 95, 128, 129, 130, 131, 132, 133, 135, 137, 138, 139, 141, 143, 145, 146, 147, 149, 151, 155, 159, 161, 163
Offset: 1

Views

Author

Gus Wiseman, Apr 22 2020

Keywords

Comments

Reversed Lyndon words are different from co-Lyndon words (A326774).
A Lyndon word is a finite sequence of positive integers that is lexicographically strictly less than all of its cyclic rotations.
The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The sequence of all reversed Lyndon words begins:
    0: ()            37: (3,2,1)         83: (2,3,1,1)
    1: (1)           39: (3,1,1,1)       85: (2,2,2,1)
    2: (2)           41: (2,3,1)         87: (2,2,1,1,1)
    4: (3)           43: (2,2,1,1)       91: (2,1,2,1,1)
    5: (2,1)         47: (2,1,1,1,1)     95: (2,1,1,1,1,1)
    8: (4)           64: (7)            128: (8)
    9: (3,1)         65: (6,1)          129: (7,1)
   11: (2,1,1)       66: (5,2)          130: (6,2)
   16: (5)           67: (5,1,1)        131: (6,1,1)
   17: (4,1)         68: (4,3)          132: (5,3)
   18: (3,2)         69: (4,2,1)        133: (5,2,1)
   19: (3,1,1)       71: (4,1,1,1)      135: (5,1,1,1)
   21: (2,2,1)       73: (3,3,1)        137: (4,3,1)
   23: (2,1,1,1)     74: (3,2,2)        138: (4,2,2)
   32: (6)           75: (3,2,1,1)      139: (4,2,1,1)
   33: (5,1)         77: (3,1,2,1)      141: (4,1,2,1)
   34: (4,2)         79: (3,1,1,1,1)    143: (4,1,1,1,1)
   35: (4,1,1)       81: (2,4,1)        145: (3,4,1)
		

Crossrefs

The non-reversed version is A275692.
The generalization to necklaces is A333943.
The dual version (reversed co-Lyndon words) is A328596.
The case that is also co-Lyndon is A334266.
Binary Lyndon words are counted by A001037.
Lyndon compositions are counted by A059966.
Normal Lyndon words are counted by A060223.
Numbers whose prime signature is a reversed Lyndon word are A334298.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Reverse is A228351 (triangle).
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Reversed Lyndon words are A334265 (this sequence).
- Co-Lyndon words are A326774.
- Reversed co-Lyndon words are A328596.
- Length of Lyndon factorization is A329312.
- Distinct rotations are counted by A333632.
- Lyndon factorizations are counted by A333940.
- Length of Lyndon factorization of reverse is A334297.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    lynQ[q_]:=Length[q]==0||Array[Union[{q,RotateRight[q,#1]}]=={q,RotateRight[q,#1]}&,Length[q]-1,1,And];
    Select[Range[0,100],lynQ[Reverse[stc[#]]]&]

A320767 Inverse Euler transform applied once to {1,-1,0,0,0,...}, twice to {1,0,0,0,0,...}, or three times to {1,1,1,1,1,...}.

Original entry on oeis.org

1, 1, -2, 1, -1, 2, -3, 4, -5, 8, -13, 18, -25, 40, -62, 90, -135, 210, -324, 492, -750, 1164, -1809, 2786, -4305, 6710, -10460, 16264, -25350, 39650, -62057, 97108, -152145, 238818, -375165, 589520, -927200, 1459960, -2300346, 3626200, -5720274, 9030450
Offset: 0

Views

Author

Gus Wiseman, Oct 20 2018

Keywords

Comments

The Euler transform of a sequence q is the sequence of coefficients of x^n, n > 0, in the expansion of Product_{n > 0} 1/(1 - x^n)^q(n). The constant term 1 is sometimes taken to be the zeroth part of the Euler transform.

Crossrefs

Programs

  • Mathematica
    EulerInvTransform[{}]={};EulerInvTransform[seq_]:=Module[{final={}},For[i=1,i<=Length[seq],i++,AppendTo[final,i*seq[[i]]-Sum[final[[d]]*seq[[i-d]],{d,i-1}]]];
    Table[Sum[MoebiusMu[i/d]*final[[d]],{d,Divisors[i]}]/i,{i,Length[seq]}]];
    Nest[EulerInvTransform,Array[DiscreteDelta,50,0],2]

A333939 Number of multisets of compositions that can be shuffled together to obtain the k-th composition in standard order.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 3, 1, 2, 2, 4, 2, 5, 4, 5, 1, 2, 2, 4, 2, 4, 5, 7, 2, 5, 4, 10, 4, 10, 7, 7, 1, 2, 2, 4, 2, 5, 5, 7, 2, 5, 3, 9, 5, 13, 11, 12, 2, 5, 5, 10, 5, 11, 13, 18, 4, 10, 9, 20, 7, 18, 12, 11, 1, 2, 2, 4, 2, 5, 5, 7, 2, 4, 4, 11, 5, 14, 11, 12, 2
Offset: 0

Views

Author

Gus Wiseman, Apr 15 2020

Keywords

Comments

Number of ways to deal out the k-th composition in standard order to form a multiset of hands.
A composition of n is a finite sequence of positive integers summing to n. The k-th composition in standard order (row k of A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			The dealings for n = 1, 3, 7, 11, 13, 23, 43:
  (1)  (11)    (111)      (211)      (121)      (2111)        (2211)
       (1)(1)  (1)(11)    (1)(21)    (1)(12)    (11)(21)      (11)(22)
               (1)(1)(1)  (2)(11)    (1)(21)    (1)(211)      (1)(221)
                          (1)(1)(2)  (2)(11)    (2)(111)      (21)(21)
                                     (1)(1)(2)  (1)(1)(21)    (2)(211)
                                                (1)(2)(11)    (1)(1)(22)
                                                (1)(1)(1)(2)  (1)(2)(21)
                                                              (2)(2)(11)
                                                              (1)(1)(2)(2)
		

Crossrefs

Multisets of compositions are counted by A034691.
Combinatory separations of normal multisets are counted by A269134.
Dealings with total sum n are counted by A292884.
Length of co-Lyndon factorization of binary expansion is A329312.
Length of Lyndon factorization of reversed binary expansion is A329313.
All of the following pertain to compositions in standard order (A066099):
- Length is A000120.
- Necklaces are A065609.
- Sum is A070939.
- Runs are counted by A124767.
- Rotational symmetries are counted by A138904.
- Strict compositions are A233564.
- Constant compositions are A272919.
- Lyndon words are A275692.
- Co-Lyndon words are A326774.
- Aperiodic compositions are A328594.
- Length of Lyndon factorization is A329312.
- Distinct rotations are counted by A333632.
- Co-Lyndon factorizations are counted by A333765.
- Lyndon factorizations are counted by A333940.
- Length of co-Lyndon factorization is A334029.
- Combinatory separations are A334030.

Programs

  • Mathematica
    nn=100;
    comps[0]:={{}};comps[n_]:=Join@@Table[Prepend[#,i]&/@comps[n-i],{i,n}];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    dealings[q_]:=Union[Function[ptn,Sort[q[[#]]&/@ptn]]/@sps[Range[Length[q]]]];
    Table[Length[dealings[stc[n]]],{n,0,nn}]

Formula

For n > 0, Sum_{k = 2^(n-1)..2^n-1} a(k) = A292884(n).

A065177 Table M(n,b) (columns: n >= 1, rows: b >= 0) gives the number of site swap juggling patterns with exact period n, using exactly b balls, where cyclic shifts are not counted as distinct.

Original entry on oeis.org

1, 0, 1, 0, 1, 1, 0, 2, 2, 1, 0, 3, 6, 3, 1, 0, 6, 15, 12, 4, 1, 0, 9, 42, 42, 20, 5, 1, 0, 18, 107, 156, 90, 30, 6, 1, 0, 30, 294, 554, 420, 165, 42, 7, 1, 0, 56, 780, 2028, 1910, 930, 273, 56, 8, 1, 0, 99, 2128, 7350, 8820, 5155, 1806, 420, 72, 9, 1, 0, 186, 5781, 26936
Offset: 0

Views

Author

Antti Karttunen, Oct 19 2001

Keywords

Examples

			Upper left corner starts as:
  1, 0,  0,   0,    0,     0,     0, ...
  1, 1,  2,   3,    6,     9,    18, ...
  1, 2,  6,  15,   42,   107,   294, ...
  1, 3, 12,  42,  156,   554,  2028, ...
  1, 4, 20,  90,  420,  1910,  8820, ...
  1, 5, 30, 165,  930,  5155, 28830, ...
  1, 6, 42, 273, 1806, 11809, 77658, ...
  ...
		

Crossrefs

Row 1: A059966, row 2: A065178, row 3: A065179, row 4: A065180.
Column 1: A002378, column 2: A059270.
Main diagonal gives A306173.
Cf. also A065167. trinv given at A054425.

Programs

  • Maple
    [seq(DistSS_table(j),j=0..119)]; DistSS_table := (n) -> DistSS((((trinv(n)-1)*(((1/2)*trinv(n))+1))-n)+1, (n-((trinv(n)*(trinv(n)-1))/2)));
    with(numtheory); DistSS := proc(n,b) local d,s; s := 0; for d in divisors(n) do s := s+mobius(n/d)*((b+1)^d - b^d); od; RETURN(s/n); end;
  • Mathematica
    trinv[n_] := Floor[(1 + Sqrt[8 n + 1])/2];
    DistSS[n_, b_] := DivisorSum[n, MoebiusMu[n/#]*((b + 1)^# - b^#)&] /n;
    a[n_] := DistSS[(((trinv[n] - 1)*(((1/2)*trinv[n]) + 1)) - n) + 1, (n - ((trinv[n]*(trinv[n] - 1))/2))];
    Table[a[n], {n, 0, 119}] (* Jean-François Alcover, Mar 06 2016, adapted from Maple *)

Formula

Row n is the inverse Euler transform of j-> n^(j-1). - Alois P. Heinz, Jun 23 2018
Previous Showing 61-70 of 144 results. Next