using namespace std;
#include<iostream>
#include<math.h>
#include<stdlib.h>
#include<fstream>
#include<vector>

int main(int argc, char* argv[]){
int N=100;
vector<double>Coord(3,0);
vector<vector< double > >Pos(N,Coord);
ofstream write;
write.open("initial_configuration.dat");
for(int n=0; n<N; n++){
Pos[n][0]=5*( sin(n*2.0*M_PI/double(N)) + 2.0*sin(2.0*n*2.0*M_PI/double(N)));
Pos[n][1]=5*( cos(n*2.0*M_PI/double(N)) - 2.0*cos(2.0*n*2.0*M_PI/double(N)));
Pos[n][2]=5*( -sin(3.0*n*2.0*M_PI/double(N)));

write << n+1 << " "<< 1 <<" "<< Pos[n][0] << " " << Pos[n][1] <<" " << Pos[n][2] <<endl; 
write << n+1 << " "<< 2  <<" "<< Pos[n][0]+1.12246 << " " << Pos[n][1] <<" " << Pos[n][2] <<endl; 
}


return 0;
}