結果

問題 No.1273 はじめのζ関数
コンテスト
ユーザー longrun
提出日時 2020-10-30 22:25:30
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,731 bytes
コンパイル時間 2,456 ms
コンパイル使用メモリ 193,672 KB
最終ジャッジ日時 2025-01-15 17:32:06
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 39 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
// #include <atcoder/all>
using namespace std;
// using namespace atcoder;
#define rep2(i,m,n) for (int i = (m); i < (n); ++i)
#define rep(i,n) rep2(i,0,n)
#define debug(x) cout << x << endl;
typedef long long int ll;
typedef long double ld;
typedef pair<int, int> P;
template<typename T> struct V : vector<T> { using vector<T>::vector; };
V() -> V<int>;
V(size_t) -> V<int>;
template<typename T> V(size_t, T) -> V<T>;
template<typename T> vector<T> make_vec(size_t n, T a) { return vector<T>(n, a); }
template<typename... Ts> auto make_vec(size_t n, Ts... ts) { return vector<decltype(make_vec(ts...))>(n, make_vec(ts...)); }
template<typename T> ostream &operator<<(ostream &os, const vector<T> &v) { for (auto &e : v) os << e << ' '; return os; }
struct fast_ios { fast_ios(){ cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;
template<typename T> inline bool chmin(T &a, T b) { if (a > b) { a = b; return true; } return false; }
template<typename T> inline bool chmax(T &a, T b) { if (a < b) { a = b; return true; } return false; }
const int INF = 1<<30;
const ll LINF = 1LL<<61;
// const ll MOD = 1000000007;
// using mint = modint1000000007;
ld eps = 1e-9;

ld power(ld x, int k)
{
    ld res = 1;
    while(k)
    {
        if(k & 1) res = (res * x);
        k >>= 1;
        x = (x * x);
    }
    return res;
}

ld sum(ll n, ld r)
{
    return (ld)r / (ld)(1 - r) - (r - power(r, n+1)) / (1 - r);
}

int main()
{
    ll n;
    cin >> n;
    ld ans = 0;
    for(ll i=2; i<1010; i++)
    {
        ld tmp =  sum(n-1, (ld)1 / (ld)i);
        ans += tmp;
    }
    ll res = (ans * 1e7) + eps;
    res /= 10;
    cout << res << endl;
    
    return 0;
}
0