結果
問題 | No.881 sin(x)/xの和 |
ユーザー | chocorusk |
提出日時 | 2019-09-07 17:55:24 |
言語 | C++11 (gcc 11.4.0) |
結果 |
AC
|
実行時間 | 385 ms / 2,000 ms |
コード長 | 1,275 bytes |
コンパイル時間 | 666 ms |
コンパイル使用メモリ | 101,944 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-05-01 06:17:37 |
合計ジャッジ時間 | 5,438 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 318 ms
6,816 KB |
testcase_01 | AC | 162 ms
6,940 KB |
testcase_02 | AC | 113 ms
6,940 KB |
testcase_03 | AC | 125 ms
6,940 KB |
testcase_04 | AC | 368 ms
6,940 KB |
testcase_05 | AC | 57 ms
6,940 KB |
testcase_06 | AC | 35 ms
6,940 KB |
testcase_07 | AC | 65 ms
6,940 KB |
testcase_08 | AC | 302 ms
6,940 KB |
testcase_09 | AC | 177 ms
6,940 KB |
testcase_10 | AC | 132 ms
6,944 KB |
testcase_11 | AC | 222 ms
6,940 KB |
testcase_12 | AC | 275 ms
6,944 KB |
testcase_13 | AC | 135 ms
6,944 KB |
testcase_14 | AC | 81 ms
6,944 KB |
testcase_15 | AC | 79 ms
6,940 KB |
testcase_16 | AC | 185 ms
6,940 KB |
testcase_17 | AC | 73 ms
6,940 KB |
testcase_18 | AC | 111 ms
6,940 KB |
testcase_19 | AC | 166 ms
6,940 KB |
testcase_20 | AC | 31 ms
6,940 KB |
testcase_21 | AC | 179 ms
6,940 KB |
testcase_22 | AC | 385 ms
6,944 KB |
ソースコード
#include <cstdio> #include <cstring> #include <iostream> #include <string> #include <cmath> #include <bitset> #include <vector> #include <map> #include <set> #include <queue> #include <deque> #include <algorithm> #include <complex> #include <unordered_map> #include <unordered_set> #include <random> #include <cassert> #include <fstream> #include <utility> #include <functional> #define popcount __builtin_popcount using namespace std; typedef long long int ll; typedef pair<int, int> P; double sinc(double x){ if(abs(x)<1e-7) return 1.0; else return sin(x)/x; } int main() { int n; cin>>n; double ans=0; double a[1010], x[1010]; for(int i=0; i<n; i++){ cin>>x[i]>>a[i]; x[i]-=floor(x[i]); } for(int i=0; i<n; i++){ ans+=a[i]*(sinc(-1+x[i])+sinc(x[i])+sinc(1+x[i])); ans+=a[i]*cos(x[i]+1.5)/(2*sin(0.5))/(2+x[i]); ans+=a[i]*sin(x[i]+2)/(4*sin(0.5)*sin(0.5))/(2+x[i])/(3+x[i]); for(int j=2; j<=10000; j++){ ans-=a[i]*(1/(x[i]+j)/(x[i]+j+2))*sinc(x[i]+j+1)/(2*sin(0.5)*sin(0.5)); } ans+=a[i]*cos(-x[i]+1.5)/(2*sin(0.5))/(2-x[i]); ans+=a[i]*sin(-x[i]+2)/(4*sin(0.5)*sin(0.5))/(2-x[i])/(3-x[i]); for(int j=2; j<=10000; j++){ ans-=a[i]*(1/(-x[i]+j)/(-x[i]+j+2))*sinc(-x[i]+j+1)/(2*sin(0.5)*sin(0.5)); } } printf("%.10lf\n", ans); return 0; }