GenQF.cc

#include <iostream>
#include <stdlib.h>

#define QDIST(XA,YA,XB,YB) ((XB-XA)*(XB-XA)+(YB-YA)*(YB-YA))

using namespace std;

int main(int argc, char **argv) {

  if(argc < 4 || argc > 5) {
    cout<<"Syntax: GenQF #variables #scenarios density [seed]"<<endl;
    return 1;
  }

  int n = atoi(argv[1]);
  int s = atoi(argv[2]);
  double d = atof(argv[3]);
  if(argc > 4) srandom(atoi(argv[4]));
  else srandom(0);

  double *posx = new double[s*n];
  double *posy = new double[s*n];

  for(int i = 0; i < s*n; i++) {
    posx[i] = (double)random()/RAND_MAX;
    posy[i] = (double)random()/RAND_MAX;
  }

  double *m = new double[s];

  for(int j = 0; j < s; j++) m[j] = 0;
  cout<<n<<" "<<(n*(n-1))/2<<" "<<s<<endl;
  for(int i1 = 0; i1 < n-1; i1++) {
    for(int i2 = i1+1; i2 < n; i2++) {
      cout<<i1+1<<" "<<i2+1<<" ";
      for(int j = 0; j < s; j++) {
        double qd = QDIST(posx[j*n+i1],posy[j*n+i1],posx[j*n+i2],posy[j*n+i2]);
        cout<<" "<<qd;
        if((double)random()/RAND_MAX < d && (double)random()/RAND_MAX < d) {
          m[j] += qd;
        }
     }
      cout<<endl;
    }
  }

  delete[] posx;
  delete[] posy;


  for(int j = 0; j < s; j++) {
    cout<<" "<<m[j];
  }
  cout<<endl;

  delete[] m;

  return 0;

}

---