HDU 1588 Gauss Fibonacci

最新推荐文章于 2020-09-05 07:25:00 发布

弱者

最新推荐文章于 2020-09-05 07:25:00 发布

阅读量789

点赞数

分类专栏：高精度

本文链接：https://blog.csdn.net/killer_in_silence/article/details/8684565

版权

高精度专栏收录该内容

2 篇文章 0 订阅

订阅专栏

这道题是一道非常有价值的高精度。通过这一个题就练习了高精度的通用算法，矩阵的运算，以及本题的关键分治法。本题让我学到了原来等比数列也可以用分治法求解。各种解释见下面的代码。

#include<iostream>
using namespace std;
struct matrix{
    _int64 a[2][2];
};
matrix A,E;
_int64    M;

//矩阵乘法，这个就不用说了。

matrix multi(matrix er,matrix et)
{
    int i,j,k;
    matrix ck;
    for(i=0;i<=1;i++)
        for(j=0;j<=1;j++)
        {
            ck.a[i][j]=0;
            for(k=0;k<=1;k++)
            {
                ck.a[i][j]=(ck.a[i][j]+(er.a[i][k]*et.a[k][j])%M)%M;
            }
        }
    return ck;
}

//分治法求等比矩阵之和，最好别偷懒直接写成return ****可能会超时，不知道原因，有空再研究再补上原因。

matrix binary_sum(_int64 k)
{
    matrix add(matrix p,matrix q);
    matrix jucheng(matrix p,_int64 k);

//如果k==1则直接返回此矩阵

    if(k==1)
        return E;

//如果k是奇数，先计算前k-1个（注意我们不算0次幂的矩阵的情况如果不懂再问）

    else if(k%2==1)
        return add(binary_sum(k-1),jucheng(E,k));
    else {
        matrix temp1=binary_sum(k/2);
        matrix temp2=multi(jucheng(E,k/2),temp1);
        return add(temp1,temp2);
    }
}
int main()
{
    _int64 k,b,n,i;
    matrix C,B,D,F;
    A.a[0][0]=1;
    A.a[0][1]=1;
    A.a[1][0]=1;
    A.a[1][1]=0;
    matrix add(matrix ap,matrix b);
    matrix jucheng(matrix q,_int64 p);
    while(cin>>k>>b>>n>>M)
    {
        B.a[0][0]=0;
        B.a[0][1]=0;
        B.a[1][0]=0;
        B.a[1][1]=0;
        E=jucheng(A,k);
        B=binary_sum(n-1);
        if(b)
        {
            F=jucheng(A,b);
            B=multi(B,F);
            B=add(B,F);
        }
        cout<<B.a[0][1]%M<<endl;
            
    }
    return 0;
}

//计算矩阵q的p次幂，也是分治法。

matrix jucheng(matrix q,_int64 p)
{
    if(p==0)
    {
        matrix g;
        g.a[0][0]=1;
        g.a[0][1]=0;
        g.a[1][0]=0;
        g.a[1][1]=1;
        return g;
    }
    else if(p==1)
    {
        return q;
    }
    else if(p%2==0)
    {
        return jucheng(multi(q,q),p/2);
    }
    else {
        return multi(jucheng(multi(q,q),p/2),q);
    }
}
//计算矩阵之和
matrix add(matrix ap,matrix b)
{
    int i,j;
    for(i=0;i<=1;i++)
        for(j=0;j<=1;j++)
            ap.a[i][j]=(ap.a[i][j]+b.a[i][j])%M;
    return ap;
}