The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 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 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.

Last modified April 5 06:13 EDT 2020. Contains 333238 sequences. (Running on oeis4.)