/* gelee.bc Nombres gelés de Saint-Exupéry http://linenn.davalan.eu/data/pdf/gelee.pdf http://linenn.davalan.eu/data/prog/gelee.bc usage : bc -q gelee.bc */ define tri() { auto i, j, k, l, m, a0, a1, a2, a3, a4, p, u, a, b, c, x; a0=1; for(i=0; i<=1; i++) { a1=1; for(j=0; j<=4; j++) { a2=1; for(k=0; k<=2; k++) { a3=1; for(l=0; l<= 1; l++) { a4=1; for(m=0; m<=1; m++) { scale=0 x = a0*a1*a2*a3*a4; b = 311850/x; a = 2*x; c = sqrt(a*a+b*b); if(c%1 == 0 && a*a+b*b-c*c==0) { print a," ",b," ",c,"\n"; } a4 *= 11; } a3 *=7; } a2 *= 5; } a1 *= 3; } a0 *= 2; } } t=tri(); quit;