結果

問題 No.456 Millions of Submits!
ユーザー rickythetarickytheta
提出日時 2016-12-07 23:06:08
言語 C++11
(gcc 11.4.0)
結果
TLE  
実行時間 -
コード長 2,959 bytes
コンパイル時間 2,517 ms
コンパイル使用メモリ 143,496 KB
実行使用メモリ 8,752 KB
最終ジャッジ日時 2023-09-05 07:16:52
合計ジャッジ時間 10,900 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 2 ms
4,376 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,380 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 3 ms
4,376 KB
testcase_07 AC 9 ms
4,376 KB
testcase_08 AC 8 ms
4,376 KB
testcase_09 AC 68 ms
4,376 KB
testcase_10 AC 68 ms
4,380 KB
testcase_11 AC 662 ms
4,380 KB
testcase_12 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

// floating
typedef double Real;
const Real EPS = 1e-11;
#define EQ0(x) (abs(x)<EPS)
#define EQ(a,b) (abs(a-b)<EPS)
typedef complex<Real> P;

// double solve1(double a,double b,double t){
//   if(a==0){
//     // ln(x)^b = t
//     t = pow(t,1.0/b);
//     return exp(t);
//   }else if(b==0){
//     // x^a = t;
//     return pow(t,1.0/a);
//   }else{
//     double n = 10.0;
//     REP(_,70){
//       double npam1 = pow(n,a-1);
//       double logn = log(n);
//       double lnpbm1 = pow(logn,b-1);
//       double fn = npam1*n*lnpbm1*logn - t;
//       double fdn = npam1*lnpbm1*(a*logn + b);
//       n -= fn/fdn;
//       DEBUG(fn);
//       DEBUG(fdn);
//       DEBUG(n);
//     }
//     return n;
//   }
// }

// double solve2(double a,double b,double t){
//   if(a==0){
//     // ln(x)^b = t
//     t = pow(t,1.0/b);
//     return exp(t);
//   }else if(b==0){
//     // x^a = t;
//     return pow(t,1.0/a);
//   }else{
//     t = log(t);
//     double lo=0.0, hi=log(25252.0);
//     REP(_,75){
//       double mid = (lo+hi)/2.0;
//       double tt = a*mid + b*log(mid);
//       if(tt<t)lo = mid;
//       else hi = mid;
//     }
//     return exp(lo);
//   }
// }

inline void solve(){
  double a,b,t;
  scanf("%lf%lf%lf",&a,&b,&t);
  if(a==0){
    // ln(x)^b = t
    t = pow(t,1.0/b);
    printf("%.12f\n",exp(t));
  }else if(b==0){
    // x^a = t;
    printf("%.12f\n",pow(t,1.0/a));
  }else{
    // t = log(t);
    // double lo=0.0, hi=log(25252.0);
    // REP(_,75){
    //   double mid = (lo+hi)/2.0;
    //   double tt = a*mid + b*log(mid);
    //   if(tt<t)lo = mid;
    //   else hi = mid;
    // }
    // printf("%.12f\n",exp(lo));
    double n = 10.0;
    REP(_,50){
      double npam1 = pow(n,a-1);
      double logn = log(n);
      double lnpbm1 = pow(logn,b-1);
      double fn = npam1*n*lnpbm1*logn - t;
      double fdn = npam1*lnpbm1*(a*logn + b);
      n -= fn/fdn;
      // DEBUG(n);
    }
    printf("%.12f\n",n);
  }
}

int main(){
  int m;
  scanf("%d",&m);
  while(m--)solve();
  // DEBUG(solve1(10.0,10.0,5.0));
  // DEBUG(solve2(10.0,10.0,5.0));
  // // for(double t = 0.1; t <= 10.0; t += 0.1){
  // //   double x = solve1(10.0,10.0,t);
  // //   double y = solve2(10.0,10.0,t);
  // //   DEBUG(abs(x-y));
  // // }
  return 0;
}
0