Queuing HDU2604

Posted bxd123

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了Queuing HDU2604相关的知识,希望对你有一定的参考价值。

一道递推题目  

得到递推关系为  f[n]=f[n-1]+f[n-3]+f[n-4]; 

用普通的枚举算法会超时

所以要用矩阵快速幂来加速

转化为矩阵即为:  

+
1 0 1 1       F(N-1)  F(N)       
1 0 0 0  *    F(N-2)  =   F(N-1)   
0 1 0 0       F(N-3)  F(N-2)
0 0 1 0       F(N-4)  F(N-3)

   

1 0 1 1(n-4)       F(4)  F(N)       
1 0 0 0            *      F(3)  =   F(N-1)   
0 1 0 0                   F(2)  F(N-2)
0 0 1 0                   F(1)  F(N-3)

 

 

所以f(n) 为 矩阵的n-4次幂  的第一行 与已知的相乘    (n-4 为  n-3-1即可   这是差值 再加一为个数)

技术图片
#include<iostream>
#include<cstdio>
#include<cstring>

using namespace std;

int n,k;

int q[5]={0,2,4,6,9};
struct matrix{
    int arr[4][4];
};


matrix multi( matrix a,matrix b )
{
         matrix c;
    for(int i=0;i<4;i++)
        for(int j=0;j<4;j++){
            c.arr[i][j]=0;
            for(int w=0;w<4;w++)
                c.arr[i][j]=(c.arr[i][j]+a.arr[i][w]*b.arr[w][j]%k)%k;
        }
    return c;
}

int fast(matrix a,int x){
    matrix ans;
    memset(ans.arr,0,sizeof(ans.arr));
    for(int i=0;i<4;i++)ans.arr[i][i]=1;

    while(x){
        if(x&1){
            ans=multi(ans,a);
               }
        x>>=1;
        a=multi(a,a);
    }
    int sum=0;
    for(int i=0;i<4;i++)
    {
        sum+=ans.arr[0][i]*q[4-i];
        sum%=k;
    }
  return sum;
}

int main(){

   while(scanf("%d%d",&n,&k)==2)
   {
       if(n<=4){printf("%d
",q[n]%k);continue;}
       else
       {
            matrix a={1,0,1,1,1,0,0,0,0,1,0,0,0,0,1,0};
            printf("%d
",fast(a,n-4)%k);
       }
   }
    return 0;
}
View Code

矩阵还是用结构体写方便。

 










以上是关于Queuing HDU2604的主要内容,如果未能解决你的问题,请参考以下文章

hdu 2604 Queuing(矩阵快速幂)

HDU 2604 Queuing,矩阵高速幂

HDU 2604Queuing

HDU-2604 Queuing(矩阵快速幂)

hdu 2604 Queuing(推推推公式+矩阵快速幂)

Queuing HDU2604