#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>

#define MAXN 60
#include "nauty22/nauty.h"

void canon(int vertices, char *s, char *canon) {
  int i,j,k,new,edges;
  char c;

  graph g[MAXN*MAXM]; 
  graph gcanon[MAXN*MAXM]; 
  nvector lab[MAXN], ptn[MAXN], orbits[MAXN];
  static DEFAULTOPTIONS(options);
  statsblk(stats);
  setword workspace[50*MAXM];
  set *gv;

  options.writemarkers=FALSE;
  options.defaultptn=FALSE;
  options.getcanon=TRUE;
  options.writeautoms=FALSE;

  for(i=0; i<MAXN; i++) {
    gv=GRAPHROW(g,i,MAXM);
    EMPTYSET(gv, MAXM);
  }
    
  edges=0;
  for(i=0; i<vertices; i++) {
    for(j=0; j<i; j++) {
      for(k=0; k<(*s-'0'); k++) {
        /*    create new edge-vertex #vertices+edges and add edges {i,new} {j,new} */
        new = vertices+edges;
        gv = GRAPHROW(g,i,MAXM);
        ADDELEMENT(gv, new);
        gv = GRAPHROW(g,j,MAXM);
        ADDELEMENT(gv, new);
        gv = GRAPHROW(g,new,MAXM);
        ADDELEMENT(gv, i);
        ADDELEMENT(gv, j);
 
        edges++;
      }
      s++;
    }
  }

  for(i=0; i<vertices+edges; i++) {
    lab[i]=i;
    ptn[i]=1;
  }

  ptn[vertices-1]=0;
  ptn[vertices+edges-1]=0;

  nauty(g, lab, ptn, NILSET, orbits, &options, &stats, workspace, 50*MAXM, MAXM, vertices+edges, gcanon);

  for(i=0; i<vertices; i++) {
    for(j=0; j<i; j++) {
      c='0';
      for(k=0; k<vertices+edges;k++) {
        if(i==k || j==k)
          continue;
        gv = GRAPHROW(gcanon,k,MAXM);
        if (ISELEMENT(gv, i) && ISELEMENT(gv, j)) {
          c++;
        }
      }
      *(canon++)=c;
    }
  }
  *canon=0;
}

int main(int argc, char**argv)
{
  char buffer[1024];
  canon(4,"201010",buffer);
  printf("%s\n",buffer);
  return 0;
}