A319630 -id:A319630 - cp's OEIS frontend

A319626 Primorial deflation of n (numerator): Let f be the completely multiplicative function over the positive rational numbers defined by f(p) = A034386(p) for any prime number p; f constitutes a permutation of the positive rational numbers; let g be the inverse of f; for any n > 0, a(n) is the numerator of g(n).

1, 2, 3, 4, 5, 3, 7, 8, 9, 10, 11, 6, 13, 14, 5, 16, 17, 9, 19, 20, 21, 22, 23, 12, 25, 26, 27, 28, 29, 5, 31, 32, 33, 34, 7, 9, 37, 38, 39, 40, 41, 21, 43, 44, 15, 46, 47, 24, 49, 50, 51, 52, 53, 27, 55, 56, 57, 58, 59, 10, 61, 62, 63, 64, 65, 33, 67, 68, 69

Offset: 1

Rémy Sigrist, Sep 25 2018

See A319627 for the corresponding denominators.
The restriction of f to the natural numbers corresponds to A108951.
The function g is completely multiplicative over the positive rational numbers with g(2) = 2 and g(q) = q/p for any pair (p, q) of consecutive prime numbers.
The ratio A319626(n)/A319627(n) can be viewed as a "primorial deflation" of n (see also A329900), with the inverse operation being n = A108951(A319626(n)) / A108951(A319627(n)), where A319627(k) = 1 for all k in A025487. - Daniel Suteu, Dec 29 2019

			f(21/5) = (2*3) * (2*3*5*7) / (2*3*5) = 42, hence g(42) = 21/5 and a(42) = 21.

Daniel Suteu, Table of n, a(n) for n = 1..10000

A left inverse of A108951. Coincides with A329900 on A025487.

Cf. A006530, A053585, A064989, A181815, A307035, A319627, A319630, A329902, A330749, A330750 (rgs-transform), A330751 (ordinal transform).

