斐波那契之通项公式

以前没有发现，斐波那契通项公式能直接KO一题。记录一下，许多的过程细节值得回味。
nyist 461 Fibonacci数列（四）

http://acm.nyist.net/JudgeOnline/problem.php?pid=461

大意：求出斐波那契f[n]的前4个数字

因为不是后几位所以不能矩阵快速幂。

利用通项公式：

#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const double f=sqrt(5.0);
int main()
{
    int n;
    while(cin>>n){
        double p1=pow((1+f)/2,n);