hdu 2604 Queuing 递推+矩阵快速幂
2015-12-11 07:28
447 查看
[code]#include<cstdio> #include<cstring> #include<algorithm> using namespace std; const int N=10; int n=4,m; struct Mat { int mat ; }; Mat operator * (Mat a, Mat b) { Mat c; memset(c.mat, 0, sizeof(c.mat)); int i, j, k; for(k = 0; k < n; ++k) { for(i = 0; i < n; ++i) { for(j = 0; j < n; ++j) { c.mat[i][j] += a.mat[i][k] * b.mat[k][j]; c.mat[i][j]=c.mat[i][j]%m; } } } return c; } Mat operator ^ (Mat a, int k) { Mat c; int i, j; for(i = 0; i < n; ++i) for(j = 0; j < n; ++j) c.mat[i][j] = (i == j); //初始化为单位矩阵 for(; k; k >>= 1) { if(k&1) c = c*a; a = a*a; } return c; } int main() { int k,i,j; while(~scanf("%d%d",&k,&m)) { if(k==0) printf("%d\n",k); else if(k==1) printf("%d\n",2%m); else if(k==2) printf("%d\n",4%m); else if(k==3) printf("%d\n",6%m); else { Mat a,b; memset(a.mat,0,sizeof(a.mat)); a.mat[0][0]=1; a.mat[0][2]=1; a.mat[0][3]=1; a.mat[1][0]=1; a.mat[2][1]=1; a.mat[3][2]=1; a=a^(k-4); /*for(i=0; i<4; i++) { for(j=0; j<4; j++) printf("%d ",a.mat[i][j]); printf("\n"); }*/ memset(b.mat,0,sizeof(b.mat)); b.mat[0][0]=9; b.mat[1][0]=6; b.mat[2][0]=4; b.mat[3][0]=2; int ans=0; for(i=0; i<4; i++) { ans=(ans+b.mat[i][0]*a.mat[0][i])%m; //printf("%d %d %d\n",b.mat[i][0],a.mat[0][i],ans); } printf("%d\n",ans%m); } } return 0; }
相关文章推荐
- oblique perspective projection
- LightOJ 1061 1061 - N Queen Again(状压dp)
- 【NGUI】当打包一个特别大的图集时候该怎么减小图集大小?
- mysql_query设置字符集为utf8,解决中文乱码的问题
- 《Elasticsearch: The Definitive Guide》读书笔记2
- UESTC 839
- 爬爬爬之路:UI(十七) 数据持久化(SQLite3数据库实现)
- UESTC 838
- IOS UIColor 自定义颜色
- IOS基础-导航栏控制器 UINavigationController
- UILongPressGestureRecognizer的selector多次调用解决方法
- 【译】UI设计基础(UI Design Basics)--自动适配与布局(Adaptivity and Layout)(四)
- 我的iOS学习历程 - UICollection详解
- 随意细解:UI -- 数据库
- IOS UISearchBar
- UI基础-数据库
- IOS UISearchBar(一)
- UI-第三讲:视图控制器
- IOS UITableView (五)
- 对web的request的全部值进行了提取