UVA - 11809 Floating-Point Numbers

#include <bits/stdc++.h>
#define DET 1e-4
using namespace std;

double A[16][64];
long long int B[16][64];

void init()
{
for(int i = 0; i <= 9; i++)
{
for(int j = 1; j <= 30; j++)
{
double E = 1 - pow(2, -(i+1));
long long int M = pow(2, j) - 1;
double X = log10(E) + M*log10(2);
long long int b = floor(X);
double a = pow(10.0, X-b);
A[i][j] = a;
B[i][j] = b;
}
}
}

int main()
{
init();
char in[32];
while(cin >> in)
{
if(strcmp(in, "0e0")==0)
break;
double a;
long long int b;
*(strchr(in, 'e')) = ' ';
sscanf(in, "%lf %lld", &a, &b);
for(int i = 0; i <= 9; i++)
{
for(int j = 1; j <= 30; j++)
{
if(fabs(A[i][j]-a)<DET && B[i][j]==b)
printf("%d %d\n", i, j);
}
}
}
return 0;
}