分块+矩阵快速幂

http://acm.hdu.edu.cn/showproblem.php?pid=6395

分块+矩阵快速幂_第1张图片

P/i会有相同的一段,所以想到分块,这个也经常用

\begin{bmatrix} F(n)\\ F(n-1) \\ 1 \end{bmatrix}= \begin{bmatrix} d &c & x\\ 1&0 &0 \\ 0& 0 &1 \end{bmatrix}\begin{bmatrix} F2\\ F1\\ 1 \end{bmatrix}

#include
#include 
#include 
#include 
#include  
#include   
#include   
#include   
#include   
#include    
#include      
#include      
using namespace std;

typedef long long ll;
const int Mod=1e9+7;
int A,B,C,D,n,P;

struct Mat
{
    int t[3][3];
    
    Mat(){memset(t,0,sizeof t);}
}I;

Mat operator * (Mat a,Mat b)
{
    Mat c;
    for(int i=0;i<3;i++)
        for(int j=0;j<3;j++)
        {
            ll t=0;
            for(int k=0;k<3;k++)
                t+=(ll)a.t[i][k]*b.t[k][j];
            c.t[i][j]=t%Mod;
        }
    return c;
}

Mat Pow(Mat a,int b)
{
    Mat c=I;
    while(b)
    {
        if(b&1)
            c=c*a;
        a=a*a;b>>=1;
    }
    return c;
}

void solve()
{
    scanf("%d%d%d%d%d%d",&A,&B,&C,&D,&P,&n);
    if(n==1)
    {
        cout<>t;
    while(t--)
        solve(); 
    return 0;
}

 

你可能感兴趣的:(数学-矩阵,比赛)