A084967 -id:A084967 - cp's OEIS frontend

A328762 Numbers n for which 2 < A257993(A276086(A276086(n))) < A257993(n), where A276086 converts the primorial base expansion of n into its prime product form, and A257993 returns the index of the least prime not present in its argument.

210, 1470, 5250, 6510, 7140, 8400, 9450, 10710, 14490, 15750, 16380, 17640, 18690, 19950, 23730, 24990, 25620, 26880, 27930, 29190, 30030, 31290, 32340, 33600, 37380, 38640, 39270, 40530, 41580, 42840, 46620, 47880, 48510, 49770, 50820, 52080, 55860, 57120, 57750, 59010, 60270, 61530, 63420, 65730, 69510, 70770, 72660, 74970

Offset: 1

Antti Karttunen, Oct 27 2019

nonn

All terms are multiples of 5 (and thus of 30), because when applied to any number which is a multiple of 6, but not of 5 (and thus not a multiple of 30, implying that the primorial expansion ends with "x00", where x <> 0, and A257993(n) = 3), A276086 will yield a number of the form 30k+5 or 30k+25 (A084967) whose primorial expansion ends either as "...021" or as "...401" (with the least significant zero either in position 2 or 3), thus A328578(n) = A257993(A276086(A276086(n))) cannot simultaneously be larger than 2 and smaller than A257993(n).

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences related to primorial base

Setwise difference A328587 \ A328632.

Cf. A257993, A276086, A328578, A328589.

A257993(n) = { for(i=1,oo,if(n%prime(i),return(i))); }
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A328578(n) = A257993(A276086(A276086(n)));
isA328762(n) = { my(u=A328578(n)); ((u > 2) && (u < A257993(n))); };

A199593 Numbers n such that 3n-2, 3n-1 and 3n are all composite.

Original entry on oeis.org

9, 12, 17, 19, 22, 26, 29, 31, 32, 39, 40, 41, 42, 45, 48, 49, 52, 54, 57, 59, 62, 63, 68, 69, 70, 72, 73, 74, 79, 82, 83, 85, 87, 89, 92, 96, 97, 99, 100, 101, 102, 107, 108, 109, 110, 112, 114, 115, 119, 121, 122, 124, 126, 129, 131, 132, 135, 136, 138, 139, 142, 143, 146, 149, 151, 152, 157, 158, 159, 161, 162, 165, 166, 169, 171, 172, 173, 176, 177, 178

Offset: 1

N. J. A. Sloane, Nov 08 2011

From Antti Karttunen, Apr 17 2015: (Start)
Other, equivalent definitions:
Numbers n such that A007310(n) is composite, from which it follows that the function c(1) = 0, c(n) = 1-A075743(n-1) is the characteristic function of this sequence.
Numbers n such that A084967(n) has at least three prime factors (when counted with bigomega, A001222).
Numbers n such that A249823(n) is composite.
(End)
There are n - pi(3n) + 1 terms in this sequence up to n; with an efficient algorithm for pi(x) this allows isolated large values to be computed. Using David Baugh and Kim Walisch's calculation that pi(10^27) = 16352460426841680446427399 one can see that a(316980872906491652886905934) = 333333333333333333333333333 (since 999999999999999999999999997 is composite). - Charles R Greathouse IV, Sep 13 2016

Ernest V. Miliauskas, letter to N. J. A. Sloane, Dec 21 1985.

Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Bogart B. Strauss, Formula Explanation, pp. 1, 2, 3.

Cf. A053726, A199595, A001222, A007310, A075743, A084967, A249823.

[n: n in [1..200] | not IsPrime(3*n) and not IsPrime(3*n-1) and not IsPrime(3*n-2)]; // Vincenzo Librandi, Apr 18 2015

remove(t -> isprime(3*t-1 - (t mod 2)),{$2..2000}); # Robert Israel, Apr 17 2015

