login
This site is supported by donations to The OEIS Foundation.

 

Logo


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: A145558 A145612 A307812 * A191463 A068941 A235943

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 | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 15 07:56 EDT 2019. Contains 328026 sequences. (Running on oeis4.)