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

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A246947 a(1)=30; for n > 1, a(n) is the least integer not occurring earlier such that a(n) shares exactly three distinct prime divisors with a(n-1). 5

%I #16 Nov 25 2015 21:43:53

%S 30,60,90,120,150,180,210,42,84,126,168,252,294,336,378,420,70,140,

%T 280,350,490,560,630,105,315,525,735,840,240,270,300,330,66,132,198,

%U 264,396,462,154,308,616,770,110,220,440,550,660,165,495,825,990,360,390,78

%N a(1)=30; for n > 1, a(n) is the least integer not occurring earlier such that a(n) shares exactly three distinct prime divisors with a(n-1).

%C All terms belong to A000977. Is this a permutation of A000977? - _Michel Marcus_, Nov 24 2015

%H Michel Lagneau, <a href="/A246947/b246947.txt">Table of n, a(n) for n = 1..2000</a>

%e 90 is in the sequence because the common prime distinct divisors between a(2)=60 and a(3)=90 are 2, 3 and 5.

%p with(numtheory):a0:={2,3,5}:lst:={}:

%p for n from 1 to 100 do:

%p ii:=0:

%p for k from 30 to 50000 while(ii=0) do:

%p y:=factorset(k):n0:=nops(y):lst1:={}:

%p for j from 1 to n0 do:

%p lst1:=lst1 union {y[j]}:

%p od:

%p a1:=a0 intersect lst1:

%p if {k} intersect lst ={} and a1 <> {} and nops(a1)=3

%p then

%p printf(`%d, `,k):lst:=lst union {k}:a0:=lst1:ii:=1:

%p else

%p fi:

%p od:

%p od:

%t f[s_List]:=Block[{m=s[[-1]],k=30},While[MemberQ[s,k]||Intersection[Transpose[FactorInteger[k]][[1]],Transpose[FactorInteger[m]][[1]]]=={}|| Length[Intersection[Transpose[FactorInteger[k]][[1]],Transpose[FactorInteger[m]][[1]]]]!=3,k++];Append[s,k]];Nest[f,{30},70]

%o (PARI) lista(nn) = {a = 30; print1(a, ", "); fa = (factor(a)[,1])~; va = [a]; k = 0; while (k!= nn, k = 1; while (!((#setintersect(fa, (factor(k)[,1])~) == 3) && (! vecsearch(va, k))), k++); a = k; print1(a, ", "); fa = (factor(a)[,1])~; va = vecsort(concat(va, k)););} \\ _Michel Marcus_, Nov 24 2015

%Y Cf. A000977, A064413, A246946.

%K nonn

%O 1,1

%A _Michel Lagneau_, Sep 08 2014

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 05:48 EDT 2024. Contains 371265 sequences. (Running on oeis4.)