A280964 a(n) is the least term that is not a power of p, the n-th prime, in the sequence of numbers whose consecutive divisors have a ratio satisfying numerator - denominator = p-1.

6, 15, 725, 91, 15851, 30589, 6977093777, 703

Offset: 1

Author

Michel Marcus, Jan 11 2017

Keywords

• nonn
• more

Comments

For primes 17 < p < 100, index n, prime p and known values: (11, 31, 1891), (12, 37, 270), (13, 41, 2893001), (14, 43, 6324702343), (16, 53, 8036549), (17, 59, 12319259), (22, 79, 12403), (23, 83, 48023219), (24, 89, 63439289), (25, 97, 18721).

For some other primes, possible values were found, but they may not be the least.

p=23, a(9) <= 78317141783;

p=29, a(10) <= 176994576151121533000319046029;

p=47, a(15) <= 7654455761751330268890575447204341560894807321851181994954582195247;

p=61, a(18) <= 9876832533549882665273701;

p=67, a(19) <= 4483224940666198270986387212125921182677251147716597066996042973310060661858239362623427;

p=71, a(20) <= 269277687648484922419868038102292093702042951;

p=73, a(21) <= 151362235513.

Examples

			a(1) = 6 is the first non-(prime power) of 2 (the 1st prime) in A140110.
a(2) = 15 is the first non-(prime power) of 3 (the 2nd prime) in A280963.
For n=3, p=5, the sequence begins: 1, 5, 25, 125, 625, 725, ... so a(3)=725.
For n=4, p=7, the sequence begins: 1, 7, 49, 91, 343, ... so a(4)=91.
		

Crossrefs

Cf. A140110, A280963.

Programs

Mathematica

Table[Function[p, SelectFirst[Range[10^5], Function[n, And[! IntegerQ@ Log[p, n], Times @@ Boole@ Map[Abs[Numerator@ # - Denominator@ #] == p - 1 &[#2/#1] & @@ # &, Partition[Divisors@ n, 2, 1]] > 0]] ]]@ Prime@ n /. k_ /; MissingQ@ k -> Nothing, {n, 6}] (* Michael De Vlieger, Jan 13 2017, Version 10.2 *)

PARI

isok(n, plus) = {my(vd = divisors(n)); for (k=1, #vd - 1, r = vd[k+1]/vd[k]; if (numerator(r) != denominator(r) + plus, return(0));); return(1);}
isokp(k) = (k!= 1) && !isprime(k) && ! isprimepower(k);
findnp(n) = {k = 1; ok = 0; nb = 0; while (! ok, okk = isok(k, n); if (okk, ok = isokp(k); if (! ok, nb++);); if (!ok, k++);); k;}
lista(nn) = for(n=1, nn, print1(findnp(prime(n)-1), ", "));

PARI

/ * alternate program; does it give the least terms? */
isokplus(n, plus) = {my(vd = divisors(n)); for (k=1, #vd - 1, r = vd[k+1]/vd[k]; if (numerator(r) != denominator(r) + plus, return(0));); return(1);}
findqq(p) = {ok = 0; ip = 1; while (!ok, ik = 1; while (!ok, if (isprime(q= ik*p^ip+(p-1)) && isokplus(p^ip*q, p-1), return([p^ip, q])); ik++; if (ik > p, break);); ip ++;); return ([]);}
listff(nn) = {for (n=1, nn, p = prime(n); v = findqq(p); pp = v[1]; q = v[2]; print1(pp*q, ", "););}