Subversion Repositories wimsdev

Rev

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

Rev 8560 Rev 12676
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
                  );\
-
 
35
               );\
34
              );\
36
          );\
35
          );\
-
 
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,\
43
  if(slib_P[slib_i]+1,slib_Q[slib_i]=1));\
43
    if(slib_P[slib_i]+1,slib_Q[slib_i]=1));\
44
print(slib_Q)}
44
    print(slib_Q)}
45
 
45
 
46
slib_out= !positionof item 1 in $slib_cc
46
slib_out= !positionof item 1 in $slib_cc
47
 
-