矩阵快速幂 HDU3483

  1 #include <iostream>
  2 #include <cstring>
  3 
  4 using namespace std;
  5 
  6 //矩阵大小上限
  7 const int SIZ=100;
  8 int MOD;
  9 
 10 //矩阵大小为n*m,初始化全部为0
 11 struct mat
 12 {
 13     int n,m;
 14     long long  ar[SIZ][SIZ];
 15     mat()
 16     {
 17         memset(ar,0,sizeof(ar));
 18         n=m=SIZ;
 19     };
 20 };
 21 
 22 //矩阵乘法
 23 mat operator *(mat a,mat b)
 24 {
 25     mat c;
 26     c=mat();
 27     c.n=a.n;
 28     c.m=b.m;
 29     for(int i=1;i<=a.n;i++)
 30         for(int j=1;j<=b.m;j++)
 31             for(int k=1;k<=a.m;k++)
 32             {
 33                 c.ar[i][j]+=(a.ar[i][k]*b.ar[k][j])%MOD;
 34                 c.ar[i][j]%=MOD;
 35             }
 36     return c;
 37 }
 38 
 39 //矩阵加法
 40 mat operator +(mat a,mat b)
 41 {
 42     mat c;
 43     c=mat();
 44     c.n=a.n;
 45     c.m=a.m;
 46     for(int i=1;i<=a.n;i++)
 47         for(int j=1;j<a.m;j++)
 48             c.ar[i][j]=a.ar[i][j]+b.ar[i][j];
 49     return c;
 50 }
 51 
 52 //矩阵快速幂
 53 mat operator ^(mat a,int k)
 54 {
 55     mat c;
 56     c=mat();
 57     c.n=a.n;
 58     c.m=a.m;
 59     for(int i=1;i<=a.n;i++)
 60         c.ar[i][i]=1;
 61     while(k)
 62     {
 63         if(k&1)
 64             c=c*a;
 65         a=a*a;
 66         k/=2;
 67     }
 68     return c;
 69 }
 70 
 71 long long  tarr[100][100];
 72 long long  C(long long  a,long long  b)
 73 {
 74     memset(tarr,0,sizeof(tarr));
 75     for(int i=0;i<=a;i++)
 76         tarr[i][0]=1;
 77     for(int i=1;i<=a;i++)
 78         for(int t=1;t<=i;t++)
 79         {
 80             tarr[i][t]=(tarr[i-1][t-1]+tarr[i-1][t])%MOD;
 81         }
 82     return tarr[a][b];
 83 }
 84 
 85 int main()
 86 {
 87     int n,x;
 88     while(cin>>n>>x>>MOD&&(n!=-1||x!=-1||MOD!=-1))
 89     {
 90         mat tt;
 91         tt=mat();
 92         tt.m=tt.n=x+2;
 93         for(int i=1;i<tt.m;i++)
 94         {
 95             for(int t=1;t<=i;t++)
 96             {
 97                 tt.ar[i][t]=C(i-1,t-1)*x%MOD;
 98             }
 99         }
100         tt.ar[tt.m][tt.m-1]=tt.ar[tt.m][tt.m]=1;
101         mat ans;
102         ans=mat();
103         ans.n=x+2;
104         ans.m=1;
105         for(int i=1;i<x+2;i++)
106         {
107             ans.ar[i][1]=x;
108         }
109         ans.ar[x+2][1]=0;
110         cout<<((tt^(n))*ans).ar[x+2][1]%MOD<<endl;
111     }
112     return 0;
113 }
View Code

 

转载于:https://www.cnblogs.com/wsruning/p/4678023.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/431561.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

哲学家就餐问题python_Python实现哲学家就餐问题实例代码

哲学家就餐问题&#xff1a; 哲学家就餐问题是典型的同步问题&#xff0c;该问题描述的是五个哲学家共用一张圆桌&#xff0c;分别坐在五张椅子上&#xff0c;在圆桌上有五个盘子和五个叉子&#xff08;如下图&#xff09;&#xff0c;他们的生活方式是交替的进行思考和进餐&am…

mongodb 库数量限制_MongoDB in 数量限制

1、查询语句本身其实是一个document&#xff0c; 最大为16MB(3.4&#xff0c;4.0 的限制&#xff0c;官方文档)2、查询语句本身&#xff0c;也就是{ : { $in : [] }}&#xff0c; 大小为 22字节3、每增加一个字段名&#xff0c;增加一位字节4、in里面每增加一个参数&#xff0…

Webx框架:Pipeline基本介绍

Pipeline。它是管道的含义。一个管道阀门可以安装非常多。有许多可能的分支。它是用来控制页处理。它需要在被定义pipeline.xml文件。该文件是为每个阀的标签。该文件可以放一些简单的控制语句。在项目。这条管线有以下配置足够好。 <services:pipeline xmlns"http://w…

fpga摄像头模块_FPGA开源项目:双目测距(一)之双目图像采集显示以及图片保存...

1.简述这个项目是大三下学期暑假(也就是2019年8份)完成的&#xff0c;当时的视频效果已发布在bilibili上&#xff0c;这是我们的省级的科研立项&#xff0c;其实就我一个人负责完成。发布bilibili后很多人比较感兴趣&#xff0c;打算年初回学校完成毕设期间开源的&#xff0c;一…

python查找公共前缀_Python实现查找字符串数组最长公共前缀示例

本文实例讲述了Python实现查找字符串数组最长公共前缀。分享给大家供大家参考&#xff0c;具体如下&#xff1a;编写一个函数来查找字符串数组中的最长公共前缀。 class Solution:def longestCommonPrefix(self, strs):""":type strs: List[str]:rtype: str&quo…

博客搬家了

