【Best Coder】#29 B GTY's birthday gift（快速幂|mod的时候记得负！）

最新推荐文章于 2018-10-04 16:37:27 发布

jason_star

最新推荐文章于 2018-10-04 16:37:27 发布

阅读量943

点赞数

分类专栏：寒假狂练 ACM_BC 方法学习文章标签：二进制拆分斐波那契

本文链接：https://blog.csdn.net/u013611908/article/details/43614307

版权

方法学习同时被 3 个专栏收录

33 篇文章 0 订阅

订阅专栏

寒假狂练

29 篇文章 0 订阅

订阅专栏

ACM_BC

9 篇文章 0 订阅

订阅专栏

题目大意：查看相关场次即可看到。
思路：推公式的题目，可以用快速幂加公式快速解决，也可以用二进制拆分运算的方法加快速度。
需要注意的一点在于：今后在mod之后有涉及到运算的都要加上一个mod之后再mod，或者统一都加一个mod
顺便复习一下二进制拆分的方法！！
二进制拆分的做法AC代码：

#define _CRT_SECURE_NO_WARNINGS
#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
long long f[1200000];
int a[100010];
#define mod 10000007

int main()
{
    f[0] = 0;
    f[1] = 1;
    f[2] = 2;
    f[3] = 4;
    f[4] = 3;
    f[5] = 5;;
    for (int i = 6; i < 1200000; i++)
    {
        f[i] = (f[i - 1] + f[i - 2]+mod) % mod;
        f[i - 2] = (f[i - 3] + f[i - 2] + mod) % mod;
    }
    int n, m;
    long long sum = 0;
    while (cin >> n >> m)
    {
        sum = 0;
        for (int i = 0; i < n; i++)
        {
            scanf("%d", &a[i]);
            sum += a[i];
            sum %= mod;
        }
        sort(a, a + n);
        long long b = a[n - 2], c = a[n - 1];
        int k = 1;
        int mid = 1;
        while (m)
        {
            k = 1;
            long long bb = b, cc = c;
            while (k <= m)
            {
                if (k>(1<<20))
                    break;
                mid = k;
                sum = (sum + (f[k] * b + (f[k + 1] - 1)*c + mod) % mod) % mod;
                bb = b; cc = c;
                /*
                f[k]-f[k-1]会出现负的情况，加上一个mod就可以解决！
                */
                c = ((((f[k] - f[k - 1] + mod)%mod)*bb) % mod + (((f[k + 1] - f[k] + mod))%mod*cc) % mod + mod) % mod;
                b = ((((f[k - 1] - f[k - 2] + mod) % mod)*bb) % mod + (((f[k] - f[k - 1] + mod)%mod)*cc) % mod) % mod;
                m -= k;
                k <<= 1;
            }
        }
        cout << sum << endl;
    }
};