#######################################################################################################################################
# #
# V I S U A L I S A T I O N D E L A C O N J E C T U R E " abc " : #
# #
# #
# Author of '$xiirv/.CABC.11.1.$U' : #
# #
# Jean-Francois COLONNA (LACTAMME, 20120919185333). #
# #
#######################################################################################################################################
$Z setParam _____NNombres 50
$Z @ NNombres2N = $_____NNombres * $_____NNombres
$Z setParam _____Zoom 0.95
$Z setParam _____RotationOX 0.0
$Z setParam _____RotationOY 0.0
$Z setParam _____RotationOZ 0.0
$Z setParam _____FRayon 0.035
$Z setParam _____Palette $xiP/trou_noir.74
$Z setParam _____Lx `GetParam $xrv/particule.10$X Lx`
$Z # Parametre introduit le 20120929092957... #
$Z setParam _____Ly `GetParam $xrv/particule.10$X Ly`
$Z # Parametre introduit le 20120929092957... #
$Z setParam _____Lz 1000
$Z # Parametre introduit le 20120929092957... #
$Z setParam _____ParticulePara1 "$K_VIDE"
$Z set PrEmIeR_NoMbRe=1
$Z FilSTmpB FNoMbReS_N
$Z FilSTmpB FNoMbReS_AN
$Z FilSTmpB FNoMbReS_BN
$Z FilSTmpB FNoMbReS_CN
$Z FilSTmpB FNoMbReS_PGCD_AB
$Z FilSTmpB FNoMbReS_PGCD_CA
$Z FilSTmpB FNoMbReS_PGCD_CB
$Z FilSTmpB FNoMbReS_A
$Z FilSTmpB FNoMbReS_B
$Z FilSTmpB FNoMbReS_C
$Z FilSTmpB FNoMbReS_ABC
$Z FilSTmpB FNoMbReS_RADICAL
$Z FilSTmpB FNoMbReS_LC
$Z FilSTmpB FNoMbReS_LRADICAL
$Z FilSTmpB FNoMbReS_Rapport
$Z FilSTmpB FNoMbReS_RAYON
$Z FilSTmpB FNoMbReS_COULEURS
$Z $xrv/accumule.01$X \
$Z ne=$_____NNombres \
$Z fichier=$PrEmIeR_NoMbRe \
$Z >>& $FNoMbReS_N
$Z $xrv/ADD2Tenso.01$X \
$Z ne=$_____NNombres \
$Z fichier1=$FNoMbReS_N \
$Z fichier2=0 \
$Z >>& $FNoMbReS_AN
$Z $xrv/ADD2Tenso.01$X \
$Z ne=$_____NNombres \
$Z fichier1=0 \
$Z fichier2=$FNoMbReS_N \
$Z >>& $FNoMbReS_BN
$Z $xrv/ADD2.11$X \
$Z ne=$NNombres2N \
$Z fichier1=$FNoMbReS_AN \
$Z fichier2=$FNoMbReS_BN \
$Z >>& $FNoMbReS_CN
$Z # Introduit le 20120922085851 uniquement pour verification... #
$Z $xrv/PGCD.11$X \
$Z ne=$NNombres2N \
$Z fichier1=$FNoMbReS_AN \
$Z fichier2=$FNoMbReS_BN | \
$Z $xrv/IFEQ.11$X \
$Z ne=0 \
$Z fichier1== \
$Z fichier2=1 \
$Z >>& $FNoMbReS_PGCD_AB
$Z # Determination des nombres {A,B} premiers entre-eux afin de ne conserver qu'eux... #
$Z $xrv/PGCD.11$X \
$Z ne=$NNombres2N \
$Z fichier1=$FNoMbReS_CN \
$Z fichier2=$FNoMbReS_AN | \
$Z $xrv/IFEQ.11$X \
$Z ne=0 \
$Z fichier1== \
$Z fichier2=1 \
$Z >>& $FNoMbReS_PGCD_CA
$Z $xrv/PGCD.11$X \
$Z ne=$NNombres2N \
$Z fichier1=$FNoMbReS_CN \
$Z fichier2=$FNoMbReS_BN | \
$Z $xrv/IFEQ.11$X \
$Z ne=0 \
$Z fichier1== \
$Z fichier2=1 \
$Z >>& $FNoMbReS_PGCD_CB
$Z $DIF -q $FNoMbReS_PGCD_CA $FNoMbReS_PGCD_CB \
$Z >& $nul
$Z if ($status != $OK) then
$Z EROR "Le calcul des PGCDs est incoherent."
$Z # Test introduit le 20120922080107 pour validation... #
$Z else
$Z endif
$Z $xrv/selection.01$X \
$Z ne=$NNombres2N \
$Z fichier=$FNoMbReS_AN \
$Z selection=$FNoMbReS_PGCD_AB \
$Z >>& $FNoMbReS_A
$Z $xrv/selection.01$X \
$Z ne=$NNombres2N \
$Z fichier=$FNoMbReS_BN \
$Z selection=$FNoMbReS_PGCD_AB \
$Z >>& $FNoMbReS_B
$Z set NNombres2A=`$CA $FNoMbReS_A | $WCl`
$Z set NNombres2B=`$CA $FNoMbReS_B | $WCl`
$Z set NNombres2=`$xcg/MIN2.01$X a=$NNombres2A b=$NNombres2B`
$Z $xrv/ADD2.11$X \
$Z ne=$NNombres2 \
$Z fichier1=$FNoMbReS_A \
$Z fichier2=$FNoMbReS_B \
$Z >>& $FNoMbReS_C
$Z # Calcul de : #
$Z # #
$Z # C = A+B #
$Z # #
$Z $xrv/MUL3.11$X \
$Z ne=$NNombres2 \
$Z fichier1=$FNoMbReS_A \
$Z fichier2=$FNoMbReS_B \
$Z fichier3=$FNoMbReS_C \
$Z >>& $FNoMbReS_ABC
$Z # Calcul de : #
$Z # #
$Z # A.B.C #
$Z # #
$Z $xrv/DivPrem.01$X \
$Z compatibilite_20220528=VRAI \
$Z ne=$NNombres2 \
$Z fichier=$FNoMbReS_ABC \
$Z editer_diviseurs_premiers=FAUX | \
$Z $SE -e 's/^.*=\([^ ]*\)$/\1/' \
$Z >>& $FNoMbReS_RADICAL
$Z # Calcul de : #
$Z # #
$Z # Radical(A.B.C) #
$Z # #
$Z # c'est-a-dire du plus grand diviseur "square free" du produit A.B.C. #
$Z # #
$Z # La "compatibilite_20220528" a ete introduite le 20220528151826 car, en effet, il semble #
$Z # qu'ici les nombres entiers auxquels commencent a 1 (et non pas a 2...). #
$Z $xrv/LOGX.01$X \
$Z ne=$NNombres2 \
$Z fichier=$FNoMbReS_C \
$Z >>& $FNoMbReS_LC
$Z $xrv/LOGX.01$X \
$Z ne=$NNombres2 \
$Z fichier=$FNoMbReS_RADICAL \
$Z >>& $FNoMbReS_LRADICAL
$Z $xrv/DIVZ.11$X \
$Z ne=$NNombres2 \
$Z fichier1=$FNoMbReS_LC \
$Z fichier2=$FNoMbReS_LRADICAL \
$Z >>& $FNoMbReS_Rapport
$Z # Calcul de : #
$Z # #
$Z # log(C) #
$Z # k(A,B,C) = --------------------- #
$Z # log(Radical(A.B.C)) #
$Z # #
$Z # La conjecture ABC dit que ce rapport est borne (par 1.6 environ...). #
$Z set MiNiMuM_Rapport=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_Rapport ms=VRAI`
$Z set MaXiMuM_Rapport=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_Rapport Ms=VRAI`
$Z echo "MIN2(k(A,B,C))=$MiNiMuM_Rapport"
$Z echo "MAX2(k(A,B,C))=$MaXiMuM_Rapport"
$Z set MiNiMuM_A=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_A ms=VRAI`
$Z set MaXiMuM_A=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_A Ms=VRAI`
$Z set AmPlItUdE_A=`calcul $MaXiMuM_A-$MiNiMuM_A`
$Z set HoMoThEtIe_A=`calcul 1/$AmPlItUdE_A`
$Z set TrAnSlAtIoN_A=`calcul -$MiNiMuM_A/$AmPlItUdE_A`
$Z set MiNiMuM_B=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_B ms=VRAI`
$Z set MaXiMuM_B=`$xrv/extrema.01$X ne=0 fichier=$FNoMbReS_B Ms=VRAI`
$Z set AmPlItUdE_B=`calcul $MaXiMuM_B-$MiNiMuM_B`
$Z set HoMoThEtIe_B=`calcul 1/$AmPlItUdE_B`
$Z set TrAnSlAtIoN_B=`calcul -$MiNiMuM_B/$AmPlItUdE_B`
$Z $xrv/RACX.01$X \
$Z ne=$NNombres2 \
$Z fichier=$FNoMbReS_Rapport \
$Z >>& $FNoMbReS_RAYON
$Z # La surface apparente des boules sera donc proportionnelle au rapport... #
$Z $xrv/normalise.01$X \
$Z ne=$NNombres2 \
$Z fichier=$FNoMbReS_Rapport \
$Z minimum=$NOIR \
$Z maximum=$BLANC \
$Z formater=VRAI entiers=VRAI \
$Z >>& $FNoMbReS_COULEURS
$Z $xrv/particule.10$X np=1 iterations=$NNombres2 \
$Z AX=$HoMoThEtIe_A BX=$TrAnSlAtIoN_A \
$Z LISTE_X=$FNoMbReS_A \
$Z AY=$HoMoThEtIe_B BY=$TrAnSlAtIoN_B \
$Z LISTE_Y=$FNoMbReS_B \
$Z LISTE_Z=0.5 \
$Z Ar=$_____FRayon \
$Z LISTE_RAYON=$FNoMbReS_RAYON \
$Z LISTE_NIVEAU=$FNoMbReS_COULEURS \
$Z paletteA=$_____Palette \
$Z Lx=$_____Lx \
$Z Ly=$_____Ly \
$Z Lz=$_____Lz \
$Z ZOOM=$_____Zoom \
$Z ROTATION_OX=$_____RotationOX \
$Z ROTATION_OY=$_____RotationOY \
$Z ROTATION_OZ=$_____RotationOZ \
$Z isoles=VRAI \
$Z RVB=VRAI \
$Z N_AU_CARRE=FAUX \
$Z editer_extrema_hors=FAUX \
$Z $_____ParticulePara1 \
$Z R=$_____ImageR \
$Z chiffres=0 \
$Z $formatI
$Z # Jusqu'au 20120920143410, il y avait : #
$Z # #
$Z # AZ=$HoMoThEtIe_C BZ=$TrAnSlAtIoN_C #
$Z # LISTE_Z=$FNoMbReS_C #
$Z # #
$Z # mais cela n'a aucun interet puisque le probleme est finalement bidimensionnel, la relation #
$Z # #
$Z # C = A+B #
$Z # #
$Z # definissant un plan... #
$Z FilSTmpE FNoMbReS_COULEURS
$Z FilSTmpE FNoMbReS_RAYON
$Z FilSTmpE FNoMbReS_Rapport
$Z FilSTmpE FNoMbReS_LRADICAL
$Z FilSTmpE FNoMbReS_LC
$Z FilSTmpE FNoMbReS_RADICAL
$Z FilSTmpE FNoMbReS_ABC
$Z FilSTmpE FNoMbReS_C
$Z FilSTmpE FNoMbReS_B
$Z FilSTmpE FNoMbReS_A
$Z FilSTmpE FNoMbReS_PGCD_CB
$Z FilSTmpE FNoMbReS_PGCD_CA
$Z FilSTmpE FNoMbReS_PGCD_AB
$Z FilSTmpE FNoMbReS_CN
$Z FilSTmpE FNoMbReS_BN
$Z FilSTmpE FNoMbReS_AN
$Z FilSTmpE FNoMbReS_N