A000004 -id:A000004 - cp's OEIS frontend

A195040 Square array read by antidiagonals with T(n,k) = kn^2/4+(k-4)((-1)^n-1)/8, n>=0, k>=0.

0, 1, 0, 0, 1, 0, 1, 1, 1, 0, 0, 3, 2, 1, 0, 1, 4, 5, 3, 1, 0, 0, 7, 8, 7, 4, 1, 0, 1, 9, 13, 12, 9, 5, 1, 0, 0, 13, 18, 19, 16, 11, 6, 1, 0, 1, 16, 25, 27, 25, 20, 13, 7, 1, 0, 0, 21, 32, 37, 36, 31, 24, 15, 8, 1, 0, 1, 25, 41, 48, 49, 45, 37, 28, 17, 9, 1, 0

Offset: 0

Omar E. Pol, Sep 27 2011

Also, if k >= 2 and m = 2*k, then column k lists the numbers of the form k*n^2 and the centered m-gonal numbers interleaved.
For k >= 3, this is also a table of concentric polygonal numbers. Column k lists the concentric k-gonal numbers.
It appears that the first differences of column k are the numbers that are congruent to {1, k-1} mod k, if k >= 3.

			Array begins:
  0,   0,   0,   0,   0,   0,   0,   0,   0,   0, ...
  1,   1,   1,   1,   1,   1,   1,   1,   1,   1, ...
  0,   1,   2,   3,   4,   5,   6,   7,   8,   9, ...
  1,   3,   5,   7,   9,  11,  13,  15,  17,  19, ...
  0,   4,   8,  12,  16,  20,  24,  28,  32,  36, ...
  1,   7,  13,  19,  25,  31,  37,  43,  49,  55, ...
  0,   9,  18,  27,  36,  45,  54,  63,  72,  81, ...
  1,  13,  25,  37,  49,  61,  73,  85,  97, 109, ...
  0,  16,  32,  48,  64,  80,  96, 112, 128, 144, ...
  1,  21,  41,  61,  81, 101, 121, 141, 161, 181, ...
  0,  25,  50,  75, 100, 125, 150, 175, 200, 225, ...
  ...

Muniru A Asiru, Rows n=0..100, flattened

Rows n: A000004 (n=0), A000012 (n=1), A001477 (n=2), A005408 (n=3), A008586 (n=4), A016921 (n=5), A008591 (n=6), A017533 (n=7), A008598 (n=8), A215145 (n=9), A008607 (n=10).

Columns k: A000035 (k=0), A004652 (k=1), A000982 (k=2), A077043 (k=3), A000290 (k=4), A032527 (k=5), A032528 (k=6), A195041 (k=7), A077221 (k=8), A195042 (k=9), A195142 (k=10), A195043 (k=11), A195143 (k=12), A195045 (k=13), A195145 (k=14), A195046 (k=15), A195146 (k=16), A195047 (k=17), A195147 (k=18), A195048 (k=19), A195148 (k=20), A195049 (k=21), A195149 (k=22), A195058 (k=23), A195158 (k=24).

nmax:=13;; T:=List([0..nmax],n->List([0..nmax],k->k*n^2/4+(k-4)*((-1)^n-1)/8));; b:=List([2..nmax],n->OrderedPartitions(n,2));;
a:=Flat(List([1..Length(b)],i->List([1..Length(b[i])],j->T[b[i][j][2]][b[i][j][1]]))); # Muniru A Asiru, Jul 19 2018

A195040 := proc(n,k)
        k*n^2/4+((-1)^n-1)*(k-4)/8 ;
end proc:
for d from 0 to 12 do
        for k from 0 to d do
                printf("%d,",A195040(d-k,k)) ;
        end do:
end do; # R. J. Mathar, Sep 28 2011

t[n_, k_] := k*n^2/4+(k-4)*((-1)^n-1)/8; Flatten[ Table[ t[n-k, k], {n, 0, 11}, {k, 0, n}]] (* Jean-François Alcover, Dec 14 2011 *)

A195151 Square array read by antidiagonals upwards: T(n,k) = n((k-2)(-1)^n+k+2)/4, n >= 0, k >= 0.

Original entry on oeis.org

