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 101-110 of 282 results. Next

A376308 Run-compression of the sequence of first differences of prime-powers.

Original entry on oeis.org

1, 2, 1, 2, 3, 1, 2, 4, 2, 1, 5, 4, 2, 4, 2, 4, 6, 2, 3, 4, 2, 6, 2, 6, 8, 4, 2, 4, 2, 4, 8, 4, 2, 1, 3, 6, 2, 10, 2, 6, 4, 2, 4, 6, 2, 10, 2, 4, 2, 12, 4, 2, 4, 6, 2, 8, 5, 1, 6, 2, 6, 4, 2, 6, 4, 14, 4, 2, 4, 14, 6, 4, 2, 4, 6, 2, 6, 4, 6, 8, 4, 8, 10, 2, 10
Offset: 1

Views

Author

Gus Wiseman, Sep 20 2024

Keywords

Comments

We define the run-compression of a sequence to be the anti-run obtained by reducing each run of repeated parts to a single part. Alternatively, run-compression removes all parts equal to the part immediately to their left. For example, (1,1,2,2,1) has run-compression (1,2,1).

Examples

			The sequence of prime-powers (A246655) is:
  2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, ...
The sequence of first differences (A057820) of prime-powers is:
  1, 1, 1, 2, 1, 1, 2, 2, 3, 1, 2, 4, 2, 2, 2, 2, 1, 5, 4, 2, 4, 2, 4, ...
The run-compression is A376308 (this sequence).
		

Crossrefs

For primes instead of prime-powers we have A037201, halved A373947.
For squarefree numbers instead of prime-powers we have A376305.
For run-lengths instead of compression we have A376309.
For run-sums instead of compression we have A376310.
For positions of first appearances we have A376341, sorted A376340.
A000040 lists the prime numbers, differences A001223.
A000961 and A246655 list prime-powers, differences A057820.
A003242 counts compressed compositions, ranks A333489.
A024619 and A361102 list non-prime-powers, differences A375708.
A116861 counts partitions by compressed sum, by compressed length A116608.
A373948 encodes compression using compositions in standard order.

Programs

  • Mathematica
    First/@Split[Differences[Select[Range[100],PrimePowerQ]]]

A378035 Greatest perfect power < prime(n).

Original entry on oeis.org

1, 1, 4, 4, 9, 9, 16, 16, 16, 27, 27, 36, 36, 36, 36, 49, 49, 49, 64, 64, 64, 64, 81, 81, 81, 100, 100, 100, 100, 100, 125, 128, 128, 128, 144, 144, 144, 144, 144, 169, 169, 169, 169, 169, 196, 196, 196, 216, 225, 225, 225, 225, 225, 243, 256, 256, 256, 256
Offset: 1

Views

Author

Gus Wiseman, Nov 23 2024

Keywords

Comments

Perfect powers (A001597) are 1 and numbers with a proper integer root, complement A007916.

Examples

			The first number line below shows the perfect powers.
The second shows each positive integer k at position prime(k).
-1-----4-------8-9------------16----------------25--27--------32------36----
===1=2===3===4=======5===6=======7===8=======9==========10==11==========12==
		

Crossrefs

Restriction of A081676 to the primes.
Positions of last appearances are also A377283.
A version for squarefree numbers is A378032.
The opposite is A378249 (run lengths A378251), restriction of A377468 to the primes.
The union is A378253.
Terms appearing exactly once are A378355.
Run lengths are A378356, first differences of A377283, complement A377436.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289.
A007916 lists the nonperfect powers, differences A375706.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A080769 counts primes between perfect powers, prime powers A067871.
A131605 lists perfect powers that are not prime powers.
A377432 counts perfect powers between primes, zeros A377436, postpositives A377466.

