A246371 -id:A246371 - cp's OEIS frontend

A246344 a(0) = 16, after which, if a(n-1) = product_{k >= 1} (p_k)^(c_k), then a(n) = (1/2) * (1 + product_{k >= 1} (p_{k+1})^(c_k)), where p_k indicates the k-th prime, A000040(k).

16, 41, 22, 20, 32, 122, 101, 52, 77, 72, 338, 434, 611, 451, 280, 1040, 4820, 7907, 3960, 30713, 15364, 22577, 12154, 9791, 4902, 8108, 9131, 5815, 4099, 2056, 3551, 2095, 1474, 1385, 984, 2903, 1455, 1768, 4361, 5869, 2940, 19058, 18845, 13227, 11053, 8707, 4357, 2182, 1640, 4064, 15917, 9432, 46238

Offset: 0

Antti Karttunen, Aug 24 2014

nonn

Iterates of A048673 starting from value 16.

Either this sequence is actually part of the cycle containing 12 (see A246342) or 16 is the smallest member of this cycle (regardless of whether this cycle is finite or infinite), which follows because all numbers 1 .. 11 are in finite cycles, and also 13 and 14 are in closed cycles and 15 is in the cycle of 12.

			Start with a(0) = 16; then after each new term is obtained by replacing each prime factor of the previous term with the next prime, to whose product is added one before it is halved:
16 = 2^4 = p_1^4 -> ((p_2^4)+1)/2 = (3^4 + 1)/2 = (81+1)/2 = 41, thus a(1) = 41.
41 = p_13 -> ((p_14)+1)/2 = (43+1)/2 = 22, thus a(2) = 22.

Antti Karttunen, Table of n, a(n) for n = 0..1001

A246345 gives the terms of the same cycle when going to the opposite direction from 16.
Cf. A048673, A064216, A246342, A246343.
Cf. also A246281, A246282, A246351, A246352, A246361, A246362, A246371, A246372.

default(primelimit, 2^30);
A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i,1] = nextprime(f[i,1]+1)); factorback(f); \\ Using code of Michel Marcus
A048673(n) = (A003961(n)+1)/2;
k = 16; for(n=0, 1001, write("b246344.txt", n, " ", k) ; k = A048673(k));
(Scheme, with memoization-macro definec)
(definec (A246344 n) (if (zero? n) 16 (A048673 (A246344 (- n 1)))))

a(0) = 16, and for n >= 1, a(n) = A048673(a(n-1)).

A246345 a(0) = 16, after which, if (2*a(n-1)) - 1 = product_{k >= 1} (p_k)^(c_k) then a(n) = product_{k >= 1} (p_{k-1})^(c_k), where p_k indicates the k-th prime, A000040(k).

Original entry on oeis.org

16, 29, 34, 61, 49, 89, 106, 199, 389, 310, 617, 524, 694, 1207, 1921, 3097, 3899, 4142, 3374, 3674, 4234, 8461, 16903, 20211, 37841, 22408, 26853, 26391, 48031, 68605, 137201, 81272, 108334, 137809, 266737, 512627, 347932, 497005, 982081, 1942279, 3855031, 5292209

Offset: 0

Antti Karttunen, Aug 24 2014

nonn

Iterates of A064216 starting from value 16.

See also the comments in A246344.

			Start with a(0) = 16; then after each new term is obtained by doubling the previous term, from which one is subtracted, after which each prime factor is replaced with the previous prime:
16 -> ((2*16)-1) = 31 = p_1, and p_10 = 29, thus a(1) = 29.
29 -> ((2*29)-1) = 57 = 3*19 = p_2 * p_8, and p_1 * p_7 = 2*17 = 34, thus a(2) = 34.

Antti Karttunen, Table of n, a(n) for n = 0..1001

A246344 gives the terms of the same cycle when going to the opposite direction from 16.
Cf. A048673, A064216, A246342, A246343.
Cf. also A246281, A246282, A246351, A246352, A246361, A246362, A246371, A246372.

nxt[n_]:=Times@@(NextPrime[#,-1]&/@(Flatten[Table[#[[1]],{#[[2]]}]&/@ FactorInteger[2 n-1]])); NestList[nxt,16,50] (* Harvey P. Dale, Apr 04 2015 *)

default(primelimit, 2^30);
A064989(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)); factorback(f)};
A064216(n) = A064989((2*n)-1);
k = 16; for(n=0, 1001, write("b246345.txt", n, " ", k); k = A064216(k));
(Scheme, with memoization-macro definec)
(definec (A246345 n) (if (zero? n) 16 (A064216 (A246345 (- n 1)))))

a(0) = 16, a(n) = A064216(a(n-1)).

A246374 Primes p such that if 2p-1 = product_{k >= 1} A000040(k)^(c_k), then p > product_{k >= 1} A000040(k-1)^(c_k).

Original entry on oeis.org

5, 11, 13, 17, 23, 41, 43, 53, 59, 61, 73, 83, 113, 131, 137, 149, 163, 167, 173, 179, 193, 233, 239, 257, 263, 281, 293, 311, 313, 347, 353, 383, 389, 401, 419, 431, 443, 449, 463, 479, 491, 503, 509, 523, 557, 563, 587, 593, 599, 613, 617, 641, 653, 677, 683, 743, 761, 773, 787, 797

Offset: 1

Antti Karttunen, Aug 25 2014

nonn

Primes p such that A064216(p) < p, or equally, A064989(2p-1) < p.

For all primes p here, 2p-1 must be composite (a necessary but not sufficient condition).

			5 is present, as 2*5 - 1 = 9 = p_2 * p_2, and p_1 * p_1 = 4, and 5 > 4.

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

Intersection of A000040 and A246371.
A246373 gives the primes not here.
Cf. A000040, A064216, A064989.

A064989(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)); factorback(f)};
n = 0; forprime(p=2,2^31, if((A064989((2*p)-1) < p), n++; write("b246374.txt", n, " ", p); if(n > 9999, break)));
(Scheme, with Antti Karttunen's IntSeq-library)
(define A246374 (MATCHING-POS 1 1 (lambda (n) (and (prime? n) (< (A064216 n) n)))))

cp's OEIS Frontend

A246344 a(0) = 16, after which, if a(n-1) = product_{k >= 1} (p_k)^(c_k), then a(n) = (1/2) * (1 + product_{k >= 1} (p_{k+1})^(c_k)), where p_k indicates the k-th prime, A000040(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A246345 a(0) = 16, after which, if (2*a(n-1)) - 1 = product_{k >= 1} (p_k)^(c_k) then a(n) = product_{k >= 1} (p_{k-1})^(c_k), where p_k indicates the k-th prime, A000040(k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A246374 Primes p such that if 2p-1 = product_{k >= 1} A000040(k)^(c_k), then p > product_{k >= 1} A000040(k-1)^(c_k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI