

A374074


Odd composite numbers k sorted by k/2^(bigomega(k)  1).


3



9, 27, 15, 81, 21, 45, 25, 243, 63, 33, 135, 35, 75, 39, 729, 189, 49, 99, 405, 51, 105, 55, 225, 57, 117, 125, 65, 2187, 69, 567, 147, 297, 1215, 153, 77, 315, 165, 675, 85, 171, 87, 175, 351, 91, 93, 375, 95, 195, 6561, 207, 1701, 441, 111, 891, 3645, 459
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,1


COMMENTS

Sorting by k/2^bigomega(k) would give the same sequence.
It appears that this sequence can be used to approximate the imaginary parts of the nontrivial zeta zeros, that is, A002410(n) is roughly equal to 2*Pi*a(n)/2^bigomega(a(n))  n/2 + sqrt(n)/2.
Calculations show that the relative error approaches 1.0+0.005 for the first 3800 zeros (z=2000 in Mathematica code). For further zeros, a better approximation may be useful, e.g. 2*Pi*a(n)/2^bigomega(a(n))  n/2 + (1/Pi) * n/log(n+1) + (...).


LINKS



FORMULA

A374022(a(n)/2^(bigomega(a(n))1)) ~ n.


EXAMPLE

The odd composite numbers (A071904) are: 9, 15, 21, 25, 27, ... .
Divide by 2^(bigomega()1): 9/2, 15/2, 21/2, 25/2, 27/4, ... .
Sort: 9/2, 27/4, 15/2, 81/8, ... .
Take numerator: this sequence = 9, 27, 15, 81, ... .


MATHEMATICA

(*This algorithm calculates and sorts all noninteger rationals of the form x = k/2^(bigomega(k)1) up to z, where the numerators of x are returned as the sequence*)
z = 100; op = Prime[Range[2, PrimePi[z]]]; (*A065091, odd primes*)
x = Select[Union[Flatten[Outer[Times, op, op/2]]], # <= z &];
For[i = 1, i < Max[1, Floor[Log[3/2, z/2]]]  1, i++, x = Select[Union[x, Flatten[Outer[Times, x, op/2]]], # <= z &]]
a = Numerator[x] (*sequence*)
zzaprx = N@Table[2 Pi a[[i]]/2^PrimeOmega[a[[i]]]  i/2 + Sqrt[i]/2, {i, 1, Length[a]}] (*approximation for zeta zeros*)


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



