結果

問題 No.456 Millions of Submits!
コンテスト
ユーザー Pachicobue
提出日時 2017-03-17 00:37:08
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
WA  
実行時間 -
コード長 1,897 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 742 ms
コンパイル使用メモリ 65,676 KB
実行使用メモリ 14,016 KB
最終ジャッジ日時 2024-06-23 18:30:49
合計ジャッジ時間 10,671 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other WA * 12 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

// clang-format off
#include <cassert>
#include <iostream>
#include <vector>
#include <algorithm>
#include <utility>
#include <cmath>

#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<int, int>
#define vi vector<int>

using namespace std;
template<class S,class T> ostream& operator<<(ostream& o,const pair<S,T>& p){return o<<"("<<p.first<<","<<p.second<<")";}
template<class T> ostream& operator<<(ostream& o,const vector<T>& vc){o<<"sz = "<<vc.size()<<endl<<"[";for(const T& v:vc) o<<v<<",";o<<"]";return o;}
typedef long long ll;

int m;

inline long double diff(int a, int b, long double x)
{
    return pow(x, 1-a)*pow(log(x), 1-b)/(a*log(x)+b);
}

inline long double val(int a, int b, long double x, long double t)
{
    return pow(x, a)*pow(log(x), b)-t;
}

long double newton(int a, int b, long double t)
{
    long double v = 10;
    constexpr ll CNT = 100000;
    for(ll i = 0; i < CNT; i++){
        //    while(true){
        const long double y = val(a,b,v,t);
        // show(y);
        if(abs(y) < 0.00000000001){
            break;
        }
        v -= y * diff(a,b,v);
        // show(v);
    }
    return v;
}

int main()
{
    cin >> m;
    rep(i, m){
        int a,b;
        long double t;
        cin >> a >> b >> t;
        if(a == 0){
            cout << exp(pow(t, 1.0/b)) << endl;
        } else if ( b == 0){
            cout << pow(t, 1.0/a) << endl;
        } else {
            cout << newton(a,b,t) << endl;            
        }
    }
    return 0;
}
0