0, 1, 0, 0, 1, 0, 3, 1, 1, 0, 0, 3, 2, 1, 0, 5, 2, 3, 3, 1, 0, 0, 5, 4, 3, 4, 1, 0, 7, 3, 5, 6, 3, 5, 1, 0, 0, 7, 6, 5, 8, 3, 6, 1, 0, 9, 4, 7, 9, 5, 10, 3, 7, 1, 0, 0, 9, 8, 7, 12, 5, 12, 3, 8, 1, 0, 11, 5, 9, 12, 7, 15, 5, 14, 3, 9, 1, 0, 0, 11, 10, 9, 16, 7

Offset: 0

Omar E. Pol, Sep 14 2011

Also square array T(n,k) read by antidiagonals in which column k lists the multiples of k and the odd numbers interleaved, n>=0, k>=0. Also square array T(n,k) read by antidiagonals in which if n is even then row n lists the multiples of (n/2), otherwise if n is odd then row n lists a constant sequence: the all n's sequence. Partial sums of the numbers of column k give the column k of A195152. Note that if k >= 1 then partial sums of the numbers of the column k give the generalized m-gonal numbers, where m = k + 4.

All columns are multiplicative. - Andrew Howroyd, Jul 23 2018

			Array begins:
.  0,   0,   0,   0,   0,   0,   0,   0,   0,   0,...
.  1,   1,   1,   1,   1,   1,   1,   1,   1,   1,...
.  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,...
.  3,   3,   3,   3,   3,   3,   3,   3,   3,   3,...
.  0,   2,   4,   6,   8,  10,  12,  14,  16,  18,...
.  5,   5,   5,   5,   5,   5,   5,   5,   5,   5,...
.  0,   3,   6,   9,  12,  15,  18,  21,  24,  27,...
.  7,   7,   7,   7,   7,   7,   7,   7,   7,   7,...
.  0,   4,   8,  12,  16,  20,  24,  28,  32,  36,...
.  9,   9,   9,   9,   9,   9,   9,   9,   9,   9,...
.  0,   5,  10,  15,  20,  25,  30,  35,  40,  45,...
...

Rows: A000004, A000012, A001477, A010701, A005843, A010716, A008585, A010727, A008586, A010734, A008587.

Columns k: A026741 (k=1), A001477 (k=2), zero together with A080512 (k=3), A022998 (k=4), A195140 (k=5), zero together with A165998 (k=6), A195159 (k=7), A195161 (k=8), A195312 k=(9), A195817 (k=10), A317311 (k=11), A317312 (k=12), A317313 (k=13), A317314 k=(14), A317315 (k=15), A317316 (k=16), A317317 (k=17), A317318 (k=18), A317319 k=(19), A317320 (k=20), A317321 (k=21), A317322 (k=22), A317323 (k=23), A317324 k=(24), A317325 (k=25), A317326 (k=26).

T(n,k) = n*((k-2)*(-1)^n+k+2)/4 \\ Andrew Howroyd, Jul 23 2018

A212957 A(n,k) is the number of moduli m such that the multiplicative order of k mod m equals n; square array A(n,k), n>=1, k>=1, read by antidiagonals.

Original entry on oeis.org

0, 1, 0, 2, 1, 0, 2, 2, 1, 0, 3, 2, 2, 2, 0, 2, 5, 4, 6, 1, 0, 4, 2, 3, 4, 4, 3, 0, 2, 6, 2, 12, 6, 10, 1, 0, 4, 4, 8, 4, 9, 16, 2, 4, 0, 3, 6, 2, 26, 4, 37, 6, 14, 2, 0, 4, 3, 12, 18, 4, 10, 3, 8, 4, 5, 0, 2, 12, 5, 14, 6, 42, 2, 28, 26, 16, 3, 0

Offset: 1

Alois P. Heinz, Jun 01 2012

			A(4,3) = 6: 3^4 = 81 == 1 (mod m) for m in {5,10,16,20,40,80}.
Square array A(n,k) begins:
  0,  1,  2,  2,  3,  2,  4,  2, ...
  0,  1,  2,  2,  5,  2,  6,  4, ...
  0,  1,  2,  4,  3,  2,  8,  2, ...
  0,  2,  6,  4, 12,  4, 26, 18, ...
  0,  1,  4,  6,  9,  4,  4,  6, ...
  0,  3, 10, 16, 37, 10, 42, 24, ...
  0,  1,  2,  6,  3,  2, 12, 10, ...
  0,  4, 14,  8, 28,  8, 48, 72, ...

