矩阵乘法求斐波那契数列

#include<iostream>
#include<cstdio>
using namespace std;
const long long mod=1926081719260817;
long long n=1926081719260817;
long long mul(long long a,long long b)
{
	long long r=0;
	while(b)
	{
		if(b&1)
			r=(r+a)%mod;
		a=(a+a)%mod;
		b>>=1;
	}
	return r;
}
struct jz
{
	long long a[3][3];
	jz operator * (const jz &b) const
	{
		jz c;
		for(int i=0;i<2;i++)
			for(int j=0;j<2;j++)
			{
				c.a[i][j]=0;
				for(int k=0;k<2;k++)
					c.a[i][j]=(c.a[i][j]+mul(a[i][k],b.a[k][j]))%mod;
			}
		return c;
	}
}a,ans;
int main()
{
	a.a[0][0]=a.a[0][1]=a.a[1][0]=1;
	ans.a[0][0]=ans.a[1][1]=1;
	while(n)
	{
		if(n&1)
			ans=ans*a;
		a=a*a;
		n>>=1;
	}cout<<ans.a[1][0]<<endl;
	printf("%lld\n",(ans.a[1][0]+1598741160551701)%mod);
	return 0;
}