The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 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. 1
 6, 15, 725, 91, 15851, 30589, 6977093777, 703 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,1 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. LINKS Table of n, a(n) for n=1..8. EXAMPLE 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. 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 *) PROG (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, ", "); ); } CROSSREFS Cf. A140110, A280963. Sequence in context: A145612 A307812 A354682 * A191463 A235943 A225304 Adjacent sequences: A280961 A280962 A280963 * A280965 A280966 A280967 KEYWORD nonn,more AUTHOR Michel Marcus, Jan 11 2017 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified March 4 23:31 EST 2024. Contains 370537 sequences. (Running on oeis4.)