Alois P. Heinz, Antidiagonals n = 1..60
Wikipedia, Multiplicative order.

Columns k=1-10 give: A000004, A059499, A059885, A059886, A059887, A059888, A059889, A059890, A059891, A059892.
Rows n=1-10 give: A000005, A059907, A059908, A059909, A059910, A059911, A218256, A218257, A218258, A218259.
Main diagonal gives A252760.
Cf. A000005, A008683.

with(numtheory):
A:= (n, k)-> add(mobius(n/d)*tau(k^d-1), d=divisors(n)):
seq(seq(A(n, 1+d-n), n=1..d), d=1..15);

a[n_, k_] := Sum[ MoebiusMu[n/d] * DivisorSigma[0, k^d - 1], {d, Divisors[n]}]; a[1, 1] = 0; Table[ a[n - k + 1, k], {n, 1, 12}, {k, n, 1, -1}] // Flatten (* Jean-François Alcover, Dec 12 2012 *)

a(n, k) = if(k == 1, 0, sumdiv(n, d, moebius(n/d) * numdiv(k^d-1))); \\ Amiram Eldar, Jan 25 2025

A(n,k) = |{m : multiplicative order of k mod m = n}|.

A(n,k) = Sum_{d|n} mu(n/d)*tau(k^d-1), mu = A008683, tau = A000005.

A031214 Initial term of sequence An.

Original entry on oeis.org

0, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 0, 1, 1, 1, 1, 1, 2, 1, 0, 1, 2, 0, 1, 0, 2, 2, 2, 1, 2, 1, 1, 2, 1, 0, 1, 1, 1, 0, 1, 2, 8, 14, 4, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 0, 4, 1, 1, 1, 1, 0, 1, 0, 1, 0, 0, 1, 0, 1, 1, 0, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1

Offset: 1

Jeff Burch

This ignores the offset and gives the first term of the actual entry.
Since the sequences in the OEIS occasionally change their initial terms (for editorial reasons), this is an especially ill-defined sequence! - N. J. A. Sloane, Jan 01 2005
Sequences like this are deprecated. - Joerg Arndt, Apr 16 2020

			A000001 begins 0,1,1,1,2,1,2,1,5,2,... so a(1) = 0 = a(31214).

PrimeFan, Esoteric Integer Sequences
PrimeFan, Esoteric Integer Sequences [Cached copy]
N. J. A. Sloane, Online Encyclopedia of Integer Sequences
Index entries for sequences whose definition involves A_n (or An).

Cf. A039928, A051070, A091967, A107357, A102288; also A000001, A000002, A000003, A000004, A000005, etc.

Data updated by Sean A. Irvine, Apr 16 2020

A212162 Triangle T(n,k), n>=1, 0<=k<=n^2, read by rows: row n gives the coefficients of the chromatic polynomial of the rhombic hexagonal square grid graph RH_(n,n), highest powers first.

Original entry on oeis.org

1, 0, 1, -5, 8, -4, 0, 1, -16, 112, -448, 1120, -1791, 1786, -1012, 248, 0, 1, -33, 510, -4898, 32703, -160859, 602408, -1749715, 3975561, -7068408, 9755858, -10265148, 7968348, -4304712, 1445104, -226720, 0, 1, -56, 1508, -25992, 321994, -3051871, 23000726, -141421592, 722137763, -3101089711

Offset: 1

Alois P. Heinz, May 02 2012

T differs from A212194 first at (n,k) = (5,10): T(5,10) = -3101089711, A212194(5,10) = -3101089710.

The rhombic hexagonal square grid graph RH_(n,n) has n^2 = A000290(n) vertices and (n-1)*(3*n-1) = A045944(n-1) edges. The chromatic polynomial of RH_(n,n) has n^2+1 = A002522(n) coefficients.

			3 example graphs:                        o--o--o
