A351227 - OEIS

%I #8 Feb 05 2022 17:20:07

%S 0,1,2,3,4,5,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,

%T 27,28,29,33,34,35,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,

%U 54,55,56,57,58,59,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87

%N Numbers k for which A276086(k) > k, where A276086 is the primorial base exp-function.

%H <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>

%t Select[Range[0, 87], Block[{i, m, n = #, p}, m = i = 1; While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m > #] &] (* _Michael De Vlieger_, Feb 05 2022 *)

%o (PARI)

%o A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };

%o isA351227(n) = (A276086(n)>n);

%Y Cf. A276086, A351226 (complement), A351229 (subsequence).

%Y Indices of positive terms in A351225.

%K nonn,base

%O 1,3

%A _Antti Karttunen_, Feb 05 2022