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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A099104 Characteristic function of badly sieved numbers (A066680). 6

%I

%S 0,1,1,0,1,0,1,1,0,0,1,1,1,0,0,0,1,1,1,0,0,0,1,0,0,0,1,0,1,1,1,0,0,0,

%T 0,0,1,0,0,0,1,0,1,0,1,0,1,0,0,1,0,0,1,0,0,0,0,0,1,0,1,0,1,0,0,0,1,0,

%U 0,1,1,0,1,0,1,0,0,0,1,1,0,0,1,0,0,0,0,0,1,0,0,0,0,0,0,0,1,1,0,0,1,0,1,0,1

%N Characteristic function of badly sieved numbers (A066680).

%C a(A066680(n)) = 1;

%C A066683(n)=Sum(a(k):1<=k<=n); A066682(n)=a(n)*A066683(n).

%H Antti Karttunen, <a href="/A099104/b099104.txt">Table of n, a(n) for n = 1..65537</a>

%H <a href="/index/Ch#char_fns">Index entries for characteristic functions</a>

%F a(n) = Prod(if n>d^2 then 1 else 1-a(d): d<n and n mod d = 0); a(1)=0.

%t a[1] = 0; a[n_] := a[n] = Product[ If[ n > d^2, 1, 1-a[d]], {d, Select[ Range[n-1], Mod[n, #] == 0 &] } ]; Table[ a[n], {n, 1, 105}](* _Jean-Fran├žois Alcover_, Feb 15 2012, after formula *)

%o (PARI) A099104(n) = if(1==n,0,my(m=1); fordiv(n,d,if((d<n)&&((d*d)>=n),m *= (1-A099104(d)))); (m)); \\ _Antti Karttunen_, Jul 02 2018

%Y Cf. A066680, A066682, A066683, A099042.

%K nonn

%O 1,1

%A _Reinhard Zumkeller_, Sep 26 2004

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 March 22 17:25 EDT 2019. Contains 321422 sequences. (Running on oeis4.)