This the guessing part of the Maple files that accompany the paper "Walks avoiding a quadrant and the reflection principl" by Mireille Bousquet-M\303\251lou and Michael Wallner. ############################################################# Abstract: We continue the enumeration of plane lattice walks with small steps avoiding the negative quadrant, initiated by the first author in 2016. We solve in detail a new case, namely the king model where all 8 nearest neighbour steps are allowed. The associated generating function is proved to be the sum of a simple, explicit D-finite series (related to the number of walks confined to the first quadrant), and an algebraic one. This was already the case for the two models solved by the first author in 2016. The principle of the approach is also the same, but challenging theoretical and computational difficulties arise as we now handle algebraic series of larger degree.
We expect a similar algebraicity phenomenon to hold for the 7 Weyl step sets, which are those for which walks confined to the first quadrant can be counted using the reflection principle. With this paper, this is now proved for 3 of them. For the remaining 4, we predict the D-finite part of the solution, and in 3 of the 4 cases, give evidence for the algebraicity of the remaining part. ############################################################# restart;
with(gfun): with(plots): gfun:-version;
with(LinearAlgebra):
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SShsaWJuYW1lR0YoNidRMEQ6XGxpYlxnZnVuLm1sYUYoUSdEOlxsaWJGKFFCL2hvbWUvbWljaGFlbHdhbGxuZXIvbGliL2dmdW4ubWxhRihROS9ob21lL21pY2hhZWx3YWxsbmVyL2xpYkYoUUBDOlxQcm9ncmFtfkZpbGVzXE1hcGxlfjIwMjBcbGliRig3J0YvRjBGMUYyRjM= JCIkJVEhIiM=
parseData := proc(data::list)
local vec,i,n;
n := nops(data):
vec := Vector(n):
for i from 1 to n do
vec[i] := op(i,data);
end:
return vec;
end:
loadData := proc(filename)
local fd,data;
#read file i
fd := fopen(filename, READ, TEXT);
data := readdata(fd,integer,1);
fclose(fd);
return parseData(data);
end:
computeParams := proc(algeq)
local degF, degC, numUndet, nn;
degF := degree(op(1,algeq),F(t));
degC := degree(subs(F(t)=F,op(1,algeq)),t);
numUndet := (degF+1)*(degC+1);
nn := nops(expand(op(1,algeq)));
return [degF,degC,numUndet,nn];
end:
infolevel[gfun] := 1;
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+JkkqaW5mb2xldmVsR0YnNiNJJWdmdW5HRiciIiI3I0Yx
fprefix := "data_King_A/":
<Text-field style="Heading 1" layout="Heading 1">R0</Text-field>
Recall that R(x) = t^2*M(x,0).
Hence, R0 = R(0) = t^2*M(0,0) = t^2*A_{-1,0}, and we will guess the GF for walks ending at (-1,0).
filename := "king-(-1,0)-2000.txt";
CC := loadData(filename):
seq(CC[i],i=1..10);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSlmaWxlbmFtZUdGKFE1a2luZy0oLTEsMCktMjAwMC50eHRGKDcjRi4= NiwiIiEiIiQiIiciI14iJFMjIiUzOyIlLSgqIiZtWCciJythVSIocjMhSA==
Divide by 3 CC := CC/3:
seq(CC[i],i=1..10);
NiwiIiEiIiIiIiMiIzwiIyEpIiRPJiIlTUsiJkE6IyInKz05IidkcCcq
Guess the algebraic relation.
As we are interested in R0 = t^2*A_{-1,0} we append to initial 0s. algeq := listtoalgeq([0,0,seq(CC[n],n=1..Dimension(CC))],F(t));
listtoalgeq: looking for equation of order 1
listtoalgeq: looking for equation of order 2
listtoalgeq: looking for equation of order 3
listtoalgeq: looking for equation of order 4
listtoalgeq: looking for equation of order 5
listtoalgeq: looking for equation of order 6
listtoalgeq: looking for equation of order 7
listtoalgeq: looking for equation of order 8
listtoalgeq: looking for equation of order 9
listtoalgeq: looking for equation of order 10
listtoalgeq: looking for equation of order 11
listtoalgeq: looking for equation of order 12
listtoalgeq: looking for equation of order 13
listtoalgeq: looking for equation of order 14
listtoalgeq: looking for equation of order 15
listtoalgeq: looking for equation of order 16
listtoalgeq: looking for equation of order 17
listtoalgeq: looking for equation of order 18
listtoalgeq: looking for equation of order 19
listtoalgeq: looking for equation of order 20
listtoalgeq: looking for equation of order 21
listtoalgeq: looking for equation of order 22
listtoalgeq: looking for equation of order 23
listtoalgeq: looking for equation of order 24
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZhbGdlcUdGKDckLGJvKiYsPComIilza1AlKSIiIilJInRHRigiI05GNEY0KiYiKEc0dyJGNClGNiIjTUY0RjQqJiIrNSgpSGJaRjQpRjYiI0xGNEY0KiYiLE1oXzQhPkY0KUY2IiNLRjRGNComIis1PyRcKCkpRjQpRjYiI0pGNEY0KiYiLG9Bdm56IkY0KUY2IiNJRjQhIiIqJiIrJW9BZyNSRjQpRjYiI0hGNEY0KiYiKyEpcDdUPEY0KUY2IiNHRjRGNComIionR0pSTUY0KUY2IiNGRjRGTComIilpViU+JkY0KUY2IiNFRjRGTComIilZOFk6RjQpRjYiI0RGNEY0KiYiKDtgOSJGNClGNiIjQ0Y0RkwqJiImT3UjRjQpRjYiI0JGNEY0RjQtSSJGR0YoNiNGNkY0RjQqJiw8KiYiKltoKjQhKUY0RjpGNEY0KiYiLCkzI1tgIj1GNEY+RjRGTComIixjTVtcYCpGNEZCRjRGNComIi13YSZSQHglRjRGRkY0RjQqJiItcV5CJkgpPUY0RkpGNEY0KiYiLVMjUkdpLSVGNEZPRjRGTComIiwzKXlbSyUpRjRGU0Y0RjQqJiIsTz4pKVFNJUY0RldGNEY0KiYiK0N4LGsqKUY0RmVuRjRGTComIispKT12QTVGNEZpbkY0RkwqJiIqbTwsQyRGNEZdb0Y0RjQqJiIpXydlRSNGNEZhb0Y0RkwqJiInI28iXEY0KUY2IiNBRjRGNEY0KUZjbyIiI0Y0RjQqJiw8KiYiLCVRPXdwP0Y0Rj5GNEZMKiYiLUdEOnZlP0Y0RkJGNEZMKiYiLjldJWVkKTQjRjRGRkY0RjQqJiIuYTlyNlEyKEY0RkpGNEY0KiYiLjRrYGcpUjxGNEZPRjRGNComIi5TeS9AVidcRjRGU0Y0RkwqJiIuJSkpPSF5Mi4iRjRGV0Y0RjQqJiItI29oKy8wJ0Y0RmVuRjRGNComIi03ZjZsJD4iRjRGaW5GNEZMKiYiLC1DQTM7IkY0Rl1vRjRGTComIit0WSY9XyRGNEZhb0Y0RjQqJiIqQXI6QSNGNEZicUY0RkwqJiIoXnI/JUY0KUY2IiNARjRGNEY0KUZjbyIiJEY0RjQqJiw8KiYiLXElPShcMklGNEZCRjRGTComIi5XU21oZ0ciRjRGRkY0RjQqJiIvJ2Y4Vk15QSRGNEZKRjRGNComIi8tVEJHOWhxRjRGT0Y0RjQqJiIudnc5aGcpeUY0RlNGNEY0KiYiL0toekckPTclRjRGV0Y0RkwqJiIuISlcYSRHbyEqRjRGZW5GNEY0KiYiLk95T0pTVCZGNEZpbkY0RjQqJiItW0IhKnAxJSpGNEZdb0Y0RkwqJiItY1cjKjRWNUY0RmFvRjRGTComIixsR1Q+ZSNGNEZicUY0RjQqJiIrd1ZgNDlGNEZic0Y0RkwqJiIpKCllakFGNClGNiIjP0Y0RjRGNClGY28iIiVGNEY0KiYsPComIi0ncHMmKnAoPUY0RkZGNEY0KiYiLylvPzRUJEdeRjRGSkY0RjQqJiIwO3pxbSl6SUxGNEZPRjRGNComIjBnKXpsb29eXkY0RlNGNEY0KiYiLil6UGkiNCNSRjRGV0Y0RjQqJiIwM010Uyw9XyNGNEZlbkY0RkwqJiIvVSZISHhEOidGNEZpbkY0RjQqJiIvJXA4WyFHMk5GNEZdb0Y0RjQqJiIuNzIzM3QjXEY0RmFvRjRGTComIi0hSHR5L1YoRjRGYnFGNEZMKiYiLV0qKilRZlIiRjRGYnNGNEY0KiYiK007bipSJ0Y0RmJ1RjRGTComIilpKSlIJilGNClGNiIjPkY0RjRGNClGY28iIiZGNEY0KiYsPComIi93JnkjKXBOayRGNEZKRjRGNComIjB3JlFvLCkqeWdGNEZPRjRGNComIjFNLm0+bDBcQ0Y0RlNGNEY0KiYiMVViJ3o3PmMqR0Y0RldGNEY0KiYiMFdIQ3Y0Syk9RjRGZW5GNEZMKiYiMV9bNnQ7YSE+IkY0RmluRjRGTComIjBDUyp5enk8TEY0Rl1vRjRGNComIjA7bldIOEh3IkY0RmFvRjRGNComIi9BbGo5RiN5IkY0RmJxRjRGTComIi5FNCQqPTg2JUY0RmJzRjRGTComIi1FamkxLmVGNEZidUY0RjQqJiIsQTk8bD0jRjRGYndGNEZMKiYiKip5ciJRI0Y0KUY2IiM9RjRGNEY0KUZjbyIiJ0Y0RjQqJiw8KiYiMEMvIlEvYCRHJUY0Rk9GNEY0KiYiMV82OS91I0hXJUY0RlNGNEY0KiYiMk0qR1hhd0FfOEY0RldGNEY0KiYiMil6XipIKSlIJSpIIkY0RmVuRjRGNComIjFqbl04JiopcFIiRjRGaW5GNEZMKiYiMTdIblUmbzxXJUY0Rl1vRjRGTComIjFXPlInUXctWCJGNEZhb0Y0RjQqJiIwUS1cXioqNDooRjRGYnFGNEY0KiYiLyEqUTNITS1VRjRGYnNGNEZMKiYiL2FkLHlHdTxGNEZidUY0RkwqJiIuJCo9JEcjKioqPUY0RmJ3RjRGNComIixTeHpzeSZGNEZieUY0RkwqJiIqSEorNSZGNClGNiIjPEY0RjRGNClGY28iIihGNEY0KiYsPComIjFWUyMzOmBFIUhGNEZTRjRGNComIjJBRlR2ajZ6SiNGNEZXRjRGNComIjIuZHNQJylbOiNlRjRGZW5GNEY0KiYiMmlRRDpHInl2WkY0RmluRjRGNComIjEzKT4+MDssImVGNEZdb0Y0RkwqJiIybSJ6Lyw3PEI4RjRGYW9GNEZMKiYiMVZjTWQheixBJkY0RmJxRjRGNComIjFJKillN0JTJVIjRjRGYnNGNEY0KiYiL1lYLU8sW1FGNEZidUY0RkwqJiIvRVQvQk5qZ0Y0RmJ3RjRGTComIi5wdzR0eChcRjRGYnlGNEY0KiYiLWV4OiMpMzdGNEZiW2xGNEZMKiYiKicqKlJ4JilGNClGNiIjO0Y0RjRGNClGY28iIilGNEY0KiYsPComIjJHYmwsc3pnUCJGNEZXRjRGNComIjJrcyJvT0lDLyMqRjRGZW5GNEY0KiYiMz80Ij50OTl1KyNGNEZpbkY0RjQqJiIzPyRvclE1UFdZIkY0Rl1vRjRGNComIjJuY0liXlBcaCJGNEZhb0Y0RkwqJiIyJUc1MSRSNE45JEY0RmJxRjRGTComIjJbV3A4NT90YyJGNEZic0Y0RjQqJiIxaXRVeDk+M25GNEZidUY0RjQqJiIwUUpPJHpaJHAiRjRGYndGNEY0KiYiMD8kemAqKUhuO0Y0RmJ5RjRGTComIi8ueltbaWI1RjRGYltsRjRGNComIi1zISp5bDw/RjRGYl1sRjRGTComIitdXW1hNkY0KUY2IiM6RjRGNEY0KUZjbyIiKkY0RjQqJiw8KiYiMidwWCY9XCZlO1xGNEZlbkY0RjQqJiIzdypbenJGR2IpR0Y0RmluRjRGNComIjM1QV83YFplYGNGNEZdb0Y0RjQqJiIzbTVFVSg+Jz4neiRGNEZhb0Y0RjQqJiIyeFdpXjcmXFZIRjRGYnFGNEZMKiYiMmdwZGdvU1smZUY0RmJzRjRGTComIjJrIlFdP0FubFJGNEZidUY0RjQqJiIydSJHKmYhKj1oZSJGNEZid0Y0RjQqJiIxYSlHViYpSCE+NUY0RmJ5RjRGNComIjBRLl1eJEhSUEY0RmJbbEY0RkwqJiIvbFcoKkhbRT1GNEZiXWxGNEY0KiYiLTFTPTE5RkY0RmJfbEY0RkwqJiIrXSV6XkUiRjQpRjYiIzlGNEY0RjQpRmNvIiM1RjRGNComLDwqJiIzS1VOPS1oaXg4RjRGaW5GNEY0KiYiM2NpNjdyRkM3dEY0Rl1vRjRGNComIjRxbXJQNUFAIj44RjRGYW9GNEY0KiYiM002c111LUAiUilGNEZicUY0RjQqJiIyenc0a184TlUjRjRGYnNGNEZMKiYiMm8yLitmdTY1KUY0RmJ1RjRGTComIjIrPSZSI1JoMl8pRjRGYndGNEY0KiYiMi1Id2tzJFIkPSRGNEZieUY0RjQqJiIxXUIpRyg0IUg/JEY0RmJbbEY0RjQqJiIwVSUqZTI/JTRwRjRGYl1sRjRGTComIi9eKlsleUohZiNGNEZiX2xGNEY0KiYiLSlSX2siZkhGNEZiYWxGNEZMKiYiK2NBblg2RjQpRjYiIzhGNEY0RjQpRmNvIiM2RjRGNComLDwqJiIzISpvMHopenRINSRGNEZdb0Y0RjQqJiI0IVFLWlN6RjJBOkY0RmFvRjRGNComIjQnbyhRTTo6cGZkI0Y0RmJxRjRGNComIjQpR2JaVDtKPiplIkY0RmJzRjRGNComIjJsI1tZWCpSKj1aRjRGYnVGNEY0KiYiMicqeTpCXl5PIW9GNEZid0Y0RkwqJiIzUXQrI3A9JHBpOkY0RmJ5RjRGNComIjItPEwjcFM5V2FGNEZiW2xGNEY0KiYiMWVCYG1rdzx0RjRGYl1sRjRGNComIjElUVcpcFonKWU1RjRGYl9sRjRGTComIi8sOEdzVjhJRjRGYmFsRjRGNComIi0jSCkpXHJpI0Y0RmJjbEY0RkwqJiIqVE1KbylGNClGNiIjN0Y0RjRGNClGY29GY2VsRjRGNComLDwqJiIzR0hVdCwkNDJyJkY0RmFvRjRGNComIjRXVGEhZU1oMEpFRjRGYnFGNEY0KiYiNE01KVtbRSplKlFVRjRGYnNGNEY0KiYiNGVlI29bUihRSGUjRjRGYnVGNEY0KiYiM1JaNUlPdSczUiNGNEZid0Y0RjQqJiIyS1NYSj04Kj46RjRGYnlGNEY0KiYiM1dRbnZTeDNgQ0Y0RmJbbEY0RjQqJiIyPSRRYy4xJzMmekY0RmJdbEY0RjQqJiIyVWhARnYtJT44RjRGYl9sRjRGNComIjFhIkc0elEoXDhGNEZiYWxGNEZMKiYiL3hcQiopNGtHRjRGYmNsRjRGNComIi0lSHk+NCI+RjRGYmVsRjRGTComIiphWHhiJkY0KUY2RmVjbEY0RjRGNClGY29GY2NsRjRGNComLDwqJiIzL2tBX1tyWHonKUY0RmJxRjRGNComIjQpMzFSK2pMSS5RRjRGYnNGNEY0KiYiNCdlZGlcJHkhKjQhZkY0RmJ1RjRGNComIjRtKCkzeSpmKmYlKmYkRjRGYndGNEY0KiYiM1hpSzFUeC8nZSZGNEZieUY0RjQqJiIzU2M0QDF0Umc8RjRGYltsRjRGNComIjN3KWV3JzNTKXlIJEY0RmJdbEY0RjQqJiIyNWQsIypvJlJBKipGNEZiX2xGNEY0KiYiMkVGZWZxZTElPkY0RmJhbEY0RjQqJiIxeSZwImY0Q0g5RjRGYmNsRjRGTComIi8uKFskMzkrQUY0RmJlbEY0RjQqJiItWU1bTV82RjRGYWdsRjRGTComIiolKj1ALSRGNClGNkZlYWxGNEY0RjQpRmNvRmNhbEY0RjQqJiw8KiYiNG9wQG1dcktoNCJGNEZic0Y0RjQqJiI0PyZcZWNXV284WUY0RmJ1RjRGNComIjQzb0o0Lyo0cmNwRjRGYndGNEY0KiYiNCdIcCNbWSxCdUclRjRGYnlGNEY0KiYiM29DSHRPSmgjPipGNEZiW2xGNEY0KiYiM1d4O3RWKCkqPnAkRjRGYl1sRjRGNComIjMlZW10QjFQMXokRjRGYl9sRjRGNComIjNrbTl5SGImcDAiRjRGYmFsRjRGNComIjJDNzZxUU5ZTyNGNEZiY2xGNEY0KiYiMSlHKVFXeicpWzdGNEZiZWxGNEZMKiYiLylbZjJLcEwiRjRGYWdsRjRGNComIix3PVlBKWVGNEZfaWxGNEZMKiYiKlVTO1MiRjQpRjZGZV9sRjRGNEY0KUZjb0ZjX2xGNEY0KiYsPComIjQmeW5NIj4yNUY6IkY0RmJ1RjRGNComIjQ1WmJsZyhHRShwJUY0RmJ3RjRGNComIjRaalliaSMzdktwRjRGYnlGNEY0KiYiNFNfUjRLMzVATSVGNEZiW2xGNEY0KiYiNEB1NzlYX21dOyJGNEZiXWxGNEY0KiYiM3knPSdleiNRQjYmRjRGYl9sRjRGNComIjMsP3A6IWUxK3IkRjRGYmFsRjRGNComIjIvaVxdPGMsZSpGNEZiY2xGNEY0KiYiMjtpaHMlcGcsQ0Y0RmJlbEY0RjQqJiIwP0lAN1JdKCkpRjRGYWdsRjRGTComIi5dYjBsRzsnRjRGX2lsRjRGNComIiw3X0pLZyNGNEZdW21GNEZMKiYiKWojKVxiRjQpRjZGZV1sRjRGNEY0KUZjb0ZjXWxGNEY0KiYsPComIjQpbyVvczwiXFEyNUY0RmJ3RjRGNComIjRzVSpldFNbcCoqUkY0RmJ5RjRGNComIjRrSjMkSHhiaDVlRjRGYltsRjRGNComIjRrUkArNyJbJSk0UEY0RmJdbEY0RjQqJiI0SiIzbUAqeV9BPCJGNEZiX2xGNEY0KiYiM0N6NC4kXFU3TiZGNEZiYWxGNEY0KiYiMyF6MDwnXCk9MTIkRjRGYmNsRjRGNComIjIrdSp5V1MmcE0oRjRGYmVsRjRGNComIjJTSEsvWHFGLiNGNEZhZ2xGNEY0KiYiMCN5Jz4hb2kmKVxGNEZfaWxGNEZMKiYiLiRcZClvcCc+RjRGXVttRjRGNComIixzNEY8KyJGNEZbXW1GNEZMKiYiKTVCdj1GNClGNkZlW2xGNEY0RjQpRmNvRmNbbEY0RjQqJiw8KiYiMyVlQGBcMiFlbnNGNEZieUY0RjQqJiI0KW9hKnklPTM3RUdGNEZiW2xGNEY0KiYiNF9Rd3ltWispZVNGNEZiXWxGNEY0KiYiNEdbJVFvMzI7V0VGNEZiX2xGNEY0KiYiM3ojZlk1ISpRP1cqRjRGYmFsRjRGNComIjNrIm8ibzxbYiRSJUY0RmJjbEY0RjQqJiIzK0MiKnlvemVFQEY0RmJlbEY0RjQqJiIyJTRUYm8uYUNaRjRGYWdsRjRGNComIjJpO3UkW19fRDlGNEZfaWxGNEY0KiYiMFdlSFwqSCYzI0Y0Rl1bbUY0RkwqJiItKkgibyd5RiRGNEZbXW1GNEY0KiYiK2FHSzJLRjRGaV5tRjRGTComIihEVFMmRjQpRjZGZXlGNEY0RjQpRmNvRmN5RjRGNComLDwqJiIzX1ooeSMqKTMoUUYlRjRGYltsRjRGNComIjRbRyNlW1YpPltqIkY0RmJdbEY0RjQqJiI0cyN5KEhdIjR0R0JGNEZiX2xGNEY0KiYiNF8+c0dNeGFwYSJGNEZiYWxGNEY0KiYiMy9ARTJjaVhlZ0Y0RmJjbEY0RjQqJiIzbyR5NFlPJT1ZR0Y0RmJlbEY0RjQqJiIzT1NIbV5DQjg3RjRGYWdsRjRGNComIjJPIylvIlEpW0heI0Y0Rl9pbEY0RjQqJiIxOzIjKj0zWCI9KUY0Rl1bbUY0RjQqJiIvJz4ucj1KXCZGNEZbXW1GNEZMKiYiLHE4MWZ0IkY0RmlebUY0RkwqJiIqPVcjeXhGNEZnYG1GNEZMKiYiKGpUSiJGNClGNkZld0Y0RjRGNClGY29GY3dGNEY0KiYsPComIjNvdmgpeUlmZysjRjRGYl1sRjRGNComIjMpb3lqd1xiZWQoRjRGYl9sRjRGNComIjQ3PG16IlxDKkgyIkY0RmJhbEY0RjQqJiIzP3prXXM3KFJEKEY0RmJjbEY0RjQqJiIzV0U2V2coXFkvJEY0RmJlbEY0RjQqJiIzS1dmMTcpcHJWIkY0RmFnbEY0RjQqJiIyLTl4IjRLISlwYkY0Rl9pbEY0RjQqJiIydTI7Qz1zRTMiRjRGXVttRjRGNComIjF1SVQvQFVqUEY0RltdbUY0RjQqJiIuIWZBYEM/PEY0RmlebUY0RkwqJiIrIVF2KWVoRjRGZ2BtRjRGTComIipXN2w4IkY0RmVibUY0RkwqJiInVDtGRjQpRjZGZXVGNEY0RjQpRmNvRmN1RjRGNComLDwqJiIyb2ZKKFFOZmdzRjRGYl9sRjRGNComIjMzPUZYcEp6OkZGNEZiYWxGNEY0KiYiMytPcmxceF5KUUY0RmJjbEY0RjQqJiIzW0MyMiIpXDtHRUY0RmJlbEY0RjQqJiIzQ2kpPiNHJXkqZjZGNEZhZ2xGNEY0KiYiMiMqNGFBeTskKVsmRjRGX2lsRjRGNComIjI7K1koeVBRJSk+RjRGXVttRjRGNComIjEpb2tVVSFIYk9GNEZbXW1GNEY0KiYiMUs1ImZVeDFNIkY0RmlebUY0RjQqJiIuOzckKkhVeSZGNEZnYG1GNEY0KiYiLDEnekIyOEY0RmVibUY0RjQqJiIodUBtJUY0RmNkbUY0RjQqJiImQl0lRjQpRjZGZXNGNEY0RjQpRmNvRmNzRjRGNComLDwqJiIya0M5dTMrLCI+RjRGYmFsRjRGNComIjJfRjcjb2x5KDQoRjRGYmNsRjRGNComIjJ3WDltclkmKSkqKkY0RmJlbEY0RjQqJiIyJVFHQjdjdEJwRjRGYWdsRjRGNComIjJPZkllJ3AoKWZKRjRGX2lsRjRGNComIjJjI29YbnQhb1wiRjRGXVttRjRGNComIjFXXHc+dS1tXkY0RltdbUY0RjQqJiIwJzRxJSoqUjM5KkY0RmlebUY0RjQqJiIwZ0pgJzN2KFskRjRGZ2BtRjRGNComIi4pM140MGxFRjRGZWJtRjRGNComIittLyFbcShGNEZjZG1GNEY0KiYiKCFILHNGNEZhZm1GNEY0KiYiJVpqRjQpRjZGZXFGNEY0RjQpRmNvRmNxRjRGNComLDwqJiIxc2lrLSoqKTRFJEY0RmJjbEY0RjQqJiIyR0hUWyhbYDI3RjRGYmVsRjRGNComIjJzdUhLXVpxcCJGNEZhZ2xGNEY0KiYiMl9OeWhtNkg9IkY0Rl9pbEY0RjQqJiIxV2hDPS5BI1wmRjRGXVttRjRGNComIjFbYzZweDguRUY0RltdbUY0RjQqJiIwTzI9eVcweClGNEZpXm1GNEY0KiYiME9aM20qcDw6RjRGZ2BtRjRGNComIi9jJT5PMyE+ZkY0RmVibUY0RjQqJiItLyZwXzBlJkY0RmNkbUY0RjQqJiIrUzUkKipwIkY0RmFmbUY0RjQqJiIoSyFvPEY0Rl9obUY0RjQqJiIkKWVGNEY2RjRGNEY0KUZjb0Zib0Y0RjQqJiw8KiYiMGNRYmUiXDxGRjRGYmVsRjRGNComIjFXeE9pIXppKyJGNEZhZ2xGNEY0KiYiMWM5cF9pPzk5RjRGX2lsRjRGNComIjAnSDo9c2ZkKSpGNEZdW21GNEY0KiYiMDcwX0Vdb2QlRjRGW11tRjRGNComIjAvalVaIkdwQEY0RmlebUY0RjQqJiIvRzw9dHkzdEY0RmdgbUY0RjQqJiIvR1A8KFxaRSJGNEZlYm1GNEY0KiYiLilHb3BdS1xGNEZjZG1GNEY0KiYiLCN6MFldWUY0RmFmbUY0RjQqJiIqPzRtVCJGNEZfaG1GNEY0KiYiJ090OUY0RjZGNEY0IiNcRjRGNClGY29GXm9GNEY0KiYiKE4qZUdGNClGNiIjT0Y0RkwqJiIpLSwpUSRGNEY1RjRGNComIio0U2tDI0Y0RjpGNEY0KiYiKmk/a2EkRjRGPkY0RjQqJiIpOkAkZiNGNEZCRjRGNComIippQVQ0I0Y0RkZGNEZMKiYiKChbR2NGNEZKRjRGNComIiknMyopKmVGNEZPRjRGNComIikzKUdUIkY0RlNGNEZMKiYiKFcvNCJGNEZXRjRGNComRmBvRjRGZW5GNEZMSSRvZ2ZHRig3I0Yu
computeParams(algeq);
NyYiI0MiI08iJEQqIiRCJA==
algR0 := primpart(subs(F(t)=R0,op(1,algeq)),R0):
algR0 := collect(algR0,R0,factor):
<Text-field style="Heading 1" layout="Heading 1">R1</Text-field>
Recall that R(x) = t^2*M(x,0).
Hence, R1 = R'(0) = t^2*M_x(0,0) = t^2*A_{-2,0}, and we will guess the GF for walks ending at (-2,0).
filename := "king-(-2,0)-2000.txt";
CC := loadData(filename):
seq(CC[i],i=1..10);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSlmaWxlbmFtZUdGKFE1a2luZy0oLTIsMCktMjAwMC50eHRGKDcjRi4= NiwiIiIiIiEiIioiI0MiJDcjIiUlMyIiJTV3IiZNdyUiJzpuSyIoazs/Iw==
Divide by 3 CC := CC/3:
seq(CC[i],i=1..10);
NiwjIiIiIiIkIiIhRiUiIikjIiQ3I0YlIyIlJTMiRiUjIiU1d0YlIiZ5ZSIiJzAqMyIiJykpUXQ=
Guess the algebraic relation.
As we are interested in R1 = t^2*A_{-2,0} we append to initial 0s. algeq := listtoalgeq([0,0,seq(CC[n],n=1..Dimension(CC))],F(t));
listtoalgeq: looking for equation of order 1
listtoalgeq: looking for equation of order 2
listtoalgeq: looking for equation of order 3
listtoalgeq: looking for equation of order 4
listtoalgeq: looking for equation of order 5
listtoalgeq: looking for equation of order 6
listtoalgeq: looking for equation of order 7
listtoalgeq: looking for equation of order 8
listtoalgeq: looking for equation of order 9
listtoalgeq: looking for equation of order 10
listtoalgeq: looking for equation of order 11
listtoalgeq: looking for equation of order 12
listtoalgeq: looking for equation of order 13
listtoalgeq: looking for equation of order 14
listtoalgeq: looking for equation of order 15
listtoalgeq: looking for equation of order 16
listtoalgeq: looking for equation of order 17
listtoalgeq: looking for equation of order 18
listtoalgeq: looking for equation of order 19
listtoalgeq: looking for equation of order 20
listtoalgeq: looking for equation of order 21
listtoalgeq: looking for equation of order 22
listtoalgeq: looking for equation of order 23
listtoalgeq: looking for equation of order 24
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZhbGdlcUdGKDckLGJyKiYiPTcoUlBbdGt2SWlRKFFuOSIiIilJInRHRigiIzpGMiEiIiomIjZBb3ZsbiJSKCoqKTNERjIpRjQiIihGMkYyKiYiP1dYOyg+VCJbMzhQQzpedlpGMilGNCIjQkYyRjIqJiI6OEJPaypcQChHXlA/JCoqRjIpRjQiIzdGMkY2KiYiLk9XNSNwQURGMilGNCIiI0YyRjIqJiI/Q3VHOENpX0hzODY0V1I5RjIpRjQiI0FGMkYyKiYiPmM1ZVQ5eW0tYjsqUnVIU0YyKUY0IiNARjJGNiomIjNdaDgnWyZ5MzUhKkYyKUY0IiImRjJGNiomIj5jQXIqKlJbKmUsciZ5NkZgKUYyKUY0IiM/RjJGNiomIj53NEIiUnhUWSFmTidcVD9sRjIpRjQiIz5GMkY2KiYiPjNZUydbbmledHJkOylSWyRGMilGNCIjPUYyRjYqJiI+d1BPNyY+PyxtKHBSTkNZIkYyKUY0IiM8RjJGNiomIjBTPXMqKVslcGZGMilGNCIiJEYyRjYqJiI9X1oyWCplJ29NKio9KFwsMCZGMilGNCIjO0YyRjYqJiI4J29QUlA8OVwyI1wpW0YyKUY0IiIqRjJGMiomIj0vKCkpKSlIVE5pWHF3UiRlbkYyKUY0IiNPRjJGNiomIj5vLlRLeCUqXDdfKT5edWIlKkYyKUY0IiNORjJGNiomIj8rVzQiUT4+Iyk+IlJTMCRwJGVGMilGNCIjTUYyRjYqJiJAS0deJDNQMz0hWy89Q1h6NCNGMilGNCIjTEYyRjYqJiJAP3Ikb2BlJT0zMVdGU1A0KlxGMilGNCIjS0YyRjYqJiJAV2REIWVrR1w+PEt2dGpnJSlGMilGNCIjSkYyRjYqJiJBKW9NO0RNYEd3LThcNzB5MSJGMilGNCIjSUYyRjYqJiJBbzJbS2koKjM/RGRhLTBSQzVGMilGNCIjSEYyRjYqJiJAWzNCJmVOZksjM2duaWlSVChGMilGNCIjR0YyRjYqJiJAR1BiU05fMUo4NW1vU0t4JEYyKUY0IiNGRjJGNiomIj8/amg5MCgqejRoSUcicF84KkYyKUY0IiNFRjJGNiomIjsrLyVbPCUzJUgoSDZUJTQoRjIpRjQiIzhGMkY2KiYiPFcoKillIlxCJ0hZUnhLKWUkRjIpRjQiIzlGMkY2KiYiNyo+OXZrZmAyYmQpZUYyKUY0IiIpRjJGMiomIj83ZDZPcjwlKWUiUihvMG1Pd0YyKUY0IiNDRjJGMiomIj8zRVRINlkiKWVpXjVPJ1w0JkYyKUY0IiNERjJGMiomLGRvKiYiPzszNFQ8cyQpZmsrLlItd1JGMkZfcEYyRjIqJiJAJT1YIWVoUyZldGVQV0s+Jz4mRjJGY3BGMkYyKiYiQScqW3VYJCpHdFZDNFonUTA1KEhGMkZncEYyRjIqJiJBTz5SI1w4YSNwZHAmNC55SCkpKkYyRltxRjJGMiomIkI7bzlxKilmNSRmV2NHUyFmJWZARjJGX3FGMkYyKiYiQj9uT01KNiY9T2skPUFLPyYpSCRGMkZjcUYyRjIqJiJCT3ojZkcyIVwpei5fWnIxVTZPRjJGZ3FGMkYyKiYiQltTQkYmUT1xSShmdG1SP3B4I0YyRltyRjJGMiomIkIzdXB6bWliYCpvTXZPcSJ6RyJGMkZfckYyRjIqJiJAWzsoXF4hZTNjbztDdC0kKT0jRjJGY3JGMkY2KiYiQUsjUklHIls5QnZAI2ZXVi5iJ0YyRmdzRjJGNiomIkFzXShROVVbKSk0JjMoZmZYZ2UnRjJGY3NGMkY2KiYiQXNJT1opKXArUGgnKnAtIil5Km8kRjJGPUYyRjYqJiJAMyUqM2YjeU0lKT5bRyd6IykpMygpRjJGSUYyRjYqJiJAdyp5TzVtQmElSExkbW9FPidGMkZNRjJGMiomIkArW1tocFwiNCxTIilmMC8/JCpGMkZVRjJGMiomIkBvVkorPyVRS0hILVZHekRwRjJGWUYyRjIqJiJAb3ptKHluU0o7QSY+VCNRYFBGMkZnbkYyRjIqJiJAb1pZSWMmSE1mJik0T3Ylb2kiRjJGW29GMkYyKiYiPysjZj1GKz9BczBwJHlucWVGMkZjb0YyRjIqJiI/KycqWy4lPl4ib2ZEUDJeKnoiRjJGM0YyRjIqJiI+YzkqSERUZ0BDOjV3diRwJUYyRltzRjJGMiomIj4nNFtVO3Y3UFF5SiRcaCwiRjJGZ3JGMkYyKiYiPUMhKVFVKSlvSHdmMmlRKnAiRjJGQUYyRjIqJiI8ayc0Q1RycjN4TFVxbjxGMilGNCIjNkYyRjIqJiI6OzlGSSFSdHJyJj0pUUpGMilGNCIjNUYyRjYqJiI6T0dPJ0dDNkUmZXkrdyZGMkZnb0YyRjYqJiI6byVlX2I2JDNCRypcSThGMkZfc0YyRjYqJiI5Wy9TWCRIeW0xJ3BgPUYyRjlGMkY2KiYiOCV5JypcIykpKVIheiZ5VjtGMilGNCIiJ0YyRjYqJiI2YytpYCwjKm94JFFqRjJGUUYyRjYqJiI1KTNHMHRRTiMpb0wkRjIpRjQiIiVGMkYyKiYiNDdXM3JbeGIkW0ZGMkZfb0YyRjIqJiIzO3VkZCwob0tUIkYyRkVGMkY2KiYiMT9iO3BZJDN6IkYyRjRGMkYyIi4zTEp3IW92RjZGMi1JIkZHRig2I0Y0RjJGMiomLGJvKiYiQTdsNDcvTT0tc2RzdFl0TjVGMkZjcEYyRjYqJiJCQ2FAVSEpKXBudyEqXCQ9Z1ViQyJGMkZncEYyRjYqJiJCc1EjZUhuIipROVdbSz03ZjpsRjJGW3FGMkY2KiYiQ1MtdC0oR2pbdyU9NFx1QVVzPkYyRl9xRjJGNiomIkNPXkopZSR5TjdSZnBaUT9QaFFGMkZjcUYyRjYqJiJDQ2VbdG1rZms6OV9cMkBNJTQmRjJGZ3FGMkY2KiYiQydwdHA5OWEyTSgpKjN1Mi5uK1dGMkZbckYyRjYqJiJDR3BaPyUpeStkSXpDWFQzWTQ+RjJGX3JGMkY2KiYiQkdwZkg6XGg6cD5FKHBsLC0kKUYyRmNyRjJGMiomIkMjKipHeWtMJW8zb1pFVjR3RTRCRjJGZ3NGMkYyKiYiQ196Tj4uPiw8RU9AJUdSLzNAI0YyRmNzRjJGMiomIkM7R1dYQEAuQl0qM0cpR2pNTzdGMkY9RjJGMiomIkJXJjNwSlEtSnZcZWwnUVpVWiNGMkZJRjJGMiomIkJXKCpSIzM6YEhXUFUvIkhZKzYkRjJGTUYyRjYqJiJCIVtbWGAnKmZYPlsoKW8zbD8uV0YyRlVGMkY2KiYiQlNpPWFGIlsmRyV5MXNTYGg8TUYyRllGMkY2KiYiQks3WjxfVUBaJio9SkI0YUsoPkYyRmduRjJGNiomIkFzOVtEXDQwbVFGbSl6Ujk8KkYyRltvRjJGNiomIkFjJSoqZT5TLVVYamMneiopXF1ORjJGY29GMkY2KiYiQTsjKilbImZQTW4neldHKCoqSGs2RjJGM0YyRjYqJiJATy5sKFwqUmpfYCFRZ2pXV0tGMkZbc0YyRjYqJiI/a1tMSERmUW9tPHdhN2N2RjJGZ3JGMkY2KiYiP2MuczxkTEFBKSlvYSpROlIiRjJGQUYyRjYqJiI+V1ZzKj1HYEEoKkhOJHlRcCJGMkZdd0YyRjYqJiI7N1w9I3BzS2hJLyR6Yj1GMkZhd0YyRjYqJiI8ZXRLdyUpbyg0TWI7RVJpRjJGZ29GMkYyKiYiPG0uS2NJbThPNF4jNDA9RjJGX3NGMkYyKiYiO1dqc2s5Y1k7eCsmNCpIRjJGOUYyRjIqJiI6b0pjK0EqZWIzaU1zSEYyRlt4RjJGMiomIjl3VV9dYygqKlJgKipbRSJGMkZRRjJGMiomIjcnSCFHTUdtTlxfS1lGMkZheEYyRjYqJiI2XSgqPStFKik+QENAJUYyRl9vRjJGNiomIjVNNCYpcDhmc246PUYyRkVGMkYyKiYiM10oUTQiKSopZXokPUYyRjRGMkY2IjBRV2E4WCoqKWZGMkYyKUZqeEZGRjJGMiomLGBvKiYiQj8qSCx4a2cmPiUpKkdreXlBQTpGMkZncEYyRjIqJiJDU1VGMnAyLWd2alFrKUg0cmoiRjJGW3FGMkYyKiYiQ0MtIykpM0pfcGcxRCsnWycqKSpSdkYyRl9xRjJGMiomIkRLdyY+LSlccydvX0wiKXokeXl0JT5GMkZjcUYyRjIqJiJEXy51by5xNSd5bDM9LF1nXjVJRjJGZ3FGMkYyKiYiRClbJlx1KVx2LCE+IkhTYmUjXHhXI0YyRltyRjJGMiomIkMpKWVvJlJ3LjgkW1IwI3A3N2VaJEYyRl9yRjJGNiomIkQ/diMpZThvIylcOz5lUEtSWGchUkYyRmNyRjJGNiomIkRDITQnR0UyIylIIVEpbz4qZjEoZihlRjJGZ3NGMkY2KiYiRDMjeUc5eG56bXIiKVtVbTZJRV9GMkZjc0YyRjYqJiJELyZIMUkwSURKTUA2MGZOPG1HRjJGPUYyRjYqJiJDWz9jTWspM1E2RWI8a0JaSVQlRjJGSUYyRjYqJiJEJXk5RVclKj1KMFYkUmRNLzk4KyJGMkZNRjJGMiomIkRTUSV5Z0ghZkcoPW5WOTpdI2ZPIkYyRlVGMkYyKiYiRENNNTc4SGAwJD5jLCV5TjlmNCJGMkZZRjJGMiomIkNnaEJtSCpwOUxqOScpZlBUJ0htRjJGZ25GMkYyKiYiQ19Wb0plJCpIbTJcTiI0MCV6S0tGMkZbb0YyRjIqJiJDSztZKlFpP1RaTV5wQCJbOTA4RjJGY29GMkYyKiYiQktPRDAjKVFLYyNcJCo0ZEE1NldGMkYzRjJGMiomIkJLMyV5RXIkXCkpb2QhZjlVRlk3RjJGW3NGMkYyKiYiQV8+MExBdjVzJypSeCE+P1EoR0YyRmdyRjJGMiomIkB3OEFKO1ZEI1wyWjNzPWBcRjJGQUYyRjIqJiI/Q3kkKj4saF1HaSM+RkVbMiVGMkZdd0YyRjIqJiI/IW8taTQwN0JgIXlPWjZNNUYyRmF3RjJGNiomIj4hXDBaYylIJTMsNyk9SC9vJkYyRmdvRjJGNiomIj5yP3c7a3krUl1YKD51WzlGMkZfc0YyRjYqJiI9KXksZDEpUmEwXzdPYGpCRjJGOUYyRjYqJiI8JVsmZUZOPyM9aDp5OyZRI0YyRlt4RjJGNiomIjtxcUtfKXlbI2YzZl4sNkYyRlFGMkY2KiYiOStERmIvIVw/NUspKVsjRjJGYXhGMkYyKiYiOGk9LT9gbTdcWk0uJEYyRl9vRjJGMiomIjcsa1tVMXQhXClRdjVGMkZFRjJGNiomIjRnVHJyJlJbVFUnKUYyRjRGMkYyIjJUamsiWy5wd0FGNkYyKUZqeEZgb0YyRjIqJixebyomIkMjKiozSjoqNDEiZVROQ3ctb2FHIkYyRltxRjJGNiomIkQvNi4nKUhRQmkiRypwRmxFaCRRNkYyRl9xRjJGNiomIkQrN0RmJUgzYVxdbTNkZFE/UFJGMkZjcUYyRjYqJiJEIUc8Xyc0JHAlNHAqUlNUU05ILmNGMkZncUYyRjYqJiJEK1NJZTlmK11pckttTXdFbl0kRjJGW3JGMkYyKiYiRXMnM3BlIXpZNCcqXEA5TyUqM2YkMyRGMkZfckYyRjIqJiJFd1Q+QS9YK15lOU1VL0JXJW9tJ0YyRmNyRjJGMiomIkVHdCN5PGBPWGMlXCFHbEFmJSlIbClGMkZnc0YyRjIqJiJFQ0FscFg8I2VZI2ZmQkpmZGozdUYyRmNzRjJGMiomIkVXcCcpZXYtbEdQdi0+NShRJWUocCRGMkY9RjJGMiomIkRbU15PRiRlcFRNXGE/aWk7I0cjRjJGSUYyRjYqJiJFWzc5LFVSU0tTYV0qKm9vZD81RUYyRk1GMkY2KiYiRW8uRC9eOl8xK1QjKSpHeik9YUhKRjJGVUYyRjYqJiJFXyJIKikqZTNYdXJeUDNEQGsjM1sjRjJGWUYyRjYqJiJFQy1hakMnKTNFUmdaMTY9KTQlPTpGMkZnbkYyRjYqJiJEW3cxZjBgKFxFbTRpKnpmLXVeKEYyRltvRjJGNiomIkQvbjhuN2lFSVRyMGYmWz8vZElGMkZjb0YyRjYqJiJEKVtWei1MZzBIOCJ5cy56XjAtIkYyRjNGMkY2KiYiQ1NBYmVNLDooZU5jcUtSM0N1I0YyRltzRjJGNiomIkJLP0Y7Sjo7KjQ7XXhyU2M9YkYyRmdyRjJGNiomIkF5Ul94Uj4sKyQqUiZmJXpBZGRGMkZBRjJGNiomIkE7WnB6KDMoZiU0YjlTcyU9JEciRjJGXXdGMkYyKiYiQF0uKGUsLkwtVXZyejBGTSYqRjJGYXdGMkYyKiYiQEs9JSo9VUgsaEJSTixPVz4kRjJGZ29GMkYyKiYiPyRlJ1tqLXEkekRtTmAlPWhzRjJGX3NGMkYyKiYiPz8qPScqUSZvKD1jKWZ3QzlcNkYyRjlGMkYyKiYiPll3TVgnUilSXiZcNCU0NzwiRjJGW3hGMkYyKiYiPFEoKipROXklKWVvNSgqZm1lRjJGUUYyRjIqJiI6JW9ZLyY9LCJ6TnFEb2lGMkZheEYyRjYqJiI6TykpSD5RSCJHeCYpeVk4RjJGX29GMkY2KiYiOCI0QTRkI1I4SyRcOFJGMkZFRjJGMiomIjZhRFsuTnRPJCpmXCNGMkY0RjJGNiIzVm4ob2RObnZgJkYyRjIpRmp4RmJ4RjJGMiomLFxvKiYiQ0d0UVclbzArZTB0NyVHI2YiW1dGMkZfcUYyRjIqJiJEJWVCQnZZWCp6KzNKSlJNPSQ+OEYyRmNxRjJGMiomIkVvLiopRzZ3SE8icCZ6IipRdzJsUDZGMkZncUYyRjYqJiJFZ0RjcXA1SyY+KmZUYCY9aTBvdClGMkZbckYyRjYqJiJGKydcYSJleXBxJHB2SjldJEdgeiJHRjJGX3JGMkY2KiYiRmtfbiw/ajc/SS0ncHBbJ29WaVwmRjJGY3JGMkY2KiYiRktjJzQrRWM7MzQ3YmslNEI5bTAoRjJGZ3NGMkY2KiYiRms7IW85cUtZSSNlQVoqPkdcYmhkRjJGY3NGMkY2KiYiRiFvOiVcNXdfWHQhR1o2el1neis+RjJGPUYyRjYqJiJGPzYhSCZmT3MyODNKMTkkKm8iNGBCRjJGSUYyRjIqJiJGcy0qUllMKD1KclpNVy1sQUQ5IVxGMkZNRjJGMiomIkZLIypSOUBMSil5Y3pmSzEpRz0iel5GMkZVRjJGMiomIkZjRWlONG1MUytOakEkPToxJFEoUkYyRllGMkYyKiYiRldkKCpRPjtBNG9jRDdSeW9FSlMjRjJGZ25GMkYyKiYiRlc0bFRQIT1bMiMqXEg2X2k7NXU2RjJGW29GMkYyKiYiRUt3PlNxW1VZOmhdSSdHby1kaSVGMkZjb0YyRjIqJiJFUyN5KmVLKVF4YyZcKDNbWSo+dEU5RjJGM0YyRjIqJiJEKVteMCUpUVhbJTRuVm13N0IyNiRGMkZbc0YyRjIqJiJDJ0h4MCE9M2ppRlUqW0xQLDFcI0YyRmdyRjJGMiomIkMzP3EyKHokKkcxKikqeURfTGNpO0YyRkFGMkY2KiYiQy94cm8hXFVAKHlsIilReSFcRjQiRjJGXXdGMkY2KiYiQjtuODpfeCNRYTclb0hxLSdHVEYyRmF3RjJGNiomIkJDVCgqZTw2Z2AnPUN5OD09YzZGMkZnb0YyRjYqJiJBLSNmOGomPmIkcCFSVSV5IzRoQ0YyRl9zRjJGNiomIkBLXTxgI2ZScSc+UGNXczokUUYyRjlGMkY2KiYiP3N3PyVRNGg8KClRJSkpZjRzUkYyRlt4RjJGNiomIj4nelwob1pWJT1dS1I+eV9ARjJGUUYyRjYqJiI7O3QlPWdKdWVgJio9NiM9RjJGYXhGMkYyKiYiO1tAN10/azd5QSVvPDclRjJGX29GMkYyKiYiOUlSSCpwNFs1NWNjJikqRjJGRUYyRjYqJiI3b11fbWMrbFtyKClcRjJGNEYyRjIiNHFZaCM0PjhXIXAqRjZGMilGanhGUkYyRjIqJixqbiomIkQpRyl5PiN5KzkjUmBcNz40bzMhR0YyRmNxRjJGMiomIkV3eD09bXlTYS5bLVVzRkRrSl5GMkZncUYyRjIqJiJGW29fW2EmRzU5byw1aUVSUWpFTEYyRltyRjJGMiomIkcjPnc+MEM6PngkMyUqUSFIciMqZWM5IkYyRl9yRjJGMiomIkczLSgpej48WidvJj1jJlJqeD5OWVMjRjJGY3JGMkYyKiYiRztTRFZaPy1SU24hSG1veTNiaDgkRjJGZ3NGMkYyKiYiR0snXCYpZSYpRzchR1xlLWBhK1JedT9GMkZjc0YyRjIqJiJGO2dzN0wnZltTXCdlaygpcDo0QiYpKUYyRj1GMkY2KiYiR186Sl9BcDZVdzpmXW5Xd2N1UVVGMkZJRjJGNiomIkc7d1smWz1aKzBNIVttVHV5MVpsaEYyRk1GMkY2KiYiR0sjUik+TkMuXko0JFt4ISo9J1x2OWdGMkZVRjJGNiomIkdLYztDPiNbbiRRcjYuN0l6Iz1BWSVGMkZZRjJGNiomIkcleVFLKFxlUCpRJyozO1IiZid5XXJoI0YyRmduRjJGNiomIkdrLyQpMysiKWYjNDNyKUc8KnomW1I5N0YyRltvRjJGNiomIkYvdigpNHNrJFs2IVtsZVFDLzdBRiVGMkZjb0YyRjYqJiJFJXkxU055Kj1hIzQlXHEjKjMiKmZzJipGMkYzRjJGNiomIkRzbWdNUHpRWi40NWp4QWdRdyJGMkZbc0YyRjIqJiJFa1s2UWVdNy9nTyoqNHdhJW9jVyJGMkZnckYyRjIqJiJEIypcJSo0cF9kcmxXME83N1AmcCcpRjJGQUYyRjIqJiJEMyk0UUJTIWVSX1osXkJEdSo0TkYyRl13RjJGMiomIkRHQGhMdSx4Wj9Qdlc3UyMzQDZGMkZhd0YyRjIqJiJDQVhNNSEpPUwoeU09M01OcyMzSEYyRmdvRjJGMiomIkJVXU9WI3lgUkhpNGkublEtZ0YyRl9zRjJGMiomIkEhM1FaUEo/O0hcQCxhJ1JMJCpGMkY5RjJGMiomIkBROk1nJlJeU2YncCc9JlwhPioqRjJGW3hGMkYyKiYiPz08YzR4alBXaj9VWzUneiZGMkZRRjJGMiomIj11XSgqM1FjO215KG9meUAlRjJGYXhGMkYyKiYiPCtuTlU+XGFdaGpWTEMqRjJGX29GMkY2KiYiO3FiU1gjeW1MYywpKlEkPUYyRkVGMkYyKiYiOG0rJlJQcnBcQFJXdEYyRjRGMkY2IjYyXEtueSF6IWYuSSJGMkYyKUZqeEZceEYyRjIqJixobiomIkUhRzB2bFEkND9JcTojW08rKSk+WyVGMkZncUYyRjYqJiJGQy0pXGctKikzaihvKlxFa2hVaHBcRjJGW3JGMkY2KiYiR1cieVlrUiRbcThOQD4lUSFRaHdqQUYyRl9yRjJGNiomIkdbayZbInpTQ24kPl8peVU5OnAqcFUmRjJGY3JGMkY2KiYiR3M5IT0lSFpoITMnM2lNZEt6KFtXUidGMkZnc0YyRjYqJiJHc1lnK0l1O0loPml4UVRJQSJ6QyJGMkZjc0YyRjIqJiJIKzspPSs8KVFMKkdWeEpmLkw2R1Y+RjJGPUYyRjIqJiJIOyc+Xihvb15oT0YzY0RvWHcoNElTRjJGSUYyRjIqJiJIP1pESzghekRpKSk9TnVJeGwjRyZ6XkYyRk1GMkYyKiYiSCVldkUmKVFXYSo0RlZBUCIqZnIqelJbRjJGVUYyRjIqJiJIa0srOGlhYUxQcylRPVBfKW8lXGJNRjJGWUYyRjIqJiJIITMpPkxPVk1DMHo8a3I5KVJTTCIqPUYyRmduRjJGMiomIkczUVksVTdsd1NZQjFmL15MISp5dUYyRltvRjJGMiomIkcpKWVnKlsmcDh0OCQzT1BvWEQxPmEiRjJGY29GMkYyKiYiRmNhOlpTO2ZmJipSKnk2Jip6WkwiKlxGMkYzRjJGNiomIkYleSUqUU0jZlpkZlkmUjlWNTRQdHJGMkZbc0YyRjYqJiJGN0hnISpISW5zWkFtPERCUjRZVSVGMkZnckYyRjYqJiJGKTNNJT52NjBEQF5EI3pKdSNmcig+RjJGQUYyRjYqJiJFL1ZHJSkpR09QKzxtTkgnem45OXJGMkZdd0YyRjYqJiJFa005MWZxXUFiIiopUlZjJXlKUUBGMkZhd0YyRjYqJiJEUXBYVlJHMWZrbyVlWEc5aXlgRjJGZ29GMkY2KiYiRFpSUShbWzc8QWFhOz07JFIxNSJGMkZfc0YyRjYqJiJDWyJvYyh6aF84KVEyVC1mSy10IkYyRjlGMkY2KiYiQkM1WTQzNUJbRjlDUy47dyo9RjJGW3hGMkY2KiYiQW01dSZvcDJUckdQekVwOT4iRjJGUUYyRjYqJiI/O20hM2pVd1hzbHU0biV6O0YyRmF4RjJGNiomIj5lLU42WyJ5RmMjZTBgdmQiRjJGX29GMkYyKiYiPCozdyVvJWVUQ2Apem92aSNGMkZFRjJGNiomIjk/JEhjcDJwZVJRWUUpRjJGNEYyRjIiNygzdHRkL0U3I1EjUiJGNkYyKUZqeEY6RjJGMiomLGZuKiYiRlsvP2NmdTo2Q0FaJClHJnBEUF5DRjJGW3JGMkYyKiYiRy8oW3E3dXM9MGdMRkEnPVlLSTQ+RjJGX3JGMkYyKiYiRzM1XFMnKj0lKVF0Z0YjUnltVD4+JVtGMkZjckYyRjIqJiJHJSl6SEtHO2hdI0cwVSl6Y3Mlb3chPUYyRmdzRjJGNiomIkhHSDBwdU5nTlkiKTRYKT4mUWlbSEYlRjJGY3NGMkY2KiYiSWNRJkhKIVx0KCl5RlNhXVksIT5dNEoiRjJGPUYyRjYqJiJJR3AmcGRaWSthSFdfQ1socHR4KkhOI0YyRklGMkY2KiYiSWshKSlHREUiPktqbDZTdkIuNkxXQUhGMkZNRjJGNiomIklfMixnJCpSRDtBRFonZTQiKilRMFZcRUYyRlVGMkY2KiYiSSM+I1JjYT8rZCR5KilSejEseCkqXEh2IkYyRllGMkY2KiYiSEM5QiFRKDQ9bVpWYiEpZjhBJVFaPXhGMkZnbkYyRjYqJiJIMyl6MTZbOTdwQCZmW1JOYTc8Ri4iRjJGW29GMkY2KiYiSFdQT14wQSk+S2AlcDhRej9cTSt6IkYyRmNvRjJGMiomIkgjekBkc3c1Jik9JUgqUVlFcXcrKVIyI0YyRjNGMkYyKiYiSHNZKEclbyFlY1RjPkdBQiEzaiU9LzlGMkZbc0YyRjIqJiJHX2okZVZCXTAjUXI4OlYmKWZOTXVyRjJGZ3JGMkYyKiYiRzBjSFIiPUJPaVFtTXkjR2kqZTgoSEYyRkFGMkYyKiYiR3E3SnR0IlxDak8jNFVVeE12dkk1RjJGXXdGMkYyKiYiRlBzYCZRazV6X0duMldzSmtXU0lGMkZhd0YyRjIqJiJFaz4qKVswWllsOllSeFwlW2AlKWYoRjJGZ29GMkYyKiYiRT1YKVxXKSlHZE95dlRNKkhbaWs6RjJGX3NGMkYyKiYiRDUib2dhJjNwM2YkXCdmSWBzKjRERjJGOUYyRjIqJiJDdE8qSERpZHpXYVkxdSlIVGNHRjJGW3hGMkYyKiYiQjlmUTRgREpCKUdLUCoqKSozQT5GMkZRRjJGMiomIkA7Pj4sVlpTSihwY2dgbVFRRjJGYXhGMkYyKiYiP3N6bz0tZWR6X3h1XmUqNCNGMkZfb0YyRjYqJiI9XjttSWlSQGBGKWUsPClIRjJGRUYyRjIqJiI6J1EmSEwlNDtvOFNbXnNGMkY0RjJGNiI4QVxSbzckKlw6KCo+QSJGMkYyKUZqeEZgc0YyRjIqJixaKiYiRiFbTzNMVlBaUSR6bE5KWDZbSWkkRjJGX3JGMkY2KiYiR3MjPVg0NS9JJnAhbzBMTjxDSHRqIkYyRmNyRjJGMiomIkhLKTM3a2JII28lKSlvK1MpeidHeWwpUiRGMkZnc0YyRjIqJiJJU00qW3kmUTYtVGlUJipcODIqKXpqTTxGMkZjc0YyRjIqJiJJKSk9bUFYIXBWLlAjZjJuLWJrI1xwJVtGMkY9RjJGMiomIkkzLSgpekYkR2EoKXprNUFHJCoqPnVXbygpRjJGSUYyRjIqJiJKJVF3cFghKTQhKkc0bD0kUUx0L2pAMjQiRjJGTUYyRjIqJiJJL050ciJSJSoqPjIzczNWKDQ1P1NCSCpGMkZVRjJGMiomIkkvTnBbVndxYkYiM2FqOUJ1Ti5EcSVGMkZZRjJGMiomIkhDQSp5YTBBbl9BTiFwQWFFLClmOzdGMkZnbkYyRjYqJiJJRyhIUnI1WHBAWVMoUksoNHR0W1cpSEYyRltvRjJGNiomIklHNCM9UiQqM3VeOil6PT9FYCJmTCNcTkYyRmNvRjJGNiomIklXQGVoTXgvem4vWiU0OV0qZitHXkZGMkYzRjJGNiomIkk7R2MhNGxIJDMhKio+IVE4YmVsaEdYO0YyRltzRjJGNiomIkgpbyd6T0EhZT02KzQzN3VLJ29HcCwpRjJGZ3JGMkY2KiYiSCUpXEY/I1FdQHZ4NjZVcmVFLEtsS0YyRkFGMkY2KiYiSD9oV2VJJz09bD5sRiR6UC01NylHNkYyRl13RjJGNiomIkcneidIJDMmKXokPjIrMXQuT3hFL1JMRjJGYXdGMkY2KiYiRid5NFg0KDNlSUZMbzVZI2UjKT06JSlGMkZnb0YyRjYqJiJGKjNqXzMwb2dVVyEzOSsjKm8meWg8RjJGX3NGMkY2KiYiRVcjSGdkdHA+XmImeSJlVk5LXyFIRjJGOUYyRjYqJiJES25mRiR6ISozM2wiR3g5TClvWU1GMkZbeEYyRjYqJiJDdSV6N2BCbio9YUp2ckNTciJbI0YyRlFGMkY2KiYiQXdWJillbjtFKWU9N0syZUVEJ0YyRmF4RjJGNiomIkBhcEJcWG5Gc1dAciZHJj1AI0YyRl9vRjJGMiomIj5YP0pAdGsmM0EpMzZnYHQjRjJGRUYyRjYqJiI7PyE9VzEmeUQ3dkldIypcRjJGNEYyRjIiODljU1UjZSIqZT5LZCopRjZGMilGanhGaG9GMkYyKiYsWComIkc3QF9xPnkpKTRuYSJbIVJrdDlBMFFGMkZjckYyRjYqJiJIWyc+LUdLeUVpIT5zLCRlOEEsKjRzJkYyRmdzRjJGNiomIklbR00oKTR5XV1nZ18rJVFiQiV5KEhKJEYyRmNzRjJGNiomIko3LGxNbypvJjRMKClmelYnKVs7OXAiWzVGMkY9RjJGNiomIko/PksnKXBvJFx3anEjPlJAISpRSWd3LyNGMkZJRjJGNiomIko/ZnknKilwc2FJI1FgR1VfenIoSD84REYyRk1GMkY2KiYiSiEpRyd6Jylbdks9SCgpeSZHLiJRJikpUTdpIkYyRlVGMkY2KiYiSUtfXF0qWyJlRkU9R1YpKmVCXWRcTldGMkZZRjJGMiomIkpDaWpKYmFRaTYxdEw2S04lXFA6NERGMkZnbkYyRjIqJiJKQ3EoNEgjW0JjQEtZKSkpUXoiNHJ1ZV0kRjJGW29GMkYyKiYiSiUpPjtdWVA8S3VvKVs+W0BJUyhcJnBLRjJGY29GMkYyKiYiSiNmJ1tmYTpnRik+PXlRcnZGM0IpKltCRjJGM0YyRjIqJiJKISlvPiUzQzlFLTldQjcuSytASGFyOEYyRltzRjJGMiomIklDJSpbR2dLMzA+J1IhKVJ0W0o3LmhuJ0YyRmdyRjJGMiomIklHTE5pIio0XilbUiVIKlIkXElrW1xZRkYyRkFGMkYyKiYiSDdgJD0/JnBLUC85eCcpRyQ0XkEtTScqRjJGXXdGMkYyKiYiSDcsUUpQRlVJS21ccih6dXJCSycqR0YyRmF3RjJGMiomIkdpIj55KClmPydbRSw6IkcsRVAmNE91RjJGZ29GMkYyKiYiRz4heS08KSlHbEEpKWZWeGZRbzddZiJGMkZfc0YyRjIqJiJGZ09XUUZvUiQ0dlVwWmBhT3Q+RkYyRjlGMkYyKiYiRVNWeSZcLHY6ciUpKVwkZmohenh4TEYyRlt4RjJGMiomIkQxKXA+QkhmKFsjNCxiVDhsOSxFRjJGUUYyRjIqJiJCJ1xCSE0lZWspKUhQKSpSc2BKIXlGMkZheEYyRjIqJiJBI1tKSC1XaygqZVsraFJLcSY9RjJGX29GMkY2KiYiP1ZnYmI5SEhqSF4rYEtoP0YyRkVGMkYyKiYiPG9eUCllTEtfJipmL1RtRUYyRjRGMkY2IjlRWEshXGQjZWx5N2ViRjJGMilGanhGYndGMkYyKiYsViomIkhvYkddIyo9cD1qbCkpKm8rIjRVRyE0SUYyRmdzRjJGMiomIkkjKipbLFNWeGczXSF5bGc6R1Y3NkpJRjJGY3NGMkYyKiYiSm8iKipwJFtXKFwoeUAoZTp2I3koRyE+NEUiRjJGPUYyRjIqJiJKP05LOF47TWZMWWJMKUdLJVE6NSVHRkYyRklGMkYyKiYiSnc8dFllWF9gTixDJXo/OiZ6aU5MbyNGMkZNRjJGMiomIkpPdlwxViFcaXh6RllROWlGcC9vVztGMkZVRjJGNiomIktTJXBReThSMnF5SiUpXCtpZi5dPnAsIkYyRllGMkY2KiYiS0coXDU2SUQmSEdbJEhCLSRRL0omZW8nPUYyRmduRjJGNiomIktDQSxEKClINHopPmJxKGZrKWZPWF8qUkFGMkZbb0YyRjYqJiJLb3Ipb0coZTdAVSwpPXFdd009R3MuLSNGMkZjb0YyRjYqJiJLXy5pXCY0XSQqeUh4M0YuaTthUz5QWCJGMkYzRjJGNiomIkpfSjZ6b04vKmYoSEhGQiI+JzNgSyc+JylGMkZbc0YyRjYqJiJKOztySk8vb10icDciKkcpeU41Y2QvSCVGMkZnckYyRjYqJiJKOyFIOlEnZmxyPl4oMywyLlk1KmU3Ij1GMkZBRjJGNiomIkknXFNwJTR0WSMpR20mXFE/JWUqMylvQ2xGMkZdd0YyRjYqJiJJW2koKmZEQnJxOjsmUmQ1NSYqKVElRywjRjJGYXdGMkY2KiYiSEVJN01iLm8pKW81TC1reUJ4W1pJJkYyRmdvRjJGNiomIkhqISkzWlAjKiozJnA0SUZBYUEpendzNkYyRl9zRjJGNiomIkd3J1JUUj4+cWFxRmtwIilcVCQzeT9GMkY5RjJGNiomIkYnZiJcNE9XKjRRcypSWm5sdyhbN0ZGMkZbeEYyRjYqJiJFcW88OkM5Rlk8KlFBPVRAel9CI0YyRlFGMkY2KiYiQzt2JEd2OEQ9VTVgTlNJZiRHeEYyRmF4RjJGNiomIkInPT0kWyZSdUR6Z0hiVEMrUzdGMkZfb0YyRjIqJiJAMkRWUEwwdENBS004Pk5IIkYyRkVGMkY2KiYiPWMzMix6PFc6J1JVNDgwIkYyRjRGMkYyIjojcDhgQ3l6aVYiKmVaSEY2RjIpRmp4Rl53RjJGMiomLFQqJiJJcyl6Q215YWtpKio9JnAyKSpING1lMTVGMkZjc0YyRjYqJiJJP1pAQCdlJVFoLjxqPiZcY08/ZHcrKEYyRj1GMkY2KiYiSiMqND1rME1LalwyZm81Lz9dVUZ4ViJGMkZJRjJGNiomIkpjYShRLjw/ZCEpKj5yb0Aob10qKVxiKCo+RjJGTUYyRjIqJiJLO2c7UXdbVmFSdHRZLTZedkYjKVtmPEYyRlVGMkYyKiYiS0NBJm9GPjk3InpKP09XcSNwI3lNempaRjJGWUYyRjIqJiJLM0k7LCgpPWRAZ3hlREVPISo9VEt4TCEpRjJGZ25GMkYyKiYiSy96WTg+MSlvXTNHOzpGZ2ctWDYiWygqRjJGW29GMkYyKiYiS1ciPVlURSV5MF9QKXA/eGFbdVhmbjQqRjJGY29GMkYyKiYiSzd0NjwhPihvQDk6Qlg1Yl9YViQ+PiFvRjJGM0YyRjIqJiJLZyhcPmhOZFU/bzNbc1pMKEdlS08hPiVGMkZbc0YyRjIqJiJLIyo0Xip5R14qUnpwTWcmPU41cSwmKmU7I0YyRmdyRjJGMiomIkpRJSkpPSdSMk0xTT9yJFxPelFRazwhXCpGMkZBRjJGMiomIkpvWTY1KillVEI5Qy85Vk5iRGl6TWEkRjJGXXdGMkYyKiYiSilIP0NBJCkpWyl6Z256PjRyQWs9XEk2RjJGYXdGMkYyKiYiSVU2YVAyQ1shKTMmPV0mZnJOTEBUeElGMkZnb0YyRjIqJiJIYjZ2Z0hLVjJqUDpsK14mcCZHby8oRjJGX3NGMkYyKiYiSF9YJDNkTSIpeWZgZSlcKTQuIylROi44RjJGOUYyRjIqJiJHYVonb3ZFL2hsLmwmKj5SKVJHcCV6IkYyRlt4RjJGMiomIkZVUSpSbyFwPnpzRyZcPCNwKSp5XmUiRjJGUUYyRjIqJiJEMyJRJHpsd0ddUz53SWEqPUsnPidGMkZheEYyRjIqJiJCQXFzWVcuJipmSHpEJkh5eiRbJ0YyRl9vRjJGNiomIkAwIlFEemcnUSEpKlFJKCozJG8lb0YyRkVGMkYyKiYiPSNSK0JaJnpBXDpRVXYyREYyRjRGMkY2IjtWJ3pFcDU0USpRbGtXOEYyRjIpRmp4RkJGMkYyKiYsUiomIkh3LCRwYWQmRzRweW9wKHplKWU6OnYpRjJGPUYyRjIqJiJJVyI9QHNMaUN2RF8oM3JqYFh5IXphJ0YyRklGMkY2KiYiS29CN2tlJSo0JD4yISp5TEYkZiFbVHg5KyJGMkZNRjJGNiomIktXKHo3TUldJUcuL1JucmVcOz82aFxaRjJGVUYyRjYqJiJMKVsielA8Xz5nOGlBZCdbTGd1Y01RI0ciRjJGWUYyRjYqJiJMNywmUVxIKkgiWyNbY2AncDwkKmU+SHA/QkYyRmduRjJGNiomIkxXNCplVGZtci5YKip5VVYrYzZubTpeSUYyRltvRjJGNiomIkxfalJWTylbOkFoaWwsVmFeQGVlLm5JRjJGY29GMkY2KiYiTGtbX2IyVDksOiN6JipcTkVbRGttaVlDRjJGM0YyRjYqJiJMNzxuNSFHXzhaNXBLLShmUyEpKik9QVNmIkYyRltzRjJGNiomIks/KFt0QnRGXGd4UWs2ZSY9cGV5eGgnKUYyRmdyRjJGNiomIktPLkR1NXZ2KmYpNDkjW04iZS5WUkN3UkYyRkFGMkY2KiYiSz94YkNqJClmLVtgWW5aI2ZpUWsiPV46RjJGXXdGMkY2KiYiSjczX1JbJXkrQT4oPUNsdGcjemxcYF5GMkZhd0YyRjYqJiJKLVxQQikqM0tNI0gnZmg8Klw6VnI5ZDlGMkZnb0YyRjYqJiJJRHJCdlpCWFpOYyhvayZIdCpROTFaJEYyRl9zRjJGNiomIkhvcTlsKD5cISpSV3hscVJ0ZVZZQm5GMkY5RjJGNiomIkdzPEUuTjwjXChcJylmZ2E9LEojeiEpKkYyRlt4RjJGNiomIkZVb180XkknZlQ3WituKlFlOiY0JCpGMkZRRjJGNiomIkU3JClwV190IW8tKD1qUWQvRV1vU0YyRmF4RjJGNiomIkNVUilHc2dpXz9iQjwza3kiUkRGMkZfb0YyRjIqJiJBSFdyIj1TaGJzJjNnJ1xddzQkRjJGRUYyRjYqJiI9M20oPS4qWyJ6JlxsdkApZiJGMkY0RjJGNiI7UVIycHohNDBuWTdzSCZGNkYyKUZqeEZockYyRjIqJixQKiYiSSlHZ0sieSgpKjMjbzE2WEQhXEYhKikzMyEpRjJGSUYyRjIqJiJLJT1kdCk0QyN6O0ZlOid6S2NeVFpAbjZGMkZNRjJGMiomIktnUCIpcDVueDAkKj1sKD5NKG96QDYiUnInRjJGVUYyRjIqJiJMKSkpKik+IVJQQWhsSnhnaikzUl16RF15QEYyRllGMkYyKiYiTGs3Km89JGZOOnp4MGlcJ3BOIjNJXSI9WUYyRmduRjJGMiomIkxbY05qOCtKeSxlJz0lKlxHIj0nem8sTnBGMkZbb0YyRjIqJiJMJT0mW1ArSCxVSEtfV00hNGYvTklwIXkoRjJGY29GMkYyKiYiTFMpPjd1K0xbIUg/TVdUUic0KTQmcGEpeSdGMkYzRjJGMiomIkwncD5VWj4sbTwoKj5gQ2tiJ1whKnpkJ2VaRjJGW3NGMkYyKiYiTCU9RCQqPnFcKHkoeUYlND9SLiIqSDopZV1GRjJGZ3JGMkYyKiYiTHM8QDg0bGZNMTZ1P18kcER3UHFJTThGMkZBRjJGMiomIkszc19jKFtaJGVkdjZdMD8kPkY8I2Z5YUYyRl13RjJGMiomIktnKillJDRKInkyXWpBL0pBIXBhREQiMz5GMkZhd0YyRjIqJiJKby0/LzJBVDRiPl0lKTRfQSFIcklNY0YyRmdvRjJGMiomIkoiM3VxIlFrTTluKEhAJj53XSMqeiE+LDlGMkZfc0YyRjIqJiJJIVEqSGBLISkpKSpSem5IO3AiXCVbT0gmR0YyRjlGMkYyKiYiSCFbIm8oPVw+Rm9WRiI+KUd6PSRwY1UlRjJGW3hGMkYyKiYiRycpMy9hIjRWJHkyOylmIkgsKHpYM2AlRjJGUUYyRjIqJiJGL2VtRDcqKipRPnF0JFwkUixeQyxBRjJGYXhGMkYyKiYiQ29bSEQ1azY5V08zSycqb1NYa0YyRl9vRjJGNiomIkJsRFBoOFg4YzRSJilIQigpUkAiRjJGRUYyRjIqJiI+S29PJzRJLz8pSDo3UE1DJkYyRjRGMkYyIjwhNGl5XHBfJypSSWQlXCE9RjJGMilGanhGXHNGMkYyKiYsTiomIkonNCcqKm8nUmZ1VCRSO3pPIlFXKVEhW0xaRjJGTUYyRjYqJiJLM0V4ZyY0ajlzWFY/X2U5PyYpPSIzalxGMkZVRjJGNiomIkxTNSRlJD4rXCU+VT4zZiZcVmpHbkFXRyNGMkZZRjJGNiomIkw/PiVRNCRHbWk6JnpBOylRWzk+KVw7NGlGMkZnbkYyRjYqJiJNWyd6Ikh6NEVcNUhdQWAnZjtkIz5mK0c2RjJGW29GMkY2KiYiTSc0W3hNKGZQd29VKm8hSGE0ckdGI1ErWiJGMkZjb0YyRjYqJiJNd0QiPkF0XkJqLG9lcXQzXGptSSJSVjlGMkYzRjJGNiomIk1TJXAnKXljbU0lKnlwI3lgOiI9bXhORi82IkYyRltzRjJGNiomIkxvJG8xPU1rIVI9M1QqekEkZTdKOHZyInBGMkZnckYyRjYqJiJMJT0lemtTb1YuXFxwcTwvTWk3cSQqcGQkRjJGQUYyRjYqJiJMN0xiOyNbO2xSTyZHNV5POnglUVp6YzpGMkZdd0YyRjYqJiJLPzAmW0ElPidHNiIqZiFmSls1P0pbU0FkRjJGYXdGMkY2KiYiS0c6R3pJWFNqZiNcNSopKVJ1OW1RZ3U8RjJGZ29GMkY2KiYiSlNyKjQxRGZJRkA/OkBvKyhSOFZdaSVGMkZfc0YyRjYqJiJJW1d3JipmQDQ1U3lMQlpmTE1sbUMqKkYyRjlGMkY2KiYiSVMnZScpM2x3UkU3KHpOa0g9MGVlVjtGMkZbeEYyRjYqJiJIT05QLkFTXj0tI2UmbyFcdCg0cU0jPUYyRlFGMkY2KiYiRmdCeHJpeCZIaV5yUmksb0VtTCkqRjJGYXhGMkY2KiYiQyFvPyY9UEZZLXFgWyxGdl1cREYyRl9vRjJGMiomIkIrOWUkZiFbbicqWzFzVmAnSG9URjJGRUYyRjYqJiI/IyoqKkhfWzt3Wj9pKGZFTjgkRjJGNEYyRjYiPGU2VygqKikpUT9faF5VOmBGNkYyKUZqeEY1RjJGMiomLEwqJiJLYyNHKFwhKW9PJlJmejJWYSg0a3kpSCopWyJGMkZVRjJGMiomIkxvOnk7aHFQSzRxQ2ZnPCoqKXAzei5LOEYyRllGMkYyKiYiTFsrLCgpUTQyWWZGcF9wRHQwYGxcIj5gRjJGZ25GMkYyKiYiTT8mPkE/RClmKCl5W2RGQVUhXFttM0dYRSJGMkZbb0YyRjIqJiJNT2Y6LGVjJClvS3BnI3BEK3FeZlw5bSwjRjJGY29GMkYyKiYiTSdweCVlPyk+bExlJilbYjx6cTMhbzUiPUojRjJGM0YyRjIqJiJNJ1xVbCNSOzMxUTsrMStHS0JOVC9yLT9GMkZbc0YyRjIqJiJNJ0hadThuYGtkImZZIUh3eCszYV0kKXBPIkYyRmdyRjJGMiomIkw+cGU7T3o9NDFvcVsnM1c2UC5IRTZ3RjJGQUYyRjIqJiJMVVNNKDRhOGV0a1oxcSdlMSllTVpCT05GMkZdd0YyRjIqJiJMXixLLDBQKVJtYU9HJ3phXjchSGk/IlEiRjJGYXdGMkYyKiYiS29vX0smPWxXSWQpeVwkNFBEV209b18lRjJGZ29GMkYyKiYiSzohPjdgQCNRXHAlUmtDKWYnPUlaKUdUN0YyRl9zRjJGMiomIko1JEdCTmFoJ0g9R2x4SmNfKzYhSDhHRjJGOUYyRjIqJiJJYFdpbV5LcCg9eS1zLz0pSGpIYiIqXEYyRlt4RjJGMiomIkgpKnBvXEQ3cTwhXCZcOSgzWS1yYk9nRjJGUUYyRjIqJiJHI1s+IkghMypIOUElPltSPGNXPjdpJEYyRmF4RjJGMiomIkRBcCYpUWM/cTt2W3VkNUZ2TipwRjJGX29GMkYyKiYiQytjdHZKIio9SXlRQyM+LmwrRSJGMkZFRjJGMiomIkAzW2JaNE5XVl9HTy8hKTNDIkYyRjRGMkYyIj1gb18zXl5XXzAieiNlWzhGMkYyKUZqeEZkb0YyRjIqJixKKiYiSz82bUFXWGNsUl0pR29DLyI0JGY2VkIkRjJGWUYyRjYqJiJMIVs/L05qJClRSDF5VDF3RS1zJVFUKG9ERjJGZ25GMkY2KiYiTE9eYlFvKyplNCRbZzU8I29HI295Zy82KkYyRltvRjJGNiomIk0hb0pJTCo+JWYxUSlbODBfNXZfPSIzJj4+RjJGY29GMkY2KiYiTUtTK14qXCMqeio+Yy1aNHUwa2dwRWYvRkYyRjNGMkY2KiYiTTsnPlBUcUZXeCk+biZRS2d4UV4qcFFNRkYyRltzRjJGNiomIk0vMjNFLEgqeSpwPlpxNlRJLnkwW1UkNCNGMkZnckYyRjYqJiJNK28iKnpXeWFGdG8xJyopZXciUTo2Q3RyN0YyRkFGMkY2KiYiTClvV2dkem4wbT8kNC1XUTV1dlIxXE4nRjJGXXdGMkY2KiYiTGtRImY6YFRQN1pCJXpZO2QoUSE0O2FsI0YyRmF3RjJGNiomIktFJio0RCdwaHNCJm9gQTJLMUFqQmVqIypGMkZnb0YyRjYqJiJLJlIzcTQ5PDlsKWVsTTVpc1tMNEImbyNGMkZfc0YyRjYqJiJKXytqKFJNaTh1Ri9kdSNlLS1oNHdWJ0YyRjlGMkY2KiYiSjtIc1FiTk8leio9Llh1NS5DOXJyQSJGMkZbeEYyRjYqJiJJUXNxSi0wI3ksTVMhM3Y6WEEoXCFIO0YyRlFGMkY2KiYiSDMoXFgiKj0sXCNcKyxKeXlPVSVHJDQiRjJGYXhGMkY2KiYiRSU+TXZ3KTQsOiEqUUVEJ1EpNEYrVEYyRl9vRjJGNiomIkNacyJwKVtjZ1otXy5QJkciR0xMRjJGRUYyRjYqJiJAIzQ4dXdFKy9MITRcVzotdCRGMkY0RjJGNiI9cUNLJXorKWUhPkMpRyQpSEhGNkYyKUZqeEZcb0YyRjIqJixIKiYiS0slZVVqIz5OTD5ncCQqeU5sT153OnNeRjJGZ25GMkYyKiYiTGslPT5XNyI0MVMhW1wvJXo6IkdZUEsocCRGMkZbb0YyRjIqJiJNW18lZUMuVWhATnIzXkVVJSlbJ0hcRnQ2RjJGY29GMkYyKiYiTTMlNG0lR1didmxIKjQ8KGVdMEImeiN6Jj4jRjJGM0YyRjIqJiJNay8yLTVyST1Oblw+blZZOXFxQWtIRkYyRltzRjJGMiomIk0zMTUoR242ISlSKnpCIlI+c0RAPmwlR0ZDRjJGZ3JGMkYyKiYiTXNxXy0zTzRidSgzTGI1VCZcNXYxSlQ7RjJGQUYyRjIqJiJMcydvR3okZlNOJEdhVWMkbyYqPWllLkkhKilGMkZdd0YyRjIqJiJMV1ZNOCZwLi13NEApPSlbJD4iPSslcFsrJUYyRmF3RjJGMiomIkwnXCE0IVxsN0kob0gwRjhvNCIpZix2KilcIkYyRmdvRjJGMiomIktINShcS3omb1s+aEwuN3VxLERRekFZRjJGX3NGMkYyKiYiS09pI3BpQFFUUm49J1JNUClcYkhFTzwiRjJGOUYyRjIqJiJKYV9GeHQpKUhsJClmajhlb2hZRTBrUyNGMkZbeEYyRjIqJiJJOzY9Pi8nUTxoIT4wPHglR1dJYFdgJEYyRlFGMkYyKiYiSG01Ymc3W0Y+VWMoSGY2cTpQIyl6RUYyRmF4RjJGMiomIkZZNXIwZi0ocDdfSUI1LDtBaTY5RjJGX29GMkYyKiYiQ3glKT5HXCk0ej1MOmxHeXRIZChGMkZFRjJGMiomIkBpKGZ5cDdQczgnZjJcanMkKSlGMkY0RjJGMiI9LjUiZmUzSDlWNkFOeFImRjJGMilGanhGaG5GMkYyKiYsRiomIkslUSkzZl1VMS0+cCJISHomeTkhemUiKj0nRjJGW29GMkY2KiYiTE8iejpCRksnXGwwUShwJFtwKEcjb1ktKyVGMkZjb0YyRjYqJiJNJT1sanc2IlslUil5JHpwMGE4VktsUmc4IkYyRjNGMkY2KiYiTUN1I3p2V2RUQ2JDLUs9LzY6OyE9XiIpPUYyRltzRjJGNiomIk1HYD5tISp6ST5qLC5rIilIdERgLyxcXT9GMkZnckYyRjYqJiJNbygpM0NVR3gsVEI/JVtEaSlvTDkpUlxmIkYyRkFGMkY2KiYiTFM9QFh6azsmKmU1ajpjSSQzKlEjPXNGJipGMkZdd0YyRjYqJiJMQzU1Mj5MSWsqUkMoXEM+PmReKUg4UllGMkZhd0YyRjYqJiJMa18vRCQqKnknSG4jKXo5UWkmeSJmQzZ3KD1GMkZnb0YyRjYqJiJLXyw9azJpYCk+OSR5VEhmT2tMU108aUYyRl9zRjJGNiomIkspKUdjRnJOWmsiM2luSkBrJHkmXCNSdTtGMkY5RjJGNiomIkpncyF5KVJQdihwZ1EheV9ANFBGVGRvJEYyRlt4RjJGNiomIklzJCllcyoqcDd0QEVYIipIWHU8Z2lSZ0YyRlFGMkY2KiYiSDt4OFAqNCIpUUNAYClvPiVSYSRlOUQmRjJGYXhGMkY2KiYiRmtLZ3pQT3cnPmA5PDwkUXUkZUZNRjJGX29GMkY2KiYiRCFSKXknKSpcKmUkNCIqKikqUXpES1lWIkYyRkVGMkY2KiYiQUtnIXBNWHlFWlV0YE4jW2c7RjJGNEYyRjYiPXQjUnJZay5fWFVLRSg0JClGNkYyKUZqeEZaRjJGMiomLEQqJiJLLyQ9N1FBPiM0emMiPiFRTyVcKSoqXHEuYkYyRmNvRjJGMiomIkwzQSUqZSp5VCc9JVI6eGsjPVBeeCQ9MDtLRjJGM0YyRjIqJiJMb1YkeXcoelxgam5GXnQrPCQ+N0FnUiIpRjJGW3NGMkYyKiYiTVN1InpqJUgjb0NjNSR5KCkpKlxoWFRmMSQ9IkYyRmdyRjJGMiomIk1zWS1sMDdxdU8ib2hnMnJBL1BCdyY9NkYyRkFGMkYyKiYiTGNpKD48PVpCODg1Nic9b2ZVcj9oVmMoRjJGXXdGMkYyKiYiTCM+KVxwem8qeTprc1NMPFp6ciNbZjcsJUYyRmF3RjJGMiomIkxPeiYzQyQ9aCllX0lsI3ByWnAmPTJtazxGMkZnb0YyRjIqJiJLb11uXmhZVWtLUDZnQF1LNTRVUmFqRjJGX3NGMkYyKiYiS18+XTUmUVs5WXJiUyEpW2xIZSc+TEM9RjJGOUYyRjIqJiJKIWY+el4iZUBSQEFELGhNO25yNDtHJUYyRlt4RjJGMiomIklhJilbInpNMjw9SkwkKilRRiZvImU9Iil5RjJGUUYyRjIqJiJIIVJhWlFxUUdIcTQhRyQ9RF9neU4hKUYyRmF4RjJGMiomIkYjR2lrXClvJ3khWykzQyo9MVBGLzInRjJGX29GMkYyKiYiRGl1IlJJMj1teSFRdicqKWZZJDQ9I0YyRkVGMkYyKiYiQSFHViJSJHl5JSkpPj1yJVE7Y1cjRjJGNEYyRjIiPkxqXS0nb0wqZmoqUThbWDVGMkYyKUZqeEZWRjJGMiomLEIqJiJLKSlRaGspNEgoKjM+cyVbcDRIJXoyIUhTTkYyRjNGMkY2KiYiTGtLWzQnNDE6ei9AJ1I9Jm9KVnRbeGs9RjJGW3NGMkY2KiYiTCV5UVNSWThWSyNSJlJjSWVwIWZLP0g8JUYyRmdyRjJGNiomIkxba2ghKT5kaFp3digzbipwIUhpTGJAYF9GMkZBRjJGNiomIkwzQVknUTBTaVkkUVBITiYpKkciKVtsMltVRjJGXXdGMkY2KiYiTGMxal1FWjk+KSlbTCVlXyNmdWxZMDJcI0YyRmF3RjJGNiomIkxfbiVlTEBiJzRaUC1PZGIqZSNHJHokKSk9IkYyRmdvRjJGNiomIksjUjxbLDU7YG0rMGE+YE5WJnl3KDNzJUYyRl9zRjJGNiomIktXOCk0V2MiZmN0b3pdYHAmSGokPVtsOUYyRjlGMkY2KiYiSksjXCRcIylbZyk+dV49Pk05IWVfMSY+T0YyRlt4RjJGNiomIklTZWtMPSFmKFxlRTQmcGBYcDRaWVsoRjJGUUYyRjYqJiJIP0YxdlE3ZlxcMi0xO0lcSz5rRCpGMkZheEYyRjYqJiJGQ14+bjgyVkIzRyVScW8kb1coUnhGMkZfb0YyRjYqJiJEaW1oKyJmKm9XJm8+PFY4XWZMREYyRkVGMkY2KiYiQStEaEpUVSZIYj12QHB1M3UjRjJGNEYyRjYiPiZ5XjtjJ3BaRSc+UiVSJVE1RjZGMilGanhGTkYyRjIqJixAKiYiSytLc2NgVUs/SCkpR0UqR0ZaKFxObGY6RjJGW3NGMkYyKiYiSyFba0hDTSdcMjMqR0MkKjNeTltuZGdPKEYyRmdyRjJGMiomIkxbYyIqKkcubmBsJG9mKD4rYChbY05yK1ciRjJGQUYyRjIqJiJMKUdrTll5YjFwUmEtazcsJFtASzwsYSJGMkZdd0YyRjIqJiJMOy91VlddRCRRP0ErVE0pPi9TXzduLyJGMkZhd0YyRjIqJiJLJWVeWytsMHU6IylbKCp5PEEjKlIrbUxQJkYyRmdvRjJGMiomIksnKm9TJzNOYSdRcmAqZkRyJXlyM2M4aUJGMkZfc0YyRjIqJiJKT1pfa0VMIVJGOE1wMiZvMWtUP3E9KUYyRjlGMkYyKiYiSmdcJCkpUUxSXDNPN2o2LlYrJipwIyopNCNGMkZbeEYyRjIqJiJJOy8qUmAkXClmaHFmKCkzMSUpM3JGM3glRjJGUUYyRjIqJiJIcylSLjl6Xzk2KiplP1JAIVxdPFNiKEYyRmF4RjJGMiomIkZfJD1iLk03ISpcVDlNaSE+OWAvdydGMkZfb0YyRjIqJiJEbU1qTXRTbUg1VUFaL3hnKCk0I0YyRkVGMkYyKiYiQXVDRVtfNEhJJCkqcGVCViYzQUYyRjRGMkYyIj1WZzBqem1sbUBhUGAnbyhGMkYyKUZqeEZKRjJGMiomLD4qJiJKLzZeUWUmKnk8dHFKMWh4KD1nWCFSQCVGMkZnckYyRjYqJiJLWyl5K1c0ZWUiW0ZLKj5Qcyk0cTo0cjxGMkZBRjJGNiomIktfNmdGIlJ6MW9rQV1jLW8qcFsyOnRIRjJGXXdGMkY2KiYiSztbPic+cicqeSN6VllHQSFcMnU2YV1pI0YyRmF3RjJGNiomIktTdU4oR0FlJHlrVE0vRF1zKSk9RSt1OUYyRmdvRjJGNiomIkpTSTkwYjhqQyJbW3FVKkhEP3kmUTdwRjJGX3NGMkY2KiYiSj9mI2VGWiZlTUcwYVNZQyhbIXpfXyJHRjJGOUYyRjYqJiJJR2hEc18vJ1xkeihlR3NmQmpyInppKEYyRlt4RjJGNiomIklvXiRcKDMqUi1HJlEsSCRmKnlMQFlYPEYyRlFGMkY2KiYiSEN5al9GQ3dLV00vOlJhKSllS2sqUUYyRmF4RjJGNiomIkZXKCoqXE5QMmVhIXA2dyFII2Uyd2kkRjJGX29GMkY2KiYiRGsrXW8zQiczayo0eWFhL0MxNSJGMkZFRjJGNiomIkFHXkI3dCZRQFdcMm54KjRWNkYyRjRGMkY2Ij0zeDpzKCoqUUdsMzo7Iip6JEY2RjIpRmp4Rj5GMkYyKiYsPComIkkpKVEiKXolcEJaVGoqRzg/WkIrMlFuX0YyRkFGMkYyKiYiSjdINVM5L2lXRTBGVidvO2hGYVxdPkYyRl13RjJGMiomIkopbykqM25ydUZ3YyMqKikqZndvP3MhPlRGRjJGYXdGMkYyKiYiSjNdMmtJKCpmZmo8MXooelxzJkdBMiI+RjJGZ29GMkYyKiYiSXdAKVFEIlwqKnpRQF5CSGRZMjhVcikpRjJGX3NGMkYyKiYiSSN6LDdKWSJldiEqXGltZigpenQ6eC9VRjJGOUYyRjIqJiJJV3QoUURYS2wxMmJuZm40Nz8hbzs5RjJGW3hGMkYyKiYiSFdMNEtJUWc4UmY+YmUnKnpXJlxeQ0YyRlFGMkYyKiYiR0MtS1Z0IUdLP3dHJikpcFFLJSp6ZyYqRjJGYXhGMkYyKiYiRTsnPkNuRSJbX1xnNiZcdSk+NTkhKkYyRl9vRjJGMiomIkNnXjtMKUcmKT4lKWVjcFFUW2V1I0YyRkVGMkYyKiYiQEcqej5XOTpULkxwJCkpXGUmR0YyRjRGMkYyIjxGJVJJJVwoNEs7eC56KFwqRjJGMilGanhGZHNGMkYyKiYiNSYqcHlISFFFMHQ3RjJGW3hGMkY2KiYiOU0oKik0aEcpMzRyci5dRjJGXXdGMkY2KiYiOWoqPXRmOHBNLCcqXG0iRjJGYXdGMkYyKiYiMjlJPS1nNWxxJUYyRmF4RjJGMkkkb2dmR0YoNyNGLg==
computeParams(algeq);
NyYiI0MiI08iJEQqIiRCJw==
algR1 := primpart(subs(F(t)=R1,op(1,algeq)),R1):
algR1 := collect(algR1,R1,factor):
<Text-field style="Heading 1" layout="Heading 1">B1 and B2</Text-field>
Recall that (1+t)^2*zeta*S'(zeta) = B1 + I*sqrt(3)*B2, where zeta = exp(2*Pi*I/3)
exp(2*Pi*I/3);
LCYjIiIiIiIjISIiKiYqJkYjRiReI0YkRiRGJCkiIiRGI0YkRiQ=
In order to ocmpute it we first generated the counting sequence for the real and imaginary parts of S'(zeta), i.e. S'(zeta) = A1 + I*A2. For this purpose recall that S(x) = t*x*M(0,x). Hence S'(x) = t*M(0,x) + t*x*M_2(0,x), and we are therefore interested in the walks A_{-1,i} for i>=0 and A_{-2,j} for j>=0 evaluated at x=zeta. In order to count them using modular arithmetic, we again multiply by 3 to have integers as initial values (recall the starting polynom at t=0: 2/3 + x^{-2}/3 + y^{-2}/3. Moreover, we use the relation 1 + zeta + zeta^2 = 0 to reduce the powers of zeta to 1 and zeta, i.e. zeta^{3*k} = 1 zeta^{3*k+1} = zeta zeta^{3*k+2} = -1 - zeta. Hence, the coefficients of S'(zeta) are elements of Q + zeta*Q. Hence, after multiplying by 2 we get that 6*S'(zeta) is part of Z + I*sqrt(3)*Z and we write 6*S'(zeta) = A1 + A2*I*sqrt(3). The integer sequences A1 and A2 have been computied in a seperate C++ programm using modular arithmetic and Chinese remaindering for the first n=2000 steps.
Let's load these numbers filename := "king-A1-2000.txt";
CCA1 := loadData(filename):
seq(CCA1[i],i=1..10);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSlmaWxlbmFtZUdGKFExa2luZy1BMS0yMDAwLnR4dEYoNyNGLg== NiwiIiFGIyEjPUYjISRNIyEkdichJW9lISZjOCQhJyNcPiMhKGUrUyI=
Divide by 6 CCA1 := CCA1/6:
seq(CCA1[i],i=1..10);
NiwiIiFGIyEiJEYjISNSIyEkRCMiIiMhJHkqISVFXyEmI2VPISdWTEI=
filename := "king-A2-2000.txt";
CCA2 := loadData(filename):
seq(CCA2[i],i=1..10);
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSlmaWxlbmFtZUdGKFExa2luZy1BMi0yMDAwLnR4dEYoNyNGLg== NiwiIiEiIichIiciI0kiI2EiJCRbIiVvRCImbW8iIichRzIiIicnNDwo
Divide by 6 CCA2 := CCA2/6:
seq(CCA2[i],i=1..10);
NiwiIiEiIiIhIiIiIiYiIiojIiRoIiIiIyIkRyUiJTZHIiYhKXkiIic7Jj4i
We are actually interested in B1 and B2.
First we compute the coefficients of zeta*S'(zeta) CCzA1 := -CCA1/2 - 3/2*CCA2:
CCzA2 := CCA1/2 - CCA2/2:
And then we compute (1+t)^2*zeta*S'(zeta) by a simple shift CCB1 := Vector(Dimension(CCzA1)):
CCB1[1] := CCzA1[1]:
CCB1[2] := 2*CCzA1[1]+CCzA1[2]:
for i from 3 to Dimension(CCB1) do CCB1[i] := CCzA1[i-2]+2*CCzA1[i-1]+CCzA1[i]: end:
seq(CCB1[i],i=1..10);
NiwiIiEjISIkIiIjRiNGJSEiJyEjZyEkdyMhJXU+ISYqKT0iISZrNyk=
This should be equal to (see B1ex(n) in king-weg-part2.mw) -3/2*t^2 - 3*t^4 - 6*t^5 - 60*t^6 - 276*t^7 - 1974*t^8 - 11889*t^9 - 81264*t^10;
LDIqJiMiIiQiIiMiIiIpSSJ0RzYiRiZGJyEiIiomRiVGJylGKSIiJUYnRisqJiIiJ0YnKUYpIiImRidGKyomIiNnRicpRilGMEYnRisqJiIkdyNGJylGKSIiKEYnRisqJiIldT5GJylGKSIiKUYnRisqJiImKik9IkYnKUYpIiIqRidGKyomIiZrNylGJylGKSIjNUYnRis=
CCB2 := Vector(Dimension(CCzA2)):
CCB2[1] := CCzA2[1]:
CCB2[2] := 2*CCzA2[1]+CCzA2[2]:
for i from 3 to Dimension(CCB2) do CCB2[i] := CCzA2[i-2]+2*CCzA2[i-1]+CCzA2[i]: end:
seq(CCB2[i],i=1..10);
NiwiIiEjISIiIiIjISIjISImISNJISRaIiEkPyohJUBiISZyZiQhJzVcQg==
This should be equal to (see B2ex(n) in king-weg-part2.mw) -1/2*t^2 - 2*t^3 - 5*t^4 - 30*t^5 - 147*t^6 - 920*t^7 - 5521*t^8 - 35971*t^9 - 234910*t^10;
LDQqJiMiIiIiIiNGJSlJInRHNiJGJkYlISIiKiZGJkYlKUYoIiIkRiVGKiomIiImRiUpRigiIiVGJUYqKiYiI0lGJSlGKEYvRiVGKiomIiRaIkYlKUYoIiInRiVGKiomIiQ/KkYlKUYoIiIoRiVGKiomIiVAYkYlKUYoIiIpRiVGKiomIiZyZiRGJSlGKCIiKkYlRioqJiInNVxCRiUpRigiIzVGJUYq
Guess the algebraic relation for B1
(wee add another 0 at the beginning, as there is another shift by t in S(x)) algeq := listtoalgeq([0,seq(CCB1[n],n=1..Dimension(CCB1))],F(t));
listtoalgeq: looking for equation of order 1
listtoalgeq: looking for equation of order 2
listtoalgeq: looking for equation of order 3
listtoalgeq: looking for equation of order 4
listtoalgeq: looking for equation of order 5
listtoalgeq: looking for equation of order 6
listtoalgeq: looking for equation of order 7
listtoalgeq: looking for equation of order 8
listtoalgeq: looking for equation of order 9
listtoalgeq: looking for equation of order 10
listtoalgeq: looking for equation of order 11
listtoalgeq: looking for equation of order 12
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZhbGdlcUdGKDckLGpuKiYsVComIjRfUnA4IzRJJWVJIyIiIilJInRHRigiI0NGNCEiIiomIjVnKHAlb2cvOiNIOiJGNClGNiIjQkY0RjgqJiI1N2hRUjNCcyMqekJGNClGNiIjQUY0RjgqJiI1Nyh6dElAZVBGZCNGNClGNiIjQEY0RjgqJiI1V0BhcCRRT2RlWyJGNClGNiIjP0Y0RjgqJiI0IyoqMz5ZeVRAT1BGNClGNiIjPkY0RjgqJiIyM2lDUGIwdSpHRjQpRjYiIz1GNEY4KiYiMyspW2pjTVdXdyRGNClGNiIjPEY0RjgqJiIzSyNmaEgzVENIJ0Y0KUY2IiM7RjRGOComIjM7a3IhenNlWTwkRjQpRjYiIzpGNEY4KiYiMjsheWk8LkopZidGNClGNiIjOUY0RjgqJiIxd2BTJVJgUCc+RjQpRjYiIzhGNEY0KiYiMVM9JHAmZT87UUY0KUY2IiM3RjRGNComIjBfej5uX3NAJ0Y0KUY2IiM2RjRGNComIi8hW2MoKSl6R0ZGNClGNiIjNUY0RjgqJiIvdzRSQjhdPUY0KUY2IiIqRjRGOComIi43eGVQQkMiRjQpRjYiIilGNEY4KiYiLXdwI3BQdyJGNClGNiIiKEY0RjQqJiIsJzRDIykqUSNGNClGNiIiJ0Y0RjQqJiIqO2M6RyNGNClGNiIiJkY0RjgqJiIqM0JhQSJGNClGNiIiJUY0RjgqJiIoZ241JEY0KUY2IiIkRjRGOComIiQpb0Y0KUY2IiIjRjRGOComIiQzIkY0RjZGNEY4RmpxRjRGNC1JIkZHRig2I0Y2RjRGNComLFQqJiI4P047eDMyTGRsIkc4RjRGNUY0RjQqJiI5Y2EoPSVHJnpyOXBmLiJGNEY7RjRGNComIjlnKls5Lyw7WTp1L20kRjRGP0Y0RjQqJiI5T0ZsRHUhXCJRdm9meEY0RkNGNEY0KiYiOkcqb3FNPnRzIypvZSw2RjRGR0Y0RjQqJiI6N1RKM0o0ZVQjUWoxNkY0RktGNEY0KiYiOXM9KnBjUUQpZik9ISk0KUY0Rk9GNEY0KiYiOVsjSCFSZVpdMTlPcFZGNEZTRjRGNComIjkjPkduZ0d2TzpIRXQiRjRGV0Y0RjQqJiI4W29PS2opR00pZU0jXEY0RmVuRjRGNComIjdfWkY2dSpwKnpbJ0cqRjRGaW5GNEY0KiYiNk86KUdPOnhFJz1eKUY0Rl1vRjRGNComIjVzOUt4em5cZz96RjRGYW9GNEY4KiYiNXd4dSc+IjRoUF9RRjRGZW9GNEY4KiYiNClbPlZ1SiUqM1dhRjRGaW9GNEY4KiYiMydcNTcmPkRlMkRGNEZdcEY0RjgqJiIybyhvejtIbCkzJEY0RmFwRjRGNComIjFLQ3deOyxJZUY0RmVwRjRGNComIjAzIyllWjk2NSVGNEZpcEY0RjQqJiIvR1RaT2ImUiJGNEZdcUY0RjQqJiItV1YuWHM9RjRGYXFGNEY0KiYiK0NNIypRNUY0RmVxRjRGOComIikjUmkxKUY0RmlxRjRGOComIigleWc/RjRGNkY0RjgiJTM+RjhGNClGXXJGanFGNEY0KiYsVComIjh3QEdyeUAoKikpPk0kR0Y0RjVGNEY0KiYiOTN1RHFIdXg2ZnRtQUY0RjtGNEY0KiYiOjdUJlxHRG4hPUo3MTYiRjRGP0Y0RjQqJiI6S0s5UTZXXXclXHRXUEY0RkNGNEY0KiYiOiV5NS5jImYzWSVmZSFmKUY0RkdGNEY0KiYiOykzWEZZdF5VRV9yPU8iRjRGS0Y0RjQqJiI7KW91OVY/eWlRdShbSTpGNEZPRjRGNComIjtXKCoqUmdga047anFZQyJGNEZTRjRGNComIjozVShIUyQqemRwR100dUY0RldGNEY0KiYiOiNmSTt4Km86UCRmU0lLRjRGZW5GNEY0KiYiOlc4MyJbUyJwJSopKUhuLCJGNEZpbkY0RjQqJiI5K19eI2V6Iz0sTzE6QUY0Rl1vRjRGNComIjg/WlRvNnUvZ1YtJ0hGNEZhb0Y0RjQqJiI3N1hlUkRcNTt0UjdGNEZlb0Y0RjQqJiI2X2JgIyo0Im8jUSd6SkY0RmlvRjRGOComIjUhWyUzUzMoKSopbydwJkY0Rl1wRjRGOComIjQleXlnN2UnZnA7IkY0RmFwRjRGOComIjMrW2NePU4xNHVGNEZlcEY0RjQqJiIzIVtPJ2ZbcXZQNkY0RmlwRjRGNComIjFnYCpbRUVOUSlGNEZdcUY0RjQqJiIwIVsjWyUqUW5nJEY0RmFxRjRGNComIi53TFJ6MD8qRjRGZXFGNEY0KiYiLWN5M2RSNkY0RmlxRjRGNComIio/KiozJGZGNEY2RjRGOCInP3ZrRjRGNClGXXJGZnFGNEY0KiYsUComIjlrMylcRjMxQj9OeSwnRjRGP0Y0RjQqJiI6Nyo9ZGhzNT9IJDNSPSVGNEZDRjRGNComIjtPZnhjZillNG4hZWtrOkY0RkdGNEY0KiYiOyU9cCFSXWV3KSoqekAnR1JGNEZLRjRGNComIjs7J1JhMlwmPiVcUndPIXBGNEZPRjRGNComIjt3YDtfOElxZWV4ZT4nKUY0RlNGNEY0KiYiO2cocCEpKnpcI3o1cm03dihGNEZXRjRGNComIjtzcUVKajFaInk2TFgyJkY0RmVuRjRGNComIjtHJD5zUkBmdilccGpIQ0Y0RmluRjRGNComIjpbS1RaO3FFZSMpKmVvJSlGNEZdb0Y0RjQqJiI6M00ybCpSPzBVUCw0QEY0RmFvRjRGNComIjlTYWNcKEg6Kls/aSplJEY0RmVvRjRGNComIjhjdWkhKT06JkgmKSo+dyRGNEZpb0Y0RjQqJiI3Zyw5OXE4MWFUMT1GNEZdcEY0RjQqJiI1R3RBMTJjMWhDVEY0RmFwRjRGNComIjVfaiQzLWlWXlpXJEY0RmVwRjRGNComIjRXVCQ+JylIU14sIipGNEZpcEY0RjQqJiI0Iyoqb1UkNEgjeik0IkY0Rl1xRjRGNComIjJXIm9PSFs8OnZGNEZhcUY0RjQqJiIxa19cYz9vN0lGNEZlcUY0RjQqJiIvT3oiWydIJnomRjRGaXFGNEY0KiYiLFcqKSo+JHAiRjRGNkY0RjgiKSMqSCgqKSlGOEY0KUZdckZicUY0RjQqJixMKiYiOkNNKTR4Im9hIT1MbikzJkY0RkdGNEY0KiYiO186dmBYVFdrWCEpXFlJRjRGS0Y0RjQqJiI7b15gISpSPE0nKj48XSVIKkY0Rk9GNEY0KiYiPCMqNCk0Oy8oeVxCYFdEJz1GNEZTRjRGNComIjxjRW1hJ1E8PT9gNkc4RUY0RldGNEY0KiYiPENhbDZqJSoqKXk/SWsjPkVGNEZlbkY0RjQqJiI8R3BuKipIR0VnLGJvayo9RjRGaW5GNEY0KiYiO2cqW3FJMjhrLnBbbSkqKkY0Rl1vRjRGNComIjsjPnNDISpSXm5EQFYjSFFGNEZhb0Y0RjQqJiI7L24mKXknKT1ZMGpoM2o1RjRGZW9GNEY0KiYiOlNNTEViTSxoMyZmNUBGNEZpb0Y0RjQqJiI5d3BnP0NZSyZHZ24nSEY0Rl1wRjRGNComIjghM1l3OEM1S3dCLkpGNEZhcEY0RjQqJiI3X2ZHayZSdDVMWEAkRjRGZXBGNEY0KiYiNiEzcTFQdlcrWDdYRjRGaXBGNEY0KiYiNSU9VEQqPmhdKG87J0Y0Rl1xRjRGNComIjRrU0lEKikpPic0JGVGNEZhcUY0RjQqJiIzN0hdTSJlcixIJEY0RmVxRjRGNComIjFDSWomNCpbMSYpRjRGaXFGNEY0KiYiLis7KSlmXDwmRjRGNkY0RjQiKyM+KWUlZiRGNEY0KUZdckZecUY0RjQqJixIKiYiOzdoYSQzKTMnek9LPC9AI0Y0Rk9GNEY0KiYiPD8+Vy4oKVwpenJjZTxPNkY0RlNGNEY0KiYiPGNlI0cuIjRkdi8lb01vR0Y0RldGNEY0KiYiPCUpKiopWzw1S0dNKSpROyJvJUY0RmVuRjRGNComIjxXRFR5KHokSC1pPkl2SyZGNEZpbkY0RjQqJiI8S0NXXkRcOmE3OGB6SiVGNEZdb0Y0RjQqJiI8LzY2WzlIU2QpZThzOkRGNEZhb0Y0RjQqJiI8NzRyY3RTJXA4KVFfKmU1RjRGZW9GNEY0KiYiOzcqKSpIVjZIW2BzZGpBJEY0RmlvRjRGNComIjo3XCh6M3lVTU1RUExyRjRGXXBGNEY0KiYiOlslKSlbLDBoZUxOP2w2RjRGYXBGNEY0KiYiOXM1bDpmRFZoSk41OkY0RmVwRjRGNComIjgnXF1xLShwc3JjZyI9RjRGaXBGNEY0KiYiNykzUDBXXSs2JClHQCNGNEZdcUY0RjQqJiI2JXk9QEIkMzpoS0sjRjRGYXFGNEY0KiYiNWdcWTJrKCk9XEY7RjRGZXFGNEY0KiYiM0MhKjM2I2VKM1UmRjRGaXFGNEY0KiYiMFNZIylcejQ0IkY0RjZGNEY0Iiw/WjB0WyZGNEY0KUZdckZqcEY0RjQqJixEKiYiOyUpekQ8dWMvIXp2TFxUJkY0RldGNEY0KiYiPGNRdnBPXz5JIXp5cy1DRjRGZW5GNEY0KiYiPDN5JypRWWkpeSJ5SGVTMCZGNEZpbkY0RjQqJiI8R2RxYDctOkA+NVNedidGNEZdb0Y0RjQqJiI8L0ojUSEpKilwKlwlW1RIW2lGNEZhb0Y0RjQqJiI8R0AyJGY7RWFgXzFYd1NGNEZlb0Y0RjQqJiI8IyoqR3krdHEzSkY3LiIqPUY0RmlvRjRGNComIjslPSRITlMjSDV1WW55RidGNEZdcEY0RjQqJiI7P1ZxUlA8JGYoM0B4MjpGNEZhcEY0RjQqJiI6azNlJ3phYCFRTnhwcCNGNEZlcEY0RjQqJiI5Q3kuJlEscEUncCVSJlFGNEZpcEY0RjQqJiI4Z0whUnhtcTRNPFVbRjRGXXFGNEY0KiYiNy9WVjlDO1w2KD1PJkY0RmFxRjRGNComIjZbKVFlJ1wzTCFva1ZGNEZlcUY0RjQqJiI1MzFVZyJRJj0oUSM9RjRGaXFGNEY0KiYiMXdUSnIhKSozeSNGNEY2RjRGNCIuZyg0QjU4U0Y4RjQpRl1yRmZwRjRGNComLEAqJiI7Ky8qPilcclUlKWZwaWR6RjRGaW5GNEY0KiYiPG9ObFMkNFVWMyQ+J1JcSUY0Rl1vRjRGNComIjwzXTZ6PFpLTyYqSDMjPmBGNEZhb0Y0RjQqJiI8ITNVeDJTdTBrRHIpR3gmRjRGZW9GNEY0KiYiPCtjWXAhNFk5IltxRV1IJUY0RmlvRjRGNComIjxnPDF2VSp6NCp5SkohPUFGNEZdcEY0RjQqJiI7KVsyaSdcRz1JPzM1SiEpRjRGYXBGNEY0KiYiOy82UipvWD47LlAqSHA/RjRGZXBGNEY0KiYiOidwOCUpW1IpNFZyR2AkUkY0RmlwRjRGNComIjlXKj1DcS9mMGIkPiMpZUY0Rl1xRjRGNComIjhjWSh6KXAhUWtBVSNHKEY0RmFxRjRGNComIjdzUWVwd0dBJ2VdJW9GNEZlcUY0RjQqJiI2N2g9O2ZXPSMpeWYkRjRGaXFGNEY0KiYiM3ciPVNbdl1tLSNGNEY2RjRGNCIvM1k/KFsrJFtGOEY0KUZdckZicEY0RjQqJiw8KiYiOyV5I3loeFVaKFJKJ3khPihGNEZhb0Y0RjQqJiI8T0Z0Lj5lInByJFstWk8jRjRGZW9GNEY0KiYiPCc0L1c8Oj4nR2Jnd2tOJEY0RmlvRjRGNComIjwlPVAlKilbPj4peWxKNSMpR0Y0Rl1wRjRGNComIjxLIW9RNGNLSyR6eit6byJGNEZhcEY0RjQqJiI7VypRQjpyNkczTE5LcSdGNEZlcEY0RjQqJiI7I2ZJY0JfYF1RJ1s5Yj1GNEZpcEY0RjQqJiI6Xz41QidbbmJmT2AmbyRGNEZdcUY0RjQqJiI5T3IneTUhKjNpJilRITRjRjRGYXFGNEY0KiYiOClHS153PHpUTiZwSSdGNEZlcUY0RjQqJiI3I3pAKHksQkljaj9WRjRGaXFGNEY0KiYiNEcsUyEpXCFvcyRvJkY0RjZGNEY0IjEhWztYYk0uVCJGNEY0KUZdckZecEY0RjQqJiw4KiYiOyNSWj49TSMqKSk0ZT5jI1JGNEZpb0Y0RjQqJiI8JXlRJ1IiKWU8PGgiUUUkNCJGNEZdcEY0RjQqJiI8RyhIOlQ6bjdMSHUkUkAiRjRGYXBGNEY0KiYiOzMtTkcqWyk9XFM1J1EieUY0RmVwRjRGNComIjtfPnUmcHhlWXRiUXNfJEY0RmlwRjRGNComIjslZT5XR1hSIT1rd3pRNUY0Rl1xRjRGNComIjpnRCsjZmxaI28kPlJAQUY0RmFxRjRGNComIjlLX3gqUnBCTGtbIik+JEY0RmVxRjRGNComIjgvJGVpJHlbSikpby02JEY0RmlxRjRGNComIjVbSzhZJXlzcVo6KEY0RjZGNEY0IjJnRCs7QSEqZUskRjRGNClGXXJGam9GNEY0KiYsNComIjtfSlYpRzV6ViI9USkzPiJGNEZhcEY0RjQqJiI7W2dFPkVeJVtwRyU9XEZGNEZlcEY0RjQqJiI7Oz91bVxNLyUpKipvaDFBRjRGaXBGNEY0KiYiOisnNC8rXypwXCUpeWpYKkY0Rl1xRjRGNComIjonSDZCKnBaVmk1PHRkJEY0RmFxRjRGNComIjkpR29tR0xPUG4/T0QoRjRGZXFGNEY0KiYiOT9uSz5AYE5CJVtLQiJGNEZpcUY0RjQqJiI2KUdDaCMzJSopemloVUY0RjZGNEY0IjM7UzApUXliKClcI0Y0RjQpRl1yRmZvRjRGNComLDAqJiI6JypHVldFdUsoXEJvVzpGNEZpcEY0RjQqJiI6X2pKeG4+LihIayYqZkdGNEZdcUY0RjQqJiI6TyZwMkdZUzE1c3VyOEY0RmFxRjRGNComIjlPPmJjcU9bbWdvPUVGNEZlcUY0RjQqJiI5d1gxTkg9YkY3Zm8/RjRGaXFGNEY0KiYiNldUPGIlMywoUjgnKSpGNEY2RjRGNCIzI3pUJ29LUkJmbEY0RjQpRl1yRmJvRjRGNComIi1nJD5GVjwjRjRGZW5GNEY0KiYiKGtvIkhGNEZlcEY0RjQqJiIsK2NdazwjRjRGYW9GNEY4KiZGZnFGNEZpcUY0RjQqJiImKW9CRjRGXXFGNEY4KiYiLzdMYFJUPjhGNEZHRjRGNComIi9fKlxSIm9vSEY0RktGNEY0KiYiL1tvJUckKSlwO0Y0Rk9GNEY0KiYiLkslKTQ5O1MkRjRGU0Y0RjgqJiIkaSJGNEZlcUY0RjgqJiIuU0EhUSQpb1BGNEZXRjRGOComIiolUWdcNkY0Rl1wRjRGOComIiwrO0NcRSNGNEZdb0Y0RjgqJiItUyU0dCoqeSRGNEZpbkY0RjQqJiIoTyoqPilGNEZhcEY0RjgqJiInP0o3RjRGaXBGNEY4KiYiK1M9TCQ9I0Y0RmVvRjRGNComIip3PGAicEY0RmlvRjRGNComIiVuS0Y0RmFxRjRGNEkkb2dmR0YoNyNGLg==
computeParams(algeq);
NyYiIzciI0MiJEQkIiRIIw==
algB1 := primpart(subs(F(t)=B1,op(1,algeq)),B1):
algB1 := collect(algB1,B1,factor):
Guess the algebraic relation for B2
(wee add another 0 at the beginning, as there is another shift by t in S(x)) algeq := listtoalgeq([0,seq(CCB2[n],n=1..Dimension(CCB2))],F(t));
listtoalgeq: looking for equation of order 1
listtoalgeq: looking for equation of order 2
listtoalgeq: looking for equation of order 3
listtoalgeq: looking for equation of order 4
listtoalgeq: looking for equation of order 5
listtoalgeq: looking for equation of order 6
listtoalgeq: looking for equation of order 7
listtoalgeq: looking for equation of order 8
listtoalgeq: looking for equation of order 9
listtoalgeq: looking for equation of order 10
listtoalgeq: looking for equation of order 11
listtoalgeq: looking for equation of order 12
listtoalgeq: looking for equation of order 13
listtoalgeq: looking for equation of order 14
listtoalgeq: looking for equation of order 15
listtoalgeq: looking for equation of order 16
listtoalgeq: looking for equation of order 17
listtoalgeq: looking for equation of order 18
listtoalgeq: looking for equation of order 19
listtoalgeq: looking for equation of order 20
listtoalgeq: looking for equation of order 21
listtoalgeq: looking for equation of order 22
listtoalgeq: looking for equation of order 23
listtoalgeq: looking for equation of order 24
LV9JLFR5cGVzZXR0aW5nRzYkJSpwcm90ZWN0ZWRHSShfc3lzbGliRzYiSSxtcHJpbnRzbGFzaEdGKDYkNyM+SSZhbGdlcUdGKDckLGZzKiYiOmMnKXo+OiozK1sneSRHdyIiIilJInRHRigiIzpGMkYyKiYiMGcvIW9aczpFRjIpRjQiIihGMkYyKiYiTV9OcTZiJCpvN3ArKHpKdk8nZSNmQThvI1xGMilGNCIjUEYyISIiKiYiTiVRU1FkSWZlSF0oKSkzPC1QPGh2eCYqZlUiRjIpRjQiI1FGMkY+KiYiTjtvZSp6RipcKmZoQSQzWmhUKlthTyRRJyo9IkYyKUY0IiNSRjJGPiomIk5DJykqPWZPWWkoZWokZkQySipIZicqPm8xUCdGMilGNCIjU0YyRjIqJiJPIVtjLj51UzVEP0k3V3M/XW5uSGJoOiRHRjIpRjQiI1RGMkYyKiYiTyFbXyEpZTtIOG9IbysjW0E5Qy0hR3RfcDYmRjIpRjQiI1VGMkYyKiYiT2dsNTIleSM+SWw7KGZQO1JoaTVldFJuPiJGMilGNCIjVkYyRjIqJiJQQ2E4b3llNUFbKWU7QmZVNm9xLkwkKXBePEYyKUY0IiNXRjJGPiomIkJjXVdlYkpzOCoqR2loKik+WCUqKUYyKUY0IiNCRjJGMiomIjc7LzVpKCo9eClmOGQkRjIpRjQiIzdGMkYyKiYiQSNSIlx5WEAwIT4mPWNaRDwpUSpGMilGNCIjQUYyRj4qJiJBS1s5bFdwRmMrM2NvJypRSjlGMilGNCIjQEYyRj4qJiIsXz1GKXk2RjIpRjQiIiZGMkYyKiYiP1ckcCI9KSpmWlZ6ZW0sNyo9KkYyKUY0IiM/RjJGMiomLF5yKiYiUE8yI1t1Y1tRaj5KOXRSKXktViYpZjwoKjRaIkYyKUY0IiNjRjJGPiomIlFbI0hZPVMuU1pvVF1xI0hoVCJIdTcoKSopb2EiRjIpRjQiI2JGMkY+KiYiUSNmIT55d0ZsVyRHZUEuOzpPSVdbMSwqeT53RjIpRjQiI2FGMkY+KiYiUilHY0g/TGYwYFIlXEUvU2g1eCIqKlwzeTBCQkYyKUY0IiNgRjJGPiomIlJvYidSKG8uZixWMFtEeSZ5REhWKmVqOnEkUVtGMilGNCIjX0YyRj4qJiJSP15XJHB2I2VYKFxlLWEtajBBIz5TbW03JnBxRjIpRjQiI15GMkY+KiYiUiVlQio0L09zZV1gNXpzPSVwMnNMXFMmKioqZXBGMilGNCIjXUYyRj4qJiJSXzYlZSVIKVxfKW9RQkhndypwUT5yazEsSlhPRjIpRjQiI1xGMkY+KiYiUnM9XlFmO1toJnApb1o6Im9WdyI0JEhFKT54KyJGMilGNCIjW0YyRjIqJiJSITMhUkgkKVFHaycpKWVMW2J3X3JcPShSPyV5K1FGMilGNCIjWkYyRjIqJiJSOykpcHhIRl16IypbIkgvJls/KjR0Jz1qPjAkUUxGMilGNCIjWUYyRjIqJiJSITMjKSo+XGVpXT4pWylwJylHKlxuP0Q7NnAiWzUiRjIpRjQiI1hGMkYyKiYiUVdEdDFcSDFfbFJjSndsayRcZzUsPTElemZGMkZZRjJGPiomIlEzdSgqKUdOeTVbR3ZAam9gUmk/LWZPJ3l5ISpGMkZVRjJGPiomIlEpR0s2dT9FQm5nSWZhczNBT1ksTTMsaVglRjJGUUYyRj4qJiJQcyl6VyxaXm0kUiUpeiMqPUd1NWtndkBoI3lCRjJGTUYyRj4qJiJRVzhPOWhHLDAqKUdiKioqZT42J29iXzAsajg1RjJGSUYyRjIqJiJQI3pYIj1dMk5mJyoqKVJcNyMzMTopSDkuZyZHakYyRkVGMkYyKiYiUFN1JlwpKnpYOHpheUF5MFo0XEglKilmWiJbOEYyRkFGMkYyKiYiT19abiY9MXJBV11wJlJfW1h4XnQnej40VyVGMkY8RjJGPiomIk9vSmkiUjYlKUdraUh0JFFHM146XV9HOU9VRjIpRjQiI09GMkY+KiYiT2szXSVwYD1YJltqMS1xKipRZl8wRCVIaUoiRjIpRjQiI05GMkY+KiYiTUNBLCQqUk1CPlNGOiUpMzQoPk5cWUNmLyhGMilGNCIjTUYyRj4qJiJNT2ZgRSVlJ1FrK2V3O2dAYidIPF80RGoqRjIpRjQiI0xGMkYyKiYiTSs3SCdRbDMkcFk+S09RNyhcTWhIdnVzJEYyKUY0IiNLRjJGMiomIkxjZScpPjpjTCQ9LTU7dXphYCE+ZUI+WiVGMilGNCIjSkYyRjIqJiJLJzQnejBeNlInRzt5PV8qW2dIYC5lYigqRjIpRjQiI0lGMkY+KiYiS0NFY0g8dWB2Q3BBRjJjZWhaR3QnUSVGMilGNCIjSEYyRj4qJiJKLygzb0J2WHhIJlwzPjgiZkRbJ2YpKUclRjIpRjQiI0dGMkY+KiYiSSFHNFkhUiJHRHZTJXAiXGZXWXV3aXgqRjIpRjQiI0ZGMkYyKiYiSVdEND1ZZ3gqRypRYVd1MkVja2ZgSEYyKUY0IiNFRjJGMiomIkdPJnB2MiYpeXQjNGMlMzNDOy9VJFt2RjIpRjQiI0RGMkYyKiYiR1snUkRqKipmIillM2shKUhjbl9TNl0oRjIpRjQiI0NGMkY+KiYiRiNmbU5nRjN5JClHdmxMY1QqZXUtKkYyRmduRjJGPiomIkYzaSlRJEd4dTR1UiMqMyE0NSRHdiwiRjJGX29GMkYyKiYiRUszTz1tJHBwIyopZilIPyl6OSRvYyNGMkZjb0YyRjIqJiJDaz9eLipvdFldWU83NnBjJD5kRjJGW3BGMkY+KiYiQytjeWghcChSP0FZLG53KTNqbiVGMilGNCIjPkYyRj4qJiJBY1Vxb3dKIj1oOGI/ZGlNXCZGMilGNCIjPUYyRj4qJiJBczEnMykpcGNfdSF5IkdOanpdJ0YyKUY0IiM8RjJGMiomIkB3Ink8KTM9dipcIz06JlFlcDlGMilGNCIjO0YyRjIqJiI/K2dkWyFlPjh5MkIiUSUpPnVGMkYzRjJGPiomIj4pM0BEc3gkKT0leUwiUjlwNkYyKUY0IiM5RjJGPiomIj1DcTxxOVtFMSRwJUdNO3BGMilGNCIjOEYyRjIqJiI7b2ZaUi8ycU8+L2snWyRGMkZbb0YyRj4qJiI7ZyZRcnpuJj0uZ1xENltGMilGNCIjNkYyRj4qJiI6U1lbLTwhUVRQVG0uPEYyKUY0IiM1RjJGMiomIjlPLlo5TytZJFsheUE+RjIpRjQiIipGMkYyKiYiOGNFbl5tNyJbM1JsOkYyKUY0IiIpRjJGPiomIjVTJ28jUi14WE8iPSJGMkY4RjJGMiomIjVnWEIicCE0Jno7cSVGMilGNCIiJ0YyRjIqJiI0KUcpbyFIRFYuXEVGMkZnb0YyRj4qJiIybypHXypmdVRFKEYyKUY0IiIlRjJGMiomIjFjXVwvN2hDNkYyKUY0IiIkRjJGPiomIi9DITNKVC4rIkYyKUY0IiIjRjJGMiomIiw7JypbRHolRjJGNEYyRj4iKXctYicqRjJGMiktSSJGR0YoNiNGNEZqW2xGMkYyKiYsZnEqJiJRayN6YlsuKj0vN2tzXyo+U3peNTY+T1B4cyRGMkZhcUYyRj4qJiJSY11nITQwImVrakUnXEQkeiU9Lyg9ME9OdmApSEYyRmVxRjJGPiomIlMvViRbP3BUIkhOYTtsOSJHPyZRKVF6WngqeWI2RjJGaXFGMkY+KiYiU1spKSoqSFRJSy1hKm87V11scVszL0NRZkVueCNGMkZdckYyRj4qJiJTV2w2XnViJ3A2KzxgaE1NZiIpRzJgMmlkJFxXRjJGYXJGMkY+KiYiU192W3ljOWpiV3hnSzNLJ3k8dHluJ29hKVtvJUYyRmVyRjJGPiomIlMncFgpPTdcQ25qUiczIj4hKXpJbkNbIj5eeituI0YyRmlyRjJGPiomIlIlUTtPYmFDMEFTYl89MUBZJVFrXzdfNi5MJUYyRl1zRjJGMiomIlN3VGJwJT1SJj09dkgxKWUrU3RGKFskPTNKPU0jRjJGWUYyRjIqJiJTbzJbN3IlPjJlTDQrTUNSJylvWChSQz9ZTyYqPkYyRlVGMkYyKiYiUkdoMHZna1dILDBPViZvcCt0b3MleUJETHIlRjJGUUYyRjIqJiJSX2pmXXkpKj1DSkA5X1x6JCl5IUdLczomPSU+Y0YyRk1GMkY+KiYiUksrIykqXFwheiQzKnkqcF9kayl6dkQ0b0w2IzNmRjJGSUYyRj4qJiJSa2tDNGgySCw3K1BfMCdwYjMtYz1NL2A4PkYyRkVGMkY+KiYiUSU9UCMpemYxKnAoKUgqbzVFZ2khPVdSYjRZO2VGMkZBRjJGMiomIlFnIjN4KXlPKCopWyJvUycqWydIYCtkJFxxTGcobyhGMkY8RjJGMiomIlE3VEYkcCRSMV95UixYYyo+IzRNdXcwc2drQkYyRmF0RjJGMiomIlAhRyhIdmUkSHhTVTlaPzNOZSQpKSlITUNUKUclRjJGZXRGMkY+KiYiUDs3M2ZhPUtTVmckMyRvbDVGNS44SSJvSWRGMkZpdEYyRj4qJiJQZ1gmZkNrSD1FdWVrWC89U1klZT9NNikpZiJGMkZddUYyRj4qJiJPLzYuXXlhOlc4Mm9UL08qKnlaSFAkbyg+OEYyRmF1RjJGMiomIk83QDcjSCopPkxcNkNfdyo0J0hyN1x3I0hnP0YyRmV1RjJGMiomIk4hWyswNVUkKUhXKmYmXF8iek9jbm0hR0dCWEYyRml1RjJGMiomIk1DWShvPzNEYlciM2Yhb1MkM2dLWXBkKVwmRjJGXXZGMkY+KiYiTSs3KlE/LilbJnkvcjo7ZWhvYWhzeSI+WkYyRmF2RjJGPiomIkxbJ1IseDRDZVIpKipSYUBdUFdGPCc0IT4oRjJGZXZGMkY+KiYiTEdcJ2YuW3YqUnlwIlFpO0VZPiFmc2g0IkYyRml2RjJGMiomIksleTEvbW5HIT5VQ1V1QSwzSmEjWzk9JkYyRl13RjJGMiomIko3TFRhJWV5YVNeQTpZOXEnZkRLekIkRjJGYXdGMkYyKiYiSnMpKlJfRFEtaS1BaFA2aSlHZ0xCMjlGMkZnbkYyRj4qJiJJWzcnUSNHJEdFNllSVzllSHo5V2NbI0YyRl9vRjJGPiomIkghR0gsaSlwKD5VXUJkJCkpeSdITDUhZiJGMkZjb0YyRjIqJiJHO0siPiJHQjtnbyczNlkoZlJ2eW5wRjJGW3BGMkYyKiYiRSNSWi45R0siKlFjXDlzbEI1MiFlRjJGXXhGMkYyKiYiRl8mZiJmQ1wyKHBAPT0sZk1XNk43RjJGYXhGMkY+KiYiRCM+JWUlZScpUlI+K3pxV0U0VmdWRjJGZXhGMkY+KiYiREtbMUZBZ2RBYD9ZZzllQ0RrIkYyRml4RjJGMiomIkIhR0h4UkhqWigpSFRzLnZEKSpmRjJGM0YyRjIqJiJCXy5NLUNxWDozKillcCFHQW48RjJGX3lGMkY+KiYiQCV5ZSZIaThMWmdQS0giW29FRjJGY3lGMkY+KiYiQCMqcCJcZShHLUlWdyVbZGdrOUYyRltvRjJGMiomIj5vKG8iUmRlYGkqPiY+LWA8JEYyRml5RjJGPiomIj0lKSpwKGZ4a1JfKHlid2E+KEYyRl16RjJGPiomIjwzVSZ5WHZLIVEuZmNCRSZGMkZhekYyRjIqJiI5OyVHeFRlVlZaXiVbTkYyRmV6RjJGPiomIjojZiMpZnNvZDFuPllaPkYyRjhGMkY+KiYiOVdQTSo0cVpiJmUpKSlbIkYyRltbbEYyRjIqJiI3R0wlW0FJQ2RgcVsnRjJGZ29GMkY+KiYiNktPMj54cidwXzQ9RjJGYVtsRjJGMiomIjQhM2FPVDZwIjQtJEYyRmVbbEYyRj4qJiIyMz1WYSFwdF9ERjJGaVtsRjJGMiomIi9bX1NbRDkpKUYyRjRGMkY+IixrIWVMciEpRjJGMilGX1xsRmJbbEYyRjIqJixecSomIlIvNkoySGhlI1wrcDMtJlxsQFN4OFFuTk5AJUYyRmFyRjJGPiomIlNTaUk+aHZDVSgzIls6aklSP19KWW43YDdoRkYyRmVyRjJGPiomIlMleSN5XD9uZTVpZmtSOzwmKm9zKUhoUjBpZF8pRjJGaXJGMkY+KiYiVHdUOlY0bnZseUElNDNYMGVjVCZRRj5cJXoiZiJGMkZdc0YyRj4qJiJUZ1RrKmZfIWU5ISlmOT0oKnknUillJGZXKj0kPW0iKj1GMkZZRjJGPiomIlQvPkgkUXJIUCtfSGZxbzpdenUrTiR5JG9ZI0ciRjJGVUYyRj4qJiJSIVtXNStLI1spW2FVMzJmemtEXDtuNzRoSicpRjJGUUYyRj4qJiJTR0wpW0xyWS5oKCopNCdIaEJmPFkrdjJAVyRveEYyRk1GMkYyKiYiUzcoPidwW0cocCF5XkAlKVtgeT1baHBVX3lVVnZGMkZJRjJGMiomIlMzeV80VVxkWTBMQjNANzlQbGZ0JlJ4IW8+PkYyRkVGMkYyKiYiU2MxMj9eJG9wKClSTjdoWzxDbklGY21vTDU+I0YyRkFGMkY+KiYiUztfPSd5XVl0NTA8Zz4td0FwYW4jKm8iXFkqRyNGMkY8RjJGPiomIlI/diFcX0lzTFlnTjxBISpvaGxyJnooZTtgNm1GMkZhdEYyRj4qJiJSN0hNZVFacykpNEopM20qUilwSWFcV0kmUjMnR0YyRmV0RjJGMiomIlIlUSVRXV1fTihbcjheJD1PQzAiNFdiSC9hP0hGMkZpdEYyRjIqJiJRNzBjJSpmIm8lNCNvdTE/dHljKTNYOnAmKillYydGMkZddUYyRjIqJiJRISkpZUtfXUcqcFdmOWlqQFRTKCl5MHE4ajNFRjJGYXVGMkY+KiYiUW9aTWxqMkxQdjl0cnlzPTFiRHdyMUFhPkYyRmV1RjJGPiomIlBzNWRoO2dfKlJhS3dDKyhlM00/Ij0mKlFoTUYyRml1RjJGPiomIk87KzYqRzNCIlJwdVsqUSpIYipcYzIhZjR1IilGMkZddkYyRjIqJiJPS3NjJj0jcD5AXm53JVwuXE1mQCJwZVA+WEYyRmF2RjJGMiomIk5vRl4zYiJbLnosYXkrLm1tS0JvO0peV0YyRmV2RjJGMiomIk53aFVYJ2UpKilbNyw1KDMqKj5XKVJxSkxhOCJGMkZpdkYyRj4qJiJNd2xcLDEkKnBWaUdDNyE+UUdlb2ZVdG4iRjJGXXdGMkY+KiYiTCsnKm8mZXZEKVx5XyUpeitJRkRxJlJ1XmJGMkZhd0YyRjIqJiJLW2siZTo6WlsjXGpaKSk0IypvJmZ5UmZGKkYyRmduRjJGMiomIktneE5tZVo9UWNGYDYybXJvM1ojUU8jRjJGX29GMkY+KiYiSls7UFRDNyo+NjkncCN5UCwwJDRSOEMmRjJGY29GMkY+KiYiSTciPXlQTS9nVW4nUmAqUWsqNGglNFYmRjJGW3BGMkYyKiYiSTdMcChHXlRFSWYpeWprSGFePzM7PUYyRl14RjJGMiomIkczMTluXzMiSCpSbmJGPjE1NChmXyZGMkZheEYyRj4qJiJHX0IqZUFDPms9ZzNAWGJwInA9JXkkRjJGZXhGMkY+KiYiRS92XiFcJ284USIqZj9jTDliT1Y9RjJGaXhGMkYyKiYiRT8iNF4qRzwxJyoqPlxEbSopPnU5OiZGMkYzRjJGMiomIkMzMUFIbD9tcz0lcHAlUnVAdUxGMkZfeUYyRj4qJiJDc2kvXDB5KD1PcXFyZSZcMDpZRjJGY3lGMkY+KiYiQiVlOlxAOj1VNSkzQkchKTMqcDdGMkZbb0YyRjIqJiJBKTMoNGJQMkdgcDJYUSIpcERARjJGaXlGMkYyKiYiQFtvUyJwOD8wS3ksVHpMRDlGMkZdekYyRj4qJiI+PyRbUWYnPVptJyp5RHYnZUtGMkZhekYyRjIqJiI8X15pSzxaTT5PeSRHRCMqRjJGZXpGMkYyKiYiPFN5MWNMak09JHpBLVtcRjJGOEYyRj4qJiI7X0pmRycpWzpUMCxwLmRGMkZbW2xGMkYyKiYiOmcwQCNbMzlKRi5LMUpGMkZnb0YyRj4qJiI4JSlSZiw0WUVvWDxRKkYyRmFbbEYyRjIqJiI3XzI8bzczYnA8WztGMkZlW2xGMkY+KiYiNSsraydbeid5PWQ4RjJGaVtsRjJGMiomIjIrSyFvOUElUnEkRjJGNEYyRjIiL29mJG9YcVYlRj5GMilGX1xsRlxbbEYyRjIqJixmcComIlNDeVZacy9sSHJ6OiJ6ZURTWjRbNE9kJj1wREYyRllGMkY+KiYiVGM5aEFBKHBqPHh3UlN6cUhOcWhCOFNrOTk5RjJGVUYyRj4qJiJUbyp6OEZ4bSd6KkhYM08hNDNVZz8hKnkuKVtBaU1GMkZRRjJGPiomIlRvcncqKVImXFAvMypwXmlGLUEhemIyIXopUjRvW0YyRk1GMkY+KiYiVHcqKVIlcCsnUT0jMzR1TldnenZxemUwSENMPFNGMkZJRjJGPiomIlRrIz5xRSlSRFMqZil5J1tKamd5dUorW0VKIWVDIkYyRkVGMkY+KiYiVEtPbFJOPillJlFGIj4peXRbQ3pJZHY6JSlHL0oiRjJGQUYyRjIqJiJUO1czUHYpUiU0dXgsZF8lNHdhNWpkOkhPIXAlPUYyRjxGMkYyKiYiUzMjPVpeLmAsUURvJG9IaEZhS1cmb2EkZkEqeSdGMkZhdEYyRjIqJiJTR0QtbzpDN05cKXBmMnB3YVcxeSRIPkd5VVhGMkZldEYyRj4qJiJTc21uLjFbPUZVayd5QlozUEIlMzoncEQpUSgqZUYyRml0RjJGPiomIlMnKm9kLXkrISpRIltPK3gzJT0xR0AtMXYyUVU9RjJGXXVGMkY+KiYiUksvPEByJz4zczNeSzAqKXoqeV94KVsuV2BSdUYyRmF1RjJGMiomIlJLbyo9MGZdUVQlUjQwOmwsJDNjTTIrMypbYihGMkZldUYyRjIqJiJST3pzJWZ1aG9VbiZwXTVGZ0xQNV1WPXRGSyJGMkZpdUYyRjIqJiJRLz44clhTUm9LKSlwXGskKnk1ISlwbG4sW0YlKkYyRl12RjJGPiomIlFDSVIiSDA3VmtnXUsvSUBzPFwiNEozJzRVJkYyRmF2RjJGPiomIlApWz4mKlxDIW96Jkd4Z1V6N19bIW9CK2xiWCVGMkZldkYyRj4qJiJQT25OKUgnNF55WSQpZndBTSNvSyZISic+NTRdRjJGaXZGMkYyKiYiUCtLbyJ5XjpcVyVwKGUoKSk9Xil6ZjElPU1ZeCJGMkZdd0YyRjIqJiJNVyplXUxLSCRbXXNqST9HRmhLeEpnUWZGMkZhd0YyRjIqJiJPU3ViTU03LTc4cXRYSjspZTBaUTQnNFc3RjJGZ25GMkY+KiYiTmMneTE1IzMocFAmM0Moenc4XVJJRXo8NENGMkZfb0YyRj4qJiJNJVFXYC47ZWdzblZwR0k1UmY1aktPLCNGMkZjb0YyRjIqJiJNS094KVtPI2ZvUClIJD1lPF1QV3pMdCM9IkYyRltwRjJGMiomIktXRCpvKVtHW2d6Ij5fT1YvTC5IO113JUYyRl14RjJGMiomIktvekExN3ZxcFYjcDUpW15yOSQpUU1WSUYyRmF4RjJGPiomIkoncFRUblwrPUBcWSMzdWpuPicqeilSSEYyRmV4RjJGPiomIkkpb1FGLzUkbzokM2YsJXl2OCQqMylmPyZGMkZpeEYyRjIqJiJIU2EvdSFHbXU7UU0pNCRcIkdBWU1mJ0YyRjNGMkYyKiYiRzcwNyRwZG53bDJMV3cnejMlUmRAKEYyRl95RjJGPiomIkZLcy9tIlw0RCVwZSlmXmE7LUpFJSlGMkZjeUYyRj4qJiJFc3FtVHEnPjxAXG9mYzAlZldJISpGMkZbb0YyRjIqJiJEK0dMb1YxdjddI0gvKkhMPS1uJkYyRml5RjJGMiomIkNfWnJbJz53Oj9qM19tdEApeSgpRjJGXXpGMkY+KiYiQHNdXiVlaFI9eUx1YjsmPXMoRjJGYXpGMkY+KiYiQUdEZXVEW2w4Mmw8W1c+ZldGMkZlekYyRjIqJiJAKzdwb1YiPU9ReWZyV1RSQUYyRjhGMkY+KiYiPDcsNCNmWE1tYSVbIioqb0hGMkZbW2xGMkYyKiYiPSNmUXk8QD9gc1J0JylSSSRGMkZnb0YyRjIqJiI8KVteekxHSlpvKSk0NmI2RjJGYVtsRjJGPiomIjo3VCJ6XnYiKj04bVMhKT1GMkZlW2xGMkYyKiYiOG9CSyZ6bDNQaEFTPkYyRmlbbEYyRj4qJiI1Myk+ZFUtWT07SChGMkY0RjJGMiIyIypIeCV6eV0+ZUY+RjIpRl9cbEZmekYyRjIqJixecComIlNnTCVwanVLdDtedlpQZTpvaiZ6IXA1ZTBsOSpGMkZJRjJGPiomIlQnXHUvKEhFO19SPz1sZiJRVjIkbzF2YnclPk5VRjJGRUYyRj4qJiJUK19CUFlRNGd0cVNJZTstOCRRXWspW0BmRzEpRjJGQUYyRj4qJiJUaydwJG9SQiV5MWY9aCNRR3ldYCtuPFpvISo0NClGMkY8RjJGPiomIlRvZipIVWchMyRlVGZKY3I+SihcYFV3aFNaZjIlRjJGYXRGMkY+KiYiU28mXGpdIzQtRzduWD10PXQ1OUdPPig9WGElUUYyRmV0RjJGMiomIlQvZmIlPiNmNl1TbFpfP0BNYSR6J1tbTiJbRHBCRjJGaXRGMkYyKiYiVClbdix0KVFrVHdKNSFbM0pNUjlzS1owSUdfOUYyRl11RjJGMiomIlNXJDR3T3AoUkV5X1IleXFQYmsjZTVtTHcsd05GMkZhdUYyRj4qJiJTJypbXSRRVndidV4qUigpMyVIR0w9bTxzSGpSQipGMkZldUYyRj4qJiJTJ1xBNExGOjhKUyE0RzN0VUtKKj1lcT05JSlmJEYyRml1RjJGPiomIlMzVTAoKik0aT1GNSR6IikqXFclZSZRYSQ+WUBhWTVGMkZddkYyRjIqJiJTJSl6UEkzN28xKD4+Z1MsWDVcI0hiI29kITRSNkYyRmF2RjJGMiomIlJfTjEkPXl3OXYtV01ncC9YZCUqKjQ6KSpwKHo1RjJGZXZGMkYyKiYiUmM1cVFDIypmQjcmPlhtL3I6TFtbclUiKlt2OkYyRml2RjJGPiomIlE/dkl6YCdRZyl6JVsjZT5bOnIieidSJ2VfenolRjJGXXdGMkY+KiYiUCUpZllKMTVWY1VkUyEpZWcxI1EzMyRSK1JwKEYyRmF3RjJGMiomIlA7Z2dGd0VMKnpoZU5lWitnUllUa11waltGMkZnbkYyRjIqJiJPXz55eSVSNmRTWU0iPUMjZnlLVzpTa002I0YyRl9vRjJGPiomIk8/WlQ3dClIJ1xHak8qPiFRUGtycFU4MiNwJEYyRmNvRjJGPiomIk43LCpbQz5sMl10TmtwSS1bSU0qMztVTCJGMkZbcEYyRj4qJiJOV0RUbSN6WHltSEJ6JnokZVsoNGtzPydcP0YyRl14RjJGMiomIk1PRlxsQFRmJSl5P1lRW2ZBSSUqPU1RJmUjRjJGYXhGMkYyKiYiTClbRnBKInkwQSRlb08jcHpEWDpDJnl2WkYyRmV4RjJGPiomIktHVFVAY2slNGFhTlRKInpjIiopW1NcbilGMkZpeEYyRj4qJiJKJypvdG1gdS5jYGBiRms5SENkcElGKUYyRjNGMkYyKiYiSjtrIipcXSQ+az41MG8rK1k/dTphPjtGMkZfeUYyRjIqJiJJJ0g6IzMsREM7MWVkLXd2cEsmKT5JOEYyRmN5RjJGPiomIkhrM1lOXDE7IVshKlxWMioqeXdKUmo8RjJGW29GMkY+KiYiR29mWiIqb05PW2pwWDFOLCVmbiIqKT1GMkZpeUYyRjIqJiJFL2ZWWGdxXV1jLWpfVDx5LSQ9KEYyRl16RjJGMiomIkVPJGU4Oy4nXFlWU19bXTs7K0c8RjJGYXpGMkY+KiYiQ092cCgqKkd3O2cjPVQqPldRLywnRjJGZXpGMkYyKiYiQjcmM21qa2dudUtfSHRNNSFcJUYyRjhGMkYyKiYiQVNpZXcoKXp4IylvL0dBW2BBZUYyRltbbEYyRj4qJiJAcydvdSgpPXN1Yk12b096JVIkRjJGZ29GMkYyKiYiPylHWzI5XDxQKWUhXFlyOEIiRjJGYVtsRjJGPiomIj1HXFNQTkZqbHFQYl0oSCNGMkZlW2xGMkYyKiYiO2srI3BYOWcnNDlKPEU3RjJGaVtsRjJGPiomIjdnUEBTcl4/LDkiKiopRjJGNEYyRj4iM0svVCh5OjAlZSkqRjJGMilGX1xsRl56RjJGMiomLGZvKiYiVC8+dHddVUAlcCUpPjtxNWhrKil6JG8xMWN1VSI+RjJGYXRGMkY+KiYiVGdIUilRSz5MRCd5KmVNJipRbEA9RksmUikzWGR0RjJGZXRGMkY+KiYiVVcwKWZsV1k0LFtXWiFwJmYtZERYRzc1J0drZzVGMkZpdEYyRj4qJiJUbzIjPmliJGVhMCRwdCcqW01wZHolZWM7YGlpKyhGMkZddUYyRj4qJiJUJzRfK1x6LTF4MSoqUTlAKm9hdm9NdEhNLld3IkYyRmF1RjJGPiomIlMlUXd4bnVRSG0hUmxyI3lGKD4jM01yWURuKjRvRjJGZXVGMkYyKiYiVFdoUzFndHMhPT9tKEgjbz1AdSkqKm8sNm1Dbi0iRjJGaXVGMkYyKiYiUykpZWtndXAxJmYmKUhPQjQiZmhfZiQ9MVEiKT47I0YyRl12RjJGMiomIlM/djkqZkVKKnlnNWczOy5EKj5JLDYmMyZRaDEnRjJGYXZGMkY+KiYiU0N1Q210biRwbWBTN0w9OjZKVCdSIm9xZmsmW0YyRmV2RjJGPiomIlEhKUclR0VGOE1lQks6PTVJbnJOYnFtTz84JkYyRml2RjJGPiomIlNPQkElKkcyRkZUKCpwZm1TTEVcdCFIWSIpXGpAIkYyRl13RjJGMiomIlJzRVxVeG9xWyVvQ05LJHp6Yjd2Kik+NW0jUUdGMkZhd0YyRjIqJiJSK18+LVhURiZIWiNIMVRkKW9da1xjPEJtVzxGMkZnbkYyRj4qJiJRK1d4VzteVFYhM1E2QVgjKnlYMjAiNHo8XiYpRjJGX29GMkY+KiYiUHNpY0cmKSlcbnEjPl85SU81NU55cG4rVFkjRjJGY29GMkYyKiYiUF82R2JVT0tZT3BrJ3B4bWNrJkd2YlI0TSlGMkZbcEYyRjIqJiJQRyJbV3g9Jnl0VldhYV8uSXVRIWVcTjY7N0YyRl14RjJGMiomIk9LKUdeNio+Syx3cWwvXztmXEU4KFx5aGsjRjJGYXhGMkY+KiYiTlsvNyVcQD84InptRTUqPTxZdzhWJTRYYChGMkZleEYyRj4qJiJMaz96P2A9LUMoejo5aVpoKkgnXGJEKT4qRjJGaXhGMkYyKiYiTUM5Rip6Zkx4XigzOi1gYlsiPUBsNTNFIkYyRjNGMkYyKiYiS3dEUi0xMWl3a1JAczJpLFVKd1NeZEYyRl95RjJGPiomIksvTiRIJT4vIVxhKGZnJD1tPiZRZU4lW3IiRjJGY3lGMkY+KiYiSiV5VSo0eF5GMSQzIkckW2htdzknemomPkYyRltvRjJGMiomIkgzIVIwJGY9ISo9ZHoxb1MmSC5wQDUkKkYyRml5RjJGMiomIkhjRVUpZVQmM0NPSi50aFgtUUgrQiRGMkZdekYyRj4qJiJHaydcNywvTUIvSTUiNGslKmYoKXBsPEYyRmF6RjJGMiomIkZfbnkmeiRIKDR6IWU9PEFNeio0IUcjRjJGZXpGMkYyKiYiRTMteiZbJCl5WV5oK0NQV2tVIlFGRjJGOEYyRj4qJiJDZ1xFO0hcKCk9JCpIZiR5KGZGcE1GMkZbW2xGMkYyKiYiQmdAXF1UJD5mYj5uXjhXUFpTRjJGZ29GMkYyKiYiQCMqNHFSbVwoKlsjenk6JDNTaihGMkZhW2xGMkY+KiYiP1MnKTMoR2slM0Mnbz4meThEQEYyRmVbbEYyRj4qJiI9S0ciKlxzNkQmPSdldVF5T0YyRmlbbEYyRjIqJiI7J0hObGpeZ1hmbSlbTURGMkY0RjJGPiI4K3dgcyFIKj04b1c3I0YyRjIpRl9cbEZcb0YyRjIqJixebyomIlRfWiZIWyo9JHpFIypvbWdvV3MsaGwwXyZmdDE+I0YyRmF1RjJGPiomIlRXbFIjcDJubGMsTVcnels5YElYPWAnRzVPSW9GMkZldUYyRj4qJiJUb2ZSJls4I1FDKltUU1xIRDRdKzA4MzcjW0hPKEYyRml1RjJGPiomIlQhM3UkKiozRFFjNEBXcD5mJipcVyJbVC4vIVw3PiRGMkZddkYyRj4qJiJTIVtPKy5WKTNXJClbQ210eTN3QG9jLXIkKXBzIylGMkZhdkYyRj4qJiJTITNpV0tYTk4hUSVRSlAqcCdcLiFSWWAsImZoRSdGMkZldkYyRj4qJiJTKVtaJW9adDMqKik9Xm8iM0InUm1fdm02JzQpeiJIRjJGaXZGMkY+KiYiUiEpKVEwJVxeRlh4QGckeiNmLjRAJ1tpJT0xTGYmRjJGXXdGMkY+KiYiUmtvSmhodGclKTRFaSQzdXdIV2p0J28qUlsicCdGMkZhd0YyRj4qJiJSK0dsJ1xKVlkuUCcpejI2a2koUUBZKmZHaF9NRjJGZ25GMkY+KiYiUCUpZkFUUDpAejY/KlJ3NydmIT1iPVRxXjgqKUYyRl9vRjJGPiomIlEpRzNQPG1WWjUhZVhpWlAxMDBZJjRcSkBYJEYyRmNvRjJGPiomIlFfalY1dHk1JUhtOjw1QCczemAnUl8lKXk8WiRGMkZbcEYyRj4qJiJQKy8ueHMmKVEkelJCZXknPSFHOSM0Tm9iSmAlRjJGXXhGMkYyKiYiUE8iUkhpVW0lXDgiZTYrKkdfTVJrcDh4JDQkKUYyRmF4RjJGMiomIk8pWzYoPjF0PiJmTHQlPkU5QV48Wl08akUoKkYyRmV4RjJGMiomIk8hM3FBNCZHQiZ6aUQkZll2Iylwb0RONnczaEYyRml4RjJGPiomIk9zLS5aO1dVKSpvRElfLW5LTzcyUjthOTlGMkYzRjJGPiomIk5XRDxrbXY7aytuJyoqM1hUYjRCMkwmXGEiRjJGX3lGMkYyKiYiTTtPJXBcTnMkUVhyTT9wWjtKXnUkSFxwJ0YyRmN5RjJGMiomIkxPaktweTJ5aXkiSG0qenIzTSlwRF5PPkYyRltvRjJGMiomIktjbTcwTTwvYSd5eT5fLk04Yix2Jj56RjJGaXlGMkY+KiYiSk86QylcSyk+JlIxUW5RITRRcV10LEAjRjJGXXpGMkY+KiYiSTtnTyYpPiMpPm09QHVRQSQqMzRdJHAiUkYyRmF6RjJGMiomIkh3eCd5OHokMyo0Qz53dSdlMzgkPikpPkYyRmV6RjJGPiomIkc7K15Wank/b1M/dChbTiFIISk+bkVGMkY4RjJGMiomIkYzUSEqKVIqKnozcDVkS05jInB6U2siRjJGW1tsRjJGMiomIkVjdVVAXiFmUUtwYHB4XG03S2clRjJGZ29GMkY+KiYiRCM+U2speidHKG8lNF1bJlJWYyFlQ0YyRmFbbEYyRjIqJiJCR3AqZnBSUGxzNU1ZME0lWyhcRjJGZVtsRjJGPiomIkBnJlFqO1dYNiVRL3YtXHVyJUYyRmlbbEYyRjIqJiI9JSk+S2xzOkkicFNiIVwpcCVGMkY0RjJGPiI7Z2hiTVp5KjRkdSpleTVGMkYyKUZfXGxGYHlGMkYyKiYsZm4qJiJUUyNRJD47bi4jcDhOJWZVJFJfVnNoW2dEI28jWzVGMkZhdkYyRj4qJiJUczEnW2lYZ2NgPCgpM1JPMyh5X2E9I0dAb2ZHRUYyRmV2RjJGPiomIlRrJypbZC5EbT4qeiRwJTRUL0wxUSdRKT1yIz05VCNGMkZpdkYyRj4qJiJUd2hdInBpaHVTWyF6R2NERUswXkFZcCVmMU44RjJGXXdGMkY+KiYiU08qPkAjeTJWVz9Dc254MSc9QU5COmtEWVNiKUYyRmF3RjJGPiomIlNzbSZIVG8rVCZIJnpAPEElZnVsLCZIYlVqRmAlRjJGZ25GMkY+KiYiUkM9LT43OEVXQUNzMSgzIjRlOSdbPi4rLihlKUYyRl9vRjJGPiomIlIhb0ZSPiZIeXghcC1XeClwZCtdRjkqRysoKjNgRjJGY29GMkYyKiYiUk9eQnEpKVwiPVdYeSNvXVc0J1w9PC05OS93XkYyRltwRjJGMiomIlJfcnRNU3JsUSEpeS0mKT1aX2dCYScqPnZGUjEiRjJGXXhGMkYyKiYiUWNpQCYzbWd5NDowQHJ4T105Yk1WIkg/USopRjJGYXhGMkY+KiYiUXNlYGYhKkhxZCI+JSpcbiJbPWZVW2tOKVtaJVtGMkZleEYyRj4qJiJPY0U1d083XjhwXSpwZzZGSGNqbDAvOVUjRjJGaXhGMkY+KiYiUC9ib2hMISlbNWtPJVx1JylHSl5YJ2ZVND9mJUYyRjNGMkYyKiYiT0NZek8qeml0JHlxenFlUEkmcCN6RzB3ayoqRjJGX3lGMkYyKiYiTktnblw0YSEpKUdeUHFxRXVgQ14iR1s/WCZGMkZjeUYyRjIqJiJNKUdjcClSLyYzQ0QqUSh5MmolcFpacUwqXCdGMkZbb0YyRjIqJiJNYy1LJW80KW9NdClROm1aTmQyMSJRInA4IkYyRml5RjJGMiomIkxPUikqZSMzaCpHSXMzencjRyZRYl9UR1MoRjJGXXpGMkY+KiYiTHdkLjR3dHQ2cF54KXAhPSY0LSY9L1l4IkYyRmF6RjJGPiomIkpbKSl6YUU6I0d1MU0mKVtbO2ZHazRZTUYyRmV6RjJGPiomIkorUyFSRylRWCcqUVUySVolKmVBQSozTTpGMkY4RjJGMiomIkhzJVxBJT5TKUc3X3JZJSpcVTNIYypRJUYyRltbbEYyRjIqJiJHYycpenguMzloc2plXDhiaSMpXFFjRjJGZ29GMkY+KiYiRkNFQ1A4PilmdypbXHVdai8mSFA3RjJGYVtsRjJGMiomIkQ7IT1gQGgnZlA4K0ZeQHNlb2coRjJGZVtsRjJGPiomIkMzJzN6QW1VJykqKTRxW2ojXHBQQEYyRmlbbEYyRjIqJiJAT2YsQiN5b1IqKkc0ZStaVzdGMkY0RjJGMiI9UzFnKClRLlNqcjpdakRDRjJGMilGX1xsRmp4RjJGMiomLFQqJiJSYyc9dCpIVlM2aVMyZV9SOzlORFddbWFWOiVGMkZhd0YyRjIqJiJSVzBJM3dCcEZyQTJFPiR5Kj10VVh3YyFILjlGMkZnbkYyRj4qJiJTJWVacHVdLHgvbWc1I29JViFwWyVmIyo9L1hZQkYyRl9vRjJGPiomIlMzKSpSUDJeK142OWRnYlosTT4xcVBjMXc7RUYyRmNvRjJGPiomIlJ3cDsmcHcrbHMoZWdFW3AycS1hME4/NCE0c0YyRltwRjJGPiomIlIlZUosWyYqeWxaJnlXMHojKnAoUlZDQnAsdkA6RjJGXXhGMkYyKiYiUUN5JGUjeilRTWQ1czskW050dyEzYVh6MFksI0YyRmF4RjJGPiomIlJPRkgnZSgqPkhPJFsmb1RDI3A6SXoraHBqakAiRjJGZXhGMkY+KiYiUSM+Q2dTSy07MjpyL2N0JSllZ0pUOkJiZmQmRjJGaXhGMkY+KiYiUWdMOSlvNSEpelNxXE1BIUdRJm9ES2A5JykpZTxGMkYzRjJGMiomIlFTMTdyUFgmelJReVRmMSZIPyY0RURiJGUzSEYyRl95RjJGMiomIlEvRnJzJTQlbykpbyxLa05FMT8kKSpRWyVIV2E1RjJGY3lGMkYyKiYiT3cieURdOTpILEwtJXpdK09bLyxDLC1jYkYyRltvRjJGMiomIk8hW19bLSc+Jyl6KFI4TXQ9YiEqKkdRU1c5OGFGMkZpeUYyRj4qJiJPO2t6Nz16dFd2J3lEc15ZbixcYGE/Y1QiRjJGXXpGMkY+KiYiTjNBUWMiZUgjPU1cZUh1W1tYKltWMUBYN0YyRmF6RjJGPiomIkxXKmUxN20jWztBInBmKClcWEAiKSpwKD1KTkYyRmV6RjJGMiomIkw/SkQsLEsnXDZhRHJBTmpgYXdqKDQ1I0YyRjhGMkYyKiYiSyN6LClvNjtXKSpbZipwZ1MqNGNjRGRoQ0YyRltbbEYyRjIqJiJKV0BxPCpIJkhqIjMyTXJLQE9ddGtnN0YyRmdvRjJGMiomIkcvMlc9JXAkKlIrKkhXbHkybSs1ISopKkYyRmFbbEYyRjIqJiJHXz4pPm9XI29bXDI/cyZwWkpHM0YjRjJGZVtsRjJGPiomIkUpRzM4bk48aF50eHNdITQvS0FJRjJGaVtsRjJGMiomIkJfenpjQFc2JUdzc0JfPiNSIlJGMkY0RjJGMiI/VyI9KD4zOVFIa3ZTISk9QkVGMkYyKUZfXGxGYnhGMkYyKiYsTComIlJ3JjNXOU96PzMyW2dLcDJFKT5WVEUkUSoqektGMkZbcEYyRjIqJiJTa2MxRWQmPlNuR0FVcGsiXEJeaFQnPUkyQkMiRjJGXXhGMkYyKiYiUy8+WG4vb1ZjZnFeJzNoRjNhJj5CaFQ3IVEnPUYyRmF4RjJGMiomIlNTKXpcSmQ1aT0mXEgyUmNxTVBmaldfMSY0VyJGMkZleEYyRjIqJiJSZyRmWGMwZVdWIypvcD9Scl8jUig0OGxNT2duRjJGaXhGMkYyKiYiUmc8RVlZJClHY1wpPiQ0PilbJTMiKltGREwncHZERjJGM0YyRjIqJiJRYz0hKWYqeldSdzxYZV1ecDt5eFhKRlVyZypGMkZfeUYyRjIqJiJRSysjZSJmN2xbV0hveGA+L0h4clBAYWpjSEYyRmN5RjJGMiomIlBnSC4oZl9OUXJKaUkzPXpwd3BPZW1DJipmRjJGW29GMkYyKiYiUFNhczRZJio0L0AneVUlXEwkUUxALC9xcDdDRjJGaXlGMkY+KiYiUFs7TG0peSpRJEdPJip6NlslSFciKmVCXC47ZCNGMkZdekYyRj4qJiJPaylHLS93LTIzbEonZTUjRyRcLnghUTVONChGMkZhekYyRj4qJiJOJVErNyVwbzVlRyRReSg+Y2opW0lzKCkqcGIkRjJGZXpGMkY+KiYiTiNmdWQ5JUcmeV5CITNTXz8+XSxPTlAmZjxGMkY4RjJGMiomIk0zVThmYnVWcih6SkQ1IUdoajc7L0AmKVJGMkZbW2xGMkYyKiYiTDdMYEYpKWU5a2VUNGB3Kil5Mld1TldtJEYyRmdvRjJGMiomIktndC8vbVJjKSlIb3MnenYiW2M/UDsiPiJGMkZhW2xGMkYyKiYiSTNFbCpIKjNqLTMjeU0pZkVoNnI2KXAiRjJGZVtsRjJGPiomIkcnXC0jMzopcEtpM0J2YlJadD14VSJGMkZpW2xGMkYyKiYiRGtbTzoiUWVUWmsiKXlgVUd5cEJGMkY0RjJGMiJBZyIzKnk4JylSRVAiXDxUZko/IkYyRjIpRl9cbEZccEYyRjIqJixEKiYiUWMtX1Y2aF5KWShwXjo8WSQqKSpSIj4xRXYnWyRGMkZpeEYyRj4qJiJSa0tPbztGIXlFXzEtcjp1OjUnPi5VdEtVRiJGMkYzRjJGPiomIlI7L00sQyRcQittKnojKnk8L3FOXzdUKnokNCI9RjJGX3lGMkY+KiYiUiNmQSMpXFgmUmw5IUhZNSl6WSUqKj5QJSo9JVJESyJGMkZjeUYyRj4qJiJRJT1YUzQpelgpUmNzMnFQV04qKnBbISpIRFtLJ0YyRltvRjJGPiomIlFTIWYuJkglSClwQEtJTm9nPWdYYlFiISl6XUZGMkZpeUYyRj4qJiJQP25PPXYkKikpUS0pSDVcUF95dGtBI3pHJzQoKkYyRl16RjJGPiomIlArXzYvRDY2Si1cSmY6a1smPkkjUnZbJjNBRjJGYXpGMkY+KiYiT1NlZkE6WShIJFw5RXhgTkckKmUjXC0jeV9bRjJGZXpGMkY+KiYiTndgL3I+KFtzQD96bzZBIz53Kzx1LCRwXUYyRjhGMkYyKiYiTm9mQjwnXF1xS0EoM1FFcVAyWjQnSFxNdUYyRltbbEYyRjIqJiJOPzp0bXV5U2o+Z0Z3ak5ALHckcE5XI1EiRjJGZ29GMkYyKiYiTC9KKWYscip6dyIpcCJINCVlbzRqJkhuInAqRjJGYVtsRjJGMiomIkk/aiRwVSVvcyIpZjwjSGhSJm9lTjNNIilGMkZlW2xGMkYyKiYiSCkzLDEpR2B4OV0rSGdhdj02VyZHTkYyRmlbbEYyRjIqJiJFX1paIz1raEowalt6WFxeUltsJkYyRjRGMkYyIkJvciNmcXMuelZzJj5gQGRiVCNGMkYyKUZfXGxGYG9GMkYyKiYsPComIk83dD4kcHJnJ3kyKlxlKlFLRXNsKm9JNCI+KSpGMkZbb0YyRjIqJiJQKSllP04oUm9BKj5OXidSKFIwQGZzMSwoKWZqJEYyRml5RjJGMiomIlA3RGpMSkcxNHAkKVxsOjEtZ1NnNygpXCo0XkYyRl16RjJGMiomIlAjPkMpUUpZJVxmeHcwczBAUDZDUytYJT1jJEYyRmF6RjJGMiomIlBDLXE2ISlRPEREbHBEO114MGhxV3pfUGwiRjJGZXpGMkYyKiYiTzMhPiMpSFBaLTwiKSopKkdTT1EjPj9cS2kjUXlGMkY4RjJGMiomIk9jIykzJCpRT2VlNyhcKW8qR2IwIjRIL0wpM2sjRjJGW1tsRjJGMiomIk5jQUY9cyFwKio+TlgjKipIZSEpXFIkUig9KnBYRjJGZ29GMkYyKiYiTnd0Tkd4IUddT1lPRDEvcHcuQSI+aUF5IkYyRmFbbEYyRjIqJiJMJXk9KlskUSk+JSo+RikzNUczVCo0U10ubyJGMkZlW2xGMkYyKiYiSVNRV2h4OE5zJypwc1M2STNRQ2o9XkYyRmlbbEYyRjIqJiJGc0FJOHVOYFFzKSopcF9EUytwQmBGMkY0RjJGMiJDW3NaXkx4Rjg+PVYrX2tecTxGMkYyKUZfXGxGYndGMkYyKiYiP2tnYDdWcSoqKVsqUSFbIj5NO0YyRl14RjJGMiomIj1DTXkmNC08WUMsKmYydSkpRjJGYXhGMkY+KiYiPWNdP1skZlVvLyozaHFrOEYyRmV4RjJGPiomIjt3QCc0I3BDLShSaWIpPiYpRjJGaXhGMkYyKiYiM29YRCRSbmNoRSZGMkZhekYyRjIqJiJMV1xCVEF1TElyI2UmPWMlUSYpKilSKnpaZCdGMkZhdEYyRj4qJiJMQzUpM0BuYndfTmBCQ3hPNClbQkZePyJGMkZldEYyRjIqJiJLU1E/WjlyWihSWmRDNGBKXyhSVi8wcEYyRml0RjJGMiomIkspKXluJ1x5ISk9Q1MjKkdNb0hEU1ZgcisiRjJGXXVGMkYyKiYiSTtjJil6Ly5sNGpwVDInXFQ+VF1YSipGMkZhdUYyRj4qJiJJY3VVUD4mKVsjKVIuKW8mUigqekZIXzJhRjJGZXVGMkY+KiYiSF9WKW9DI1xhWE94NEFYWGFjNHAlXEYyRml1RjJGPiomIkc7Y0R5LXVTImYhKVwlWy86ZiRcJioqKSpGMkZddkYyRjIqJiJHK186OidIbkIramEkWzVSNS5tQkNGMkZhdkYyRjIqJiJFP3YpSC8vVFVVWVJWZyoqXEBpaCJGMkZldkYyRj4qJiJFJ1whSHQ4W09MNXgtdidwRHBzXCZGMkZpdkYyRj4qJiI4czU4emkiSEthZEo9RjJGY3lGMkY+KiYiOT9SSiJ6NnZxJj1nU29GMkZfeUYyRj4qJiIyPFNkcTwpKSkqZiJGMkZlekYyRj4qJiJDR3RpeUVgJWYlNDNOSiU0NHdEKUYyRmF3RjJGMiomIkRvWkVOLkRnZThDVHZ2WkhlRiRGMkZdd0YyRj4qJiJQd2xUbkExOUJvOiZIQSMqb3BNXXo3WWxUWkYyRl1zRjJGPiomIlBbSzw4JUhqQ1wyKSopKjNqbWZNWmB0ZU1fY0YyRmlyRjJGPiomIk9nMEQkSHBMd1pzcUlNYSRcJDNmVDYjNCxGRjJGZXJGMkYyKiYiUSdIOl1xLFQyazRyLlpYaUduaFkpNE8vdjhGMkZhckYyRjIqJiJRP3poNDExcXgyKnkrX0BBL01VbCM0ZSlmbCNGMkZdckYyRjIqJiJRO3dDQlAhKW9ARCxPV3klKmU5K0A1Jyp5LFVERjJGaXFGMkYyKiYiUE9yZTVfW292JikqZU5aLCs7c0g4Pj8zWlhGMkZlcUYyRjIqJiJRW187VT89R0RaIik9VGRFMyJIKXlLcSR6P2cjRjJGYXFGMkY+KiYiUW9iXzdPXmREdWRQJEhGPzBaIUd6Si4walpGMkZdcUYyRj4qJiJRJ0gkKXBaRyg9MGghXFdqKltMelVSTGhLN0pcRjJGaXBGMkY+KiYiUWNVKSoqPThJcCgqPjh5JXoiSHIlZjZTIT1zJ1tORjJGZXBGMkY+KiYiUUd0UVc7RG04b3ZDQzIhb088TUMyXlMoXD1GMkZhcEYyRj4qJiJQRyopZWU4cFUlKnAscj92P05GeTR5Mzk3I3BGMilGNCIjZEYyRj4qJiJQY21XYD9QMk1ibi8mWyIqb2ItbCtsRSdSeCJGMilGNCIjZUYyRj4qJiJPI3pgSkxoTjg8Z2xZL2tyJXpTVD4iUSUpeiNGMilGNCIjZkYyRj4qJiJOVyo9R2BRQjdmXlIiW3lPSG8tUyg9IjQwI0YyKUY0IiNnRjJGPiomIi5lLURyZVMjRjJGW1tsRjJGPiomIjU7b21TJW8jM0hJd0YyRml5RjJGPiomIjRPNk1JV3U9IWVxRjJGXXpGMkY+KiYiKXB2OENGMkZhW2xGMkY+SSRvZ2ZHRig3I0Yu
computeParams(algeq);
NyYiI0MiI2ciJUQ6IiR4JQ==
algB2 := primpart(subs(F(t)=B2,op(1,algeq)),B2):
algB2 := collect(algB2,B2,factor):
Save all computed equations save(algR0,algR1,algB1,algB2,guessedeqs):
LUklbXJvd0c2Iy9JK21vZHVsZW5hbWVHNiJJLFR5cGVzZXR0aW5nR0koX3N5c2xpYkdGJzYjLUkjbWlHRiQ2I1EhRic=