login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
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). 14
1, 2, 3, 25, 26, 33, 93, 1034, 970225, 8550146, 325422273, 414690595, 1864797542, 2438037206 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

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.

LINKS

Table of n, a(n) for n=1..14.

EXAMPLE

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.

MAPLE

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;

MATHEMATICA

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 *)

CROSSREFS

Fixed points of permutation pair A048673/A064216.

Subsequence of A245449.

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.

Sequence in context: A326226 A291262 A307922 * A320223 A296275 A295328

Adjacent sequences:  A048671 A048672 A048673 * A048675 A048676 A048677

KEYWORD

nonn,hard,more

AUTHOR

Antti Karttunen, Jul 14 1999

EXTENSIONS

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

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 21 14:18 EDT 2019. Contains 328301 sequences. (Running on oeis4.)