Hi, Thanks for your posts. No good today, travel again, just time for a few remarks:
Li function indeed better behaved than pi; continuous and différentiable, just that singularity to deal with. Have seen two versions, with different lower limits of integral; one with 0 the other with 2. Don’t suppose it makes much difference as
x—>00 but as matter of interest which are you using?
Despite the change from p/2 to 2*p I got correct data; strange.
Might be better for oeis if code for data can be decoupled from code for plot. Latter could then be a link (they have a separate section for links in any submission).
Thought: Consider triple (a,b,c) where a,b,c are all primes <prime p and such that a*b*c==1mod p. Presumably your code could be adapted to find such triples for arbitrary p, just as it did for the pairs (a,b) ?
eg: 2*3*5==1mod 29.
This might lead to an asymptotic formula, residues etc, similar to the case for prime pairs. Likewise, generalisation to r-tupples:
(a,b,c,d.....r) might be worth looking at. Does this sound interesting?
Dashing now catch up later