FJNU 2011 华丽的A+B

华丽的A+B

Time Limit: 1 Seconds     Memory Limit: 32768 K Total Submit:94     Accepted:10

Description 由于没人辅导，wekooo学的东西总是那么地散，最近他又迷上了方程求解问题. 他发现一个方程的的系数总会牵连到数学里面重要的结论或者性质.为此不时有那么多神奇的东西可以挖掘. 下面给你一个方程如下： (a^2-ab-b^2)^2=1;(其中a,b均为正整数) 我们规定将解的情况按a+b的大小从小到大排序. 输入一个正整数N，求满足排序的第N小a+b并输出这个结果. 记住结果是很美丽的.o(∩_∩)o... Input 有多组测试数据，每行一个输入一个正整数N（1<=N<=10^9). Output 对应输入每行输出一个解，a+b的值.为了避免高精度我们将结果对2^64取模. Sample Input 6 28 496 8128 Sample Output 21 832040 17929816685963478168 6075550357158932568 Source wekooo@fjnu

 

 

生成几组小数据，马上得到规律，a+b满足
f[1]=2;
f[2]=3;
f[3]=5;
....
f[n]=f[n-1]+f[n-2](n>=3)
由于题目需要对2^64取mod
所以直接使用矩阵乘法就可以快速的解决
f[n]      1 1     f[n-1]
       =       *
f[n-1]    1 0     f[n-2]
一直化下去可以得到实际上:
f[n]      1 1                  f[2]
       为       的n-2次方   * 
f[n-1]    1 0                  f[1]
矩阵乘法使用二分

同时其他的类似的递推也可以使用矩阵-__________________-

以下是代码:

 

1. #include <iostream>
2. using namespace std;
3. struct mc
4. {
5.     unsigned long long a[2][2];
6.     mc operator*(const mc& m) 
7.     {
8.         mc n;
9.         int i,j,k;
10.         for(i=0;i<2;i++)
11.             for(j=0;j<2;j++)
12.                 n.a[i][j]=0;
13.         for(i=0;i<2;i++)
14.             for(j=0;j<2;j++)
15.                 for(k=0;k<2;k++)
16.                 {
17.                     n.a[i][j]+=a[i][k]*m.a[k][j];
18.                 }
19.         return n;
20.     }
21. }A;
22. mc pow(int k,mc* mat)
23.     {
24.         mc temp;
25.         if(k==1) return (*mat);
26.         temp=(*mat)*(*mat);
27.         if(k&1)
28.             return pow(k/2,&temp)*(*mat);
29.         else
30.             return pow(k/2,&temp);
31.     }
32. int main()
33. {
34.     int n;
35.     mc kk;
36.     A.a[0][0]=1;
37.     A.a[0][1]=1;
38.     A.a[1][0]=1;
39.     A.a[1][1]=0;
40.     while(cin>>n)
41.         {
42.             if(n==1)
43.                 {
44.                     cout<<2<<endl;
45.                     continue;
46.                 }
47.             if(n==2)
48.                 {
49.                     cout<<3<<endl;
50.                     continue;
51.                 }
52.             kk=pow(n-2,&A);
53.             cout<<kk.a[0][0]*3+kk.a[0][1]*2<<endl;
54.         }
55. }