#+begin_src javascript
  console.log("Me ejecuto en node!");
  #+end_src

(org-babel-do-load-languages
 'org-babel-load-languages '((C . t)))

#include <iostream>
#include <fstream>
#include <climits>
#include <cstdlib>
#include <string>

using namespace std;


ofstream log("saving-the-universe-again.log");

long value( string p ){
  long ret = 0;
  long strength = 1;
  for( string::iterator i = p.begin() ; i < p.end(); i++ ){
    if( *i == 'C' ){
      strength *= 2;
    }
    else{
      ret += strength;
    }
  }
  return ret;
}

void swap( string &p, int pos ){
  //log << "1swap " << p << endl;
  char c = p[pos];
  p[pos] = p[pos-1];
  p[pos-1] = c;
  //log << "2swap " << p << endl;
}

bool nextBestChange( string &p ){
  
  for( int i = p.length()-1 ; i > 0 ; i-- ){
    log << i << endl;
    if( p[i] == 'S' && p[i-1] == 'C' ){
      swap(p,i);
      return true;
    }
  }
  return false;
}

long saveTheUniverse( string p, long d ){
  long changes = 0;
  log << "saveTheUniverse:" << d << ":" << value(p) << endl;
  while(  value(p) > d ){
    if( !nextBestChange(p) ){
      log << "IMPOSIBLE" << endl;
      return -1;
    }
    changes++;
  }
  log << "solucion " << p << " " << d << ":" << changes;
  return changes;
}

string output(long l){
  if( l == -1)
    return "IMPOSSIBLE";
  else{
    char b[100];
    sprintf( b, "%ld", l );
    return b;
  }
}

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

  int T;
  cin >> T;

  for( int i = 0 ; i < T ; i++ ){
    long D;
    string P;
    cin >> D >> P;
    log << D << " " << P << endl;
    cout << "Case #" << (i+1) << ": " << output(saveTheUniverse(P,D)) << endl;
  }
}

#include <iostream>
#include <fstream>
#include <climits>
#include <cstdlib>
#include <string>
#include <cmath>

using namespace std;

int main(){
  printf( "%lf\n", atan( 1 / (2*sqrt(2)) ) );
  printf( "%lf\n", acos( (2*sqrt(2))/3 ) );
}

#include <iostream>
#include <fstream>
#include <climits>
#include <cstdlib>
#include <string>
#include <cmath>

using namespace std;


ofstream logger("cubic-ufo.log");


class Mat{
  public:
  double v[3][3];
  Mat(double _v[][3]){
    for( int f = 0 ; f < 3 ; f++ ){
      for( int c = 0 ; c < 3 ; c++ ){
        v[f][c] = _v[f][c];
      }
    }
  }
};

Mat rotX(double a){
  double ret[][3] = {
    { 1, 0, 0},
    { 0, cos(a), -sin(a)},
    { 0, sin(a), cos(a) }
  };
  return Mat(ret);
}

Mat rotZ(double a){
  double ret[][3] = {
    {  cos(a), -sin(a), 0},
    { sin(a), cos(a),0 },
    {0,0,1}
  };
  return Mat(ret);
}

class Point{
  public:
  double x,y,z;
  Point(double _x, double _y, double _z):x(_x),y(_y),z(_z){}
  
  Point rotate(double aroundX, double aroundZ){
    // https://es.mathworks.com/help/phased/ref/rotx.html?requestedDomain=true
    Mat rx = rotX(aroundX);
    Mat rz = rotZ(aroundZ);
    return times(rx).times(rz);
  }

  Point times(Mat m){
    //logger << "\n\n TIMES *******" << endl;

    //logger << "this " << toString() << endl;

    //logger << m.v[0][0] << "\t" << m.v[0][1] << "\t" << m.v[0][2] << endl;
    //logger << m.v[1][0] << "\t" << m.v[1][1] << "\t" << m.v[1][2] << endl;
    //logger << m.v[2][0] << "\t" << m.v[2][1] << "\t" << m.v[2][2] << endl;

    double retx = x*m.v[0][0] + y*m.v[0][1] + z*m.v[0][2];
    double rety = x*m.v[1][0] + y*m.v[1][1] + z*m.v[1][2];
    double retz = x*m.v[2][0] + y*m.v[2][1] + z*m.v[2][2];

    //logger << "ret " << retx << " " <<rety << " " << retz;

    return Point(retx,rety,retz);
  }

  string toString(){
     char b[1000];
     snprintf(b, sizeof(b)/sizeof(*b), "%.20lf %.20lf %.20lf", x, y,z );
     //////logger << "toString" << " " << x << " " <<y << " " <<z << " " <<b << endl;
     return b;
  }
};

class Cube{
  public:
  Point a, b, c;

  Cube():a(0.5,0,0), b(0,0.5,0), c(0,0,0.5){}

  Cube rotate(double aroundX, double aroundZ ){
    Cube ret;
    ret.a = a.rotate(aroundX,aroundZ);
    ret.b = b.rotate(aroundX,aroundZ);
    ret.c = c.rotate(aroundX,aroundZ);
    //logger << "rotate " << toString() << endl;
    return ret;
  }

  string toString(){
     char buf[1000];
     snprintf(buf, sizeof(buf)/sizeof(*buf), "%s \n%s \n%s", a.toString().c_str(), b.toString().c_str(), c.toString().c_str() );
     return string(buf);
  }

};


const double maxAroundZ = acos( sqrt(2)/sqrt(3) );
const double SQRT2 = sqrt(2);
Cube onlyAroundX( double a ){
  double aroundX = M_PI/4 + acos( a/SQRT2 );
  return Cube().rotate(aroundX,0);
}

double maxAreaForAroundZ(double aroundz){
  return SQRT2*cos(aroundz) + sin(aroundz);
}

double findAroundZ(double a){
  double minz = 0;
  double maxz = maxAroundZ;

  double ret = (maxz + minz)/2;
  double area = maxAreaForAroundZ(ret);
  while( fabs( area - a) > 0.00000001 ){
    
    if( area > a )
      maxz = ret;
    else
      minz = ret;

    ret = (maxz + minz)/2;
    area = maxAreaForAroundZ(ret);
  }

  return ret;
}




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

  int T;
  cin >> T;

  for( int i = 0 ; i < T ; i++ ){
    double A;
    cin >> A;
    if( A <= SQRT2){
      Cube c = onlyAroundX(A);
      printf( "Case #%d:\n%s\n", i+1, c.toString().c_str() );
    }
    else{
      Cube c = onlyAroundX(SQRT2).rotate(0,findAroundZ(A));
      printf( "Case #%d:\n%s\n", i+1, c.toString().c_str() );
    }
  }
}