This site is supported by donations to The OEIS Foundation. Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A005171 0 if n is prime else 1. 33
 1, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 COMMENTS Number of orbits of length n in map whose periodic points are A023890. - Thomas Ward (t.ward(AT)uea.ac.uk) Characteristic function of nonprimes A018252. - Jonathan Vos Post, Dec 30 2007 Triangle A157423 = A005171 in every column. A052284 = INVERT transform of A005171, and the eigensequence of triangle A157423. - Gary W. Adamson, Feb 28 2009 REFERENCES Douglas Hofstadter, Fluid Concepts and Creative Analogies: Computer Models of the Fundamental Mechanisms of Thought. LINKS Antti Karttunen, Table of n, a(n) for n = 1..65537 R. K. Guy, The strong law of small numbers. Amer. Math. Monthly 95 (1988), no. 8, 697-712. [Annotated scanned copy] Y. Puri and T. Ward, Arithmetic and growth of periodic orbits, J. Integer Seqs., Vol. 4 (2001), #01.2.1. Yash Puri and Thomas Ward, A dynamical property unique to the Lucas sequence, Fibonacci Quarterly, Volume 39, Number 5 (November 2001), pp. 398-402. FORMULA If b(n) is the n-th term of A023890, then a(n)=(1/n)* Sum_{ d divides n } mu(d)*a(n/d). E.g., a(6) = 1 since the 6th term of A023890 is 7 and the first term is 1. a(n) = 1-[(n-1)!^2 mod n], with n>=1. - Paolo P. Lava, Jun 11 2007 a(n) = NOT(A010051(n)) = 1 - A010051(n). - Jonathan Vos Post, Dec 30 2007 a(n) equals the first column in a table T defined by the recurrence: If n = k then T(n,k) = 1 else if k = 1 then T(n,k) = 1 - Product_{k divides n} of T(n,k), else if k divides n then T(n,k) = T(n/k,1). This is true since T(n,k) = 0 when k divides n and n/k is prime which results in Product_{k divides n} = 0 for the composite numbers and where k ranges from 2 to n. Therefore there is a remaining 1 in the expression 1-Product_{k divides n}, in the first column. Provided below is a Mathematica program as an illustration. - Mats Granvik, Sep 21 2013 a(n) = A057427(A239968(n)). - Reinhard Zumkeller, Mar 30 2014 MAPLE A005171 := proc(n)     if isprime(n) then         0 ;     else         1 ;     end if; end proc: # R. J. Mathar, May 26 2017 MATHEMATICA f[n_] := If[PrimeQ@ n, 0, 1]; Array[f, 105] (* Robert G. Wilson v, Jun 20 2011 *) Clear[t, n, k, i, nn]; nn = 105; t[n_, k_] :=  t[n, k] = If[n == k, 1, If[k == 1, 1 - Product[t[n, k + i], {i, 1, n - 1}], If[Mod[n, k] == 0, t[n/k, 1], 1], 1]]; Table[t[n, 1], {n, 1, nn}] (* Mats Granvik, Sep 21 2013 *) PROG (PARI) a(n)=if(n<1, 0, !isprime(n)) /* Michael Somos, Jun 08 2005 */ (Haskell) a005171 = (1 -) . a010051  -- Reinhard Zumkeller, Mar 30 2014 CROSSREFS Cf. A010051, A018252, A023890. Cf. A157423, A157424, A052284, A050374. Sequence in context: A242252 A100810 A174889 * A283265 A181406 A285252 Adjacent sequences:  A005168 A005169 A005170 * A005172 A005173 A005174 KEYWORD nonn,easy AUTHOR 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 March 24 00:12 EDT 2019. Contains 321444 sequences. (Running on oeis4.)