fork(1) download
copy 
  1. #include <iostream>
  2. #include <cmath>
  3. #include <iostream>
  4. #include <cmath>
  5.  
  6. template<class T>
  7. // F(1) == F(2) == 1 なので X != 1 とする
  8. T Fi(const T& X) {
  9. 	T Phi = (1 + std::sqrt(5)) / 2;
  10. 	T A01 = std::pow(-1, std::fmod(X + 1, 3));
  11. 	return (std::log(std::sqrt(5) * X + std::sqrt(5 * (X * X) - 4 * A01)) - std::log(2)) / std::log(Phi);
  12. }
  13. template<class T>
  14. T F(const T& N) {
  15. 	T Phi = (1 + std::sqrt(5)) / 2;
  16.  
  17. 	return (std::pow(Phi, N) - std::pow(-Phi, -N)) / std::sqrt(5);
  18. }
  19.  
  20. int main() {
  21. 	for (std::size_t i = 2; i < 128; i++) {
  22. 		double T = F<double>(i);
  23. 		std::cout << T << ',' << Fi<double>(T) << std::endl;
  24.  
  25. 	}
  26.  
  27. 	return 0;
  28.  
  29. }
  30.  
Success #stdin #stdout 0s 4192KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
1,1
2,3
3,4
5,-nan
8,-nan
13,-nan
21,-nan
34,-nan
55,-nan
89,-nan
144,-nan
233,-nan
377,-nan
610,-nan
987,-nan
1597,-nan
2584,-nan
4181,-nan
6765,-nan
10946,-nan
17711,-nan
28657,-nan
46368,-nan
75025,-nan
121393,-nan
196418,-nan
317811,-nan
514229,-nan
832040,-nan
1.34627e+06,-nan
2.17831e+06,-nan
3.52458e+06,-nan
5.70289e+06,-nan
9.22747e+06,-nan
1.49304e+07,-nan
2.41578e+07,-nan
3.90882e+07,-nan
6.3246e+07,-nan
1.02334e+08,-nan
1.6558e+08,-nan
2.67914e+08,-nan
4.33494e+08,-nan
7.01409e+08,-nan
1.1349e+09,-nan
1.83631e+09,-nan
2.97122e+09,-nan
4.80753e+09,-nan
7.77874e+09,-nan
1.25863e+10,-nan
2.0365e+10,-nan
3.29513e+10,-nan
5.33163e+10,-nan
8.62676e+10,-nan
1.39584e+11,-nan
2.25851e+11,-nan
3.65435e+11,-nan
5.91287e+11,-nan
9.56722e+11,-nan
1.54801e+12,-nan
2.50473e+12,-nan
4.05274e+12,-nan
6.55747e+12,-nan
1.06102e+13,-nan
1.71677e+13,-nan
2.77779e+13,-nan
4.49456e+13,-nan
7.27235e+13,-nan
1.17669e+14,-nan
1.90392e+14,-nan
3.08062e+14,-nan
4.98454e+14,-nan
8.06516e+14,73
1.30497e+15,74
2.11149e+15,-nan
3.41645e+15,-nan
5.52794e+15,77
8.94439e+15,78
1.44723e+16,79
2.34167e+16,80
3.78891e+16,81
6.13058e+16,82
9.91949e+16,83
1.60501e+17,84
2.59695e+17,85
4.20196e+17,86
6.79892e+17,87
1.10009e+18,88
1.77998e+18,89
2.88007e+18,90
4.66005e+18,91
7.54011e+18,92
1.22002e+19,93
1.97403e+19,94
3.19404e+19,95
5.16807e+19,96
8.36211e+19,97
1.35302e+20,98
2.18923e+20,99
3.54225e+20,100
5.73148e+20,101
9.27373e+20,102
1.50052e+21,103
2.42789e+21,104
3.92841e+21,105
6.35631e+21,106
1.02847e+22,107
1.6641e+22,108
2.69257e+22,109
4.35668e+22,110
7.04925e+22,111
1.14059e+23,112
1.84552e+23,113
2.98611e+23,114
4.83163e+23,115
7.81774e+23,116
1.26494e+24,117
2.04671e+24,118
3.31165e+24,119
5.35836e+24,120
8.67001e+24,121
1.40284e+25,122
2.26984e+25,123
3.67267e+25,124
5.94251e+25,125
9.61519e+25,126
1.55577e+26,127
https://ideone.com/qN8CCx
language:
C++14 (gcc 8.3)
created:
4 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!