Project-Euler-002矩阵快速幂+二分

题意：

​ 斐波那契数列中每一项都是前两项的和。由1和2开始生成的斐波那契数列前10项为：

​ 1, 2, 3, 5, 8, 13, 21, 34, 55, 89, …

​ 考虑该斐波那契数列中不超过四百万的项，求其中为偶数的项之和

思路：

​ 本题也想到了两个思路：

​ 思路一：由于斐波那契数列增长速度非常快，所以不超过四百万的斐波那契数顶多有100个，极限了。因此我们可以通过while(Fib[i] <= 4e6) 循环查找出为偶数的项并累加求和

​ 思路二（有点脱了裤子放屁的感觉～）：我们通过斐波那契数列的性质可知，对于该斐波那契数列：

​ 1 2 3 5 8 13 21 34…

​ 有以下两个有用的性质：

​ 1. S(n) = Fib[n + 2] - 2

​ //S(n)代表斐波那契数列前n项和 ，Fib[i]代表下标为i的斐波那契数

​ 2. 对于前下角标为3的倍数 + 2个斐波那契数，项数为偶数的和等于项数为奇数的和 + 1

​

​ 因此我们可以通过二分法套矩阵快速幂来求取4e6以内的最大斐波那契数，并返回该数在斐波那契数中的下角标n，然后用矩阵快速幂求Fib[n + 2]，求的n + 2的下角标为3的倍数 + 2，所以偶数项的和为（Fib[n+2] - 1） / 2;

代码如下：

#include <stdio.h>
#include <string.h>
#define ll long long
int N = 2;
//用于矩阵相乘
void multi(ll *a, ll *b) {		
    ll c[N][N];
    memset(c, 0, sizeof(c));
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            for (int k = 0; k < N; k++) {
                c[i][j] += *(a + i * N + k) * *(b + k * N + j);
            }
        }
    }
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            *(a + i * N + j) = c[i][j];
        }
    }
}
//矩阵快速幂
void quick(ll *a, ll n) {	
    ll ans[N][N];
    memset(ans, 0, sizeof(ans));
    for (int i = 0; i < N; i++) {
        ans[i][i] = 1;
    }
    while (n) {
        if (n & 1) {
            multi(ans[0], a);
        }
        n >>= 1;
        multi(a, a);
    }
    for (int i = 0; i < N; i++) {
        for (int j = 0; j < N; j++) {
            *(a + i * N + j) = ans[i][j];
        }
    }
}
//二分法找4e6以内最大斐波那契数，并返回该下角标
int erfen(int l, int r, ll val) {	
    while (l <= r) {
        int mid = (l + r) >> 1;
        ll root[N][N], a[N][N];
        memset(root, 0, sizeof(root));
        memset(a, 0, sizeof(a));
        root[1][0] = 1;
        a[0][1] = 1; a[1][0] = 1; a[1][1] = 1;
        quick(a[0], mid);
        multi(a[0], root[0]);
        ll res = a[1][0];
        if(res <= val) {
            l = mid + 1;
        }else {
            r = mid - 1;
        }
    }
    return r;
}

int main() {
	//感觉不会超过100个啦，就二分的下标1到100
    int res = erfen(1, 100, 4000000);	
    ll root[N][N], a[N][N], sum;
    memset(root, 0, sizeof(root));
    memset(a, 0, sizeof(a));
    root[1][0] = 1;
    a[0][1] = 1; a[1][0] = 1; a[1][1] = 1;
    if (res % 3 == 2) {
        quick(a[0], res + 2);
        multi(a[0], root[0]);
        sum = (a[1][0] - 1 )/ 2;
    }else if(res % 3 == 1) {
        //由于满足res % 3 == 2 就不再考虑
    }else {
        //由于满足res % 3 == 2 就不再考虑
    }
    printf("%lld\n", sum);	
    return 0;
}

​ 最后求出结果为：4613732

​ 第二个想法虽好，但由于最后发现下角标的值才40多！！可能二分+矩阵快速幂比纯暴力还慢。。。就酱紫～

转载请注明出处！！！

如果有写的不对或者不全面的地方 可通过主页的联系方式进行指正，谢谢！