/*************************************************************************************************************************************/
/* */
/* E N C O D A G E D E S D E C I M A L E S D E ' P I ' G R A C E A U X N O M B R E S P R E M I E R S */
/* " A L A G O D E L " : */
/* */
/* */
/* Author of '$xtb/SuitesGodel.01.SubmissionSloane.pi' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20220320113103). */
/* */
/*************************************************************************************************************************************/
scale = 1000;
print "\n";
pn[0]=2;pn[1]=3;pn[2]=5;pn[3]=7;pn[4]=11;pn[5]=13;pn[6]=17;pn[7]=19;pn[8]=23;pn[9]=29;pn[10]=31;pn[11]=37;pn[12]=41;
pn[13]=-1;
/* Definition des premiers nombres premiers ("Prime Numbers"). */
pi[0]=1;pi[1]=4;pi[2]=1;pi[3]=5;pi[4]=9;pi[5]=2;pi[6]=6;pi[7]=5;
pi[8]=-1;
/* Definition des premieres decimales de 'PI'... */
define pow(x,n) {
produit = 1;
i = 1;
while (i <= n)
{
produit = produit*x;
i=i+1;
/* L'utilisation de 'i++' provoquerait l'edition de 'i', d'ou 'i=i+1'... */
}
return(produit);
}
iterer = 0;
i1 = 0;
i2 = 0;
while (iterer == 0)
{
if ((pn[i1] > 0) && (pi[i1] >= 0))
{
term = 1;
for (i2 = 0 ; i2 <= i1 ; i2 = i2+1)
{
term = term*pow(pn[i2],pi[i2]);
if (i2 == i1) print "a(",i2+1,")=",term,"\n";
}
} else
{
if (pn[i1] <= 0) print "OVERFLOW\n";
iterer = 1;
}
i1 = i1+1;
/* L'utilisation de 'i1++' provoquerait l'edition de 'i1', d'ou 'i1=i1+1'... */
}
quit;
Copyright © Jean-François COLONNA, 2022-2024.
Copyright © CMAP (Centre de Mathématiques APpliquées) UMR CNRS 7641 / École polytechnique, Institut Polytechnique de Paris, 2022-2024.