Oh, cool, I forgot this is essentially solvable and I'm correct- ish:

2657/prime(7^2-7/4) ~= 12

prime(7^2-7/4) ~= 2657/12

7^2-7/4 ~= pi function(2657/12)

And

7^2-2 = pi function(2657/12)

prime(23)*x = 2657

--> x~=32

prime(24)*x = 2657

--> x~=30

The prior was based on:

(3^2-1) = pi function(2657/120)
(7^2-2) = pi function(2657/12)

And made me curious about

15 = 3*5
31 ~= 3^(3/2)*6
63 = 3^2*7

sqrt(pi function({1,2,3,7,11,19,43,67,163}))

-->

{0, 1, sqrt(2), 2, sqrt(5), sqrt(8), sqrt(14), sqrt(19), sqrt(38)}

{0, 1, 2}

sqrt({5, 8}) --> 4 + 1, 9 - 1; (8 - 5) + 1 = 4, (8 + 5) + 3 = 16

sqrt({2, 19}) --> ?

sqrt({14, 38}) --> 16 - 2, 36 + 2; (38 - 14) + 1 = 25, (38 + 14) - 3 = 49

2+19 = 22 - 1
2*19 = 37 + 1
2<*>19 ?= 58 // can't find an operator

19*2^e ~= 5^3
19*2^pi ~= 13^2

(19*2^pi - 13^2)/(19*2^e - 5^3) ~= -36

This made me curious about:

(7)^(23/9) ~= 12^2

Which has interesting complex roots as well.

I just realized

19*2^((e+pi)/2) ~= 144.79...
19*2^((e*pi)^.5) ~= 144.02...

So it checks out pretty well with the previous idea.

sqrt(144.5*2) = 17

(19^2 + 2)/3 = 11^2
(12^2 + 1/2)*2 = 17^2

Which I'm gonna try to relate back to 12 & 19 in music theory.

((15.5^2)*(2/3))-14^2 ~= -36

(36 - (19*2^pi - 13^2)/(19*2^e - 5^3))/6 ~= 12

2657 = 2000*2^(-e + π) - 5^2

19*2^2.6854520010653064453097 ~= 122.2236...
19*2^3.2758229187218111597876 ~= 184.0240...

https://en.m.wikipedia.org/wiki/Khinchin%27s_constant
https://en.m.wikipedia.org/wiki/L%C3%A9vy%27s_constant

For the above, the values are essentially 11^2 + epsilon & (17/2)^2 + epsilon, which sort of relates back to:

(19^2 + 2)/3 = 11^2
(12^2 + 1/2)*2 = 17^2

28 = 3^4-prime(2^4)

-163 = 4^4-prime(3^4)

factor 14^4-prime(3^4) --> 37997
factor 15^4-prime(6^4) --> 39998 = 2*7*2857

I'm curious about:

28,30,34,43,67,163,953...

Which is based on

floor(sqrt(k*28)), k=163

Exploration is like:

floor(sqrt(k*28)), k=32437
floor(sqrt(k*28)), k=32504

primes near 32470 (7 candidates)

a(n+1) = floor(sqrt(a(n)*163)), a(0)=28

n | a(n)
0 | 28
1 | 67
2 | 104
3 | 130
4 | 145
5 | 153
6 | 157
7 | 159
8 | 160
9 | 161

Where 67, 130, 145, 157 are of interest to me, particularly in the context of

https://math.stackexchange.com/questions/1309541/relatives-of-heegner-numbers

I'm thinking 145 & 34 are the relevant ones here:

(6^2 - 2)/2 = 17
(12^2 + 1/2)*2 = 17^2

Which may imply 17 is relevant to music theory.

https://en.m.wikipedia.org/wiki/34_equal_temperament

So I'm thinking 72, 77, 80 is relevant here, for:

https://en.m.wikipedia.org/wiki/Music_and_mathematics
https://en.m.wikipedia.org/wiki/34_equal_temperament

9^2, 80 for 12 notes
5^2, 24/"72" for 17 notes

