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.
k = 1; t = Table[0, {50}]; While[k < 500000001, a = PrimeOmega[k] + PrimeOmega[3 k - 1] - 1; If[ t[[a + 1]] == 0, t[[a + 1]] = k; Print[{k, a}]]; k++]; t
a(3) = 28, because 28 is a triangular number with 3 prime factors (counted with multiplicity) {2, 2, 7} and this is the smallest such number.
g:= proc(s) local n, p, F;
for n from 1 to 10^7 do
p:= (s-2)*n*(n-1)/2 + n;
if numtheory:-bigomega(p) = s then return p fi;
od
end proc:
map(g, [$2..30]); # Robert Israel, Jan 15 2023
sng[n_]:=Module[{k=1},While[PrimeOmega[PolygonalNumber[n,k]]!=n,k++];PolygonalNumber[ n,k]]; Array[sng,21,2] (* The program generates the first 20 terms of the sequence. *) (* Harvey P. Dale, Feb 19 2023 *)
a(n) = if(n<3, return()); for(k=1, oo, my(t=(k*(n*k - n - 2*k + 4))\2); if(bigomega(t) == n, return(t))); \\ Daniel Suteu, Dec 04 2022
bigomega_polygonals(A, B, n, k) = A=max(A, 2^n); (f(m, p, n) = my(list=List()); if(n==1, forprime(q=max(p,ceil(A/m)), B\m, my(t=m*q); if(ispolygonal(t,k), listput(list, t))), forprime(q = p, sqrtnint(B\m, n), my(t=m*q); if(ceil(A/t) <= B\t, list=concat(list, f(t, q, n-1))))); list); vecsort(Vec(f(1, 2, n))); a(n, k=n) = if(k < 3, return()); my(x=2^n, y=2*x); while(1, my(v=bigomega_polygonals(x, y, n, k)); if(#v >= 1, return(v[1])); x=y+1; y=2*x); \\ Daniel Suteu, Dec 04 2022
Comments