|
|
A333967
|
|
Subsequence of A071395. The extra constraint is m is not a term if m*q/p is abundant where prime p|m and q is the least prime larger than p.
|
|
2
|
|
|
70, 2002, 3230, 4030, 5830, 8415, 8925, 20482, 32445, 45885, 51765, 83265, 107198, 131054, 133042, 178486, 206770, 253270, 253946, 258970, 270470, 310930, 330310, 334305, 362710, 442365, 474045, 497835, 513890, 544310, 567765, 589095, 592670, 602175, 617265, 631670, 654675
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
EXAMPLE
|
70 is in the sequence as it's abundant. Its prime factorization is 2 * 5 * 7. Each of 3 * 5 * 7, 2 * 7 * 7 and 2 * 5 * 11 are deficient and no divisor of 70 is in this sequence.
|
|
MATHEMATICA
|
primabQ[n_] := DivisorSigma[1, n] > 2n && AllTrue[Most @ Divisors[n], DivisorSigma[1, #] < 2# &]; seqQ[n_] := Module[{f = FactorInteger[n]}, p = f[[;; , 1]]; e = f[[;; , 2]]; q = NextPrime[p]; AllTrue[n*(q/p), DivisorSigma[1, #] <= 2# &]]; Select[Range[10^5], primabQ[#] && seqQ[#] &] (* Amiram Eldar, Jul 05 2020 *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|