{VERSION 3 0 "DEC ALPHA UNIX" "3.0" }
{USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 
1 0 0 0 0 }{CSTYLE "2D Math" -1 2 "Times" 0 1 0 0 0 0 0 0 2 0 0 0 0 0 
0 }{CSTYLE "2D Output" 2 20 "" 0 1 0 0 255 1 0 0 0 0 0 0 0 0 0 }
{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 
0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Text Output" -1 2 1 
{CSTYLE "" -1 -1 "Courier" 1 10 0 0 255 1 0 0 0 0 0 1 3 0 3 }1 0 0 -1 
-1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 2 6 1 {CSTYLE "" -1 -1 "" 0 1 0 0 
0 0 0 0 0 0 0 0 2 0 0 }0 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "Maple
 Output" 0 11 1 {CSTYLE "" -1 -1 "" 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }3 
3 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }{PSTYLE "" 11 12 1 {CSTYLE "" -1 -1 "" 
0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 }1 0 0 -1 -1 -1 0 0 0 0 0 0 -1 0 }}
{SECT 0 {EXCHG {PARA 0 "" 0 "" {TEXT -1 70 "On reprend la version de l
'exercice precedent, a titre de comparaison." }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 97 "L2:=[1, 2, 3, 7, 12, 21, 29, 38, 43, 47, 48, 49,
 51, 52, 53, 57, 62, 71, 79, 88, 93, 97, 98, 99];" }}{PARA 12 "" 1 "" 
{XPPMATH 20 "6#>%#L2G7:\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#H\"#Q\"#V\"#Z\"
#[\"#\\\"#^\"#_\"#`\"#d\"#i\"#r\"#z\"#))\"#$*\"#(*\"#)*\"#**" }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 235 "is149_2:=proc(n)\n  local r
,s;\n  r:=n;\n  while r<>0 do\n    s:=irem(r,100,'r');\n    if not (me
mber(s,\{11,14,19,41,44,49,91,94,99\})\n            or \n            (
member(s,\{1,4,9\}) and r=0)) then\n      RETURN(false)\n    fi\n  od;
\n  true\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 232 "tricol_2
:=proc(N)\n  local k,res,q,r,n,i;\n  global L2;\n  k:=0;\n  for q from
 0 by 100 while k<N do\n    for r in L2 do\n      if is149_2((q+r)^2) \+
then\n        k:=k+1;\n        res[k]:=q+r;      fi\n    od\n  od;\n  \+
[seq(res[i],i=1..N)]\nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
42 "start:=time():\ntricol_2(12);\ntime()-start;" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#7.\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J\"&
2,(\"')G(Q" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"&81*!\"$" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 25 "Ensuite, c'est trop long." }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 58 "On ameliore. On met la largeur de la tran
che en parametre." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 13 "slices
ize:=4;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*slicesizeG\"\"%" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 89 "On construit les mots de longueur \+
inferieure ou egale a slicesize sur l'alphabet \{1,4,9\}." }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 17 "with(combstruct):" }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 89 "word:=\{W=Union(Epsilon,Prod(A,W)),\n    \+
   A=Union(a,b,c),\n       a=Atom, b=Atom, c=Atom\}:" }}}{EXCHG {PARA 
0 "" 0 "" {TEXT -1 125 "QQ contient les mots de longueur strictement p
lus petite que slicesize.\nQ contient les mots de longueur exactement \+
slicesize." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 292 "for s to sli
cesize do     W.s:=allstructs([W,word],size=s);\n  WW.s:=\{seq([eval(s
ubs(Prod=proc() args end,Epsilon=NULL,a=1,b=4,c=9,i))],i=W.s)\};\n  LL
.s:=map(proc(l) local i;add(l[i]*10^(i-1),i=1..nops(l)) end,WW.s)\nod:
\nQQ.slicesize:=\{seq(op(LL.s),s=1..slicesize-1)\};\nQ.slicesize:=LL.s
licesize;" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%$QQ4G<I\"\"\"\"\"%\"\"*
\"#6\"#9\"#>\"#T\"#W\"#\\\"#\"*\"#%*\"#**\"$6\"\"$9\"\"$>\"\"$T\"\"$W
\"\"$\\\"\"$*\\\"$W%\"$>%\"$%>\"$\\%\"$*>\"$6*\"$6%\"$W*\"$T%\"$9*\"$
\\*\"$\"**\"$9%\"$%**\"$***\"$>*\"$T*\"$\">\"$%\\\"$\"\\" }}{PARA 12 "
" 1 "" {XPPMATH 20 "6#>%#Q4G<]p\"%W6\"%\\**\"%9>\"%\">\"\"%\">*\"%%\\%
\"%>>\"%>\"*\"%TW\"%%*\\\"%****\"%6**\"%%*>\"%6>\"%*\\\"\"%*>\"\"%W**
\"%*\\*\"%\\9\"%>T\"%96\"%>6\"%%***\"%6T\"%W\"*\"%W>\"%\"\\%\"%6%*\"%
\\6\"%\"\\*\"%**\\\"%9W\"%66\"%\\\"*\"%\"\\\"\"%T%*\"%TT\"%WW\"%T9\"%9
**\"%9T\"%T\"*\"%%>\"\"%99\"%\\%*\"%W\\\"%T**\"%T6\"%6\\\"%%\\*\"%*>*
\"%6W\"%\\W\"%%>%\"%**>\"%%>*\"%%\\\"\"%\"***\"%\\\\\"%\">%\"%6\"*\"%>
\\\"%9\"*\"%69\"%\"*\\\"%>**\"%>W\"%>9\"%9\\\"%T\\\"%W9\"%\\>\"%\"*>\"
%T>\"%W%*\"%WT\"%*\\%\"%>%*\"%*>%\"%\\T\"%9%*" }}}{EXCHG {PARA 0 "" 0 
"" {TEXT -1 104 "On cherche les residus modulo 10^slicesize des motifs
 precedents; on les ordonne et on les range dans L." }}}{EXCHG {PARA 
0 "> " 0 "" {MPLTEXT 1 0 116 "L.slicesize:=sort(map2(op,2,map(op,\n  [
seq(msolve(r^2=d,10^slicesize),\n       d=QQ.slicesize union Q.slicesi
ze)])));" }}{PARA 12 "" 1 "" {XPPMATH 20 "6#>%#L4G7dt\"\"\"\"\"#\"\"$
\"\"(\"#7\"#@\"#Q\"$2\"\"$V\"\"$7#\"$$R\"$z%\"$@&\"$Q&\"$H(\"$r(\"$)y
\"$d)\"%75\"%z6\"%H7\"%\\7\"%`7\"%d7\"%d8\"%$R\"\"%)[\"\"%V;\"%7<\"%H<
\"%i>\"%@?\"%2@\"%)G#\"%HC\"%iC\"%)[#\"%)\\#\"%-D\"%7D\"%QD\"%rD\"%7F
\"%$*G\"%zH\"%QI\"%rK\"%)G$\"%dL\"%7N\"%2O\"%VO\"%VP\"%ZP\"%^P\"%rP\"%
@Q\"%))R\"%VT\"%7U\"%HU\"%rU\"%iW\"%zW\"%@X\"%2Y\"%)y%\"%d[\"%$*[\"%i
\\\"%z\\\"%))\\\"%$*\\\"%(*\\\"%)*\\\"%**\\\"%,]\"%-]\"%.]\"%2]\"%7]\"
%@]\"%Q]\"%2^\"%V^\"%7_\"%$R&\"%za\"%@b\"%Qb\"%Hd\"%rd\"%)y&\"%de\"%7g
\"%zh\"%Hi\"%\\i\"%`i\"%di\"%dj\"%$R'\"%)['\"%Vm\"%7n\"%Hn\"%ip\"%@q\"
%2r\"%)G(\"%Hu\"%iu\"%)[(\"%)\\(\"%-v\"%7v\"%Qv\"%rv\"%7x\"%$*y\"%zz\"
%Q!)\"%r#)\"%)G)\"%d$)\"%7&)\"%2')\"%V')\"%V()\"%Z()\"%^()\"%r()\"%@))
\"%))*)\"%V\"*\"%7#*\"%H#*\"%r#*\"%i%*\"%z%*\"%@&*\"%2'*\"%)y*\"%d)*\"
%$*)*\"%i**\"%z**\"%))**\"%$***\"%(***\"%)***\"%****" }}}{EXCHG {PARA 
0 "" 0 "" {TEXT -1 69 "On ecrit une procedure de selection dependant d
u parametre slicesize." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 310 "
is149_.slicesize:=eval(subs(_Q=Q.slicesize,_QQ=QQ.slicesize,_slicesize
=slicesize,\nproc(n::posint)\n  local r,s;\n  r:=n;\n  while r<>0 do\n
    s:=irem(r,10^_slicesize,'r');\n    if not (member(s,_Q)\n         \+
   or \n            (member(s,_QQ) and r=0)\n           ) then\n      \+
RETURN(false)\n    fi\n  od;\n  true\nend));" }}{PARA 12 "" 1 "" 
{XPPMATH 20 "6#>%(is149_4GR6#'%\"nG%'posintG6$%\"rG%\"sG6\"F-C%>8$9$?(
F-\"\"\"F3F-0F0\"\"!C$>8%-%%iremG6%F0\"&++\".F0@$45-%'memberG6$F8<]p\"
%W6\"%\\**\"%9>\"%\">\"\"%\">*\"%%\\%\"%>>\"%>\"*\"%TW\"%%*\\\"%****\"
%6**\"%%*>\"%6>\"%*\\\"\"%*>\"\"%W**\"%*\\*\"%\\9\"%>T\"%96\"%>6\"%%**
*\"%6T\"%W\"*\"%W>\"%\"\\%\"%6%*\"%\\6\"%\"\\*\"%**\\\"%9W\"%66\"%\\\"
*\"%\"\\\"\"%T%*\"%TT\"%WW\"%T9\"%9**\"%9T\"%T\"*\"%%>\"\"%99\"%\\%*\"
%W\\\"%T**\"%T6\"%6\\\"%%\\*\"%*>*\"%6W\"%\\W\"%%>%\"%**>\"%%>*\"%%\\
\"\"%\"***\"%\\\\\"%\">%\"%6\"*\"%>\\\"%9\"*\"%69\"%\"*\\\"%>**\"%>W\"
%>9\"%9\\\"%T\\\"%W9\"%\\>\"%\"*>\"%T>\"%W%*\"%WT\"%*\\%\"%>%*\"%*>%\"
%\\T\"%9%*3-FB6$F8<IF3\"\"%\"\"*\"#6\"#9\"#>\"#T\"#W\"#\\\"#\"*\"#%*\"
#**\"$6\"\"$9\"\"$>\"\"$T\"\"$W\"\"$\\\"\"$*\\\"$W%\"$>%\"$%>\"$\\%\"$
*>\"$6*\"$6%\"$W*\"$T%\"$9*\"$\\*\"$\"**\"$9%\"$%**\"$***\"$>*\"$T*\"$
\">\"$%\\\"$\"\\/F0F5-%'RETURNG6#%&falseG%%trueGF-F-F-" }}}{EXCHG 
{PARA 0 "" 0 "" {TEXT -1 134 "On ecrit la procedure de generation des \+
mots tricolores ou le test est base sur les motifs de taille slicesize
 pour les poids faibles." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
440 "tricol_.slicesize:=subs(_L=L.slicesize,\n                        \+
_slicesize=slicesize,\n                        _is149=is149_.slicesize
,\nproc(N::posint)\n  local k,res,q,r,n,i;\n  k:=0;\n  for q from 0 by
 10^_slicesize while k<N do\n    userinfo(3,procname,\"testing q=\",q)
;\n    for r in _L do\n      if _is149((q+r)^2) then\n        k:=k+1;
\n        res[k]:=q+r;\n        userinfo(2,procname,res[k])\n      fi
\n    od\n  od;\n  [seq(res[i],i=1..N)]\nend);" }}{PARA 12 "" 1 "" 
{XPPMATH 20 "6#>%)tricol_4GR6#'%\"NG%'posintG6(%\"kG%$resG%\"qG%\"rG%
\"nG%\"iG6\"F1C%>8$\"\"!?(8&F5\"&++\"F12F49$C$-%)userinfoG6&\"\"$9!Q+t
esting~q=6\"F7?&8'7dt\"\"\"\"\"#F?\"\"(\"#7\"#@\"#Q\"$2\"\"$V\"\"$7#\"
$$R\"$z%\"$@&\"$Q&\"$H(\"$r(\"$)y\"$d)\"%75\"%z6\"%H7\"%\\7\"%`7\"%d7
\"%d8\"%$R\"\"%)[\"\"%V;\"%7<\"%H<\"%i>\"%@?\"%2@\"%)G#\"%HC\"%iC\"%)[
#\"%)\\#\"%-D\"%7D\"%QD\"%rD\"%7F\"%$*G\"%zH\"%QI\"%rK\"%)G$\"%dL\"%7N
\"%2O\"%VO\"%VP\"%ZP\"%^P\"%rP\"%@Q\"%))R\"%VT\"%7U\"%HU\"%rU\"%iW\"%z
W\"%@X\"%2Y\"%)y%\"%d[\"%$*[\"%i\\\"%z\\\"%))\\\"%$*\\\"%(*\\\"%)*\\\"
%**\\\"%,]\"%-]\"%.]\"%2]\"%7]\"%@]\"%Q]\"%2^\"%V^\"%7_\"%$R&\"%za\"%@
b\"%Qb\"%Hd\"%rd\"%)y&\"%de\"%7g\"%zh\"%Hi\"%\\i\"%`i\"%di\"%dj\"%$R'
\"%)['\"%Vm\"%7n\"%Hn\"%ip\"%@q\"%2r\"%)G(\"%Hu\"%iu\"%)[(\"%)\\(\"%-v
\"%7v\"%Qv\"%rv\"%7x\"%$*y\"%zz\"%Q!)\"%r#)\"%)G)\"%d$)\"%7&)\"%2')\"%
V')\"%V()\"%Z()\"%^()\"%r()\"%@))\"%))*)\"%V\"*\"%7#*\"%H#*\"%r#*\"%i%
*\"%z%*\"%@&*\"%2'*\"%)y*\"%d)*\"%$*)*\"%i**\"%z**\"%))**\"%$***\"%(**
*\"%)***\"%****%%trueG@$-%(is149_4G6#*$),&F7FFFDFFFG\"\"\"C%>F4,&F4FFF
FFF>&8%6#F4F^w-F=6%FGF@Fdw7#-%$seqG6$&Few6#8)/F_x;FFF:F1F1F1" }}}
{EXCHG {PARA 0 "" 0 "" {TEXT -1 9 "On teste." }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 32 "infolevel[tricol_.slicesize]:=3:" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "start:=time():\ntricol_.slicesize(1
0);\ntime()-start;" }}{PARA 6 "" 1 "" {TEXT -1 34 "`\\`tricol_4:\\`   \+
\"testing q=\"   0`" }}{PARA 6 "" 1 "" {TEXT -1 19 "`\\`tricol_4:\\`  \+
 1`" }}{PARA 6 "" 1 "" {TEXT -1 19 "`\\`tricol_4:\\`   2`" }}{PARA 6 "
" 1 "" {TEXT -1 19 "`\\`tricol_4:\\`   3`" }}{PARA 6 "" 1 "" {TEXT -1 
19 "`\\`tricol_4:\\`   7`" }}{PARA 6 "" 1 "" {TEXT -1 20 "`\\`tricol_4
:\\`   12`" }}{PARA 6 "" 1 "" {TEXT -1 20 "`\\`tricol_4:\\`   21`" }}
{PARA 6 "" 1 "" {TEXT -1 20 "`\\`tricol_4:\\`   38`" }}{PARA 6 "" 1 "
" {TEXT -1 21 "`\\`tricol_4:\\`   107`" }}{PARA 6 "" 1 "" {TEXT -1 21 
"`\\`tricol_4:\\`   212`" }}{PARA 6 "" 1 "" {TEXT -1 38 "`\\`tricol_4:
\\`   \"testing q=\"   10000`" }}{PARA 6 "" 1 "" {TEXT -1 38 "`\\`tric
ol_4:\\`   \"testing q=\"   20000`" }}{PARA 6 "" 1 "" {TEXT -1 38 "`\\
`tricol_4:\\`   \"testing q=\"   30000`" }}{PARA 6 "" 1 "" {TEXT -1 
23 "`\\`tricol_4:\\`   31488`" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7,\"
\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"$7(!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 
32 "infolevel[tricol_.slicesize]:=1:" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 13 "slicesize:=6;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#>%*s
licesizeG\"\"'" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 292 "for s to
 slicesize do     W.s:=allstructs([W,word],size=s);\n  WW.s:=\{seq([ev
al(subs(Prod=proc() args end,Epsilon=NULL,a=1,b=4,c=9,i))],i=W.s)\};\n
  LL.s:=map(proc(l) local i;add(l[i]*10^(i-1),i=1..nops(l)) end,WW.s)
\nod:\nQQ.slicesize:=\{seq(op(LL.s),s=1..slicesize-1)\}:\nQ.slicesize:
=LL.slicesize:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 120 "L.slices
ize:=sort(map2(op,2,map(op,\n    [seq(msolve(r^2=d,10^slicesize),\n   \+
      d=QQ.slicesize union Q.slicesize)]))):" }}}{EXCHG {PARA 0 "> " 
0 "" {MPLTEXT 1 0 368 "is149_.slicesize:=eval(subs(_Q=Q.slicesize,\n  \+
                          _QQ=QQ.slicesize,\n                         \+
   _slicesize=slicesize,\nproc(n::posint)\n  local r,s;\n  r:=n;\n  wh
ile r<>0 do\n    s:=irem(r,10^_slicesize,'r');\n    if not (member(s,_
Q)\n            or \n            (member(s,_QQ) and r=0)\n           )
 then\n      RETURN(false)\n    fi\n  od;\n  true\nend)):" }}}{EXCHG 
{PARA 0 "> " 0 "" {MPLTEXT 1 0 440 "tricol_.slicesize:=subs(_L=L.slice
size,\n                        _slicesize=slicesize,\n                \+
        _is149=is149_.slicesize,\nproc(N::posint)\n  local k,res,q,r,n
,i;\n  k:=0;\n  for q from 0 by 10^_slicesize while k<N do\n    userin
fo(3,procname,\"testing q=\",q);\n    for r in _L do\n      if _is149(
(q+r)^2) then\n        k:=k+1;\n        res[k]:=q+r;\n        userinfo
(2,procname,res[k])\n      fi\n    od\n  od;\n  [seq(res[i],i=1..N)]\n
end):" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 51 "start:=time():\ntr
icol_.slicesize(10);\ntime()-start;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6
#7,\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J" }}{PARA 11 "" 
1 "" {XPPMATH 20 "6#$\"%?G!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 
1 0 51 "start:=time():\ntricol_.slicesize(12);\ntime()-start;" }}
{PARA 11 "" 1 "" {XPPMATH 20 "6#7.\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$
2\"\"$7#\"&)[J\"&2,(\"')G(Q" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\"%)>$
!\"$" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 20 "Revenons en arriere." }}}
{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 42 "start:=time():\ntricol_4(12)
;\ntime()-start;" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7.\"\"\"\"\"#\"\"$
\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J\"&2,(\"')G(Q" }}{PARA 11 "" 1 "" 
{XPPMATH 20 "6#$\"%&H(!\"$" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 150 "Bi
lan: pour les tailles 2, 4, 6 on a respectivement 90, 7 et 3 secondes \+
pour N=12. \nLe gain est clair, mais le pretraitement est de plus en p
lus long." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 113 "Pour la taille 6,  \+
on n'atteint pas N=14. Si on passe a la taille 8, le pretraitement dev
ient excessivement long." }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 33 "On pa
sse a une version dynamique." }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 
0 1603 "tricol:=proc(N::posint)\n  local res,k,A,AA,oldQ,newQ,newL,\n \+
       slicesize,oldL,is149,qq,q,l,r,n,n2,\n        tail,newT,i,j,st;
\n  st:=time();\n  res[1]:=1;\n  res[2]:=2;\n  res[3]:=3;\n  res[4]:=7
;\n  res[5]:=12;\n  res[6]:=21;\n  res[7]:=38;\n  k:=7;\n  A:=\{1,4,9
\};\n  AA:=\{11,14,19,41,44,49,91,94,99\};\n  oldQ:=A;\n  newQ:=AA;\n \+
 newL:=[1,2,3,7,12,21,29,38,43,47,48,49,51,52,53,\n          57,62,71,
79,88,93,97,98,99];\n  for slicesize from 2 by 2 while k<N do\n    use
rinfo(2,procname,\"slicesize=\",slicesize);\n    oldL:=newL;\n    user
info(3,procname,\"oldL=\",oldL);\n\n    is149:=subs(_slicesize=slicesi
ze,\n                _oldQ=oldQ,_newQ=newQ,\n    proc(n)\n      local \+
r,s;\n      option remember;\n      r:=n;\n      while r<>0 do\n      \+
  s:=irem(r,10^_slicesize,'r');\n        if not (member(s,_newQ)\n    \+
          or \n              (member(s,_oldQ) and r=0)\n             )
 then\n        RETURN(false)\n        fi\n      od; \n      true\n    \+
end);\n    for qq to 99 while k<N do\n      q:=qq*10^slicesize;\n     \+
 l:=0;\n      for r in oldL do\n        n:=q+r;\n        n2:=n^2;\n   \+
     tail:=irem(n2,10^(slicesize+2));\n        if is149(tail) then\n  \+
        l:=l+1;\n          newT[qq,l]:=n;\n           if is149(n2) the
n\n            k:=k+1;\n            res[k]:=n;\n            userinfo(2
,procname,k,res[k]);\n            userinfo(2,procname,\"time=\",time()
-st);\n            if k=N then break fi;\n          fi\n        fi\n  \+
    od\n    od;\n    newL:=sort(convert(newT,list));\n    oldQ:=oldQ \+
\n          union \{seq(seq(10*i+j,i=newQ),j=A)\}\n          union new
Q;\n    newQ:=\{seq(seq(100*i+j,i=newQ),j=AA)\};\n  od;\n  [seq(res[i]
,i=1..N)] \nend:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "infolev
el[tricol]:=2:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "start:=ti
me():\ntricol(12);\ntime()-start;" }}{PARA 6 "" 1 "" {TEXT -1 26 "tric
ol:   \"slicesize=\"   2" }}{PARA 6 "" 1 "" {TEXT -1 17 "tricol:   8  \+
 107" }}{PARA 6 "" 1 "" {TEXT -1 26 "tricol:   \"time=\"   .14e-1" }}
{PARA 6 "" 1 "" {TEXT -1 17 "tricol:   9   212" }}{PARA 6 "" 1 "" 
{TEXT -1 26 "tricol:   \"time=\"   .49e-1" }}{PARA 6 "" 1 "" {TEXT -1 
26 "tricol:   \"slicesize=\"   4" }}{PARA 6 "" 1 "" {TEXT -1 20 "trico
l:   10   31488" }}{PARA 6 "" 1 "" {TEXT -1 25 "tricol:   \"time=\"   \+
1.895" }}{PARA 6 "" 1 "" {TEXT -1 20 "tricol:   11   70107" }}{PARA 6 
"" 1 "" {TEXT -1 25 "tricol:   \"time=\"   2.596" }}{PARA 6 "" 1 "" 
{TEXT -1 21 "tricol:   12   387288" }}{PARA 6 "" 1 "" {TEXT -1 25 "tri
col:   \"time=\"   6.234" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#7.\"\"\"\"
\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J\"&2,(\"')G(Q" }}{PARA 11 "
" 1 "" {XPPMATH 20 "6#$\"%`j!\"$" }}}{EXCHG {PARA 0 "> " 0 "" 
{MPLTEXT 1 0 40 "start:=time():\ntricol(14);\ntime()-start;" }}{PARA 
6 "" 1 "" {TEXT -1 26 "tricol:   \"slicesize=\"   2" }}{PARA 6 "" 1 "
" {TEXT -1 17 "tricol:   8   107" }}{PARA 6 "" 1 "" {TEXT -1 26 "trico
l:   \"time=\"   .13e-1" }}{PARA 6 "" 1 "" {TEXT -1 17 "tricol:   9   \+
212" }}{PARA 6 "" 1 "" {TEXT -1 26 "tricol:   \"time=\"   .44e-1" }}
{PARA 6 "" 1 "" {TEXT -1 26 "tricol:   \"slicesize=\"   4" }}{PARA 6 "
" 1 "" {TEXT -1 20 "tricol:   10   31488" }}{PARA 6 "" 1 "" {TEXT -1 
25 "tricol:   \"time=\"   1.645" }}{PARA 6 "" 1 "" {TEXT -1 20 "tricol
:   11   70107" }}{PARA 6 "" 1 "" {TEXT -1 25 "tricol:   \"time=\"   2
.258" }}{PARA 6 "" 1 "" {TEXT -1 21 "tricol:   12   387288" }}{PARA 6 
"" 1 "" {TEXT -1 25 "tricol:   \"time=\"   5.447" }}{PARA 6 "" 1 "" 
{TEXT -1 26 "tricol:   \"slicesize=\"   6" }}{PARA 6 "" 1 "" {TEXT -1 
23 "tricol:   13   95610729" }}{PARA 6 "" 1 "" {TEXT -1 27 "tricol:   \+
\"time=\"   107.229" }}{PARA 6 "" 1 "" {TEXT -1 26 "tricol:   \"slices
ize=\"   8" }}{PARA 6 "" 1 "" {TEXT -1 24 "tricol:   14   446653271" }
}{PARA 6 "" 1 "" {TEXT -1 27 "tricol:   \"time=\"   681.336" }}{PARA 
11 "" 1 "" {XPPMATH 20 "6#70\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$
7#\"&)[J\"&2,(\"')G(Q\")H2h&*\"*rKlY%" }}{PARA 11 "" 1 "" {XPPMATH 20 
"6#$\"'g[p!\"$" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 21 "infolevel
[tricol]:=1:" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 40 "start:=time
():\ntricol(16);\ntime()-start;" }}}{PARA 12 "" 1 "" {XPPMATH 20 "6#72
\"\"\"\"\"#\"\"$\"\"(\"#7\"#@\"#Q\"$2\"\"$7#\"&)[J\"&2,(\"')G(Q\")H2h&
*\"*rKlY%\"+2rr[J\",z%z53@" }}{PARA 11 "" 1 "" {XPPMATH 20 "6#$\")+&H!
>!\"$" }}}{MARK "40 0 0" 40 }{VIEWOPTS 1 1 0 1 1 1803 }