cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A109616 Numbers k > 0 such that (10's complement factorial of k) + 1 is prime.

This page as a plain text file.
%I A109616 #26 Aug 13 2025 19:28:01
%S A109616 2,5,10,25,36,44,65,67,138,149,176,212,279,1293,2367,2463,2707,3130,
%T A109616 4150,4635,6070,6355,10111,10560
%N A109616 Numbers k > 0 such that (10's complement factorial of k) + 1 is prime.
%C A109616 Larger values not certified.
%C A109616 Some of the larger entries may only correspond to probable primes.
%C A109616 a(25) > 16000. - _Michael S. Branicky_, Apr 27 2025
%e A109616 5 is a term because (10-5)*(10-4)*(10-3)*(10-2)*(10-1) + 1 = 15121 is prime.
%e A109616 10 is a term because (100-10)*(10-9)*..*(10-1) + 1 = 32659201 is prime.
%p A109616 b:= proc(n) option remember; `if`(n=0, 1,
%p A109616       (10^length(n)-n)*b(n-1))
%p A109616     end:
%p A109616 q:= n-> isprime(1+b(n)):
%p A109616 select(q, [$1..300])[];  # _Alois P. Heinz_, Aug 13 2025
%t A109616 f[n_] := 10^Length[IntegerDigits[n]] - n; p = 1; Do[p *= f[n]; If[PrimeQ[p + 1], Print[n]], {n, 4635}] (* _Ryan Propper_, May 20 2006 *)
%Y A109616 Cf. A110396, A109617.
%K A109616 nonn,base,hard,more
%O A109616 1,1
%A A109616 _Jason Earls_, Aug 01 2005
%E A109616 a(14)-a(20) from _Ryan Propper_, May 20 2006
%E A109616 a(21)-a(24) from _Michael S. Branicky_, Apr 25 2025