結果

問題 No.3072 Sum of sqrt(x)
ユーザー AtamaokaCAtamaokaC
提出日時 2021-10-25 22:50:52
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 1,146 bytes
コンパイル時間 8,955 ms
コンパイル使用メモリ 364,316 KB
実行使用メモリ 10,624 KB
最終ジャッジ日時 2024-10-03 21:27:56
合計ジャッジ時間 24,370 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
10,624 KB
testcase_01 AC 3 ms
5,248 KB
testcase_02 AC 3 ms
5,248 KB
testcase_03 AC 3 ms
5,248 KB
testcase_04 AC 3 ms
5,248 KB
testcase_05 AC 3 ms
5,248 KB
testcase_06 AC 3 ms
5,248 KB
testcase_07 TLE -
testcase_08 -- -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
testcase_19 -- -
testcase_20 -- -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <iomanip>
#include <boost/multiprecision/cpp_int.hpp>
using namespace std;
using namespace boost::multiprecision;
typedef cpp_int Bint;
typedef long long int llint;
typedef long double real_t;

Bint sqrt_floor(Bint x)
{
    Bint ceil = x, floor = 1;
    while (ceil > floor + 1) {
        auto m = (ceil+floor)/2;
        if (m*m <= x)
            floor = m;
        else
            ceil = m;
    }
    return floor;
}

int main()
{
    Bint unit = pow(Bint(10),20);
    cout << fixed << setprecision(17);

    int N;
    cin >> N;
    Bint ans = 0;
    for (auto i=0; i<N; ++i) {
        Bint x;
        cin >> x;

        auto p = sqrt_floor(x);
        if (p*p == x) {
            ans += p*unit;
        } else {
            Bint R = 1;
            Bint U = 1, V = p;
            Bint a = 2*p, b = x-p*p;
            while (R < unit) {
                R *= 2*p;
                auto prev_V = V;
                V = a*V + b*U;
                U = prev_V;
            }
            ans -= p * unit;
            ans += (V * unit)/U;
        }
        cout << real_t(ans)/real_t(unit) << endl;
    }

    return 0;
}
0