The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A249909 Smallest prime factor of A241601(n), or 1 if A241601(n) = 1. 3

%I

%S 1,1,1,1,1,1,1,61,1,277,1,19,691,43,1,47,3617,228135437,43867,79,283,

%T 41737,131,31,103,2137,657931,67,9349,71,1721,15669721,37,930157,

%U 151628697551,4153,26315271553053477373,9257,154210205991661,23489580527043108252017828576198947741,137616929,763601

%N Smallest prime factor of A241601(n), or 1 if A241601(n) = 1.

%C Also the smallest prime factor of A246006(n) that is >= n+2.

%C a(n) = A020639(A241601(n)).

%C a(n) = 1 iff n is in the set {0, 1, 2, 3, 4, 5, 6, 8, 10, 14}.

%C a(189) is currently unknown; a(190)..a(199) = {5101, 559570609330768709, 40833790860803270336710504624737304862569304959957, 311, 467, 34110029, 461, 26034939865747697437451558982836040663625026070193, 34470847, 1879}.

%C All terms are Bernoulli or Euler irregular primes.

%H Eric Chen, <a href="/A249909/b249909.txt">Table of n, a(n) for n = 0..188</a>

%H Samuel Wagstaff, <a href="http://homes.cerias.purdue.edu/~ssw/bernoulli/index.html">Factorization of Bernoulli and Euler numbers</a>

%t a246006[n_] := If[EvenQ[n], Abs[Numerator[BernoulliB[n]]], Abs[EulerE[n-1]]];

%t a241601[n_] := a246006[n]/GCD[a246006[n], n!]

%t a = {}; Do[a = Append[a, FactorInteger[a241601[n]][[1, 1]]], {n, 0, 99} ]; a

%K nonn

%O 0,8

%A _Eric Chen_, Dec 15 2014

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

Last modified May 26 04:26 EDT 2022. Contains 354074 sequences. (Running on oeis4.)