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.

A070893 Let r, s, t be three permutations of the set {1,2,3,..,n}; a(n) = value of Sum_{i=1..n} r(i)*s(i)*t(i), with r={1,2,3,..,n}; s={n,n-1,..,1} and t={n,n-2,n-4,...,1,...,n-3,n-1}.

Original entry on oeis.org

1, 6, 19, 46, 94, 172, 290, 460, 695, 1010, 1421, 1946, 2604, 3416, 4404, 5592, 7005, 8670, 10615, 12870, 15466, 18436, 21814, 25636, 29939, 34762, 40145, 46130, 52760, 60080, 68136, 76976, 86649, 97206, 108699, 121182, 134710, 149340
Offset: 1

Views

Author

Wouter Meeussen, May 22 2002

Keywords

Comments

See A070735 for the minimal values for these products. This sequence is an upper bound. The third permutation 't'= ceiling(abs(range(n-1/2,-n,-2))) is such that it associates its smallest factor with the largest factor of the product 'r'*'s'.
We observe that is the transform of A002717 by the following transform T: T(u_0,u_1,u_2,u_3,...) = (u_0,u_0+u_1, u_0+u_1+u_2, u_0+u_1+u_2+u_3+u_4,...). In other words, v_p = Sum_{k=0..p} u_k and the g.f. phi_v of v is given by phi_v = phi_u/(1-z). - Richard Choulet, Jan 28 2010

Examples

			{1,2,3,4,5,6,7}*{7,6,5,4,3,2,1}*{7,5,3,1,2,4,6} gives {49,60,45,16,30,48,42}, with sum 290, so a(7)=290.
		

Crossrefs

Cf. A070735, A082289. a(n)=A082290(2n-2).
Cf. A002717 (first differences). - Bruno Berselli, Aug 26 2011
Column k=3 of A166278. - Alois P. Heinz, Nov 02 2012

Programs

  • Magma
    [(1/96)*(2*n*(n+2)*(3*n^2+10*n+4)+3*(-1)^n-3): n in [1..40]]; // Vincenzo Librandi, Aug 26 2011
  • Mathematica
    Table[Plus@@(Range[n]*Range[n, 1, -1]*Ceiling[Abs[Range[n-1/2, -n, -2]]]), {n, 49}];
    (* or *)
    CoefficientList[Series[ -(1+2x)/(-1+x)^5/(1+x), {x, 0, 48}], x]//Flatten
  • PARI
    a(n)=sum(i=1,n,i*(n+1-i)*ceil(abs(n+3/2-2*i)))
    
  • PARI
    a(n)=polcoeff(if(n<0,x^4*(2+x)/((1+x)*(1-x)^5),x*(1+2*x)/((1+x)*(1-x)^5))+x*O(x^abs(n)),abs(n))
    

Formula

G.f.: x*(1+2*x)/((1+x)*(1-x)^5). - Michael Somos, Apr 07 2003
a(n) = 3*a(n-1) - 2*a(n-2) - 2*a(n-3) + 3*a(n-4) - a(n-5) + 3. If sequence is also defined for n <= 0 by this equation, then a(n)=0 for -3 <= n <= 0 and a(n)=A082289(-n) for n <= -4. - Michael Somos, Apr 07 2003
a(n) = (1/96)*(2*n*(n+2)*(3*n^2+10*n+4)+3*(-1)^n-3). a(n) - a(n-2) = A002411(n). - Bruno Berselli, Aug 26 2011

A261490 The total element sum of the n-fold f transform applied to the length n sequence of 1's. And f returns a sorted result after multiplying the elements in its input sequence with 1, 2, 3,... in descending size order.

Original entry on oeis.org

0, 1, 4, 19, 100, 633, 4626, 37878, 348224, 3542952, 39339852, 478962252, 6289532928, 89038853856, 1346224983936, 21729308136720, 371924399416320, 6740200653419520, 128878557725067264, 2598800542616444724, 54986036469506668800, 1217069235297874269792
Offset: 0

Views

Author

Alois P. Heinz, Aug 21 2015

Keywords

Crossrefs

Main diagonal of A166278.
Cf. A031971 (for ascending sort), A036740 (when sum is replaced by product).

Programs

  • Maple
    f:= l-> sort([seq(sort(l, `>`)[i]*i, i=1..nops(l))]):
    a:= n-> add(i, i=(f@@n)([1$n])):
    seq(a(n), n=0..35);
Showing 1-2 of 2 results.