Select[Range[200], Union[PrimeQ[{3# - 2, 3# - 1, 3#}]] == {False} &] (* Alonso del Arte, Jul 06 2013 *)

is(n)=!isprime(bitor(3*n-2,1)) && n>1 \\ Charles R Greathouse IV, Oct 27 2013
(Scheme, after Greathouse's PARI-program above, requiring also Antti Karttunen's IntSeq-library)
(define A199593 (MATCHING-POS 1 2 (lambda (n) (not (prime? (A003986bi (+ n n n -2) 1)))))) ;; A003986bi implements binary inclusive or (A003986).
;; Antti Karttunen, Apr 17 2015

from sympy import isprime
def ok(n): return n > 0 and not any(isprime(3*n-i) for i in [2, 1, 0])
print([k for k in range(179) if ok(k)]) # Michael S. Branicky, Apr 16 2022

((1+(-1)^k)((-1)^n)(2n+3)+2k(6n+9+(-1)^n)+((-1)^k)+(12n^2)+36n+29)/4 n,k are all natural numbers and zero. - Bogart B. Strauss, Jul 10 2013

a(n) = n + 3n/log n + o(n/log n). - Charles R Greathouse IV, Oct 27 2013, corrected Aug 07 2016

A249827 Row 3 of A246278: replace in 2n each prime factor p(k) with prime p(k+2).

Original entry on oeis.org

5, 25, 35, 125, 55, 175, 65, 625, 245, 275, 85, 875, 95, 325, 385, 3125, 115, 1225, 145, 1375, 455, 425, 155, 4375, 605, 475, 1715, 1625, 185, 1925, 205, 15625, 595, 575, 715, 6125, 215, 725, 665, 6875, 235, 2275, 265, 2125, 2695, 775, 295, 21875, 845, 3025, 805, 2375, 305, 8575, 935, 8125, 1015, 925, 335, 9625, 355, 1025, 3185, 78125, 1045, 2975, 365

Offset: 1

Antti Karttunen, Nov 23 2014

nonn

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

Row 3 of A246278.

Cf. A003961, A016945, A048673, A084967, A249824, A249734.

(define (A249827 n) (A003961 (A249734 n)))
(define (A249827 n) (A246278bi 3 n)) ;; Alternative implementation. Code for A246278bi given in A246278.

a(n) = A003961(A249734(n)).
a(n) = A003961(A016945(A048673(n)-1)).
a(n) = A084967(A249824(n)). [Permutation of A084967.]

A251758 Let n>=2 be a positive integer with divisors 1 = d_1 < d_2 < ... < d_k = n, and s = d_1d_2 + d_2d_3 + ... + d_(k-1)*d_k. The sequence lists the values a(n) = floor(n^2/s).

Original entry on oeis.org

2, 3, 1, 5, 1, 7, 1, 2, 1, 11, 1, 13, 1, 2, 1, 17, 1, 19, 1, 2, 1, 23, 1, 4, 1, 2, 1, 29, 1, 31, 1, 2, 1, 4, 1, 37, 1, 2, 1, 41, 1, 43, 1, 2, 1, 47, 1, 6, 1, 2, 1, 53, 1, 4, 1, 2, 1, 59, 1, 61, 1, 2, 1, 4, 1, 67, 1, 2, 1, 71, 1, 73, 1, 2, 1, 6, 1, 79, 1, 2, 1

Offset: 2

Michel Lagneau, Dec 08 2014

nonn

s is always less than n^2 and if n is a prime number then s divides n^2.
For n >= 2, the sequence has the following properties:
a(n) = n if n is prime.
a(n) = 1 if n is in A005843 and > 2;
a(n) <= 2 if n is in A016945 and > 3;
a(n) <= 4 if n is in A084967 and > 5;
a(n) <= 6 if n is in A084968 and > 7;
a(n) = 8: <= 35336848261, ...;
a(n) <= 10 if n is in A084969 and > 11;
a(n) <= 12 if n is in A084970 and > 13;
a(n) = 14: 6678671, ...;
This is different from A250480 (a(n) = n for all prime n, and a(n) = A020639(n) - 1 for all composite n), which thus satisfies the above conditions exactly, while with this sequence A020639(n)-1 gives only the guaranteed upper limit for a(n) at composite n. Note that the first different term does not occur until at n = 2431 = 11*13*17, for which a(n) = 9. (See the example below.)
Conjecture: Terms x, where a(x)=n, x=p#k/p#j, p#i is the i-th primorial, k>j is suitable large k and j is the number of primes less than n. As an example, n=9, x = p#7/p#4 = 2431. For n=10, x = p#6/p#4 = 143 although 121 = 11^2 is the least x where a(x)=10 (see formula section). For n=8, x = p#12/p#4, p#13/p#4, p#14/p#4, p#15/p#4, p#16/p#4, etc. But is p#12/p#4 the least such x? - Robert G. Wilson v, Dec 18 2014
n^2/s is only an integer iff n is prime. - Robert G. Wilson v, Dec 18 2014
First occurrence of n >= 1: 4, 2, 3, 25, 5, 49, 7, ??? <= 35336848261, 2431, 121, 11, 169, 13, 6678671, 7429, 289, 17, 361, 19, 31367009, 20677, 529, 23, ..., . - Robert G. Wilson v, Dec 18 2014

			For n = 2431 = 11*13*17, we have (as the eight divisors of 2431 are [1, 11, 13, 17, 143, 187, 221, 2431]) a(n) = floor((2431*2431) / ((1*11)+(11*13)+(13*17)+(17*143)+(143*187)+(187*221)+(221*2431))) = floor(5909761/608125) = floor(9.718) = 9.

Michel Lagneau, Table of n, a(n) for n = 2..10000
International Mathematical Olympiad, Problems, IMO-2002, Problem 4.

Cf. A000040 (prime numbers), A005843 (even numbers), A016945 (6n+3), A084967 (GCD( 5k, 6) =1), A084968 (GCD( 7k, 30) =1), A084969 (GCD( 11k, 30) =1), A084970 (Numbers whose smallest prime factor is 13).
Cf. also A020639 (the smallest prime divisor), A055396 (its index) and arrays A083140 and A083221 (Sieve of Eratosthenes).
Differs from A250480 for the first time at n = 2431, where a(2431) = 9, while A250480(2431) = 10.
Cf. A078730 (sum of products of two successive divisors of n).

with(numtheory):nn:=100:
for n from 2 to nn do:
   x:=divisors(n):n0:=nops(x):s:=sum('x[i]*x[i+1]','i'=1..n0-1):
   z:=floor(n^2/s):printf(`%d, `,z):
od:

f[n_] := Floor[ n^2/Plus @@ Times @@@ Partition[ Divisors@ n, 2, 1]]; Array[f, 81, 2] (* Robert G. Wilson v, Dec 18 2014 *)

a(n) <= A250480(n), and especially, for all composite n, a(n) < A020639(n). [Cf. the Comments section above.] - Antti Karttunen, Dec 09 2014
From Robert G. Wilson v, Dec 18 2014: (Start)
a(n) = floor(n^2/A078730(n));
a(n) = n iff n is prime. (End)

Comments section edited by Antti Karttunen, Dec 09 2014

Instances of n for which a(n) = 8 and 14 found by Robert G. Wilson v, Dec 18 2014

A376839 a(1) = 1. For n > 1 if A007947(a(n-1)) is in A002110, a(n) is the smallest prime not already a term. Otherwise a(n) is the least novel multiple of the smallest non divisor prime of a(n-1).

Original entry on oeis.org

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

Offset: 1

David James Sycamore, Oct 06 2024

nonn

A non divisor prime of a(n-1) is any prime p < Gpf(a(n-1)) which does not divide a(n-1). A007947(a(n-1)) is in A002110 iff a(n-1) is a term in A055932. Sequence is conjectured to be a permutation of the natural numbers (A000027) with primes in order.

Scatterplot shows trajectories of numbers whose smallest prime factor is prime p, e.g., for p = 5, numbers in A084967, p = 7, those in A084968, p = 11 those in A084969, etc. - Michael De Vlieger, Oct 09 2024

			a(1) = 1 = A002110(0) so a(2) = 2 (smallest prime not already a term).
a(2) = 2 = A002110(1) so a(3) = 3.
a(3) = 3 not a term in A002110 so a(4) is least novel multiple of 2, the least non divisor prime of 3. Therefore a(4) = 4 since 2 has occurred earlier.
a(39) = 42, not a term in A002110 so a(40) = 25, the least novel multiple of 5, the smallest non divisor prime of 42.

Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^20.
Michael De Vlieger, Log log scatterplot of a(n), n = 1..2^14, showing primes in red, perfect prime powers in gold, squarefree composites in green, and numbers neither squarefree nor prime powers in blue or purple, with the latter additionally representing powerful numbers that are not prime powers.

Cf. A000040, A000027, A002110, A007947, A055932, A084967, A084968, A084969, A372368, A376865.

nn = 120; c[] := False; m[] := 1; f[x_] := FactorInteger[x][[All, 1]];
  Array[Set[{a[#], c[#], m[#]}, {#, True, 2}] &, 2]; j = 2; v = 3;
  Do[If[Or[IntegerQ@ Log2[j], And[EvenQ[j], Union@ Differences@ PrimePi[#] == {1}]],
     k = v; While[c[k*m[k]], m[k]++]; k *= m[k],
     k = 2; While[Divisible[j, k], k = NextPrime[k]];
     While[c[k*m[k]], m[k]++]; k *= m[k]] &[f[j]];
  Set[{a[n], c[k], j}, {k, True, k}];
  If[k == v, While[c[v], v = NextPrime[v]]], {n, 3, nn}];
Array[a, nn] (* Michael De Vlieger, Oct 09 2024 *)

More terms from Michael De Vlieger, Oct 09 2024

A120320 RF(5): refactorable numbers with smallest prime factor 5.

Original entry on oeis.org

625, 1500625, 9150625, 17850625, 37515625, 52200625, 73530625, 81450625, 174900625, 442050625, 577200625, 1171350625, 1766100625, 1838265625, 2136750625, 3049800625, 4931550625, 7573350625, 8653650625, 12594450625, 15882300625, 17748900625, 21970650625, 24343800625

Offset: 1

Walter Kehowski, Jun 20 2006

nonn

Numbers that are odd squares, 5 is their smallest prime factor, and are refactorable.
See A033950 for references. For any prime p, p^(p-1) is the smallest element of RF(p), the refactorable numbers whose smallest prime factor is p. Thus 5^(5-1) = 625 is the first element. Other elements would also be 5^4*17^4 or 5^16*17^4.
All the terms are of the form 5^2 * A084967(k)^2 = 5^4 * A007310(k)^2. - Amiram Eldar, Aug 01 2024

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

Intersection of A033950 and A084967.

Cf. A007310, A036896, A036897.

with(numtheory); RF5:=[]: p:=5: for w to 1 do for j from 1 to 12^5 do k:=2*j+1; if k mod 3 <> 0 and k mod p = 0 then n:=k^2; t:=tau(n); if (n mod t = 0) then RF5:=[op(RF5),n]; print(ifactor(n)); fi fi; od od;

lista(kmax) = {my(m); for(k = 1, kmax, m = 25*(k\2*6-(-1)^k)^2; if(!(m % numdiv(m)), print1(m, ", ")));} \\ Amiram Eldar, Aug 01 2024

a(37)-a(40) from Amiram Eldar, Aug 01 2024

A138045 Triangle read by rows: largest proper divisor of n as a table, ones excluded.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 0

Offset: 1

Mats Granvik, Mar 02 2008

The numbers in the triangle form lines that begin at T(A001248,A000040). The first line of numbers from the right, is T(A005843,A000027). The second line is T(A016945,A005408). The third line is T(A084967,A007310).

			The first few terms of the table are:
  0
  0,0
  0,0,0
  0,2,0,0
  0,0,0,0,0
  0,0,3,0,0,0
  0,0,0,0,0,0,0
  0,0,0,4,0,0,0,0
  0,0,3,0,0,0,0,0,0

Antti Karttunen, Table of n, a(n) for n = 1..23220 (the first 215 rows of the triangle).
Eric Weisstein's World of Mathematics, Proper Divisor.

Cf. A001248, A016945, A084967, A007310, A127093, A032742, A007775.

up_to = 23220; \\ binomial(215+1,2)
A032742(n) = if(1==n,n,n/vecmin(factor(n)[,1]));
A138045tr(n, k) = if((k>1) && (A032742(n)==k), k, 0);
A138045list(up_to) = { my(v = vector(up_to), i=0); for(n=1,oo, for(k=1,n, i++; if(i > up_to, return(v)); v[i] = A138045tr(n,k))); (v); };
v138045 = A138045list(up_to);
A138045(n) = v138045[n]; \\ Antti Karttunen, Dec 24 2018

T(n,k) = if k==A032742(n) and n(T(n,k))==n(A032742(n)) and k>1 then k else 0 (1<=k<=n), T(1,1)=0.

cp's OEIS Frontend

A328762 Numbers n for which 2 < A257993(A276086(A276086(n))) < A257993(n), where A276086 converts the primorial base expansion of n into its prime product form, and A257993 returns the index of the least prime not present in its argument.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A199593 Numbers n such that 3n-2, 3n-1 and 3n are all composite.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

Formula

A249827 Row 3 of A246278: replace in 2n each prime factor p(k) with prime p(k+2).

Views

Author

Keywords

Links

Crossrefs

Programs

Scheme

Formula

A251758 Let n>=2 be a positive integer with divisors 1 = d_1 < d_2 < ... < d_k = n, and s = d_1*d_2 + d_2*d_3 + ... + d_(k-1)*d_k. The sequence lists the values a(n) = floor(n^2/s).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A376839 a(1) = 1. For n > 1 if A007947(a(n-1)) is in A002110, a(n) is the smallest prime not already a term. Otherwise a(n) is the least novel multiple of the smallest non divisor prime of a(n-1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Extensions

A120320 RF(5): refactorable numbers with smallest prime factor 5.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

PARI

Extensions

A138045 Triangle read by rows: largest proper divisor of n as a table, ones excluded.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A251758 Let n>=2 be a positive integer with divisors 1 = d_1 < d_2 < ... < d_k = n, and s = d_1d_2 + d_2d_3 + ... + d_(k-1)*d_k. The sequence lists the values a(n) = floor(n^2/s).