A023894 -id:A023894 - cp's OEIS frontend

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

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

Gus Wiseman, Oct 05 2024

nonn

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.

			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, ...

Gus Wiseman, Points of nonzero curvature in the prime-powers inclusive.

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).
Cf. A025475, A053707, A069623, A093555, A174965, A251092, A333254, A336416, A361102.

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

A051613 a(n) = partitions of n into powers of distinct primes (1 not considered a power).

Original entry on oeis.org

1, 0, 1, 1, 1, 2, 0, 3, 2, 3, 2, 4, 3, 4, 4, 4, 8, 4, 8, 6, 9, 8, 10, 10, 13, 12, 13, 16, 16, 19, 17, 21, 23, 23, 25, 29, 31, 31, 31, 37, 40, 42, 44, 48, 49, 54, 55, 64, 67, 68, 70, 77, 84, 90, 92, 99, 102, 108, 115, 127, 133, 135, 138, 150, 165, 171, 183, 186, 198, 201, 220

Offset: 0

Vladeta Jovovic

			a(16) = 8 because we can write 16 = 2^4 = 3+13 = 5+11 = 3^2+7 = 2+3+11 = 2+3^2+5 = 2^3+3+5 = 2^2+5+7.

T. D. Noe, Table of n, a(n) for n = 0..1000
J. Bamberg, G. Cairns and D. Kilminster, The crystallographic restriction, permutations and Goldbach's conjecture, Amer. Math. Monthly, 110 (March 2003), 202-209.

Cf. A023894, A009490, A054685, A008475, A106245, A000961, A025473, A106244.

import Data.MemoCombinators (memo3, integral)
a051613' = p 1 2 where
   p x _ 0 = 1
   p x k m | m < qq       = 0
           | mod x q == 0 = p x (k + 1) m
           | otherwise    = p (q * x) (k + 1) (m - qq) + p x (k + 1) m
           where q = a025473 k; qq = a000961 k
-- Reinhard Zumkeller, Nov 23 2015

b:= proc(n,i) option remember; local p;
      p:= `if`(i<1, 1, ithprime(i));
      `if`(n=0, 1, `if`(i<1, 0, b(n,i-1)+
      add(b(n-p^j, i-1), j=1..ilog[p](n))))
    end:
a:= n-> b(n, numtheory[pi](n)):
seq(a(n), n=0..100);  # Alois P. Heinz, Feb 15 2013

max = 70; f[x_] := Product[ 1 + Sum[x^(Prime[n]^k), {k, 1, If[n > 4, 1, 6]}], {n, 1, PrimePi[max]}]; CoefficientList[ Series[f[x], {x, 0, max}] , x](* Jean-François Alcover, Sep 12 2012 *)

first(n)=my(x='x,pr=O(x^(n+1))+1); forprime(p=sqrtint(n)+1,n, pr*=1+x^p); forprime(p=2,sqrtint(n), pr*=1+sum(e=1,logint(n,2), x^p^e)); Vec(pr) \\ Charles R Greathouse IV, Jun 25 2017

a(n) = number of m such that A008475(m) = n.

G.f.: Product_{p prime} (1 + Sum_{k >= 1} x^(p^k)).

Better description from David W. Wilson, Apr 19 2000

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

Marc LeBrun, Jan 08 2005

All entries above main diagonal are = 1.

			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)

Alois P. Heinz, Rows n = 2..500, flattened

Cf. A102431, A102432, A102433, A102434, A001700, A018819, A062051, A008645, A008650, A008652, A008648.
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.
Cf. A000961, A001597, A007916, A023894, A052410, A112344.

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

