A322361 -id:A322361 - cp's OEIS frontend

A372572 Lexicographically earliest infinite sequence such that a(i) = a(j) => A009194(i) = A009194(j), A322361(i) = A322361(j) and A342671(i) = A342671(j), for all i, j >= 1.

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

Offset: 1

Antti Karttunen, May 24 2024

nonn

Restricted growth sequence transform of the triple [A009194(n), A322361(n), A342671(n)].
For all i, j:
  a(i) = a(j) => A349167(i) = A349167(j),
  a(i) = a(j) => A353666(i) = A353666(j),
  a(i) = a(j) => A372565(i) = A372565(j).

Cf. A009194, A322361, A342671, A349167, A353666, A372565.

up_to = 65537;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
Aux372572(n) = [gcd(n, sigma(n)), gcd(n, A003961(n)), gcd(sigma(n), A003961(n))];
v372572 = rgs_transform(vector(up_to, n, Aux372572(n)));
A372572(n) = v372572[n];

A372570 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = [A009194(n), A009195(n), A009223(n), A322361(n), A342671(n)], for all i, j >= 1.

Original entry on oeis.org

1, 2, 3, 4, 3, 5, 3, 6, 7, 8, 3, 9, 3, 10, 11, 12, 3, 13, 3, 14, 15, 8, 3, 16, 17, 10, 18, 19, 3, 20, 3, 21, 22, 8, 23, 24, 3, 10, 25, 26, 3, 27, 3, 28, 29, 8, 3, 30, 31, 32, 33, 34, 3, 35, 36, 37, 38, 8, 3, 39, 3, 10, 40, 41, 42, 43, 3, 44, 22, 45, 3, 46, 3, 10, 47, 48, 49, 50, 3, 51, 52, 8, 3, 53, 54, 10, 55, 56, 3, 57, 58, 28, 15, 8, 59, 60, 3, 61, 62, 63, 3

Offset: 1

Antti Karttunen, May 25 2024

nonn

Restricted growth sequence transform of the quintuple [A009194(n), A009195(n), A009223(n), A322361(n), A342671(n)].
For all i, j:
  A305801(i) = A305801(j) => a(i) = a(j),
  a(i) = a(j) => A372569(i) = A372569(j),
  a(i) = a(j) => A372572(i) = A372572(j).

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

Cf. A009194, A009195, A009223, A322361, A342671.

Cf. also A305801, A372569, A372572.

up_to = 65537;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
Aux372570(n) = [gcd(n, sigma(n)), gcd(n, eulerphi(n)), gcd(eulerphi(n), sigma(n)), gcd(n, A003961(n)), gcd(sigma(n), A003961(n))];
v372570 = rgs_transform(vector(up_to, n, Aux372570(n)));
A372570(n) = v372570[n];

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Original entry on oeis.org

1, 3, 5, 9, 7, 15, 11, 27, 25, 21, 13, 45, 17, 33, 35, 81, 19, 75, 23, 63, 55, 39, 29, 135, 49, 51, 125, 99, 31, 105, 37, 243, 65, 57, 77, 225, 41, 69, 85, 189, 43, 165, 47, 117, 175, 87, 53, 405, 121, 147, 95, 153, 59, 375, 91, 297, 115, 93, 61, 315, 67, 111, 275, 729, 119

Offset: 1

Marc LeBrun

Meyers (see Guy reference) conjectures that for all r >= 1, the least odd number not in the set {a(i): i < prime(r)} is prime(r+1). - N. J. A. Sloane, Jan 08 2021
Meyers' conjecture would be refuted if and only if for some r there were such a large gap between prime(r) and prime(r+1) that there existed a composite c for which prime(r) < c < a(c) < prime(r+1), in which case (by Bertrand's postulate) c would necessarily be a term of A246281. - Antti Karttunen, Mar 29 2021
a(n) is odd for all n and for each odd m there exists a k with a(k) = m (see A064216). a(n) > n for n > 1: bijection between the odd and all numbers. - Reinhard Zumkeller, Sep 26 2001
a(n) and n have the same number of distinct primes with (A001222) and without multiplicity (A001221). - Michel Marcus, Jun 13 2014
From Antti Karttunen, Nov 01 2019: (Start)
More generally, a(n) has the same prime signature as n, A046523(a(n)) = A046523(n). Also A246277(a(n)) = A246277(n) and A287170(a(n)) = A287170(n).
Many permutations and other sequences that employ prime factorization of n to encode either polynomials, partitions (via Heinz numbers) or multisets in general can be easily defined by using this sequence as one of their constituent functions. See the last line in the Crossrefs section for examples.
(End)

			a(12) = a(2^2 * 3) = a(prime(1)^2 * prime(2)) = prime(2)^2 * prime(3) = 3^2 * 5 = 45.
