ljw的博客

这个人既想要学软件又想要搞硬件,吃枣药丸

矩阵快速幂,longlong初始值

一道堵了我很久的题,被别人一说,就开窍了再见

题目描述

    有一群一共 n 个妖精和兽排成一队,其中这一队列可以全是妖精或者全是兽,当然也可以是妖精和兽并存。

    排队时还规定:兽不能落单(落单了就要被妖精砍死了),再明确一些,就是要么这个队列里面没有兽,如果有兽的话,必须 至少 两个兽肩并肩排在一起。

    请问有多少种排法呢?由于数可能很大,请你对这个结果对 1e9 + 7 取模。

输入

多组输入,请处理到文件结束。

数据组数不超过 20000。

一行一个数,代表 (n <= 10^10)。

输出

输出方案数,对 1e9 + 7 取模。


输入样例

1
2
3

输出样例

1
2
4

提示

假设妖精为 G,兽为 B。

n = 1 的情况: G n = 2 的情况: BB GG n = 3 的情况: BBG GGG GBB BBB

注意 n = 1 中没有 B 这种情况,因为这样算落单。


--------------------------------我是分界线-------------------------------------------


从它给的范围((n <= 10^10)),很显然这是一个找规律或者是快速幂的题目。

题目说B(兽,借用题目的提示)不能落单,必须连续的2个或者两个以上才合法,举例的话GBBGBB是合法的。

那么开始找规律:

1. 假设Y(n)表示为,输入为n的输出。

2. 假设T(n)表示为:合法的排列,且最后一个为G。例如:GGGG。

3. 假设S(n)表示为:“缺B合法排列”的排列数。表示在序列的最后,额外添加一个B,能够使这个序列变成合法的,例如:GBBGGGB。

4. 已知:Y(n-1),T(n-1),S(n-1)。

5. 求Y(n):那么我在这Y(n-1)个合法的序列中最后面添加一个G,这个新序列也是合法的,共有Y(n-1)中;我在这Y(n-1)个合法的序列中最后面添加一个B,只有[Y(n-1)-T(n-1)]种新序列是合法的;然后我在“缺B合法排列”后面加一个B,就会有新的S(n-1)中合法排列。Y(n)=Y(n-1)+[Y(n-1)-T(n-1)]+S(n-1)。

6. 求T(n):显然,T(n)=Y(n-1)。在n-1长的合法序列末加G。

7. 求S(n):显然,S(n)=T(n-1)。在n-1长的合法序列且以G结尾后面加个B。

8. 得到Y(n)=2*Y(n-1)-Y(n-2)+Y(n-3)。

得到上面那个式子后我就束手无策了.......在同学的点拨下,想起了矩阵快速幂大笑

其他的没什么说的,一个点被坑了一下,就是要使用long long。初始化的时候不知道为什么,初始化写 0L 也自动转化为int型数据了,最后我尝试初始化为 1000000007000L 就没事了,这也算是一点收获。

最后AC代码:

#include <iostream>

using namespace std;

#define M (long long)(1e9+7)

typedef struct Mat{
    long long c[3][3];
}Mat;

Mat calone(Mat a, Mat b);
long long speedcal(Mat s, long long n);

int main()
{
    long long N;
    while(cin>>N){
        Mat sta;//矩阵初始化操作很low...
        sta.c[0][0]=4L;sta.c[0][1]=2L;sta.c[0][2]=1L;
        sta.c[1][0]=0L;sta.c[1][1]=0L;sta.c[1][2]=0L;
        sta.c[2][0]=0L;sta.c[2][1]=0L;sta.c[2][2]=0L;
        long long ans=speedcal(sta,N);
        cout<<ans<<endl;
    }
    return 0;
}

long long speedcal(Mat s, long long n)
{
    if(n==1)
        return s.c[0][2];
    else if(n==2)
        return s.c[0][1];
    else if(n==3)
        return s.c[0][0];
    else{
        Mat m,t=s;
        m.c[0][0]=2L;m.c[0][1]=1L;m.c[0][2]=0L;
        m.c[1][0]=-1L;m.c[1][1]=0L;m.c[1][2]=1L;
        m.c[2][0]=1L;m.c[2][1]=0L;m.c[2][2]=0L;
        n-=3;
        while(n>0){
            if(n%2==1)
                t=calone(t,m);
            m=calone(m,m);
            n/=2;
        }
        return t.c[0][0];
    }
    return -1L;
}

Mat calone(Mat a, Mat b)
{
    Mat c;
    for(int i=0;i<3;i++){
        for(int j=0;j<3;j++){
            long long t=1000000007000L;//就是这里
            for(int k=0;k<3;k++){
                t=(t+a.c[i][k]%M*b.c[k][j]%M)%M;
            }
            c.c[i][j]=t;
            if(c.c[0][0]<0)
                cout<<"error"<<endl;
        }
    }
    return c;
} 
希望每天进步一点点~~







阅读更多
版权声明: https://blog.csdn.net/wxwd14388/article/details/79962878
文章标签: 编程 C
个人分类: 软件
上一篇基于FPGA的SSRAM读写速度测试实验
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