I'm thinking about https://oeis.org/A107360

{3, 5, 7, 13, 17, 19, 31, 61, 127 }

2^n - 3 = p --> {3, 4, 6} | {5, 13, 61}
2^n - 1 = p --> {2, 3, 5, 7} | {3, 7, 31, 127}
2^n + 1 = p --> {1, 2, 4} | {3, 5, 17}
2^n + 3 = p --> {-inf, 1, 2, 4} | {3, 5, 7, 19}

3 occurs 3x
5 occurs 3x
7 occurs 2x
All other terms occur 1x

Another aggregation:

2 happens 3x
4 happens 3x
1 happens 2x
3 happens 2x
All others -inf, 5-7 happen 1x

I'm looking at

2^({1, 2, 4, 6, 8, 12, 16, 18, 20, 24, 28, 30, 32, 36, 40, 42, 48, 54, 56, 60, 64, 66, 72, 78, 80, 84, 88, 90, 96, 100, 104, 108, 112, 120, 126, 128, 132, 140, 144, 150}+2) - 3 (practical numbers)

Example:

factor {5, 13, 61, 253, 1021, 16381, 262141, 1048573, 4194301, 67108861, 1073741821, 4294967293, 17179869181, 274877906941, 4398046511101, 17592186044413, 1125899906842621, 72057594037927933}

Prime for

{5, 13, 61, 1021, 16381, 1048573, 4194301}

It looks like practical numbers may not be relevant, based on

Table[(n-2)*Boole[isprime(2^n - 3)]], n= 2 to 1000

I was looking at

"3, 4, 6 -seq:7 -seq:9 -seq:11 -seq:13 -seq:15" on oeis.org

And realized the 3, 4, 6... may just be prime(n) + 1

There's two features in this picture that interest me:

https://grantjenks.com/wiki/_detail/projects/automata/primediffabs500x1.jpg%3Fid=projects:automata:absolute_differences_in_consecutive_prime_numbers.html

One is at coordinate (165, 363), the other is at (318, 477), which may essentially be related to the Heegner number 163 (& 2*163).

You can draw a 45 degree angle and find another feature of interest at (250, 450) -- that's my estimate based on the features at (240, 422) & (254, 437):

This makes me interested in prime(54), and prime(38) / prime(70).

I'm leaning towards 163, 251, 342 now.

Now I'm leaning towards 163, 251, 339:

This image shows a line that essentially is supposed to represent 339 at location (338, 367):

For the above image, you can't see the bottom right of the line segment, but it seems to have 5/20 pixels to the left, unlike the other one which is a 50:50 split.

I'm a little interested in:

prime(20)*e ~= 193
prime(30)*pi ~= 355 = 5*prime(20)

(339+71)/(5/2) = 163+1 = 41*4

So if sigma(n) < 2*n - 1 (powers of 2), I think n cannot be a practical number. I'm curious if there is an upper limit where n must be a practical number.

Creating code isn't that hard, although sigma and practical numbers are obscure-ish.

from itertools import chain, cycle, accumulate, combinations
from typing import List, Tuple

# %% Factors

def factors5(n: int) -> List[int]:
"""Factors of n, (but not n)."""
def prime_powers(n):
# c goes through 2, 3, 5, then the infinite (6n+1, 6n+5) series
for c in accumulate(chain([2, 1, 2], cycle([2,4]))):
if c*c > n: break
if n%c: continue
d,p = (), c
...

...

while not n%c:
n,p,d = n//c, p*c, d + (p,)
yield(d)
if n > 1: yield((n,))

r = [1]
for e in prime_powers(n):
r += [a*b for a in r for b in e]
return r[:-1]

# %% Powerset

def powerset(s: List[int]) -> List[Tuple[int, ...]]:
"""powerset([1,2,3]) --> () (1,) (2,) (3,) (1,2) (1,3) (2,3) (1,2,3) ."""
return chain.from_iterable(combinations(s, r) for r in range(1, len(s)+1))
...