Subversion Repositories wimsdev

Rev

Rev 7692 | Rev 8560 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 7692 Rev 8321
Line 18... Line 18...
18
$   to ; in $slib_G
18
$   to ; in $slib_G
19
 
19
 
20
slib_n=!exec pari slib_N=Mat([$slib_G]);matsize(slib_N)[1]
20
slib_n=!exec pari slib_N=Mat([$slib_G]);matsize(slib_N)[1]
21
 
21
 
22
slib_dstx=!exec pari {slib_N=Mat([$slib_G]);\
22
slib_dstx=!exec pari {slib_N=Mat([$slib_G]);\
23
	slib_P=vector($slib_n,slib_i,-1);\
23
    slib_P=vector($slib_n,slib_i,-1);\
24
	slib_P[$slib_point]=0;\
24
    slib_P[$slib_point]=0;\
25
	slib_modif=1;\
25
    slib_modif=1;\
26
	while(slib_modif,\
26
    while(slib_modif,\
27
  		slib_modif=0;\
27
          slib_modif=0;\
28
 		 for(slib_i=1,$slib_n,\
28
          for(slib_i=1,$slib_n,\
29
    		for(slib_j=1,$slib_n,\
29
            for(slib_j=1,$slib_n,\
30
     			 if(slib_N[slib_i,slib_j]*(slib_P[slib_i]+1),\
30
                  if(slib_N[slib_i,slib_j]*(slib_P[slib_i]+1),\
31
       				 if((slib_P[slib_j]+1),\
31
                        if((slib_P[slib_j]+1),\
32
         				 if(slib_P[slib_j]>slib_P[slib_i]+1, slib_P[slib_j]=slib_P[slib_i]+1;slib_modif=1),\
32
                          if(slib_P[slib_j]>slib_P[slib_i]+1, slib_P[slib_j]=slib_P[slib_i]+1;slib_modif=1),\
33
          				slib_P[slib_j]=slib_P[slib_i]+1; slib_modif=1);\
33
                          slib_P[slib_j]=slib_P[slib_i]+1; slib_modif=1);\
34
          		);\
34
                  );\
35
   			);\
35
               );\
36
 		 );\
36
          );\
37
	);\
37
    );\
38
print(slib_P)}
38
print(slib_P)}
39
 
39
 
40
slib_cc=!exec pari {slib_P=Vec([$slib_dstx]);\
40
slib_cc=!exec pari {slib_P=Vec([$slib_dstx]);\
41
slib_Q=vector($slib_n);\
41
slib_Q=vector($slib_n);\
42
for(slib_i=1,$slib_n,\
42
for(slib_i=1,$slib_n,\