.                                        | /| /|
.                                        |/ |/ |
.                            o--o        o--o--o
.                            | /|        | /| /|
.                            |/ |        |/ |/ |
.               o            o--o        o--o--o
Graph:       RH_(1,1)      RH_(2,2)      RH_(3,3)
Vertices:       1             4             9
Edges:          0             5            16
The rhombic hexagonal square grid graph RH_(2,2) has chromatic polynomial q*(q-1)*(q-2)^2 = q^4 -5*q^3 +8*q^2 -4*q => row 2 = [1, -5, 8, -4, 0].
Triangle T(n,k) begins:
1,    0;
1,   -5,     8,      -4,        0;
1,  -16,   112,    -448,     1120,      -1791, ...
1,  -33,   510,   -4898,    32703,    -160859, ...
1,  -56,  1508,  -25992,   321994,   -3051871, ... , -3101089711, ...
1,  -85,  3520,  -94620,  1855860,  -28306676, ...
1, -120,  7068, -272344,  7720110, -171656543, ...
1, -161, 12782, -667058, 25738055, -783003395, ...

Alois P. Heinz, Rows n = 1..8, flattened
Wikipedia, Chromatic Polynomial

Columns 1-2 give: A000012, (-1)*A045944(n-1).
Row sums (for n>1) and last elements of rows give: A000004, row lengths give: A002522.
Cf. A000290, A212163, A212194.

A228275 A(n,k) = Sum_{i=1..k} n^i; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

0, 0, 0, 0, 1, 0, 0, 2, 2, 0, 0, 3, 6, 3, 0, 0, 4, 14, 12, 4, 0, 0, 5, 30, 39, 20, 5, 0, 0, 6, 62, 120, 84, 30, 6, 0, 0, 7, 126, 363, 340, 155, 42, 7, 0, 0, 8, 254, 1092, 1364, 780, 258, 56, 8, 0, 0, 9, 510, 3279, 5460, 3905, 1554, 399, 72, 9, 0

Offset: 0

Alois P. Heinz, Aug 19 2013

A(n,k) is the total sum of lengths of longest ending contiguous subsequences with the same value over all s in {1,...,n}^k:
A(4,1) = 4 = 1+1+1+1: [1], [2], [3], [4].
A(1,4) = 4: [1,1,1,1].
A(3,2) = 12 = 2+1+1+1+2+1+1+1+2: [1,1], [1,2], [1,3], [2,1], [2,2], [2,3], [3,1], [3,2], [3,3].
A(2,3) = 14 = 3+1+1+2+2+1+1+3: [1,1,1], [1,1,2], [1,2,1], [1,2,2], [2,1,1], [2,1,2], [2,2,1], [2,2,2].

			Square array A(n,k) begins:
  0, 0,  0,   0,    0,     0,      0,      0, ...
  0, 1,  2,   3,    4,     5,      6,      7, ...
  0, 2,  6,  14,   30,    62,    126,    254, ...
  0, 3, 12,  39,  120,   363,   1092,   3279, ...
  0, 4, 20,  84,  340,  1364,   5460,  21844, ...
  0, 5, 30, 155,  780,  3905,  19530,  97655, ...
  0, 6, 42, 258, 1554,  9330,  55986, 335922, ...
  0, 7, 56, 399, 2800, 19607, 137256, 960799, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0-10 give: A000004, A001477, A002378, A027444, A027445, A152031, A228290, A228291, A228292, A228293, A228294.
Rows n=0-11 give: A000004, A001477, A000918(k+1), A029858(k+1), A080674, A104891, A105281, A104896, A052379(k-1), A052386, A105279, A105280.
Main diagonal gives A031972.
Lower diagonal gives A226238.
Cf. A228250.

A:= (n, k)-> `if`(n=1, k, (n/(n-1))*(n^k-1)):
seq(seq(A(n, d-n), n=0..d), d=0..12);

a[0, 0] = 0; a[1, k_] := k; a[n_, k_] := n*(n^k-1)/(n-1); Table[a[n-k, k], {n, 0, 12}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 16 2013 *)

A(1,k) = k, else A(n,k) = n/(n-1)*(n^k-1).
A(n,k) = Sum_{i=1..k} n^i.
A(n,k) = Sum_{i=1..k+1} binomial(k+1,i)*A(n-i,k)*(-1)^(i+1) for n>k, given values A(0,k), A(1,k),..., A(k,k). - Yosu Yurramendi, Sep 03 2013

