ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

hdu 2604 Queuing (矩阵快速幂)

2021-04-03 15:01:48  阅读:155  来源: 互联网

标签:2604 tmp hdu mat int Queuing ++ MXM base


Queuing

分析

  1. 一般递推式
    在这里插入图片描述

  2. 向量递推式
    在这里插入图片描述

代码

#include<bits/stdc++.h>
using namespace std;
#define MXM 4
#define mod(x) ((x)%M)
int L, M, f[5]={0, 2, 4, 6, 9}, A[MXM][MXM]={{1,0,1,1},{1,0,0,0},{0,1,0,0},{0,0,1,0}};
struct mat{ 
    int d[MXM][MXM];
    mat operator*(const mat x){
        mat ret;
        int tmp;
        for(int i = 0; i < MXM; i++){
            for(int j = 0; j < MXM; j++){
                tmp = 0;
                for(int k = 0; k < MXM; k++){
                    tmp = mod(tmp + d[i][k]* x.d[k][j]);
                }
                ret.d[i][j] = tmp;
            }
        }
        return ret;
    }
    void init_unit(){
        for(int i = 0; i < MXM; i++)
            for(int j = 0; j < MXM; j++)
                d[i][j] = i == j ? 1 : 0;
    }
    void init(){
        for(int i = 0; i < MXM; i++){
            for(int j = 0; j < MXM; j++){
                this->d[i][j] = A[i][j];
            }
        }
    }
}ma;
mat matrixPow(mat base, int pow){
    mat res;
    res.init_unit();
    while(pow){
        if(pow & 1) res = res * base;
        base = base * base;
        pow >>= 1;
    }
    return res;
}
int main(){
	int tmp;
    while(scanf("%d%d", &L, &M) == 2){
        if(L <= MXM){
			printf("%d\n", f[L]%M);
			continue;
		}
        ma.init();
        ma = matrixPow(ma, L-4);
        int ans = 0;
        for(int i = 0; i < MXM; i++)
            ans = mod(ans+ma.d[0][i]*f[MXM-i]);
        printf("%d\n", ans);
    }
	return 0;
}

标签:2604,tmp,hdu,mat,int,Queuing,++,MXM,base
来源: https://blog.csdn.net/jpphy0/article/details/115412339

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有