|
Search: id:A129733
|
|
|
| A129733 |
|
List of primitive prime divisors of the numbers (3^n-1)/2 (A003462) in their order of occurrence. |
|
+0
6
|
|
| 2, 13, 5, 11, 7, 1093, 41, 757, 61, 23, 3851, 73, 797161, 547, 4561, 17, 193, 34511, 1871, 19, 37, 1597, 363889, 1181, 368089, 67, 661, 47, 1001523179, 6481, 391151, 8951, 398581, 109, 433, 8209, 29, 16493, 59, 20381027, 28537, 31, 271, 683, 4404047, 102673, 21523361
(list; graph; listen)
|
|
|
OFFSET
|
1,1
|
|
|
COMMENT
|
Read A003462 term-by-term, factorize each term, write down any primes not seen before.
|
|
REFERENCES
|
G. Everest et al., Primes generated by recurrence sequences, Amer. Math. Monthly, 11 4 (No. 5, 2007), 417-431.
K. Zsigmondy, Zur Theorie der Potenreste, Monatsh. Math., 3 (1892), 265-284.
|
|
MAPLE
|
# produce sequence
s1:=(a, b, M)->[seq( (a^n-b^n)/(a-b), n=0..M)];
# find primes and their indices
s2:=proc(s) local t1, t2, i; t1:=[]; t2:=[];
for i from 1 to nops(s) do if isprime(s[i]) then
t1:=[op(t1), s[i]];
t2:=[op(t2), i-1]; fi; od; RETURN(t1, t2); end;
# get primitive prime divisors in order
s3:=proc(s) local t2, t3, i, j, k, np; t2:=[]; np:=0;
for i from 1 to nops(s) do t3:=ifactors(s[i])[2];
for j from 1 to nops(t3) do p := t3[j][1]; new:=1;
for k from 1 to np do if p = t2[k] then new:= -1; break; fi; od;
if new = 1 then np:=np+1; t2:=[op(t2), p]; fi; od; od;
RETURN(t2); end;
|
|
CROSSREFS
|
Cf. A003462, A076481, A028491. If 3 is replaced by 2 we get A000225, A004668, A000043, A108974 respectively.
Sequence in context: A120863 A093079 A095417 this_sequence A084160 A128155 A124869
Adjacent sequences: A129730 A129731 A129732 this_sequence A129734 A129735 A129736
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
njas, May 13 2007
|
|
|
Search completed in 0.002 seconds
|