A269129 Number A(n,k) of sequences with k copies each of 1,2,...,n avoiding the pattern 12...n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 1, 5, 1, 0, 0, 1, 43, 23, 1, 0, 0, 1, 374, 1879, 119, 1, 0, 0, 1, 3199, 173891, 102011, 719, 1, 0, 0, 1, 26945, 16140983, 117392909, 7235651, 5039, 1, 0, 0, 1, 224296, 1474050783, 142951955371, 117108036719, 674641325, 40319, 1

Offset: 0

Alois P. Heinz, Feb 19 2016

			Square array A(n,k) begins:
  0,   0,      0,         0,            0,               0, ...
  1,   0,      0,         0,            0,               0, ...
  1,   1,      1,         1,            1,               1, ...
  1,   5,     43,       374,         3199,           26945, ...
  1,  23,   1879,    173891,     16140983,      1474050783, ...
  1, 119, 102011, 117392909, 142951955371, 173996758190594, ...

Alois P. Heinz, Antidiagonals n = 0..50, flattened
J. D. Horton and A. Kurn, Counting sequences with complete increasing subsequences, Congressus Numerantium, 33 (1981), 75-80. MR 681905

Columns k=0-10 give: A057427, A033312, A267532, A269113, A269114, A269115, A269116, A269117, A269118, A269119, A269120.
Rows n=0-10 give: A000004, A000007, A000012, A269121, A269122, A269123, A269124, A269125, A269126, A269127, A269128.
Main diagonal gives: A268751.
Cf. A047909, A089759, A187783, A331562.

g:= proc(l) option remember; (n-> f(l[1..nops(l)-1])*
      binomial(n-1, l[-1]-1)+add(f(sort(subsop(j=l[j]
      -1, l))), j=1..nops(l)-1))(add(i, i=l))
    end:
f:= l->(n->`if`(n=0, 1, `if`(l[1]=0, 0, `if`(n=1 or l[-1]=1, 1,
    `if`(n=2, binomial(l[1]+l[2], l[1])-1, g(l))))))(nops(l)):
A:= (n, k)-> (k*n)!/k!^n - f([k$n]):
seq(seq(A(n, d-n), n=0..d), d=0..12);
# second Maple program:
b:= proc(k, p, j, l, t) option remember;
      `if`(k=0, (-1)^t/l!, `if`(p<0, 0, add(b(k-i, p-1,
       j+1, l+i*j, irem(t+i*j, 2))/(i!*p!^i), i=0..k)))
    end:
A:= (n, k)-> (n*k)!*(1/k!^n-b(n, k-1, 1, 0, irem(n, 2))*n!):
seq(seq(A(n, d-n), n=0..d), d=0..12);  # Alois P. Heinz, Mar 03 2016

b[k_, p_, j_, l_, t_] := b[k, p, j, l, t] = If[k == 0, (-1)^t/l!, If[p < 0, 0, Sum[b[k-i, p-1, j+1, l + i j, Mod[t + i j, 2]]/(i! p!^i), {i, 0, k}]] ];
A[n_, k_] := (n k)! (1/k!^n - b[n, k-1, 1, 0, Mod[n, 2]] n!); Table[ Table[ A[n, d-n], {n, 0, d}], {d, 0, 12}] // Flatten (* Jean-François Alcover, Apr 07 2016, after Alois P. Heinz *)

A(n,k) = A089759(k,n) - A047909(k,n) = A187783(n,k) - A047909(k,n).

A173787 Triangle read by rows: T(n,k) = 2^n - 2^k, 0 <= k <= n.

Original entry on oeis.org

0, 1, 0, 3, 2, 0, 7, 6, 4, 0, 15, 14, 12, 8, 0, 31, 30, 28, 24, 16, 0, 63, 62, 60, 56, 48, 32, 0, 127, 126, 124, 120, 112, 96, 64, 0, 255, 254, 252, 248, 240, 224, 192, 128, 0, 511, 510, 508, 504, 496, 480, 448, 384, 256, 0, 1023, 1022, 1020, 1016, 1008, 992, 960, 896, 768, 512, 0

Offset: 0

Reinhard Zumkeller, Feb 28 2010

			Triangle begins as:
   0;
   1,  0;
   3,  2,  0;
   7,  6,  4,  0;
  15, 14, 12,  8,  0;
  31, 30, 28, 24, 16, 0;

G. C. Greubel, Rows n = 0..50 of the triangle, flattened