a(A002110(n)) = A002110(n + 1) / 2.

Richard K. Guy, editor, Problems From Western Number Theory Conferences, Labor Day, 1983, Problem 367 (Proposed by Leroy F. Meyers, The Ohio State U.).

Indranil Ghosh, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)
Index entries for sequences computed from indices in prime factorization.
Index entries for sequences related to Heinz numbers.

See A045965 for another version.
Row 1 of table A242378 (which gives the "k-th powers" of this sequence), row 3 of A297845 and of A306697. See also arrays A066117, A246278, A255483, A308503, A329050.
Cf. A064989 (a left inverse), A064216, A000040, A002110, A000265, A027746, A046523, A048673 (= (a(n)+1)/2), A108228 (= (a(n)-1)/2), A191002 (= a(n)*n), A252748 (= a(n)-2n), A286385 (= a(n)-sigma(n)), A283980 (= a(n)*A006519(n)), A341529 (= a(n)*sigma(n)), A326042, A049084, A001221, A001222, A122111, A225546, A260443, A245606, A244319, A246269 (= A065338(a(n))), A322361 (= gcd(n, a(n))), A305293.
Cf. A191555, A252738.
Cf. A249734, A249735 (bisections).
Cf. A246261 (a(n) is of the form 4k+1), A246263 (of the form 4k+3), A246271, A246272, A246259, A246281 (n such that a(n) < 2n), A246282 (n such that a(n) > 2n), A252742.
Cf. A275717 (a(n) > a(n-1)), A275718 (a(n) < a(n-1)).
Cf. A003972 (Möbius transform), A003973 (Inverse Möbius transform), A318321.
Cf. A300841, A305421, A322991, A250469, A269379 for analogous shift-operators in other factorization and quasi-factorization systems.
Cf. also following permutations and other sequences that can be defined with the help of this sequence: A005940, A163511, A122111, A260443, A206296, A265408, A265750, A275733, A275735, A297845, A091202 & A091203, A250245 & A250246, A302023 & A302024, A302025 & A302026.
A version for partition numbers is A003964, strict A357853.
A permutation of A005408.
Applying the same transformation again gives A357852.
Other multiplicative sequences: A064988, A357977, A357978, A357980, A357983.
A056239 adds up prime indices, row-sums of A112798.
Cf. A000720, A076610, A296150.

a003961 1 = 1
a003961 n = product $ map (a000040 . (+ 1) . a049084) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012, Oct 09 2011
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A003961 n) (apply * (map A000040 (map 1+ (map A049084 (factor n))))))
;; Antti Karttunen, May 20 2014