新博客&#xff1a;http://lindawhite.gq转载于:https://www.cnblogs.com/linda586586/p/4682027.html

github ssh 配置_Github远程仓库克隆更新本机,SSH协议免密操作配置和注意事项

Github远程仓库-克隆远程仓库到本机 【重点】目标如何克隆远程仓库到本机呢&#xff1f;小结下载项目的命令&#xff1a;git clone https://gitee.com/kekesam/sassweb777.git注意&#xff1a;1&#xff1a;它会自动创建本地仓库&#xff1b;2&#xff1a;它也会自动和远程仓库…

mysql查询库表变更信息_PythonMySQL进行数据库表变更和查询

这篇文章主要介绍了Python MySQL进行数据库表变更和查询的相关资料,需要的朋友可以参考下Python连接MySQL&#xff0c;进行数据库表变更和查询&#xff1a;python mysql insert delete query&#xff1a;#!/usr/bin/pythonimport MySQLdbdef doInsert(cursor,db):#insert# Prep…

CSS实现不固定宽度和高度的自动居中

有时候我们需要实现下面这种效果&#xff1a;嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#xff0c;撑大高度不让你剧中嘎嘎&#…

golang中文文档_【译】Go 语言源码贡献官方指导文档

以前给 Go 语言项目源码提交过一些 commits&#xff0c;期间阅读他们的官方指导文档的时候觉得这篇指导文档可以作为绝佳的关于大型软件项目的规范管理的参考&#xff0c;因为最近又提交了几个 commits&#xff0c;就又把这篇文档再看了一遍&#xff0c;有感于 Go 团队在项目管…

mysql增加最大连接数_mysql最大连接数怎么设置

设置mysql最大连接数的方法&#xff1a;首先打开mysql的控制台&#xff1b;然后输入语句【set GLOBAL max_connections1000;】即可直接设置最大连接数。通常&#xff0c;mysql的最大连接数默认是100, 最大可以达到16384&#xff1b;如果我们想修改mysql的最大连接数&#xff0c…

Hbase Interface HConnection

HTablePool 在Hbase 0.94、0.95、0.97被废弃&#xff0c;在0.98中被清除&#xff08; HTablePool 对比HConnection.getTable&#xff09;&#xff0c;hbase0.98 HTablePool会被逐渐废弃&#xff0c;直到0.98的发行版本中会被清楚。使用新的API&#xff0c;HConnection.getTable…

让一个动画一直执行的属性是_iOS 动画 一

View AnimationsAnimatable properties• bounds: 改变 bounds 属性可以在当前 view 内改变子视图等的相对位置。• frame: 改变 frame 可以移动或者缩放 view 。• center: 当你想移动 view 到屏幕的新位置时&#xff0c;可以改变此属性。• backgroundColor: 背景颜色。• al…

阿里云java mysql环境_阿里云搭建centos java mysql nginx环境

1.上传下载yum install lrzszsz下载rz上传2.安装压缩命令yum install -y unzip zip;3.安装javarpm -ivh jdk-8u161-linux-x64.rpm4.安装mysqlyum install libaioyum -y remove mariadb-libs-5.5*rpm -ivh mysql-community-common-5.7.21-1.el7.x86_64.rpmrpm -ivh mysql-commun…

ide快捷键_一款好用的IDE怎么可以没有代码提示?

我在使用过IntelliJ IDEA(一款Java开发工具)之后&#xff0c;感觉Delphi的IDE在编写代码方面尤其是代码提示方面和IntelliJ IDEA相比真的天差地别&#xff0c;所以决定认真的研究一下Delphi的代码提示功能一般情况下我们使用Delphi的开发工具都会安装cnpack专家包&#xff0c;但…

IOS开发学习记录第4天之C语言学习

&#xff08;一&#xff09;、今天我们要学习的主要包括一下内容&#xff1a; 1、标示符概念及其命名原则 在C语言中&#xff0c;符号常量&#xff0c;变量&#xff0c;数组&#xff0c;函数等都需要一定的名称&#xff0c;我们把这种名称称之为标识符。 标识符划分&#xff1a…

mysql pmm 布署_给 mysql 安装 pmm 监控

PMM 说明PMM(Percona Monitoring and Management) 是一款监控和分析 MySQL 服务的一套工具&#xff0c;可以从图形化的方式看到 MySQL 服务的各种性能指标&#xff0c;慢查询分析、连接数、线程状态、查询信息、缓存信息等等&#xff0c;对分析 MySQL 运行时问题很有帮助。PMM …

java 自定义注解_Java注解

前言近日在阅读开源项目&#xff0c;发现项目里好多奇奇怪怪的注解(DataScope、Log...)看得我一脸懵&#xff0c;不知道大家是否也有过这样的经历&#xff0c;回想了一下&#xff0c;发现自己对于注解的知识&#xff0c;好像只停留在Override。。。异常尴尬&#xff0c;所以今天…

JS 判断滚动底部并加载更多效果。。。。。。。。。

JS 判断滚动底部并加载更多效果。。。。。。。。。 <html lang"zh-cn"> <head> <meta http-equiv"Content-Type" content"text/html; charsetutf-8"/> <title>js网页滚动条滚动事件 </title> <style type"…

mysql 格林时间转换_oracle中将格林威治时间转化为一般时间

下面是编程之家 jb51.cc 通过网络收集整理的代码片段。编程之家小编现在分享给大家&#xff0c;也给大家做个参考。to_char(trunc(b.starttime/1000/60/60/24)TO_DATE(1970-01-01 00:00:00,yyyy-MM-dd hh24:mi:ss),yyyy-mm-dd)|| ||lpad(mod(trunc(b.starttime/1000/60/60),24)…