Array[#1/GCD[#1, #2] & @@ {#, Apply[Times, Map[If[#1 <= 2, 1, NextPrime[#1, -1]]^#2 & @@ # &, FactorInteger[#]]]} &, 120] (* Michael De Vlieger, Aug 27 2020 *)

a(n) = my (f=factor(n)); numerator(prod(i=1, #f~, my (p=f[i,1]); (p/if (p>2, precprime(p-1), 1))^f[i,2]))

a(n) = n / gcd(n, A064989(n)) = n / A330749(n).
a(n) <= n with equality iff n belongs to A319630.
A006530(a(n)) = A006530(n).
A053585(a(n)) = A053585(n).
From Antti Karttunen, Dec 29 2019: (Start)
a(A108951(n)) = n.
a(A025487(n)) = A329900(A025487(n)) = A181815(n).
Many of the formulas given in A329900 apply here as well:
a(n!) = A307035(n), a(A002182(n)) = A329902(n), and so on.
(End)

"Primorial deflation" prefixed to the name by Antti Karttunen, Dec 29 2019

A350842 Number of integer partitions of n with no difference -2.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 16, 24, 30, 40, 54, 69, 89, 118, 146, 187, 239, 297, 372, 468, 575, 711, 880, 1075, 1314, 1610, 1947, 2359, 2864, 3438, 4135, 4973, 5936, 7090, 8466, 10044, 11922, 14144, 16698, 19704, 23249, 27306, 32071, 37639, 44019, 51457, 60113

Offset: 0

Gus Wiseman, Jan 20 2022

nonn

			The a(1) = 1 through a(7) = 12 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (21)   (22)    (32)     (33)      (43)
             (111)  (211)   (41)     (51)      (52)
                    (1111)  (221)    (222)     (61)
                            (2111)   (321)     (322)
                            (11111)  (411)     (511)
                                     (2211)    (2221)
                                     (21111)   (3211)
                                     (111111)  (4111)
                                               (22111)
                                               (211111)
                                               (1111111)

Gus Wiseman, Sequences counting and ranking partitions and compositions by their differences and quotients.

Heinz number rankings are in parentheses below.
The version for no difference 0 is A000009.
The version for subsets of prescribed maximum is A005314.
The version for all differences < -2 is A025157, non-strict A116932.
The version for all differences > -2 is A034296, strict A001227.
The opposite version is A072670.
The version for no difference -1 is A116931 (A319630), strict A003114.
The multiplicative version is A350837 (A350838), strict A350840.
The strict case is A350844.
The complement for quotients is counted by A350846 (A350845).
A000041 = integer partitions.
A027187 = partitions of even length.
A027193 = partitions of odd length (A026424).
A323092 = double-free partitions (A320340), strict A120641.
A325534 = separable partitions (A335433).
A325535 = inseparable partitions (A335448).
A350839 = partitions with a gap and conjugate gap (A350841).
Cf. A000070, A000929, A001511, A003242, A007359, A018819, A040039, A045690, A045691, A101417, A154402, A323093.

Table[Length[Select[IntegerPartitions[n],FreeQ[Differences[#],-2]&]],{n,0,30}]

A048674 Fixed points of A048673 and A064216: Numbers n such that if n = product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) = (2*n)-1, where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

1, 2, 3, 25, 26, 33, 93, 1034, 970225, 8550146, 325422273, 414690595, 1864797542, 2438037206

Offset: 1

Antti Karttunen, Jul 14 1999

Equally: after 1, numbers n such that, if the prime factorization of 2n-1 = Product_{k >= 1} (p_k)^(c_k) then Product_{k >= 1} (p_{k-1})^(c_k) = n.
Factorization of the initial terms: 1, 2, 3, 5^2, 2*13, 3*11, 3*31, 2*11*47, 5^2*197^2, 2*11*47*8269, 3*11*797*12373, 5*11^2*433*1583, 2*23*59*101*6803, 2*11*53*1201*1741.
The only 3-cycle of permutation A048673 in range 1 .. 402653184 is (2821 3460 5639).
For 2-cycles, take setwise difference of A245449 and this sequence.
Numbers k for which A336853(k) = k-1. - Antti Karttunen, Nov 26 2021

			25 is present, as 2*25 - 1 = 49 = p_4^2, and p_3^2 = 5*5 = 25.
26 is present, as 2*26 - 1 = 51 = 3*17 = p_2 * p_8, and p_1 * p_7 = 2*13 = 26.
Alternatively, as 26 = 2*13 = p_1 * p_7, and ((p_2 * p_8)+1)/2 = ((3*17)+1)/2 = 26 also, thus 26 is present.

Fixed points of permutation pair A048673/A064216.
Positions of zeros in A349573.
Subsequence of the following sequences: A245449, A269860, A319630, A349622, A378980 (see also A379216).
This sequence is also obtained as a setwise difference of the following pairs of sequences: A246281 \ A246351, A246352 \ A246282, A246361 \ A246371, A246372 \ A246362.
Cf. A000040, A003961 (A045965), A064989, A285701, A336853.
Cf. also A348514 (fixed points of map A108228, similar to A048673).

A048673 := n -> (A003961(n)+1)/2;
A048674list := proc(upto_n) local b,i; b := [ ]; for i from 1 to upto_n do if(A048673(i) = i) then b := [ op(b), i ]; fi; od: RETURN(b); end;

Join[{1}, Reap[For[n = 1, n < 10^7, n++, ff = FactorInteger[n]; If[Times @@ Power @@@ (NextPrime[ff[[All, 1]]]^ff[[All, 2]]) == 2 n - 1, Print[n]; Sow[n]]]][[2, 1]]] (* Jean-François Alcover, Mar 04 2016 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA048674(n) = ((n+n)==(1+A003961(n))); \\ Antti Karttunen, Nov 26 2021

Entry revised and the names in Maple-code cleaned by Antti Karttunen, Aug 25 2014

Terms a(11) - a(14) added by Antti Karttunen, Sep 11-13 2014

A325160 Products of distinct, non-consecutive primes. Squarefree numbers not divisible by any two consecutive primes.

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 31, 33, 34, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 67, 69, 71, 73, 74, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 110, 111, 113, 115, 118, 119

Offset: 1

Gus Wiseman, Apr 05 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions into distinct non-consecutive parts (counted by A003114). The nonsquarefree case is A319630, which gives the Heinz numbers of integer partitions with no consecutive parts (counted by A116931).

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 6, 52, 515, 5146, 51435, 514416, 5144232, 51442384, ... . Apparently, the asymptotic density of this sequence exists and equals 0.51442... . - Amiram Eldar, Sep 24 2022

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   5: {3}
   7: {4}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  17: {7}
  19: {8}
  21: {2,4}
  22: {1,5}
  23: {9}
  26: {1,6}
  29: {10}
  31: {11}
  33: {2,5}
  34: {1,7}
  37: {12}

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001227, A003114, A005117, A025157, A034296, A056239, A073485, A073491, A089995, A112798, A116931, A319630, A325161, A325162.

Select[Range[100],Min@@Differences[Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]>1&]

isok(k) = {if (issquarefree(k), my(v = apply(primepi, factor(k)[,1])); ! #select(x->(v[x+1]-v[x] == 1), [1..#v-1]));} \\ Michel Marcus, Jan 09 2021

A104210 Positive integers divisible by at least 2 consecutive primes.

Original entry on oeis.org

6, 12, 15, 18, 24, 30, 35, 36, 42, 45, 48, 54, 60, 66, 70, 72, 75, 77, 78, 84, 90, 96, 102, 105, 108, 114, 120, 126, 132, 135, 138, 140, 143, 144, 150, 154, 156, 162, 165, 168, 174, 175, 180, 186, 192, 195, 198, 204, 210, 216, 221, 222, 225, 228, 231, 234, 240

Offset: 1

Leroy Quet, Mar 13 2005

nonn

If a perfect square is in this sequence, then so is its square root (e.g., 144 and 12). - Alonso del Arte, May 07 2012

The numbers of terms not exceeding 10^k, for k=1,2,..., are 1, 22, 242, 2456, 24632, 246414, 2464272, 24643281, 246433426, ... Apparently, the asymptotic density of this sequence is 0.24643... - Amiram Eldar, Apr 10 2021

			35 is divisible by both 5 and 7, and 5 and 7 are consecutive primes.
77 is divisible by both 7 and 11, and 7 and 11 are consecutive primes.
110 is not in the sequence because, although it is divisible by 2, 5 and 11, it is not divisible by 3 or 7.

Cf. A003961, A296210 (characteristic function), A319630 (complement), A379230 [= A252748(a(n))].
Positions of terms larger than 1 in A300820 and in A322361.
Subsequences: A006094, A349169 (conjectured, after its initial 1), A349176, A355527 (squarefree terms), A372566, A378884, A379232.

N:= 1000: # for terms <= N
R:= {}:
p:= 2:
do
  q:= p; p:= nextprime(p);
  if p*q > N then break fi;
  R:= R union {seq(i,i=p*q..N,p*q)}
od:
sort(convert(R,list)); # Robert Israel, Apr 13 2020

fQ[n_] := Block[{lst = PrimePi /@ Flatten[ Table[ #[[1]], {1}] & /@ FactorInteger[n]]}, Count[ Drop[lst, 1] - Drop[lst, -1], 1] > 0]; Select[ Range[244], fQ[ # ] &] (* Robert G. Wilson v, Mar 16 2005 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
is_A104210(n) = (gcd(n,A003961(n))>1); \\ Antti Karttunen, Dec 24 2024

{k such that gcd(k, A003961(k)) > 1}. - Antti Karttunen, Dec 24 2024

More terms from Robert G. Wilson v, Mar 16 2005

A349162 a(n) = sigma(n) / gcd(sigma(n), A003961(n)), where A003961 shifts the prime factorization of n one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 1, 4, 7, 6, 4, 8, 5, 13, 6, 12, 28, 14, 8, 24, 31, 18, 13, 20, 2, 32, 12, 24, 4, 31, 14, 8, 56, 30, 24, 32, 7, 48, 18, 48, 91, 38, 20, 56, 10, 42, 32, 44, 28, 78, 24, 48, 124, 57, 31, 72, 98, 54, 8, 72, 40, 16, 30, 60, 8, 62, 32, 104, 127, 12, 48, 68, 14, 96, 48, 72, 13, 74, 38, 124, 140, 96, 56, 80, 62, 121, 42

Offset: 1

Antti Karttunen, Nov 09 2021

Denominator of ratio A003961(n) / A000203(n).
Small values are rare, but are not limited to the beginning. For example in range 1 .. 2^25, a(n) = 4 at n = 3, 6, 24, 792, 2720, 122944, 31307472.
Question: Would it be possible to prove that a(n) > 1 for all n > 2?
Obviously, 1's may occur only on squares & twice squares (A028982). See also comments in A350072. - Antti Karttunen, Feb 16 2022

Cf. A000203, A003961, A028982 (positions of odd terms), A319630, A336702, A342671, A348992 (the odd part), A348993, A349161 (numerators), A349163, A349164, A349627, A349628, A350072 [= a(n^2)].

Cf. also A349745, A351551, A351554.

Array[#1/GCD[##] & @@ {DivisorSigma[1, #], If[# == 1, 1, Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]]} &, 82] (* Michael De Vlieger, Nov 11 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A349162(n) = { my(s=sigma(n)); (s/gcd(s,A003961(n))); };

a(n) = A000203(n) / A342671(n) = A000203(n) / gcd(A000203(n), A003961(n)).

A350839 Number of integer partitions of n with a difference < -1 and a conjugate difference < -1.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 2, 3, 7, 11, 17, 26, 39, 54, 81, 108, 148, 201, 269, 353, 467, 601, 779, 995, 1272, 1605, 2029, 2538, 3171, 3941, 4881, 6012, 7405, 9058, 11077, 13478, 16373, 19817, 23953, 28850, 34692, 41599, 49802, 59461, 70905, 84321, 100155, 118694

Offset: 0

Gus Wiseman, Jan 24 2022

nonn

We define a difference of a partition to be a difference of two adjacent parts.

			The a(5) = 1 through a(10) = 17 partitions:
  (311)  (411)   (511)    (422)     (522)      (622)
         (3111)  (4111)   (611)     (711)      (811)
                 (31111)  (3311)    (4221)     (4222)
                          (4211)    (4311)     (4411)
                          (5111)    (5211)     (5221)
                          (41111)   (6111)     (5311)
                          (311111)  (33111)    (6211)
                                    (42111)    (7111)
                                    (51111)    (42211)
                                    (411111)   (43111)
                                    (3111111)  (52111)
                                               (61111)
                                               (331111)
                                               (421111)
                                               (511111)
                                               (4111111)
                                               (31111111)

Allowing -1 gives A144300 = non-constant partitions.
Taking one of the two conditions gives A239955, ranked by A073492, A065201.
These partitions are ranked by A350841.
A000041 = integer partitions, strict A000009.
A034296 = flat (contiguous) partitions, strict A001227.
A073491 = numbers whose prime indices have no gaps, strict A137793.
A090858 = partitions with a single hole, ranked by A325284.
A116931 = partitions with differences != -1, strict A003114.
A116932 = partitions with differences != -1 or -2, strict A025157.
A277103 = partitions with the same number of odd parts as their conjugate.
A350837 = partitions with no adjacent doublings, strict A350840.
A350842 = partitions with differences != -2, strict A350844, sets A005314.
Cf. A000070, A000929, A008284, A040039, A183558, A319630, A321440, A350838.

conj[y_]:=If[Length[y]==0,y,Table[Length[Select[y,#>=k&]],{k,1,Max[y]}]];
Table[Length[Select[IntegerPartitions[n],(Min@@Differences[#]<-1)&&(Min@@Differences[conj[#]]<-1)&]],{n,0,30}]

A350844 Number of strict integer partitions of n with no difference -2.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 3, 4, 4, 7, 7, 8, 11, 12, 15, 18, 21, 23, 31, 32, 40, 45, 54, 59, 73, 78, 94, 106, 122, 136, 161, 177, 203, 231, 259, 293, 334, 372, 417, 476, 525, 592, 663, 742, 821, 931, 1020, 1147, 1271, 1416, 1558, 1752, 1916, 2137, 2357, 2613, 2867

Offset: 0

Gus Wiseman, Jan 21 2022

nonn

			The a(1) = 1 through a(12) = 11 partitions (A..C = 10..12):
  1   2   3    4   5    6     7    8     9     A      B     C
          21       32   51    43   62    54    73     65    84
                   41   321   52   71    63    82     74    93
                              61   521   72    91     83    A2
                                         81    541    92    B1
                                         432   721    A1    543
                                         621   4321   632   651
                                                      821   732
                                                            741
                                                            921
                                                            6321

Gus Wiseman, Sequences counting and ranking partitions and compositions by their differences and quotients.

The version for no difference 0 is A000009.
The version for no difference > -2 is A001227, non-strict A034296.
The version for no difference -1 is A003114 (A325160).
The version for subsets of prescribed maximum is A005314.
The version for all differences < -2 is A025157, non-strict A116932.
The opposite version is A072670.
The multiplicative version is A350840, non-strict A350837 (A350838).
The non-strict version is A350842.
A000041 counts integer partitions.
A027187 counts partitions of even length.
A027193 counts partitions of odd length (A026424).
A116931 counts partitions with no difference -1 (A319630).
A323092 counts double-free integer partitions (A320340) strict A120641.
A325534 counts separable partitions (A335433).
A325535 counts inseparable partitions (A335448).
Cf. A000929, A003000, A018819, A040039, A045690, A045691, A154402, A303362, A323094, A342095, A342097.

Table[Length[Select[IntegerPartitions[n],FreeQ[Differences[#],0|-2]&]],{n,0,30}]

A356236 Number of integer partitions of n with a neighborless part.

Original entry on oeis.org

0, 1, 2, 2, 4, 4, 8, 9, 16, 20, 31, 40, 59, 76, 105, 138, 184, 238, 311, 400, 515, 656, 831, 1052, 1322, 1659, 2064, 2572, 3182, 3934, 4837, 5942, 7264, 8872, 10789, 13109, 15865, 19174, 23105, 27796, 33361, 39956, 47766, 56985, 67871, 80675, 95750, 113416

Offset: 0

Gus Wiseman, Aug 24 2022

nonn

A part x of a partition is neighborless if neither x - 1 nor x + 1 are parts.

			The a(1) = 1 through a(8) = 9 partitions:
  (1)  (2)   (3)    (4)     (5)      (6)       (7)
       (11)  (111)  (22)    (41)     (33)      (52)
                    (31)    (311)    (42)      (61)
                    (1111)  (11111)  (51)      (331)
                                     (222)     (421)
                                     (411)     (511)
                                     (3111)    (4111)
                                     (111111)  (31111)
                                               (1111111)

The complement is counted by A355394, singleton case A355393.
The singleton case is A356235, ranked by A356237.
The strict case is A356607, complement A356606.
These partitions are ranked by the complement of A356736.
A000041 counts integer partitions, strict A000009.
A000837 counts relatively prime partitions, ranked by A289509.
A007690 counts partitions with no singletons, complement A183558.
Cf. A066205, A112798, A319630, A325160, A328171, A328172, A328187, A328221.

Table[Length[Select[IntegerPartitions[n],Function[ptn,Or@@Table[!MemberQ[ptn,x-1]&&!MemberQ[ptn,x+1],{x,Union[ptn]}]]]],{n,0,30}]

a(n) = A000041(n) - A355394(n).

A349164 a(n) = A064989(A003961(n) / gcd(sigma(n), A003961(n))), where A003961 shifts the prime factorization of n one step towards larger primes, while A064989 shifts it back towards smaller primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 1, 3, 4, 5, 3, 7, 4, 9, 5, 11, 12, 13, 7, 15, 16, 17, 9, 19, 2, 21, 11, 23, 4, 25, 13, 9, 28, 29, 15, 31, 8, 33, 17, 35, 36, 37, 19, 39, 10, 41, 21, 43, 22, 45, 23, 47, 48, 49, 25, 51, 52, 53, 9, 55, 28, 19, 29, 59, 6, 61, 31, 63, 64, 13, 33, 67, 17, 69, 35, 71, 12, 73, 37, 75, 76, 77, 39, 79, 40, 81, 41, 83, 84

Offset: 1

Antti Karttunen, Nov 09 2021

nonn

Cf. A000203, A003961, A319630, A326042, A336702, A342671, A348993, A349161, A349162, A349169, A349174.

Cf. A349144 and A349168 [positions where a(n) is / is not relatively prime with A349163(n) = n/a(n)].

Array[Times @@ Map[If[#1 <= 2, 1, NextPrime[#1, -1]]^#2 & @@ # &, FactorInteger[#2/GCD[##]]] & @@ {DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &, 84] (* Michael De Vlieger, Nov 11 2021 *)

A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A064989(n) = { my(f=factor(n)); if((n>1 && f[1,1]==2), f[1,2] = 0); for(i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); };
A349164(n) = { my(u=A003961(n)); A064989(u/gcd(u,sigma(n))); };

a(n) = A064989(A349161(n)).

a(n) = n / A349163(n).

cp's OEIS Frontend

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A350842 Number of integer partitions of n with no difference -2.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

A048674 Fixed points of A048673 and A064216: Numbers n such that if n = product_{k >= 1} (p_k)^(c_k), then Product_{k >= 1} (p_{k+1})^(c_k) = (2*n)-1, where p_k indicates the k-th prime, A000040(k).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

PARI

Extensions

A325160 Products of distinct, non-consecutive primes. Squarefree numbers not divisible by any two consecutive primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A104210 Positive integers divisible by at least 2 consecutive primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A349162 a(n) = sigma(n) / gcd(sigma(n), A003961(n)), where A003961 shifts the prime factorization of n one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A350839 Number of integer partitions of n with a difference < -1 and a conjugate difference < -1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A350844 Number of strict integer partitions of n with no difference -2.