a:= n-> mul(nextprime(i[1])^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017

a[p_?PrimeQ] := a[p] = Prime[ PrimePi[p] + 1]; a[1] = 1; a[n_] := a[n] = Times @@ (a[#1]^#2& @@@ FactorInteger[n]); Table[a[n], {n, 1, 65}] (* Jean-François Alcover, Dec 01 2011, updated Sep 20 2019 *)
Table[Times @@ Map[#1^#2 & @@ # &, FactorInteger[n] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[n == 1], {n, 65}] (* Michael De Vlieger, Mar 24 2017 *)

a(n)=local(f); if(n<1,0,f=factor(n); prod(k=1,matsize(f)[1],nextprime(1+f[k,1])^f[k,2]))

a(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ Michel Marcus, May 17 2014

use ntheory ":all";  sub a003961 { vecprod(map { next_prime($) } factor(shift)); }  # _Dana Jacobsen, Mar 06 2016

from sympy import factorint, prime, primepi, prod
def a(n):
    f=factorint(n)
    return 1 if n==1 else prod(prime(primepi(i) + 1)**f[i] for i in f)
[a(n) for n in range(1, 11)] # Indranil Ghosh, May 13 2017

If n = Product p(k)^e(k) then a(n) = Product p(k+1)^e(k).
Multiplicative with a(p^e) = A000040(A000720(p)+1)^e. - David W. Wilson, Aug 01 2001
a(n) = Product_{k=1..A001221(n)} A000040(A049084(A027748(n,k))+1)^A124010(n,k). - Reinhard Zumkeller, Oct 09 2011 [Corrected by Peter Munn, Nov 11 2019]
A064989(a(n)) = n and a(A064989(n)) = A000265(n). - Antti Karttunen, May 20 2014 & Nov 01 2019
A001221(a(n)) = A001221(n) and A001222(a(n)) = A001222(n). - Michel Marcus, Jun 13 2014
From Peter Munn, Oct 31 2019: (Start)
a(n) = A225546((A225546(n))^2).
a(A225546(n)) = A225546(n^2).
(End)
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Product_{p prime} ((p^2-p)/(p^2-nextprime(p))) = 2.06399637... . - Amiram Eldar, Nov 18 2022

A319630 Positive numbers that are not divisible by two consecutive prime numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 13, 14, 16, 17, 19, 20, 21, 22, 23, 25, 26, 27, 28, 29, 31, 32, 33, 34, 37, 38, 39, 40, 41, 43, 44, 46, 47, 49, 50, 51, 52, 53, 55, 56, 57, 58, 59, 61, 62, 63, 64, 65, 67, 68, 69, 71, 73, 74, 76, 79, 80, 81, 82, 83, 85, 86, 87

Offset: 1

Rémy Sigrist, Sep 25 2018

nonn

This sequence is the complement of A104210.
Equivalently, this sequence corresponds to the positive numbers k such that:
- A300820(k) <= 1,
- A087207(k) is a Fibbinary number (A003714).
For any n > 0 and k >= 0, a(n)^k belongs to the sequence.
The numbers of terms not exceeding 10^k, for k=1,2,..., are 9, 78, 758, 7544, 75368, 753586, 7535728, 75356719, 753566574, ... Apparently, the asymptotic density of this sequence is 0.75356... - Amiram Eldar, Apr 10 2021
Numbers not divisible by any term of A006094. - Antti Karttunen, Jul 29 2022

			The number 10 is only divisible by 2 and 5, hence 10 appears in the sequence.
The number 42 is divisible by 2 and 3, hence 42 does not appear in the sequence.

Robert Israel, Table of n, a(n) for n = 1..10000

Cf. A003714, A006094, A087207, A104210, A300820, A356171 (odd terms only).

Positions of 1's in A322361 and in A356173 (characteristic function).

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({$1..N} minus R,list)); # Robert Israel, Apr 13 2020

q[n_] := SequenceCount[FactorInteger[n][[;; , 1]], {p1_, p2_} /; p2 == NextPrime[p1]] ==  0; Select[Range[100], q] (* Amiram Eldar, Apr 10 2021 *)

is(n) = my (f=factor(n)); for (i=1, #f~-1, if (nextprime(f[i,1]+1)==f[i+1,1], return (0))); return (1)

A300820(a(n)) <= 1.

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

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 15, 105, 3003, 3465, 13923, 45045, 264537, 459459, 745875, 1541475, 5221125, 8729721, 10790325, 14171625, 29288025, 34563375, 57034575, 71430975, 99201375, 109643625, 144729585, 205016175, 255835125, 295708875, 356080725, 399242025, 419159475, 449323875, 928602675, 939495375, 1083656925, 1941623775, 1962350685, 2083228875

Offset: 1

Antti Karttunen, Nov 10 2021

Numbers k such that A348990(k) [= k/gcd(k, A003961(k))] is equal to A348992(k), which is the odd part of A349162(k), thus all terms must be odd, as A348990 preserves the parity of its argument.
Equally, numbers k for which gcd(A064987(k), A191002(k)) is equal to A000265(gcd(A064987(k), A341529(k))).
Also odd numbers k for which A348993(k) = A319627(k).
Odd terms of A336702 are given by the intersection of this sequence and A349174.
Conjectures:
  (1) After 1, all terms are multiples of 3. (Why?)
  (2) After 1, all terms are in A104210, in other words, for all n > 1, gcd(a(n), A003961(a(n))) > 1. Note that if we encountered a term k with gcd(k, A003961(k)) = 1, then we would have discovered an odd multiperfect number.
  (3) Apart from 1, 15, 105, 3003, 13923, 264537, all other terms are abundant. [These apparently are also the only terms that are not Zumkeller, A083207. Note added Dec 05 2024]
  (4) After 1, all terms are in A248150. (Cf. also A386430).
  (5) After 1, all terms are in A348748.
  (6) Apart from 1, there are no common terms with A349753.
Note: If any of the last four conjectures could be proved, it would refute the existence of odd perfect numbers at once. Note that it seems that gcd(sigma(k), A003961(k)) < k, for all k except these four: 1, 2, 20, 160.
Questions:
  (1) For any term x here, can 2*x be in A349745? (Partial answer: at least x should be in A191218 and should not be a multiple of 3). Would this then imply that x is an odd perfect number? (Which could explain the points (1) and (4) in above, assuming the nonexistence of opn's).

Cf. A000203, A003961, A007949, A064989, A083207, A104210, A161942, A191002, A191218, A228058, A319627, A322361, A326134, A329963, A342671, A348748, A348990, A348992, A348993, A349162, A349164, A349174.
Cf. also A349745, A349746, A349753.
Subsequence of A349749, and of A386430.

Select[Range[10^6], #1/GCD[#1, #3] == #2/(2^IntegerExponent[#2, 2]*GCD[#2, #3]) & @@ {#, DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &] (* Michael De Vlieger, Nov 11 2021 *)

A000265(n) = (n >> valuation(n, 2));
A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
isA349169(n) = { my(s=sigma(n),u=A003961(n)); (n*gcd(s,u) == A000265(s)*gcd(n,u)); }; \\ (Program simplified Nov 30 2021)

For all n >= 1, A007949(A000203(a(n))) = A007949(a(n)). [sigma preserves the 3-adic valuation of the terms of this sequence] - Antti Karttunen, Nov 29 2021

Name changed and comment section rewritten by Antti Karttunen, Nov 29 2021

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 120, 216, 672, 2464, 22176, 228480, 523776, 640640, 837760, 5581440, 5765760, 7539840, 12999168, 19603584, 33860736, 38342304, 71344000, 95472000, 102136320, 197308800, 220093440, 345080736, 459818240, 807009280, 975576960, 1476304896, 1510831360, 1773584640

Offset: 1

Antti Karttunen, Nov 29 2021

nonn

Numbers k for which k * A342671(k) = A000203(k) * A322361(k).
Numbers k such that gcd(A064987(k), A191002(k)) = gcd(A064987(k), A341529(k)).
Obviously, all odd terms in this sequence must be squares.
All the terms k of A005820 that satisfy A007949(k) < A007814(k) [i.e., whose 3-adic valuation is strictly less than their 2-adic valuation] are also terms of this sequence. Incidentally, the first six known terms of A005820 satisfy this condition, while on the other hand, any hypothetical odd 3-perfect number would be excluded from this sequence. Also, as a corollary, any hypothetical 3-perfect numbers of the form 4u+2 must not be multiples of 3 if they are to appear here. Similarly for any k which occurs in A349169, for 2*k to occur in this sequence, it shouldn't be a multiple of 3 and k should also be a term of A191218. See question 2 and its partial answer in A349169.
From Antti Karttunen, Feb 13-20 2022: (Start)
Question: Are all terms/2 (A351548) abundant, from n > 1 onward?
Note that of the 65 known 5-multiperfect numbers, all others except these three 1245087725796543283200, 1940351499647188992000, 4010059765937523916800 are also included in this sequence. The three exceptions are distinguished by the fact that their 3 and 5-adic valuations are equal. In 62 others the former is larger.
If k satisfying the condition were of the form 4u+2, then it should be one of the terms of A191218 doubled as only then both k and sigma(k) are of the form 4u+2, with equal 2-adic valuations for both. More precisely, one of the terms of A351538.
(End)

Amiram Eldar, Table of n, a(n) for n = 1..52 (terms below 10^11)
Antti Karttunen, Ratio A342671(n)/A322361(n) plotted with OEIS Plot2 tool
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to sigma(n)

Cf. A000203, A003961, A005101, A005820, A007814, A007949, A046060, A064987, A191002, A191218, A248150, A322361, A341529, A342671, A347870, A351534, A351536, A351538.
Cf. also A349169, A349746, A351458, A351549 for other variants.
Subsequence of A351554 and also of its subsequence A351551.
Cf. A351459 (subsequence, intersection with A351458), A351548 (terms halved).

f1[p_, e_] := (p^(e + 1) - 1)/(p - 1); f2[p_, e_] := NextPrime[p]^e; q[1] = True; q[n_] := n * GCD[(s = Times @@ f1 @@@ (f = FactorInteger[n])), (r = Times @@ f2 @@@ f)] == s*GCD[n, r]; Select[Range[10^6], q] (* Amiram Eldar, Nov 29 2021 *)

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

For all n >= 1, A007814(A000203(a(n))) = A007814(a(n)). [sigma preserves the 2-adic valuation of the terms of this sequence]

A355442 a(n) = gcd(A003961(n), A276086(n)), where A003961 is fully multiplicative with a(p) = nextprime(p), and A276086 is primorial base exp-function.

Original entry on oeis.org

1, 3, 1, 9, 1, 5, 1, 3, 5, 3, 1, 5, 1, 3, 5, 9, 1, 25, 1, 3, 5, 3, 1, 5, 1, 3, 125, 9, 1, 7, 1, 3, 1, 3, 7, 5, 1, 3, 5, 63, 1, 5, 1, 3, 175, 3, 1, 5, 1, 21, 5, 9, 1, 125, 7, 3, 5, 3, 1, 7, 1, 3, 1, 9, 7, 5, 1, 3, 5, 21, 1, 25, 1, 3, 245, 9, 1, 5, 1, 21, 125, 3, 1, 5, 7, 3, 5, 9, 1, 7, 1, 3, 1, 3, 7, 5, 1, 3, 5, 441

Offset: 1

Antti Karttunen, Jul 13 2022

nonn

Antti Karttunen, Table of n, a(n) for n = 1..11550
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Index entries for sequences computed from indices in prime factorization
Index entries for sequences related to primorial base

Cf. A003961, A020639, A276086, A355001 [smallest prime factor of a(n)], A355456 [= gcd(sigma(n), a(n))], A355692 (Dirichlet inverse), A355820, A355821 (positions of 1's).

Cf. also A322361, A324198, A351459.

A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A355442(n) = gcd(A003961(n), A276086(n));

a(n) = gcd(A003961(n), A276086(n)).

A330749 a(n) = gcd(n, A064989(n)), where A064989 is fully multiplicative with a(2) = 1 and a(prime(k)) = prime(k-1) for odd primes.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 5, 4, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 2, 1, 1, 1, 5, 1, 4, 1, 1, 3, 1, 7, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 6, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 15

Offset: 1

Antti Karttunen, Dec 29 2019

nonn

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

Cf. A064989, A319626, A322361.

A330749(n) = {my(f); 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)); gcd(n,factorback(f)); };

a(n) = gcd(n, A064989(n)).
a(n) = n / A319626(n).
a(n) = A064989(A322361(n)).

A349174 Odd numbers k for which gcd(k, A003961(k)) is equal to gcd(sigma(k), A003961(k)), where A003961(n) is fully multiplicative with a(prime(k)) = prime(k+1), and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 5, 7, 9, 11, 13, 17, 19, 21, 23, 25, 29, 31, 33, 37, 39, 41, 43, 47, 49, 51, 53, 55, 59, 61, 63, 67, 69, 71, 73, 79, 81, 83, 85, 89, 91, 93, 95, 97, 101, 103, 107, 109, 111, 113, 115, 117, 119, 121, 123, 125, 127, 129, 131, 133, 135, 137, 139, 141, 145, 147, 149, 151, 153, 155, 157, 159, 161, 163, 167, 169

Offset: 1

Antti Karttunen, Nov 10 2021

nonn

Odd numbers k for which A322361(k) = A342671(k).
Odd numbers k for which A348994(k) = A349161(k).
Odd numbers k such that A319626(k) = A349164(k).
Odd terms of A336702 form a subsequence of this sequence. See also A349169.
Ratio of odd numbers residing in this sequence, vs. in A349175 seems to slowly decrease, but still apparently stays > 2 for a long time. E.g., for range 2 .. 2^28, it is 95302074/38915653 = 2.4489...

Cf. A000203, A003961, A319626, A322361, A336702, A342671, A348994, A349161, A349164, A349169.
Cf. A349175 (complement among the odd numbers).
Union of A349176 and A349177.

Select[Range[1, 169, 2], GCD[#1, #3] == GCD[#2, #3] & @@ {#, DivisorSigma[1, #], Times @@ Map[NextPrime[#1]^#2 & @@ # &, FactorInteger[#]]} &] (* 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); };
isA349174(n) = if(!(n%2),0,my(u=A003961(n)); gcd(u,sigma(n))==gcd(u,n));

cp's OEIS Frontend

A372572 Lexicographically earliest infinite sequence such that a(i) = a(j) => A009194(i) = A009194(j), A322361(i) = A322361(j) and A342671(i) = A342671(j), for all i, j >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A372570 Lexicographically earliest infinite sequence such that a(i) = a(j) => f(i) = f(j), where f(n) = [A009194(n), A009195(n), A009223(n), A322361(n), A342671(n)], for all i, j >= 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A003961 Completely multiplicative with a(prime(k)) = prime(k+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Perl

Python

Formula

A319630 Positive numbers that are not divisible by two consecutive prime numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A104210 Positive integers divisible by at least 2 consecutive primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A349169 Numbers k such that k * gcd(sigma(k), A003961(k)) is equal to the odd part of {sigma(k) * gcd(k, A003961(k))}, where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A349745 Numbers k for which k * gcd(sigma(k), A003961(k)) is equal to sigma(k) * gcd(k, A003961(k)), where A003961 shifts the prime factorization one step towards larger primes, and sigma is the sum of divisors function.

Views

Author

Keywords