editing
proposed
editing
proposed
2, 4, 17, 130, 1283, 6889, 40037, 638521, 10126943, 186814849, 2092495862
a(12) <= 8167862431, a(13) <= 1052676193433, a(14) <= 30964627320559. - Daniel Suteu, Jan 07 2025
a(n)*2^n - 1 >= A360438(n). - Daniel Suteu, Jan 07 2025
a(11) from Daniel Suteu, Jan 07 2025
approved
editing
editing
proposed
(PARI)
almost_primes(A, B, n) = A=max(A, 2^n); (f(m, p, n) = my(list=List()); if(n==1, forprime(q=max(p, ceil(A/m)), B\m, listput(list, m*q)), forprime(q=p, sqrtnint(B\m, n), list=concat(list, f(m*q, q, n-1)))); list); vecsort(Vec(f(1, 2, n)));
a(n) = my(x=2^n, y=2*x); while(1, my(v=almost_primes(x, y, n)); for(k=1, #v, if(isprime(v[k]-1), return(v[k]-1))); x=y+1; y=2*x); \\ Daniel Suteu, Jan 07 2025
approved
editing
editing
proposed
(PARI)
almost_primes(A, B, n) = A=max(A, 2^n); (f(m, p, n) = my(list=List()); if(n==1, forprime(q=max(p, ceil(A/m)), B\m, listput(list, m*q)), forprime(q=p, sqrtnint(B\m, n), list=concat(list, f(m*q, q, n-1)))); list); vecsort(Vec(f(1, 2, n)));
a(n) = if(n==0, return(2)); my(x=2^n, y=2*x); while(1, my(v=almost_primes(x, y, n)); for(k=1, #v, if(isprime(v[k]+1), return(v[k]+1))); x=y+1; y=2*x); \\ Daniel Suteu, Jan 07 2025
approved
editing
editing
proposed
2 * 10^6 < a(16) <= 206874667; a(18) = 33577; a(20) <= 3258569.
a(3) = 5 is the smallest prime number of the set {p(i)} = {5, 11, 13, 19, 23, ...} where omega(p(i)^3 + 1) = 3.
a[n_] := Module[{p = 2}, While[PrimeNu[p^n + 1] != n, p = NextPrime[p]]; p]; Print[Array[a, 11]]
allocated for Daniel Suteu
a(n) is the smallest prime p such that omega(p^n + 1) = n.
2, 3, 5, 43, 17, 47, 151, 1697, 59, 2153, 521, 13183, 30089, 66569, 761
1,1
2 * 10^6 < a(16) <= 206874667; a(18) = 33577;
a(3) = 5 is the smallest prime number of the set {p(i)} = {5, 11, 13, 19, 23, ...} where omega(p(i)^3 + 1) = 3.
(PARI) a(n) = forprime(p=2, oo, if(omega(p^n+1) == n, return(p)));
allocated
nonn,more,hard
Daniel Suteu, Jan 06 2025
approved
editing
allocated for Daniel Suteu
recycled
allocated