// clang-format off #pragma GCC optimize ("O3") #pragma GCC target ("avx") #include #include #include #include #include #include #include #define FOR(i, a, b) for(int i = (a); i < (b); i++) #define RFOR(i, a, b) for(int i = (b)-1; i >= (a); i--) #define rep(i, n) for(int i = 0; i < (n); i++) #define rep1(i,n) for(int i = 1; i <= (n); i++) #define rrep(i, n) for(int i = (n)-1; i >= 0; i--) #define pb push_back #define mp make_pair #define fst first #define snd second #define show(x) cout << #x << " = " << x << endl #define chmin(x,y) x=min(x,y) #define chmax(x,y) x=max(x,y) #define pii pair #define vi vector using namespace std; template ostream& operator<<(ostream& o,const pair& p){return o<<"("< ostream& operator<<(ostream& o,const vector& vc){o<<"sz = "<= 0){ long double p = 1; rep(i, a){ p*=x; } return p; } else { const long double y = 1.0/x; long double p = 1; rep(i, -a){ p*=y; } return p; } } inline long double func(int a, int b, long double x, long double t) { return exp(a*x)*power(x,b)-t; } inline long double diff(int a, int b, long double x) { return (a*x+b)*exp(a*x)*power(x, b-1); } inline long double ddiff(int a, int b, long double x) { return ((a*x+b)*(a*x+b)-b)*exp(a*x)*power(x, b-2); } long double newton(int a, int b, long double t) { long double v = 10; long double prev = 0; while(true){ const long double y = func(a,b,v,t); prev = v; v -= y / diff(a,b,v); if(abs(v-prev) < eps){ // if(abs(exp(v)-exp(prev)) < eps){ break; } } return exp(v); } long double halley(int a, int b, long double t) { long double v = 10; long double prev = 0; while(true){ const long double y = func(a,b,v,t); prev = v; const long double d = diff(a,b,v); v -= y / (d-(y*ddiff(a,b,v)/(2*d))); if(abs(v-prev) < eps){ // if(abs(exp(v)-exp(prev)) < eps){ break; } } return exp(v); } int main() { cin >> m; rep(i, m){ int a,b; long double t; cin >> a >> b >> t; if(a == 0){ cout << fixed << setprecision(12) << exp(pow(t, 1.0/b)) << endl; } else if ( b == 0){ cout << fixed << setprecision(12) << pow(t, 1.0/a) << endl; } else if ( a == 1 or b == 1){ cout << fixed << setprecision(12) << newton(a,b,t) << endl; } else { cout << fixed << setprecision(12) << halley(a,b,t) << endl; } } return 0; }