Table[Length[Select[IntegerPartitions[n],And@@(IntegerQ[Log[k,#]]&/@#)&]],{n,2,10},{k,2,n}] (* Gus Wiseman, Jun 07 2019 *)

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)

Corrected and rewritten by Gus Wiseman, Jun 07 2019

A356065 Squarefree numbers whose prime indices are all prime-powers.

Original entry on oeis.org

1, 3, 5, 7, 11, 15, 17, 19, 21, 23, 31, 33, 35, 41, 51, 53, 55, 57, 59, 67, 69, 77, 83, 85, 93, 95, 97, 103, 105, 109, 115, 119, 123, 127, 131, 133, 155, 157, 159, 161, 165, 177, 179, 187, 191, 201, 205, 209, 211, 217, 227, 231, 241, 249, 253, 255, 265, 277

Offset: 1

Gus Wiseman, Jul 25 2022

nonn

			105 has prime indices {2,3,4}, all three of which are prime-powers, so 105 is in the sequence.

The multiplicative version (factorizations) is A050361, non-strict A000688.
Heinz numbers of the partitions counted by A054685, with 1's A106244, non-strict A023894, non-strict with 1's A023893.
Counting twice-partitions of this type gives A279786, non-strict A279784.
Counting twice-factorizations gives A295935, non-strict A296131.
These are the odd products of distinct elements of A302493.
Allowing prime index 1 gives A302496, non-strict A302492.
The case of primes (instead of prime-powers) is A302590, non-strict A076610.
These are the squarefree positions of 1's in A355741.
This is the squarefree case of A355743, complement A356066.
A001222 counts prime-power divisors.
A005117 lists the squarefree numbers.
A034699 gives maximal prime-power divisor.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
A355742 chooses a prime-power divisor of each prime index.
Cf. A001970, A055887, A063834, A302601, A355731, A355744, A356064.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],SquareFreeQ[#]&&And@@PrimePowerQ/@primeMS[#]&]

Intersection of A005117 and A355743.

A377052 Antidiagonal-sums of the array A377051(n,k) = n-th term of k-th differences of powers of primes.

Original entry on oeis.org

1, 3, 4, 5, 6, 13, -6, 45, -50, 113, -98, 73, 274, -1159, 3563, -8707, 19024, -36977, 64582, -98401, 121436, -81961, -147383, 860871, -2709964, 7110655, -17077217, 38873213, -85085216, 179965720, -367884935, 725051361, -1372311916, 2481473639, -4257624155

Offset: 0

Gus Wiseman, Oct 22 2024

sign

These are the row-sums of the triangle-version of A377051.

			The sixth antidiagonal of A377051 is (8, 1, -1, -2, -3, -4, -5), so a(6) = -6.

The version for primes is A140119, noncomposites A376683, composites A377034.
For squarefree numbers we have A377039, nonsquarefree A377047.
These are the antidiagonal-sums of A377051.
The unsigned version is A377053.
For leaders we have A377054, for primes A007442 or A030016.
For first zero-positions we have A377055.
A version for partitions is A377056, cf. A175804, A053445, A281425, A320590.
A000040 lists the primes, differences A001223, seconds A036263.
A001597 lists perfect-powers, complement A007916.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.
Cf. A000961, A025475, A053707, A057820, A093555, A174965, A246655, A361102, A376340, A376596, A376598.

nn=20;
t=Table[Differences[NestList[NestWhile[#+1&, #+1,!PrimePowerQ[#]&]&,1,2*nn],k],{k,0,nn}];
Total/@Table[t[[j,i-j+1]],{i,nn},{j,i}]

A280195 Expansion of 1/(1 - Sum_{k>=2} floor(1/omega(k))*x^k), where omega(k) is the number of distinct prime factors (A001221).

Original entry on oeis.org

1, 0, 1, 1, 2, 3, 4, 8, 11, 19, 28, 47, 72, 116, 182, 289, 460, 724, 1153, 1820, 2891, 4572, 7249, 11482, 18190, 28821, 45651, 72338, 114582, 181549, 287596, 455647, 721847, 1143588, 1811748, 2870239, 4547232, 7203907, 11412882, 18080833, 28644680, 45380392, 71894054, 113898439, 180443915, 285869028, 452888824, 717490903, 1136687237

Offset: 0

Ilya Gutkovskiy, Dec 28 2016

nonn

Number of compositions (ordered partitions) into prime powers (1 excluded).

			a(6) = 4 because we have [4, 2], [3, 3], [2, 4] and [2, 2, 2].

G. C. Greubel, Table of n, a(n) for n = 0..1000
Eric Weisstein's World of Mathematics, Prime Power
Index entries for sequences related to compositions

Cf. A001221, A023360, A023894, A246655.

nmax = 48; CoefficientList[Series[1/(1 - Sum[Floor[1/PrimeNu[k]] x^k, {k, 2, nmax}]), {x, 0, nmax}], x]

G.f.: 1/(1 - Sum_{k>=2} floor(1/omega(k))*x^k).

A321346 Number of integer partitions of n containing no prime powers > 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 5, 5, 6, 7, 8, 8, 11, 11, 14, 16, 19, 19, 25, 26, 31, 34, 40, 41, 52, 54, 63, 69, 81, 86, 105, 109, 126, 137, 160, 169, 201, 211, 242, 264, 303, 320, 375, 396, 453, 490, 557, 590, 682, 726, 823, 888, 1002, 1065, 1219

Offset: 0

Gus Wiseman, Dec 11 2018

nonn

First differs from A285798 at a(30) = 52, A285798(30) = 51.

			The a(20) = 14 integer partitions:
  (20)
  (10,10)
  (14,6)
  (18,1,1)
  (12,6,1,1)
  (6,6,6,1,1)
  (10,6,1,1,1,1)
  (15,1,1,1,1,1)
  (14,1,1,1,1,1,1)
  (12,1,1,1,1,1,1,1,1)
  (6,6,1,1,1,1,1,1,1,1)
  (10,1,1,1,1,1,1,1,1,1,1)
  (6,1,1,1,1,1,1,1,1,1,1,1,1,1,1)
  (1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1)

Fausto A. C. Cariboni, Table of n, a(n) for n = 0..350

Cf. A000607, A000961, A001597, A002095, A023893, A023894, A096258, A246655, A320322, A321347, A321378, A321665, A322452, A322454.

nn=100;
ser=Product[If[PrimePowerQ[n],1,1/(1-x^n)],{n,nn}];
CoefficientList[Series[ser,{x,0,nn}],x]

A356068 Number of integers ranging from 1 to n that are not prime-powers (1 is not a prime-power).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 2, 2, 3, 3, 4, 4, 5, 6, 6, 6, 7, 7, 8, 9, 10, 10, 11, 11, 12, 12, 13, 13, 14, 14, 14, 15, 16, 17, 18, 18, 19, 20, 21, 21, 22, 22, 23, 24, 25, 25, 26, 26, 27, 28, 29, 29, 30, 31, 32, 33, 34, 34, 35, 35, 36, 37, 37, 38, 39, 39, 40, 41, 42

Offset: 1

Gus Wiseman, Jul 31 2022

nonn

			The a(30) = 14 numbers: 1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30.

The complement is counted by A025528, with 1's A065515.
For primes instead of prime-powers we have A062298, with 1's A065855.
The version treating 1 as a prime-power is A085970.
One more than the partial sums of A143731.
A000688 counts factorizations into prime-powers.
A001222 counts prime-power divisors.
A246655 lists the prime-powers (A000961 includes 1), towers A164336.
Cf. A000720, A023894, A036234, A050361, A054685, A069513, A355743, A356064, A356065, A356066.

Table[Length[Select[Range[n],!PrimePowerQ[#]&]],{n,100}]

a(n) = A085970(n) + 1.

A376653 Sorted positions of first appearances in the second differences of consecutive prime-powers inclusive (A000961).

Original entry on oeis.org

1, 4, 5, 10, 12, 18, 25, 45, 47, 48, 60, 68, 69, 71, 80, 118, 121, 178, 179, 199, 206, 207, 216, 244, 245, 304, 325, 327, 402, 466, 484, 605, 801, 880, 939, 1033, 1055, 1077, 1234, 1281, 1721, 1890, 1891, 1906, 1940, 1960, 1962, 2257, 2290, 2410, 2880, 3150

Offset: 1

Gus Wiseman, Oct 06 2024

nonn

			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, 3, ...
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 first appearances (A376653):
  1, 4, 5, 10, 12, 18, 25, 45, 47, 48, 60, 68, 69, 71, 80, 118, 121, 178, 179, 199, ...

For first differences we had A057820, sorted firsts A376340(n)+1 (except first term).
These are the sorted positions of first appearances in A376596.
The exclusive version is a(n) - 1 = A376654(n), except first term.
For squarefree instead of prime-power we have A376655.
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 (inflections and undulations), A376598 (nonzero curvature).
For second differences: A036263 (prime), A073445 (composite), A376559 (perfect-power), A376562 (non-perfect-power), A376590 (squarefree), A376593 (nonsquarefree), A376599 (non-prime-power).
Cf. A053707, A069623, A093555, A174965, A251092, A333254, A361102.

q=Differences[Select[Range[100],#==1||PrimePowerQ[#]&],2];
Select[Range[Length[q]],!MemberQ[Take[q,#-1],q[[#]]]&]

A377053 Antidiagonal-sums of the absolute value of the array A377051(n,k) = n-th term of k-th differences of powers of primes.

Original entry on oeis.org

1, 3, 4, 5, 6, 13, 24, 45, 80, 123, 174, 229, 382, 1219, 3591, 8849, 19288, 37899, 67442, 108323, 156054, 206733, 311525, 860955, 2710374, 7111657, 17080759, 38884849, 85124764, 180097856, 368321633, 726482493, 1377039690, 2496856437, 4306569569, 7016267449

Offset: 0

Gus Wiseman, Oct 22 2024

nonn

These are the row-sums of the absolute value of the triangle-version of A377051.

			The sixth antidiagonal of A377051 is (8, 1, -1, -2, -3, -4, -5), so a(6) = 24.

The version for primes is A376681, noncomposites A376684, composites A377035.
For squarefree numbers we have A377040, nonsquarefree A377048.
This is the antidiagonal-sums of the absolute value of A377051.
The signed version is A377052.
For leaders we have A377054, for primes A007442 or A030016.
For first zero-positions we have A377055.
A version for partitions is A377056, cf. A175804, A053445, A281425, A320590.
A000040 lists the primes, differences A001223, seconds A036263.
A008578 lists the noncomposites, differences A075526.
A023893 and A023894 count integer partitions into prime-powers, factorizations A000688.
Cf. A000961, A025475, A053707, A057820, A069623, A093555, A174965, A246655, A361102, A376340, A376596.

nn=20;
t=Table[Differences[NestList[NestWhile[#+1&, #+1,!PrimePowerQ[#]&]&,1,2*nn],k],{k,0,nn}];
Total/@Abs[Table[t[[j,i-j+1]],{i,nn},{j,i}]]

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A051613 a(n) = partitions of n into powers of distinct primes (1 not considered a power).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A356065 Squarefree numbers whose prime indices are all prime-powers.

Views

Author

Keywords

Examples

Crossrefs

Programs

Mathematica

Formula

A377052 Antidiagonal-sums of the array A377051(n,k) = n-th term of k-th differences of powers of primes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A280195 Expansion of 1/(1 - Sum_{k>=2} floor(1/omega(k))*x^k), where omega(k) is the number of distinct prime factors (A001221).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A321346 Number of integer partitions of n containing no prime powers > 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A356068 Number of integers ranging from 1 to n that are not prime-powers (1 is not a prime-power).

Views

Author