NEUOJ 1040 Count (快速幂求第 n 个斐波那契数)

最新推荐文章于 2020-09-09 15:26:35 发布

RuHua27

最新推荐文章于 2020-09-09 15:26:35 发布

阅读量457

点赞数

分类专栏： HDOJ 快速幂

本文链接：https://blog.csdn.net/u013112097/article/details/20400625

版权

HDOJ 同时被 2 个专栏收录

32 篇文章 0 订阅

订阅专栏

快速幂

1 篇文章 0 订阅

订阅专栏

RT。。

#include <cstdio>
#include <cstring>
using namespace std;

int m, n, p;

int a[2] = {2, 1}, b[4] = {1, 1, 1, 0};

void po2()
{
    int t[4];
    t[0] = b[0]*b[0] + b[1]*b[2];
    t[1] = b[0]*b[1] + b[1]*b[3];
    t[2] = b[2]*b[0] + b[3]*b[2];
    t[3] = b[2]*b[1] + b[3]*b[3];
    for(int i = 0; i < 4; i++) b[i] = t[i]%p;
}

void mul()
{
    int temp = a[0]*b[0] + a[1]*b[1];
    a[1] = (a[0]*b[2] + a[1]*b[3]) % p;
    a[0] = temp % p;
}

int main()
{
//    freopen("1","w",stdout);
    scanf("%d", &m);
    while(m--) {
        scanf("%d%d", &n, &p);
//    p = 2009;
//    for(int n1 = 1; n1 <= 1000; n1++) {
//        n = n1;
        a[0] = 2; a[1] = 1;
        b[0] = b[1] = b[2] = 1;
        b[3] = 0;
        n--;
        while(n) {
            if(n & 1) mul();
            n /= 2;
            po2();
        }
        printf("%d\n",a[1]);
    }

    return 0;
}