use Array: all;
use StdIO : all;
use Numerical : all;
use CommandLine: all;
use String: {to_string,tochar,sscanf};
use ArrayFormat: all;
use Bits: all;



double[.,.] floydXDD_slow(double[.,.] D)
{
  for( k=0; k<shape(D)[0]; k++) {
    for( i=0; i<shape(D)[0]; i++) {
      for( j=0; j<shape(D)[1]; j++) {
        D[i,j] = min(D[i,k]+D[k,j], D[i,j]);
      }
    }
  }
  return( D);
}

double[.,.] floydXDD_fast(double[.,.] D)
{
  siz0 = shape(D)[0];
  siz1 = shape(D)[1];
  
  for( k=0; k<siz0; k++) {
    for( i=0; i<siz0; i++) {
      for( j=0; j<siz1; j++) {
        D[i,j] = min(D[i,k]+D[k,j], D[i,j]);
      }
    }
  }
  return( D);
}

int main()
{
  D = tod(reshape([10,10], reverse(iota(100))));
  print(D);
  
  S = floydXDD_slow(D);
  F = floydXDD_fast(D);
  
  print(S);
  print(F);
  return 0;
}