|
|
A067665
|
|
The start of a record-setting run of consecutive integers i with distinct A001222(i).
|
|
4
|
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The list of indices of record terms in A068796;
n is in the sequence if A068796(n) is larger than A068796(m) when n is larger than m. For the known terms, f(a(n)) = n+1. Is that true for all n? In other words, is the monotonic subsequence of A068796 identical to A020725?
|
|
LINKS
|
|
|
EXAMPLE
|
The values of f(n) for n=1 to 15 are 2,1,2,2,2,3,3,2,1,3,2,3,2,1,4. Records occur at f(1)=2, f(6)=3 and f(15)=4.
|
|
MATHEMATICA
|
bigomega[n_] := Plus@@Last/@FactorInteger[n]; f[n_] := For[k=1; s={bigomega[n]}, True, k++, If[MemberQ[s, z=bigomega[n+k]], Return[k], AppendTo[s, z]]]; For[n=1; max=0, True, n++, If[f[n]>max, Print[n, " ", max=f[n]]]]
|
|
PROG
|
(PARI) a(n, lim=1e12, startAt=1)={
forstep(i=startAt-1, lim, 10^6-n,
my(v=vectorsmall(min(10^6, lim\1-i), j, bigomega(j+i)));
for(j=n, #v, if(#vecsort(v[j-n+1..j], , 8)==n, return(j+i-n+1)))
)
|
|
CROSSREFS
|
|
|
KEYWORD
|
more,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|