結果
問題 |
No.2971 無理積分
|
ユーザー |
![]() |
提出日時 | 2024-11-29 22:51:32 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 86 ms / 1,000 ms |
コード長 | 757 bytes |
コンパイル時間 | 2,594 ms |
コンパイル使用メモリ | 243,616 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-29 22:51:39 |
合計ジャッジ時間 | 7,189 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 36 |
ソースコード
#define rep(i,n) for (int i=0;i < (int)(n);i++) #define ALL(v) v.begin(),v.end() typedef long long ll; #include <bits/stdc++.h> using namespace std; double n; double n3; //被積分関数 (acos(x)でx>1 や sqrt(x)でx<0 だと nan になるので端っこ注意) double f(double x) { return sqrt(x*x*x+n3); } int main() { ios::sync_with_stdio(false); std::cin.tie(nullptr); cin>>n; n3=n*n*n; double A=0; //下端 double B=1; //上端 int N=10000000; const double delta=(B-A)/N; double value=0.0; //i=1からでも大して影響ない for(int i=0;i<N;i++){ double x=A+i*delta; value+=(f(x)+f(x+delta)+4.0*f(x+delta/2.0))*delta/6.0; } cout<<fixed<<setprecision(15)<<value<<endl; return 0; }