Fibonacci Numbers（斐波那契数列 前四位 后四位）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <math.h>
using namespace std;
#define mod 1000000000
struct mat
{
	__int64 a[2][2];
};
mat mul(mat a,mat b)
{
	mat ret;
	for(int i=0;i<2;i++)
	{
		for(int j=0;j<2;j++)
		{
			ret.a[i][j] = (a.a[i][0]*b.a[0][j]+a.a[i][1]*b.a[1][j])%mod;
		}
	}
	return ret;
}
mat mpower(mat a,int x)
{
	mat I;
	I.a[0][0] = I.a[1][1] = 1;
	I.a[0][1] = I.a[1][0] = 0;
	while(x)
	{
		if(x&1) I = mul(I,a);
		x>>=1;
		a = mul(a,a);
	}
	return I;
}
int main()
{
	int n;
	while(cin>>n)
	{
		mat tmp;
		tmp.a[0][0] = tmp.a[0][1] = tmp.a[1][0]=1;
		tmp.a[1][1] = 0;
		if(n==0||n==1) 
		{
			cout<<n<<endl;
			continue;
		}
		tmp = mpower(tmp,n-1);
		if(n<=39)
		{
			int sum = tmp.a[0][0];
			cout<<sum<<endl;
		}
		else
		{
			__int64 a1 = tmp.a[0][0];
			double ai = log10((1*1.0)/sqrt(5*1.0));  
        	double bi = n*(log10((1+sqrt(5*1.0))*1.0/2))+ai;  
        	bi = bi-floor(bi);  
        	bi = pow(10,bi);  
        	int a2 = (int)(bi*1000);  
			cout<<a2<<"...";
			int p[4];
			for(int i=0;i<4;i++)
			{
				p[i] = a1%10;
				a1/=10;
			}
			for(int i=3;i>=0;i--) cout<<p[i];
			cout<<endl;
		}
	}
	return 0;
}