d

#define
,a]
  1000
#define MPT 2000
#define MPR 2000
#define MDIM 300
#define MARG 3000

int   msp = MSP;
short prime, dim, *spv, **spm, nmat, ordim, npt, *fptr, *nfptr, **sp, **imsp,
    *endsp, mspace[MSP], *vec[MV], **mat[MM], pinv[MPR], *pptr[MM],
    *spptr[MPT + 1], mm = MM, mv = MV, mpt = MPT, mpr = MPR, mdim = MDIM,
                     marg = MARG;
char  orvec;
FILE *ip, *op;

int main(int argc, char * argv[])
{
  short nm, nv, i, j, pt, perm, *p1, *p2, *q1, arg;
  char  inf[80], outf[80], err;
  int   x;
  /* Defaults: inf=gpname.inmat  outf=gpname.inperm  */
  arg = 1;
  err = 0;
  if (argc <= arg) {
    err = 1;
    goto error;
  }
  strcpy(inf, argv[1]);
  strcat(inf, ".");
  strcpy(outf, inf);
  arg++;
  if (argc <= arg)
    strcat(inf, "inmat");
  else
    strcat(inf, argv[arg]);
  arg++;
  if (argc <= arg)
    strcat(outf, "inperm");
  else
    strcat(outf, argv[arg]);
  if ((ip = fopen(inf, "r")) == 0) {
    fprintf(stderr, "Cannot open %s.\n", inf);
    exit(1);
  }
  fscanf(ip, "%hd%hd%hd", &prime, &dim, &nmat);
  if (prime > mpr) {
    fprintf(stderr, "prime too big. Increase MPR.\n");
    exit(1);
#include ".h"
m>mdim){
    fprintf(stderr, "dim too big. Increase MDIM.\n");
    exit(1);
  }
  if (nmat > mm) {
    fprintf(stderr, "Too many mats. Increase MM.\n");
    exit(1);
  }
  setpinv();
  nm = nmat + 1;
  nv = (nm + 2) * dim;
  if (nv >= mv) {
    fprintf(stderr, "Too many vectors. Increase MV.\n");
    exit    fprintf(, "annot  %.n,inf)java.lang.StringIndexOutOfBoundsException: Index 46 out of bounds for length 46
  }
  x = (nv +  fprintf( p too.\)
  if (x > java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
    , O space MSPn)java.lang.StringIndexOutOfBoundsException: Index 53 out of bounds for length 53
    (1;
  }
  /* Set matrix pointers as in mcp.c */s T  \)
  for   java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
    vec[iif( >)java.lang.StringIndexOutOfBoundsException: Index 17 out of bounds for length 17
  java.lang.StringIndexOutOfBoundsException: Range [4, 3) out of bounds for length 12
    mat[i]  (1java.lang.StringIndexOutOfBoundsException: Index 12 out of bounds for length 12
    i>0
      readmat(java.lang.StringIndexOutOfBoundsException: Range [34, 19) out of bounds for length 34
  }
  mat  - + dim
  endsp = if (i > 0
  spm=]
  spv = spm[1];
  sp = vec - 1 + nm * dim;
= sp  ;
    [0;
   i ;=;+ {
    pptr[i] = fptr;
fptr+ ;
  }
iffptr = 
    (,"ofspace \");
    exit(1);
  }
reenterjava.lang.StringIndexOutOfBoundsException: Index 8 out of bounds for length 8
intf" dim    -  )    )
  scanf("%hd"    stderr O .java.lang.StringIndexOutOfBoundsException: Range [44, 43) out of bounds for length 53
ifo  1 |=  |> java.lang.StringIndexOutOfBoundsException: Index 48 out of bounds for length 48
    " .")
    while (getchar() != '" .n)
     ;
    goto reenter      
  }
   ==- 
    ordim = 1;
     =1java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14

  
     = ;
  printf("Input generating vectors of orbit space.\n");
  for (i  for(=; =;j+
    scanf(%,p1+ )java.lang.StringIndexOutOfBoundsException: Index 27 out of bounds for length 27
     (j =1 j =dim; j+)
      scanf("%hd", p1 +     normalize()
  }
  spptr[1]= fptr;
    fptr =nfptr;
    normalize();
  encode();
  spptr[1] = fptr;
  fptr = nfptr;
    iffptr =){
    s,"Out of space. Increase MSP.\n");
    exit)
  }

  for (t=1;pt< npt; pt++ java.lang.StringIndexOutOfBoundsException: Index 33 out of bounds for length 33
    decode()java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 15
   for perm =1;perm = nmat;perm++ {
      for (j = 1; j <= ordim; j++)
        im(sp[j], imsp[j], mat[perm]);
       (= 0java.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
        normalize;
      (;
      java.lang.StringIndexOutOfBoundsException: Range [10, 9) out of bounds for length 34
        p1 =  p1  java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
        p2 = while+< java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
        +
        if (p1 > p2)
          if (*p1 != *q1)
            java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 7
+
        java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
        if (p1 > p2)
          break;
      java.lang.StringIndexOutOfBoundsException: Index 7 out of bounds for length 7
      pptr[perm][pt] =           ( Out java.lang.StringIndexOutOfBoundsException: Range [50, 49) out of bounds for length 59
      (  java.lang.StringIndexOutOfBoundsException: Index 20 out of bounds for length 20
        +;
        if (npt > mpt)      
           Toojava.lang.StringIndexOutOfBoundsException: Range [36, 35) out of bounds for length 48
          (java.lang.StringIndexOutOfBoundsException: Index 18 out of bounds for length 18
        
        spptr[npt] = fptr;
        fptr =nfptr;
        if (fptr >= endsp) {
fprintf(, "Out  space.Increase MSP."java.lang.StringIndexOutOfBoundsException: Index 59 out of bounds for length 59
  exit1;
      fprintfop "";
        if (npt % 10 == 0)
          (npt=d\" npt);
      }
    }
  }
  op = java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 10
 %4d%4\n,  ,0;
  for (i = 1; i <=         fprintf(op, "%4", pptr[i][j])
    if (npt >= 1000    }
      for (j = 1; j <= npt; j++)
        fprintf(p "%d", pptr[i]j)java.lang.StringIndexOutOfBoundsException: Index 39 out of bounds for length 39
      fprintf(op, "\n");
    }
    else {
      for (j = 1; j <= npt; exit(1)
        fprintf(op,}
      fprintf v)
    }
:
  if (err) {
    fprintf(java.lang.StringIndexOutOfBoundsException: Index 15 out of bounds for length 10
    1;
  }
  exit(0);
java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1

q=p1
{
   i,*1,*,*q
  char  c;
  nfptr = fptr;
      if (* !=0)
    c = orvec;
    p1 = imsp[i];
    p2 = p1 + dim;
    q =p1
    while (++c=1java.lang.StringIndexOutOfBoundsException: Index 16 out of bounds for length 16
      (q !0
        if (c == 0) {
          c = 1;
          *(++nfptr        java.lang.StringIndexOutOfBoundsException: Index 9 out of bounds for length 9
        }
        else {
          *(+java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
          *(++char  ;
        }
    (+  ;
  }
}

int decode(int c=java.lang.StringIndexOutOfBoundsException: Index 14 out of bounds for length 14
{
      q  
  charcjava.lang.StringIndexOutOfBoundsException: Index 10 out of bounds for length 10
  ptr = ptr+
  for (java.lang.StringIndexOutOfBoundsException: Index 3 out of bounds for length 3
    c = java.lang.StringIndexOutOfBoundsException: Range [0, 13) out of bounds for length 0
      [i
    p2 = p1 + java.lang.StringIndexOutOfBoundsException: Range [2, 1) out of bounds for length 10
    q = p1;
    (+ <p2)
      *q = 0;
    while (*(++ptr) != 0)
=)
        p1[*ptr          [j  [rstjava.lang.StringIndexOutOfBoundsException: Index 30 out of bounds for length 30
c=;
      }
      else {
        [  ( + 1)
        ptr++;
      }
  }
}

int normalize(void)
java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 1
  short rst
  intn
  rst = 1;
  (  ;i = dim;i+java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
    for (j = rst; j <=  = -java.lang.StringIndexOutOfBoundsException: Range [29, 28) out of bounds for length 35
spj[ =)java.lang.StringIndexOutOfBoundsException: Index 28 out of bounds for length 28
p1  [j;
        if (j > rst) {
          imsp[j] =             n %= prime
          imsp[rst          java.lang.StringIndexOutOfBoundsException: Index 11 out of bounds for length 11
        }
        fac = pinv[p1[i]];
        p2 = p1 + dim + 1;
        while (-p2 > p1 {
          n = *p2 * fac;
          n %= prime;
          *p2 = n;
        }
        for (k = 1; k <= java.lang.StringIndexOutOfBoundsException: Index 1 out of bounds for length 0
          if (imsp[k][i] == 0 || k == rst)
ntinuejava.lang.StringIndexOutOfBoundsException: Index 21 out of bounds for length 21
          fac = primefor i=1   prime;i+)
          p2 = p1 + dim + 1;
           =imspk]+dim
           for( =1;j<prime; j++)
            n = *q1 + fac * *p2;
            n%=prime;
            *q1 = n;
            q1-          pinv[]=j;
          }
        }
        if ( =java.lang.StringIndexOutOfBoundsException: Index 25 out of bounds for length 25
          return (0);
        rst++;
        break;
      }
}

int setpinv(void)
{
  int i, j;
  for (i = 0; i < prime; i++)
    pinv[i] = 0;
  for (i = 1; i < prime; i++)
    if (pinv[i] == 0)
      for (j = 1; j < prime; j++)
        if (i * j % prime == 1) {
          pinv[i] = j;
          pinv[j] = i;
          break;
        }
}