DATA nombre_magique;
INPUT a;
CARDS;
1
;
RUN;


DATA nombre_magique2;
SET nombre_magique;
DO a=1 TO 9; 
DO b=1 TO 9;
DO c=1 TO 9;
DO d=1 TO 9;
e=5;
DO f=1 TO 9;
DO g=1 TO 9;
DO h=1 TO 9;
i=45-a-b-c-d-e-f-g-h;                         /*Si chaque chiffre n'apparait qu'une fois alors leur somme vaut 45 */
if a*b*c*d*e*f*g*h*i= 362880           /*Et leur produit vaut 362880 */                                         
and mod(a*10+b,2)=0                      /* Puis on verifie les conditions demandées*/
and mod(a*10**2+b*10+c,3)=0 
and mod(a*10**3+b*10**2+c*10+d,4)=0 
and mod(a*10**5+b*10**4+c*10**3+d*10**2+e*10+f,6)=0
and mod(a*10**6+b*10**5+c*10**4+d*10**3+e*10**2+f*10+g,7)=0
and mod(a*10**7+b*10**6+c*10**5+d*10**4+e*10**3+f*10**2+g*10+h,8)=0 
then output;
END;
END;
END;
END;
END;
END;
END;
RUN;

proc print; run;

*Je crée une table avec tout les nombres de 1000 à 9999 -> 9000 obs
    éviter un proc print sur la totalité ^^;

data nombre4;
do mil=1 to 9;
    do cent=0 to 9;
        do diz=0 to 9;
            do unit=0 to 9;
            output;
            end;
        end;
    end;
end;
run;
proc contents; run;

*Je crée une première table dans laquelle seuls les nombres fraternels de la forme 1112, 3234, 9899
(nombres consécutifs mais avec les 2 premiers plus petits que les deux derniers);

data nombre4_2;
set nombre4;
                                                                       *pour obtenir les nombres dans la configuration voulue, j'élimine étape par étape tout les autres;
if cent=unit then delete;                                   *le chiffre des centaines est forcément différent de celui des unités etc...;
if diz>(mil+1) then delete; 
if unit>(cent+1) and unit^=0 then delete;
if (mil+cent)>(diz+unit) and unit^=0 then delete;
if unit=0 and diz ne (mil+1) then delete;
if unit^=0 and mil ne diz then delete;
if unit=0 and cent^=9 then delete;                    *là par exemple c'est pour gérer les cas 1920, 2930, 3940...;
run;
proc print data=nombre4_2; run;                         *à la fin j'ai une table qui contient 89 obs de nombres fraternels à quatre chiffres
                                           de la forme 1112, 1213..
                                            il suffira ensuite de permuter pour obtenir les autres (1211 1312...;

*Je concatène les variables (elles seront transformées en variables caractères);

data nombre4_3 (drop= mil cent diz unit);
set nombre4_2;
a1=strip(mil)||strip(cent)||strip(diz)!!strip(unit);
run;

*Je regroupe le chiffre des milliers et des centaines puis celui des dizaines et des unités
    pour ensuite pouvoir permuter;

data nombre4_4 (drop= mil cent diz unit b1 b2);
set nombre4_2;
b1=strip(mil)!!strip(cent);
b2=strip(diz)!!strip(unit);
a1=strip(b2)!!strip(b1);                             *(1112 -> 1211, 6768 -> 6867 etc...);
run;

*J'utilise la procédure APPEND pour ajouter les observations de la table que je viens d'obtenir.
    une étape data avec instruction set sur les deux tables aurait pu tout aussi bien faire l'affaire;

proc append base=nombre4_3 data=nombre4_4;
run;
proc print data=nombre4_3; run;                   *j'ai maintenant tout les nombres à quatre chiffres fraternels;

data nombre_fraternel (drop=a1);
set nombre4_3;
a2=input(a1,8.);                                              *j'utilise INPUT pour transformer la variable caractère en variable numérique;
a3=sqrt(a2);
if round(a3)^= a3 then delete;                         *l'arrondi d'un nombre entier est égal à lui même,
                                          pour conserver uniquement le/les carrés parfaits, je demande à supprimer toutes les obs qui
                                          n’auraient pas un nombre entier pour racine carrée;
run;
proc print; run;

data nf;
do a=10 to 99999;
l=int(log10(a));
if l^=log10(a) then do;    /*on enlève les cas où a=10, 100, 1000 ou 10000*/
b=a-1;
N1=a*10**(l+1)+b;        /* si a=112, alors b=111 et l+1=3 : N1=112*10^3+111=112111*/
X1=N1**(1/2);            /* ainsi les N1 et N2 que l'on crée sont fraternelle*/
N2=b*10**(l+1)+a;
X2=N2**(1/2);
if int(X1)=X1 then do;   /*on vérifie que N1 est un carré parfait*/
N=N1;X=X1;
output;
end;
if int(X2)=X2 then do;
N=N2;X=X2;
output;
end;
end;
end;
Keep N X;
run;

proc print; run;

data toto;
   do i=10 to 98;
      X=i*100+i+1;
      output;
   end;
   do i=100 to 998;
      X=i*1000+i+1;
      output;
   end;
   do i=1000 TO 9998;
      X=i*10000+i+1;
      OUTPUT;
   END;
   Do i=10000 to 99998;
      X=i*100000+i+1;
      OUTPUT;
   END;
run;

data toto;set toto;
   where x**0.5=int(x**.5);
   Y=x**.5;
run;

proc print;run;

Obs      i               X      Y

 1       183        183184      428
 2       328        328329      573
 3       528        528529      727
 4       715        715716      846
 5      6099      60996100     7810
 6     13224    1322413225    36365
 7     40495    4049540496    63636