xmuC语言程序实践week 1 大作业
时间:2022-07-25
本文章向大家介绍xmuC语言程序实践week 1 大作业,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
算法提高 矩阵乘方
描述
给定一个矩阵A,一个非负整数b和一个正整数m,求A的b次方除m的余数。 其中一个nxn的矩阵除m的余数得到的仍是一个nxn的矩阵,这个矩阵的每一个元素是原矩阵对应位置上的数除m的余数。 要计算这个问题,可以将A连乘b次,每次都对m求余,但这种方法特别慢,当b较大时无法使用。下面给出一种较快的算法(用A^b表示A的b次方): 若b=0,则A^b%m=I%m。其中I表示单位矩阵。 若b为偶数,则A^b%m=(A^(b/2)%m)^2%m,即先把A乘b/2次方对m求余,然后再平方后对m求余。 若b为奇数,则A^b%m=(A^(b-1)%m)*a%m,即先求A乘b-1次方对m求余,然后再乘A后对m求余。 这种方法速度较快,请使用这种方法计算A^b%m,其中A是一个2x2的矩阵,m不大于10000。
输入
输入描述: 输入第一行包含两个整数b, m,第二行和第三行每行两个整数,为矩阵A。 输入样例: 2 2 1 1 0 1
输出
输出描述: 输出两行,每行两个整数,表示A^b%m的值。 输出样例: 1 0 0 1
输入样例 1
参考上文
输出样例 1
参考上文
提示
HINT:时间限制:1.0s 内存限制:512.0MB
矩阵快速幂板子题~~~
#include<bits/stdc++.h>
#define LL long long
using namespace std;
inline int read()//快读
{
char ch=getchar();int s=0,w=1;
while(ch<48||ch>57){if(ch=='-')w=-1;ch=getchar();}
while(ch>=48&&ch<=57){s=(s<<1)+(s<<3)+(ch^48);ch=getchar();}
return s*w;
}
inline void write(int x)//快写
{
if(x<0)putchar('-'),x=-x;
if(x>9)write(x/10);
putchar(x%10+48);
}
struct Mat //定义矩阵结构体
{
int m[101][101];
};
int n,b,m; //n*n方阵 A,求矩阵A^b%m
Mat a,e; //a是输入的矩阵,e是输出的矩阵
Mat Mul(Mat x,Mat y)
{
Mat c;
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
c.m[i][j] = 0;
}
}
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j){
for(int k=1;k<=n;++k){
c.m[i][j] += (x.m[i][k]*y.m[k][j]);
c.m[i][j]%=m;
}
}
}
return c;
}
Mat pow(Mat x,int y)
{
Mat ans;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
i==j?ans.m[i][j]=1:ans.m[i][j]=0;//单位矩阵*任何矩阵=任何矩阵本身 单位矩阵定义:对角线上元素为1,其他为0
}
}
/* for(int i=1;i<=n;i++) //注释掉部分可供调试
{
for(int j=1;j<=n;j++)
{
cout<<ans.m[i][j]<<" ";
}cout<<endl;
}*/
if(!y)
{
ans=Mul(ans,x);
return ans;
}
while(y) //矩阵快速幂模板
{
if(y&1)
{
ans=Mul(ans,x);
/*for(int i=1;i<=n;i++) //注释掉部分可供调试
{
for(int j=1;j<=n;j++)
{
cout<<ans.m[i][j]<<" ";
}cout<<endl;
}*/
}
x = Mul(x,x);
y>>=1;
}
return ans;
}
int main()
{
n=read(),b=read(),m=read();
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
a.m[i][j]=read();
}
}
/* for(int i=1;i<=2;i++)
{
for(int j=1;j<=2;j++)
{
cout<<a.m[i][j]<<" ";
}
cout<<endl;
}*/
Mat k=pow(a,b);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
{
j==n?cout<<k.m[i][j]%m:cout<<k.m[i][j]%m<<" ";//注意格式,最后一个没有空格输出
}
cout<<endl;
}
return 0;
}
复杂度分析:暴力解法o(n^3*m)->矩阵快速幂o(n^3logm)
矩阵相乘有个o(n^2.7)的听说
ac截图,没办法,作业要这个
- UE4新手之编程指南
- redis 学习笔记(5)-Spring与Jedis的集成
- Linq to SQL 查询Tips
- SQLite 预写式日志
- java调用.net asmx / wcf
- mybatis3.2.8 与 hibernate4.3.6 混用
- mybatis的物理分页:mybatis-paginator
- 使用 WMI 进行诊断WCF
- java:快速文件分割及合并
- 暴涨210倍的一个数字货币正悄无声息崛起
- QT Creator 快速入门教程 读书笔记(一)
- .NET程序优化(GCServer )
- redis 学习笔记(4)-HA高可用方案Sentinel配置
- oracle: job使用
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- Linux模拟网络丢包与延迟的方法
- centos6.5通过yum安装nginx
- Linux系统利用cp命令实现强制覆盖功能的方法
- leetcode队列之设计循环双端队列
- Centos7学习之添加用户和用户组的方法
- Linux静态库与动态库实例详解
- Linux字符终端如何用鼠标移动一个红色矩形详解
- Linux中出现“No space left on device”错误的排查与解决方法
- 浏览器是如何调度进程和线程的?
- Linux shell利用sed如何批量更改文件名详解
- linux下通过xinetd服务管理 rsync 实现开机自启动
- linux实现自动删除最旧的几个文件详解
- 基于Ubuntu 16.04设置固定IP的方法教程
- CentOS添加和删除用户以及用户组的方法
- Linux下使用inode删除指定文件方法示例