[2^n -2^k: k in [0..n], n in [0..15]]; // G. C. Greubel, Jul 13 2021

Table[2^n -2^k, {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Jul 13 2021 *)

flatten([[2^n -2^k for k in (0..n)] for n in (0..15)]) # G. C. Greubel, Jul 13 2021

A000120(T(n,k)) = A025581(n,k).
Row sums give A000337.
Central terms give A020522.
T(2*n+1, n) = A006516(n+1).
T(2*n+3, n+2) = A059153(n).
T(n, k) = A140513(n,k) - A173786(n,k), 0 <= k <= n.
T(n, k) = A173786(n,k) - A059268(n+1,k+1), 0 < k <= n.
T(2*n, 2*k) = T(n,k) * A173786(n,k), 0 <= k <= n.
T(n, 0) = A000225(n).
T(n, 1) = A000918(n) for n>0.
T(n, 2) = A028399(n) for n>1.
T(n, 3) = A159741(n-3) for n>3.
T(n, 4) = A175164(n-4) for n>4.
T(n, 5) = A175165(n-5) for n>5.
T(n, 6) = A175166(n-6) for n>6.
T(n, n-4) = A110286(n-4) for n>3.
T(n, n-3) = A005009(n-3) for n>2.
T(n, n-2) = A007283(n-2) for n>1.
T(n, n-1) = A000079(n-1) for n>0.
T(n, n) = A000004(n).

A213191 Total sum A(n,k) of k-th powers of parts in all partitions of n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Original entry on oeis.org

0, 0, 1, 0, 1, 3, 0, 1, 4, 6, 0, 1, 6, 9, 12, 0, 1, 10, 17, 20, 20, 0, 1, 18, 39, 44, 35, 35, 0, 1, 34, 101, 122, 87, 66, 54, 0, 1, 66, 279, 392, 287, 180, 105, 86, 0, 1, 130, 797, 1370, 1119, 660, 311, 176, 128, 0, 1, 258, 2319, 5024, 4775, 2904, 1281, 558, 270, 192

Offset: 0

Alois P. Heinz, Feb 28 2013

In general, if k > 0 then column k is asymptotic to 2^((k-3)/2) * 3^(k/2) * k! * Zeta(k+1) / Pi^(k+1) * exp(Pi*sqrt(2*n/3)) * n^((k-1)/2). - Vaclav Kotesovec, May 27 2018

			Square array A(n,k) begins:
:   0,  0,   0,   0,    0,     0,     0, ...
:   1,  1,   1,   1,    1,     1,     1, ...
:   3,  4,   6,  10,   18,    34,    66, ...
:   6,  9,  17,  39,  101,   279,   797, ...
:  12, 20,  44, 122,  392,  1370,  5024, ...
:  20, 35,  87, 287, 1119,  4775, 21447, ...
:  35, 66, 180, 660, 2904, 14196, 73920, ...

Alois P. Heinz, Antidiagonals n = 0..140, flattened

Columns k=0-10 give: A006128, A066186, A066183, A229325, A229326, A229327, A229328, A229329, A229330, A229331, A229332.
Rows n=0-10 give: A000004, A000012, A052548, A229354, A229355, A229356, A229357, A229358, A229359, A229360, A229361.
Main diagonal gives A252761.
Cf. A213180.

b:= proc(n, p, k) option remember; `if`(n=0, [1, 0], `if`(p<1, [0, 0],
      add((l->`if`(m=0, l, l+[0, l[1]*p^k*m]))
          (b(n-p*m, p-1, k)), m=0..n/p)))
    end:
A:= (n, k)-> b(n, n, k)[2]:
seq(seq(A(n, d-n), n=0..d), d=0..10);

b[n_, p_, k_] := b[n, p, k] = If[n == 0, {1, 0}, If[p < 1, {0, 0}, Sum[Function[l, If[m == 0, l, l + {0, First[l]*p^k*m}]][b[n - p*m, p - 1, k]], { m, 0, n/p}]]] ; a[n_, k_] := b[n, n, k][[2]]; Table[Table[a[n, d - n], {n, 0, d}], {d, 0, 10}] // Flatten (* Jean-François Alcover, Dec 12 2013, translated from Maple *)
(* T = A066633 *) T[n_, n_] = 1; T[n_, k_] /; k, ] = 0; A[n_, k_] := Sum[T[n, j]*j^k, {j, 1, n}]; Table[A[n-k, k], {n, 0, 10}, {k, n, 0, -1}] // Flatten (* Jean-François Alcover, Dec 15 2016 *)

