C++の練習を兼ねて, 第10回 アルゴリズム実技検定 過去問 の 問題G (方程式) を解いてみた.
■感想.
1. 問題Gは, 方針を絞り込めたので, AC版に到達できたと思う
2. 個人的には, ニュートン法について, 学習できたので, 良かったと思う.
3. 引き続き, 時間を見つけて, 過去問の学習を進めていきたいと思う.
本家のサイト 第10回 アルゴリズム実技検定 過去問 の 解説 の 各リンク を ご覧下さい.
■C++版プログラム(問題G/AC版).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 |
// C++(GCC 9.2.1) #include <bits/stdc++.h> using namespace std; #define repex(i, a, b, c) for(int i = a; i < b; i += c) #define repx(i, a, b) repex(i, a, b, 1) #define rep(i, n) repx(i, 0, n) #define repr(i, a, b) for(int i = a; i >= b; i--) int main(){ // 1. 入力情報. double a, b, c; scanf("%lf %lf %lf", &a, &b, &c); // 2. 方程式. auto f = [&](double x){ return a * x * x * x * x * x + b * x + c; }; // 3. 導関数. auto g = [&](double x){ return 5.0 * a * x * x * x * x + b; }; // 4. ニュートン法. // x[n + 1] = x[n] - f[x[n]] / f'[x[n]] double ans = 1.5; rep(i, 32) ans -= f(ans) / g(ans); // 5. 出力. printf("%.18lf\n", ans); return 0; } |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 |
[入力例] 32 2 -246 [出力例] 1.500000000000000000 ※AtCoderのテストケースより [入力例] 12 3 -45 [出力例] 1.279562760087743278 ※AtCoderのテストケースより [入力例] 1 1 -3 [出力例] 1.132997565885065283 [入力例] 1 2 -5 [出力例] 1.208917813386894968 [入力例] 12345678 23456789 -200000000 [出力例] 1.670925787919859573 |
■参照サイト
第10回 アルゴリズム実技検定 過去問