|
|
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
PROG
|
(Python)
from bisect import bisect
from sympy import totient
k = 10**n
plist, qlist, c = tuple(-totient(i) for i in range(1, k+1)), [0]*(k+1), 0
for i in range(k):
qlist[a:=bisect(qlist, plist[i], lo=1, hi=c+1, key=lambda x:plist[x])]=i
c = max(c, a)
return c
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,hard,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|