A(n,k) = Sum_{j=1..n} A066633(n,j) * j^k.

A349136 Möbius transform of Kimberling's paraphrases, A003602.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 3, 0, 3, 0, 5, 0, 6, 0, 4, 0, 8, 0, 9, 0, 6, 0, 11, 0, 10, 0, 9, 0, 14, 0, 15, 0, 10, 0, 12, 0, 18, 0, 12, 0, 20, 0, 21, 0, 12, 0, 23, 0, 21, 0, 16, 0, 26, 0, 20, 0, 18, 0, 29, 0, 30, 0, 18, 0, 24, 0, 33, 0, 22, 0, 35, 0, 36, 0, 20, 0, 30, 0, 39, 0, 27, 0, 41, 0, 32, 0, 28, 0, 44, 0, 36, 0, 30, 0, 36

Offset: 1

Antti Karttunen, Nov 13 2021

nonn

Antti Karttunen, Table of n, a(n) for n = 1..20000

Agrees with A055034 on odd arguments.
Cf. A000010, A003602, A008683, A349134.
Cf. A000004, A072451 (even and odd bisection).
Cf. also A347233, A349127, A349137.

with(numtheory): a:=proc(n) if n=1 then 1; elif n mod 2 = 0 then 0; else phi(n)/2; fi: end proc: seq(a(n), n=1..60); # Ridouane Oudra, Jul 13 2023

k[n_] := (n/2^IntegerExponent[n, 2] + 1)/2; a[n_] := DivisorSum[n, MoebiusMu[#] * k[n/#] &]; Array[a, 100] (* Amiram Eldar, Nov 13 2021 *)

A349136(n) = if(1==n,1, if(n%2, eulerphi(n)/2, 0));

A003602(n) = (1+(n>>valuation(n,2)))/2;
A349136(n) = sumdiv(n,d,moebius(d)*A003602(n/d));

from sympy import totient
def A349136(n): return totient(n)+1>>1 if n&1 else 0 # Chai Wah Wu, Nov 24 2023

a(n) = Sum_{d|n} A008683(d) * A003602(n/d).
a(1) = 1, a(n) = A000010(n)/2 for odd n > 1, a(n) = 0 for even n.
For all n >= 1, a(2*n-1) = A055034(2*n-1) = A072451(n).
a(n) = phi(n) - (1/2)*phi(2n), for n>1. - Ridouane Oudra, Jul 13 2023
Sum_{k=1..n} a(k) ~ (1/Pi^2)*n^2. - Amiram Eldar, Jul 15 2023

cp's OEIS Frontend

A195040 Square array read by antidiagonals with T(n,k) = k*n^2/4+(k-4)*((-1)^n-1)/8, n>=0, k>=0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Maple

Mathematica

A195151 Square array read by antidiagonals upwards: T(n,k) = n*((k-2)*(-1)^n+k+2)/4, n >= 0, k >= 0.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

A212957 A(n,k) is the number of moduli m such that the multiplicative order of k mod m equals n; square array A(n,k), n>=1, k>=1, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A031214 Initial term of sequence An.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Extensions

A212162 Triangle T(n,k), n>=1, 0<=k<=n^2, read by rows: row n gives the coefficients of the chromatic polynomial of the rhombic hexagonal square grid graph RH_(n,n), highest powers first.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

A228275 A(n,k) = Sum_{i=1..k} n^i; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A269129 Number A(n,k) of sequences with k copies each of 1,2,...,n avoiding the pattern 12...n; square array A(n,k), n>=0, k>=0, read by antidiagonals.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A173787 Triangle read by rows: T(n,k) = 2^n - 2^k, 0 <= k <= n.

Views

Author

Keywords

Examples

Links

A195040 Square array read by antidiagonals with T(n,k) = kn^2/4+(k-4)((-1)^n-1)/8, n>=0, k>=0.

A195151 Square array read by antidiagonals upwards: T(n,k) = n((k-2)(-1)^n+k+2)/4, n >= 0, k >= 0.