結果
問題 | No.454 逆2乗和 |
ユーザー |
![]() |
提出日時 | 2016-12-05 02:02:39 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 2 ms / 2,000 ms |
コード長 | 1,641 bytes |
コンパイル時間 | 596 ms |
コンパイル使用メモリ | 88,664 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-30 02:40:57 |
合計ジャッジ時間 | 1,495 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 31 |
ソースコード
// #includes {{{#include <algorithm>#include <numeric>#include <iostream>#include <string>#include <vector>#include <queue>#include <list>#include <deque>#include <stack>#include <set>#include <map>#include <cstdio>#include <cstdlib>#include <cassert>#include <cstring>#include <cmath>using namespace std;// }}}// pre-written code {{{#define REP(i,n) for(int i=0;i<(int)(n);++i)#define RREP(i,a,b) for(int i=(int)(a);i<(int)(b);++i)#define FOR(i,c) for(__typeof((c).begin()) i=(c).begin();i!=(c).end();++i)#define LET(x,a) __typeof(a) x(a)//#define IFOR(i,it,c) for(__typeof((c).begin())it=(c).begin();it!=(c).end();++it,++i)#define ALL(c) (c).begin(), (c).end()#define MP make_pair#define EXIST(e,s) ((s).find(e)!=(s).end())#define RESET(a) memset((a),0,sizeof(a))#define SET(a) memset((a),-1,sizeof(a))#define PB push_back#define DEC(it,command) __typeof(command) it=command//debug#define dump(x) cerr << #x << " = " << (x) << endl;#define debug(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl;#define debug2(x) cerr << #x << " = [";REP(__ind,(x).size()){cerr << (x)[__ind] << ", ";}cerr << "] (L" << __LINE__ << ")" << endl;const int INF=0x3f3f3f3f;typedef long long Int;typedef unsigned long long uInt;typedef long double rn;// }}}int main(){rn x;cin>>x;rn ans = 1.0L/(x+1.0L)+1.0L/2.0L*(1.0L/(x+1.0L)-1.0L/(x+2.0L))+1.0L/((x+1.0L)*(x+1.0L)*(x+2.0L));rn d = x;for(int n=2;;n++){rn b = x+(rn)n;rn ans2 = ans-1.0L/b/b/(b+1.0L)/(b-1.0L);if(abs(ans-ans2)<1e-15L){break;}ans = ans2;}printf("%.10Lf\n",ans);return 0;}