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

 

Logo

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A019565 If n = Sum 2^e_i, e_i distinct, then a(n) = Product prime_{e_i+1}. 26
1, 2, 3, 6, 5, 10, 15, 30, 7, 14, 21, 42, 35, 70, 105, 210, 11, 22, 33, 66, 55, 110, 165, 330, 77, 154, 231, 462, 385, 770, 1155, 2310, 13, 26, 39, 78, 65, 130, 195, 390, 91, 182, 273, 546, 455, 910, 1365, 2730, 143, 286, 429, 858, 715, 1430, 2145, 4290 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

A permutation of the squarefree numbers A005117.  The missing positive numbers are in A013929. - Alois P. Heinz, Sep 06 2014

LINKS

R. Zumkeller, Table of n, a(n) for n = 0..8191

FORMULA

G.f. prod(k>=0, 1 + prime(k+1)*x^2^k), where prime(k)=A000040(k). - Ralf Stephan, Jun 20 2003

a(n) = f(n, 1, 1) with f(x, y, z) = if x > 0 then f(floor(x/2), y*prime(z)^(x mod 2), z+1) else y. - Reinhard Zumkeller, Mar 13 2010

EXAMPLE

5 = 2^2+2^0, e_1 = 2, e_2 = 0, prime(2+1) = prime(3) = 5, prime(0+1) = prime(1) = 2, so a(5) = 5*2 = 10.

MAPLE

a:= proc(n) local i, m, r; m:=n; r:=1;

      for i while m>0 do if irem(m, 2, 'm')=1

        then r:=r*ithprime(i) fi od; r

    end:

seq(a(n), n=0..60);  # Alois P. Heinz, Sep 06 2014

MATHEMATICA

Do[m=1; o=1; k1=k; While[ k1>0, k2=Mod[k1, 2]; If[k2\[Equal]1, m=m*Prime[o]]; k1=(k1-k2)/ 2; o=o+1]; Print[m], {k, 0, 55}] (Lei Zhou, Feb 15 2005)

PROG

(PARI) a(n)=factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))  \\ M. F. Hasler, Mar 26 2011, updated Aug 22 2014

(Haskell)

a019565 n = product $ zipWith (^) a000040_list (a030308_row n)

-- Reinhard Zumkeller, Apr 27 2013

CROSSREFS

Cf. A101278, A054842, A007088, A030308, A000040.

Cf. A110765.

Sequence in context: A073740 A239956 A077320 * A133477 A039653 A106379

Adjacent sequences:  A019562 A019563 A019564 * A019566 A019567 A019568

KEYWORD

nonn

AUTHOR

Marc LeBrun

EXTENSIONS

Definition corrected by Klaus-R. Löffler, Aug 20 2014.

STATUS

approved

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

Content is available under The OEIS End-User License Agreement .

Last modified December 21 20:27 EST 2014. Contains 252326 sequences.