Programs

  • Mathematica
    radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
    Table[NestWhile[#-1&,Prime[n],radQ[#]&],{n,100}]
  • PARI
    a(n) = my(k=prime(n)-1); while (!(ispower(k) || (k==1)), k--); k; \\ Michel Marcus, Nov 25 2024
    
  • Python
    from sympy import mobius, integer_nthroot, prime
    def A378035(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(x-1+sum(mobius(k)*(integer_nthroot(x,k)[0]-1) for k in range(2,x.bit_length())))
        m = (p:=prime(n)-1)-f(p)
        return bisection(lambda x:f(x)+m,m,m) # Chai Wah Wu, Nov 25 2024

A378249 Least perfect power > prime(n).

Original entry on oeis.org

4, 4, 8, 8, 16, 16, 25, 25, 25, 32, 32, 49, 49, 49, 49, 64, 64, 64, 81, 81, 81, 81, 100, 100, 100, 121, 121, 121, 121, 121, 128, 144, 144, 144, 169, 169, 169, 169, 169, 196, 196, 196, 196, 196, 216, 216, 216, 225, 243, 243, 243, 243, 243, 256, 289, 289, 289
Offset: 1

Views

Author

Gus Wiseman, Nov 21 2024

Keywords

Comments

Perfect-powers (A001597) are numbers with a proper integer root, complement A007916.
Which terms appear only once? Just 128, 225, 256, 64009, 1295044?

Examples

			The first number line below shows the perfect powers. The second shows each prime.
-1-----4-------8-9------------16----------------25--27--------32------36------------------------49--
===2=3===5===7======11==13======17==19======23==========29==31==========37======41==43======47======
		

Crossrefs

A version for prime powers (but starting with prime(k) + 1) is A345531.
Positions of last appearances are A377283, complement A377436.
Restriction of A377468 to the primes, for prime powers A000015.
The opposite is A378035, restriction of A081676.
The union is A378250.
Run lengths are A378251.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289, seconds A376559.
A007916 lists numbers that are not perfect powers, differences A375706, seconds A376562.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A131605 lists perfect powers that are not prime powers.
A377432 counts perfect powers between primes, zeros A377436, postpositives A377466.

Programs

  • Mathematica
    radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
    Table[NestWhile[#+1&,Prime[n],radQ[#]&],{n,100}]
  • PARI
    f(p) = p++; while(!ispower(p), p++); p;
    lista(nn) = apply(f, primes(nn)); \\ Michel Marcus, Dec 19 2024

A320813 Number of non-isomorphic multiset partitions of an aperiodic multiset of weight n such that there are no singletons and all parts are themselves aperiodic multisets.

Original entry on oeis.org

1, 0, 1, 2, 5, 13, 33, 104, 293, 938, 2892
Offset: 0

Views

Author

Gus Wiseman, Nov 08 2018

Keywords

Comments

Also the number of nonnegative integer matrices up to row and column permutations with sum of elements equal to n and no zero rows or columns, in which (1) the row sums are all > 1, (2) the positive entries in each row are relatively prime, and (3) the column-sums are relatively prime.
A multiset is aperiodic if its multiplicities are relatively prime.
The weight of a multiset partition is the sum of sizes of its parts. Weight is generally not the same as number of vertices.

Examples

			Non-isomorphic representatives of the a(2) = 1 through a(5) = 13 multiset partitions:
  {{1,2}}  {{1,2,2}}  {{1,2,2,2}}    {{1,1,2,2,2}}
           {{1,2,3}}  {{1,2,3,3}}    {{1,2,2,2,2}}
                      {{1,2,3,4}}    {{1,2,2,3,3}}
                      {{1,2},{3,4}}  {{1,2,3,3,3}}
                      {{1,3},{2,3}}  {{1,2,3,4,4}}
                                     {{1,2,3,4,5}}
                                     {{1,2},{1,2,2}}
                                     {{1,2},{2,3,3}}
                                     {{1,2},{3,4,4}}
                                     {{1,2},{3,4,5}}
                                     {{1,3},{2,3,3}}
                                     {{1,4},{2,3,4}}
                                     {{2,3},{1,2,3}}
		

Crossrefs

This is the case of A320804 where the underlying multiset is aperiodic.

Programs

  • Mathematica
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mpm[n_]:=Join@@Table[Union[Sort[Sort /@ (#/.x_Integer:>s[[x]])]&/@sps[Range[n]]],{s,Flatten[MapIndexed[Table[#2,{#1}]&,#]]& /@ IntegerPartitions[n]}];
    brute[m_]:=First[Sort[Table[Sort[Sort /@ (m/.Rule@@@Table[{i,p[[i]]},{i,Length[p]}])], {p,Permutations[Union@@m]}]]];
    aperQ[m_]:=Length[m]==0||GCD@@Length/@Split[Sort[m]]==1;
    Table[Length[Union[brute /@ Select[mpm[n],And[Min@@Length/@#>1,aperQ[Join@@#]&&And@@aperQ /@ #]&]]],{n,0,7}] (* Gus Wiseman, Jan 19 2024 *)

Extensions

Definition corrected by Gus Wiseman, Jan 19 2024

A367584 Least number whose multiset multiplicity kernel (in which each prime exponent becomes the least prime factor with that exponent) is n. First position of n in A367580.

Original entry on oeis.org

1, 2, 3, 6, 5, 12, 7, 30, 15, 20, 11, 90, 13, 28, 45, 210, 17, 60, 19, 150, 63, 44, 23, 630, 35, 52, 105, 252, 29, 360, 31, 2310, 99, 68, 175, 2100, 37, 76, 117, 1050, 41, 504, 43, 396, 525, 92, 47, 6930, 77, 140, 153, 468, 53, 420, 275, 1470, 171, 116, 59
Offset: 1

Views

Author

Gus Wiseman, Nov 29 2023

Keywords

Comments

We define the multiset multiplicity kernel (MMK) of a positive integer n to be the product of (least prime factor with exponent k)^(number of prime factors with exponent k) over all distinct exponents k appearing in the prime factorization of n. For example, 90 has prime factorization 2^1 * 3^2 * 5^1, so for k = 1 we have 2^2, and for k = 2 we have 3^1, so MMK(90) = 12. As an operation on multisets, MMK is represented by the triangle A367579, and as an operation on their ranks it is represented by A367580.

Examples

			The least number with multiset multiplicity kernel 9 is 15, so a(9) = 15.
The terms together with their prime indices begin:
   1 ->  1: {}
   2 ->  2: {1}
   3 ->  3: {2}
   4 ->  6: {1,2}
   5 ->  5: {3}
   6 -> 12: {1,1,2}
   7 ->  7: {4}
   8 -> 30: {1,2,3}
   9 -> 15: {2,3}
  10 -> 20: {1,1,3}
  11 -> 11: {5}
  12 -> 90: {1,2,2,3}
  13 -> 13: {6}
  14 -> 28: {1,1,4}
  15 -> 45: {2,2,3}
  16 ->210: {1,2,3,4}
		

Crossrefs

Positions of primes are A000040.
Positions of squarefree numbers are A000961.
All terms are rootless A007916.
Contains no nonprime prime powers A246547.
The MMK triangle is A367579, sum A367581, min A055396, max A367583.
Positions of first appearances in A367580.
The sorted version is A367585.
The complement is A367768.
A007947 gives squarefree kernel.
A027746 lists prime factors, length A001222, indices A112798.
A027748 lists distinct prime factors, length A001221, indices A304038.
A071625 counts distinct prime exponents.
A124010 gives prime signature, sorted A118914.

Programs

  • Mathematica
    nn=1000;
    mmk[q_]:=With[{mts=Length/@Split[q]}, Sort[Table[Min@@Select[q,Count[q,#]==i&], {i,mts}]]];
    spnm[y_]:=Max@@NestWhile[Most, Sort[y], Union[#]!=Range[Max@@#]&];
    qq=Table[Times@@mmk[Join@@ConstantArray@@@FactorInteger[n]], {n,nn}];
    Table[Position[qq,i][[1,1]], {i,spnm[qq]}]

Formula

a(p) = p for all primes p.

A375705 Sum of the n-th maximal run of adjacent (increasing by one at a time) non-perfect-powers.

Original entry on oeis.org

5, 18, 75, 164, 26, 118, 102, 510, 791, 1160, 1629, 2210, 369, 253, 2040, 3756, 4745, 3914, 1764, 3978, 2994, 8720, 10421, 6003, 5984, 14459, 16820, 19425, 13446, 8328, 25415, 28824, 32525, 36530, 40851, 45500, 50489, 55830, 37259, 23276, 67616, 74085, 80954
Offset: 1

Views

Author

Gus Wiseman, Aug 29 2024

Keywords

Comments

Non-perfect-powers (A007916) are numbers without a proper integer root.

Examples

			The list of all non-perfect-powers, split into runs, begins:
   2   3
   5   6   7
  10  11  12  13  14  15
  17  18  19  20  21  22  23  24
  26
  28  29  30  31
  33  34  35
  37  38  39  40  41  42  43  44  45  46  47  48
Row n begins with A375703(n), ends with A375704(n), adds up to a(n), and has length A375702(n).
		

Crossrefs

For nonprime numbers we have A054265, anti-runs A373404.
For nonsquarefree numbers we have A373414, anti-runs A373412.
For squarefree numbers we have A373413, anti-runs A373411.
For prime-powers we have A373675, anti-runs A373576.
For non-prime-powers we have A373678, anti-runs A373679.
The anti-run version is A375737, sums of A375736.
A001597 lists perfect-powers, differences A053289.
A007916 lists non-perfect-powers, differences A375706.
A046933 counts composite numbers between primes.
For runs of non-perfect-powers:
- length: A375702 = A053289(n+1) - 1
- first: A375703 (same as A216765 with 2 exceptions)
- last: A375704 (same as A045542 with 8 removed)
- sum: A375705 (this)

Programs

  • Mathematica
    radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
    Total/@Split[Select[Range[100],radQ],#1+1==#2&]//Most

A376598 Points of nonzero curvature in the sequence of prime-powers inclusive (A000961).

Original entry on oeis.org

4, 5, 7, 9, 10, 11, 12, 13, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 62, 63, 64, 65, 66, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80
Offset: 1

Views

Author

Gus Wiseman, Oct 05 2024

Keywords

Comments

These are points at which the second differences (A376596) are nonzero.
Inclusive means 1 is a prime-power. For the exclusive version, subtract 1 from all terms.

Examples

			The prime-powers inclusive (A000961) are:
  1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, ...
with first differences (A057820):
  1, 1, 1, 1, 2, 1, 1, 2, 2, 3, 1, 2, 4, 2, 2, 2, 2, 1, 5, 4, 2, 4, 2, 4, 6, 2, ...
with first differences (A376596):
  0, 0, 0, 1, -1, 0, 1, 0, 1, -2, 1, 2, -2, 0, 0, 0, -1, 4, -1, -2, 2, -2, 2, 2, ...
with nonzeros at (A376598):
  4, 5, 7, 9, 10, 11, 12, 13, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, ...
		

Crossrefs

The first differences were A057820, see also A376340.
First differences are A376309.
These are the nonzeros of A376596 (sorted firsts A376653, exclusive A376654).
The complement is A376597.
A000961 lists prime-powers inclusive, exclusive A246655.
A001597 lists perfect-powers, complement A007916.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.
`A064113 lists positions of adjacent equal prime gaps.
For prime-powers inclusive: A057820 (first differences), A376597 (second differences), A376597 (inflections and undulations), A376653 (sorted firsts in second differences).
For points of nonzero curvature: A333214 (prime), A376603 (composite), A376589 (non-perfect-power), A376592 (squarefree), A376595 (nonsquarefree), A376601 (non-prime-power).

Programs

  • Mathematica
    Join@@Position[Sign[Differences[Select[Range[1000], #==1||PrimePowerQ[#]&],2]],1|-1]

A376652 Points of downward concavity in the sequence of composite numbers (A002808).

Original entry on oeis.org

2, 6, 10, 13, 19, 24, 28, 31, 36, 42, 47, 51, 56, 59, 64, 71, 75, 79, 82, 95, 98, 104, 114, 119, 124, 127, 132, 138, 148, 152, 163, 174, 178, 181, 187, 196, 201, 206, 212, 217, 221, 230, 243, 247, 250, 263, 268, 278, 281, 286, 293, 298, 303, 306, 311, 318, 321
Offset: 1

Views

Author

Gus Wiseman, Oct 06 2024

Keywords

Comments

These are points at which the second differences (A073445) are negative.
Also positions of strict descents in the first differences (A073783) of composite numbers (A002808).

Examples

			The composite numbers are (A002808):
  4, 6, 8, 9, 10, 12, 14, 15, 16, 18, 20, 21, 22, 24, 25, 26, 27, 28, 30, 32, 33, ...
with first differences (A073783):
  2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, ...
with second differences (A073445):
  0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 1, 0, -1, 0, 0, 0, 1, -1, 0, ...
with negative terms at (A376651):
  2, 6, 10, 13, 19, 24, 28, 31, 36, 42, 47, 51, 56, 59, 64, 71, 75, 79, 82, 95, 98, ...
		

Crossrefs

The version for A000002 is A156242, positive A022297.
Partitions into composite numbers are counted by A023895, factorizations A050370.
For first differences we had A065310 or A073783, ones A375929.
These are the positions of negative terms in A073445, positive A376651.
For prime instead of composite we have A258026, positive A258025.
For zero second differences instead of negative we have A376602.
For composite numbers: A002808 (terms), A073783 (first differences), A073445 (second differences), A376602 (inflections and undulations), A376603 (nonzero curvature), A376651 (concave-up).

Programs

  • Maple
    Comps:= remove(isprime, [seq(i,i=4..1000)]):
    D1:= Comps[2..-1]-Comps[1..-2]:
    D2:= D1[2..-1]-D1[1..-2]:
    select(t -> D2[t] < 0, [$1..nops(D2)]); # Robert Israel, Nov 06 2024
  • Mathematica
    Join@@Position[Sign[Differences[Select[Range[1000],CompositeQ],2]],-1]

A378251 Number of primes between consecutive perfect powers, zeros omitted.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 23 2024

Keywords

Comments

First differences of A377283 and A378365. Run-lengths of A378035 and A378249.
Perfect powers (A001597) are 1 and numbers with a proper integer root, complement A007916.

Examples

			The first number line below shows the perfect powers. The second shows each prime. To get a(n) we count the primes between consecutive perfect powers, skipping the cases where there are none.
-1-----4-------8-9------------16----------------25--27--------32------36----
===2=3===5===7======11==13======17==19======23==========29==31==========37==
		

Crossrefs

Same as A080769 with 0's removed (which were at positions A274605).
First differences of A377283 and A378365 (union of A378356).
Run-lengths of A378035 (union A378253) and A378249 (union A378250).
The version for nonprime prime powers is A378373, with zeros A067871.
A000040 lists the primes, differences A001223.
A000961 lists the powers of primes, differences A057820.
A001597 lists the perfect powers, differences A053289, run-lengths of A377468.
A007916 lists the non-perfect powers, differences A375706.
A069623 counts perfect powers <= n.
A076411 counts perfect powers < n.
A131605 lists perfect powers that are not prime powers.
A377432 counts perfect powers between primes, see A377434, A377436, A377466.

Programs

  • Maple
    N:= 10^6: # to use perfect powers up to N
    PP:= {1,seq(seq(i^j,j=2..ilog[i](N)),i=2..isqrt(N))}:
    PP:= sort(convert(PP,list)):
    M:= map(numtheory:-pi, PP):
    subs(0=NULL, M[2..-1]-M[1..-2]): # Robert Israel, Jan 23 2025
  • Mathematica
    radQ[n_]:=n>1&&GCD@@Last/@FactorInteger[n]==1;
    Length/@Split[Table[NestWhile[#+1&,Prime[n],radQ[#]&],{n,100}]]

A102430 Triangle read by rows where T(n,k) is the number of integer partitions of n > 1 into powers of k > 1.

Original entry on oeis.org

2, 2, 2, 4, 2, 2, 4, 2, 2, 2, 6, 3, 2, 2, 2, 6, 3, 2, 2, 2, 2, 10, 3, 3, 2, 2, 2, 2, 10, 5, 3, 2, 2, 2, 2, 2, 14, 5, 3, 3, 2, 2, 2, 2, 2, 14, 5, 3, 3, 2, 2, 2, 2, 2, 2, 20, 7, 4, 3, 3, 2, 2, 2, 2, 2, 2, 20, 7, 4, 3, 3, 2, 2, 2, 2, 2, 2, 2, 26, 7, 4, 3, 3, 3, 2, 2, 2, 2, 2, 2, 2
Offset: 2

Views

Author

Marc LeBrun, Jan 08 2005

Keywords

Comments

All entries above main diagonal are = 1.

Examples

			The T(9,3)=5 partitions of 9 into powers of 3: 111111111, 1111113, 11133, 333, 9.
From _Gus Wiseman_, Jun 07 2019: (Start)
Triangle begins:
   2
   2  2
   4  2  2
   4  2  2  2
   6  3  2  2  2
   6  3  2  2  2  2
  10  3  3  2  2  2  2
  10  5  3  2  2  2  2  2
  14  5  3  3  2  2  2  2  2
  14  5  3  3  2  2  2  2  2  2
  20  7  4  3  3  2  2  2  2  2  2
  20  7  4  3  3  2  2  2  2  2  2  2
  26  7  4  3  3  3  2  2  2  2  2  2  2
  26  9  4  4  3  3  2  2  2  2  2  2  2  2
  36  9  6  4  3  3  3  2  2  2  2  2  2  2  2
  36  9  6  4  3  3  3  2  2  2  2  2  2  2  2  2
  46 12  6  4  4  3  3  3  2  2  2  2  2  2  2  2  2
Row n = 8 counts the following partitions:
  8          3311       44         5111       611        71         8
  44         311111     41111      11111111   11111111   11111111   11111111
  422        11111111   11111111
  2222
  4211
  22211
  41111
  221111
  2111111
  11111111
(End)
		

Crossrefs

Same as A308558 except for the k = 1 column.
Row sums are A102431.
First column (k = 2) is A018819.
Second column (k = 3) is A062051.

Programs

  • Maple
    b:= proc(n, i, k) option remember; `if`(n=0, 1, `if`(i<0, 0,
          b(n, i-1, k)+(p-> `if`(p>n, 0, b(n-p, i, k)))(k^i)))
        end:
    T:= (n, k)-> b(n, ilog[k](n), k):
    seq(seq(T(n, k), k=2..n), n=2..20);  # Alois P. Heinz, Oct 12 2019
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And@@(IntegerQ[Log[k,#]]&/@#)&]],{n,2,10},{k,2,n}] (* Gus Wiseman, Jun 07 2019 *)

Formula

T(1, k) = 1, T(n, 1) = choose(2n-1, n), T(n>1, k>1) = T(n-1, k) + (T(n/k, k) if k divides n, else 0)

Extensions

Corrected and rewritten by Gus Wiseman, Jun 07 2019
Previous Showing 101-110 of 282 results. Next