結果

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

ソースコード

diff #
raw source code 

// clang-format off
#pragma GCC optimize ("O3")
#pragma GCC target ("avx")

#include <cassert>
#include <iostream>
#include <vector>
#include <algorithm>
#include <utility>
#include <cmath>
#include <iomanip>

#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;
constexpr long double eps = 0.000000001;

inline long double power(long double x, int a)
{
    if(a >= 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;
}
0