WUSTOJ 1232: 矩阵乘法(C)
1232: 矩阵乘法
Time Limit: 1 Sec Memory Limit: 128 MB 64bit IO Format: %lld
Description
小明明正在学习线性代数,老师布置了好多矩阵相乘的题目,手算累死了。最近刚学习了c语言,小明明灵机一动,决定让小伙伴帮他编一个程序用来计算矩阵相乘。亲爱的小伙伴们,不要让小明明童鞋失望哦哟嚯嚯嚯嚯v
Input
首先输入一个n,然后n个数据块。
每个数据块先给3个正整数N,M,P(N,M,P均不大于13)。
然后给一个NM的矩阵A,再接着给一个MP的矩阵B。
矩阵A和B的元素都是一个不大于37的非负整数。
详情请见sample。
Output
输出n个矩阵,分别表示A*B的结果。注意以%6d的格式输出,行末不能有空格。每个矩阵输出后再输出一个空行,具体见sample
Sample Input
2
2 3 2
1 2 3
4 5 6
1 2
1 2
1 2
1 2 3
1 2
1 2 3
4 5 6
Sample Output
6 12
15 30
9 12 15
HINT
注意格式:每个数字宽度为6,即%6d的格式输出每个数字,每行末尾不能有多余的空格。
每个矩阵输出后再输出一个空行以示分隔。看清样例。
代码
原理比较简单,不单独分析,看注释就可以。
网页评判是404ms,但我看到了不到200ms的大神,如果愿意分享代码,在此感谢。可以评论附上关键代码或者邮箱1。
// wustoj ac 404ms
#include<stdio.h>
#include<string.h>
#define MAX 14 // 矩阵最大为13*13的
int main()
{
// A[N][M], B[M][P]
int n, N, M, P, i, j, k;
int matrixA[MAX][MAX], matrixB[MAX][MAX], matrixAns[MAX][MAX];
scanf("%d", &n); // 输入数据组数
while (n--)
{
scanf("%d%d%d", &N, &M, &P);
// 输入A矩阵
for (i = 0; i < N; i++)
{
for (j = 0; j < M; j++)
{
scanf("%d", &matrixA[i][j]);
}
}
// 输入B矩阵
for (i = 0; i < M; i++)
{
for (j = 0; j < P; j++)
{
scanf("%d", &matrixB[i][j]);
}
}
// 重置为0
memset(matrixAns, 0, sizeof(matrixAns));
// 计算乘法
for (i = 0; i < P; i++)
{
for (j = 0; j < M; j++)
{
for (k = 0; k < N; k++)
{
matrixAns[k][i] += matrixA[k][j] * matrixB[j][i];
}
}
}
// 输出结果
for (i = 0; i < N; i++)
{
for (j = 0; j < P; j++)
{
printf("%6d", matrixAns[i][j]);
}
printf("\n");
}
// 每组数据空一行
printf("\n");
}
return 0;
}
Java版(提交超时,暂未解决,希望有人指点)
代码在Eclipse运行结果正确,OJ提交一直超时(Time Limit Exceed)
import java.util.Scanner;
public class Main {
public Main() {
Scanner sc = new Scanner(System.in);
int n;
// 输入数据组数
n = sc.nextInt();
while(n > 0) {
// 输入矩阵大小
int N, M, P;
N = sc.nextInt();
M = sc.nextInt();
P = sc.nextInt();
// 申请空间
int[][] matrixA = new int[N][M];
int[][] matrixB = new int[M][P];
// 输入矩阵A
for(int i = 0; i < N; i++) {
for(int j = 0; j < M; j++) {
matrixA[i][j] = sc.nextInt();
}
}
// 输入矩阵B
for(int i = 0; i < M; i++) {
for(int j = 0; j < P; j++) {
matrixB[i][j] = sc.nextInt();
}
}
// 申请保存结果的数组
int[][] matrixAns = new int[N][P];
// 计算乘法
for(int i = 0; i < N; i++) {
for(int j = 0; j < P; j++) {
for(int k = 0; k < M; k++) {
matrixAns[i][j] += matrixA[i][k] * matrixB[k][j];
}
}
}
// 输出结果
for(int i = 0; i < N; i++) {
for(int j = 0; j < P; j++) {
System.out.printf("%6d", matrixAns[i][j]);
}
System.out.println();
}
// 每组数据空一行
System.out.println();
n--;
}
sc.close();
}
public static void main(String[] args) {
new Main();
}
}
【End】
QQ邮箱:929450073@qq.com ↩︎
WUSTOJ 1232: 矩阵乘法(C)的更多相关文章
- *HDU2254 矩阵乘法
奥运 Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submissi ...
- *HDU 1757 矩阵乘法
A Simple Math Problem Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Ot ...
- CH Round #30 摆花[矩阵乘法]
摆花 CH Round #30 - 清明欢乐赛 背景及描述 艺术馆门前将摆出许多花,一共有n个位置排成一排,每个位置可以摆花也可以不摆花.有些花如果摆在相邻的位置(隔着一个空的位置不算相邻),就不好看 ...
- POJ3070 Fibonacci[矩阵乘法]
Fibonacci Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 13677 Accepted: 9697 Descri ...
- bzoj 2738 矩阵乘法
其实这题跟矩阵乘法没有任何卵关系,直接整体二分,用二维树状数组维护(刚刚学会>_<),复杂度好像有点爆炸(好像有十几亿不知道是不是算错了),但我们不能怂啊23333. #include&l ...
- 【BZOJ-2476】战场的数目 矩阵乘法 + 递推
2476: 战场的数目 Time Limit: 1 Sec Memory Limit: 128 MBSubmit: 58 Solved: 38[Submit][Status][Discuss] D ...
- 【BZOJ-1898】Swamp 沼泽鳄鱼 矩阵乘法
1898: [Zjoi2005]Swamp 沼泽鳄鱼 Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 1012 Solved: 566[Submit][S ...
- 【Codeforces718C】Sasha and Array 线段树 + 矩阵乘法
C. Sasha and Array time limit per test:5 seconds memory limit per test:256 megabytes input:standard ...
- 矩阵乘法的MapReduce实现
对于任意矩阵M和N,若矩阵M的列数等于矩阵N的行数,则记M和N的乘积为P=M*N,其中mik 记做矩阵M的第i行和第k列,nkj记做矩阵N的第k行和第j列,则矩阵P中,第i行第j列的元素可表示为公式( ...
随机推荐
- WPF中打开网页的两种方法
1.浏览器打开 Process proc = new System.Diagnostics.Process(); proc.StartInfo.FileName = "http://www. ...
- meshing-八分之一圆球
原视频下载地址:https://yunpan.cn/cqwiFDCg6PbFj 访问密码 d1c8
- Vue学习手记01-安装和项目创建
1.安装Vue 注:node版本必须大于等于8.9 vue-cli3.x:npm install -g @vue/cli vue-cli2.x:npm install -g @vue/cli-i ...
- STM32F429中LTDC的DMA2D加速
液晶屏的时序问题?每个液晶屏的时序都不一样,但总体上是类似的.如下图: VDEN: 数据使能信号.HSYNC: 每一行扫描的起始点, 在扫描过程中, 不会管上一行扫描有没有结束, 当出现 ...
- OpenJudge计算概论-最长单词2
/*======================================================================== 最长单词2 总时间限制: 1000ms 内存限制: ...
- 猎豹网校C++ Primer学习笔记2
14.数组 数组定义时的长度必须是在编译时就能确定的值. 全局数组会自动初始化为0. size_t 15.指针 其指向类型要相同. 指针和引用: 指针可以先不初始化,可以修改指向.有指针的指针. 16 ...
- VGG网络-ILSVRC-2014亚军
用于大尺度图片识别的非常深的卷积网络 使用一个带有非常小的(3*3)的卷积核的结构去加深深度,该论文的一个十分重要的改进就是它将卷机网络的深度增加到了16-19层,且可以用于比较大224*224的图片 ...
- 40 Flutter仿京东商城项目签名验证原理、签名验证算法
加群452892873 下载对应40课文件,运行方法,建好项目,直接替换lib目录 pubspec.yaml crypto: ^ SignServices.dart import 'dart:conv ...
- Redis键通知机制
Redis键通知机制 一.概念 自从redis2.8.0以后出了一个新特性,Keyspace Notifications 称为“键空间通知”. 这个特性大概是,凡是实现了Redis的Pub/Sub的客 ...
- Treeview控件如何获得子节点的所有父节点的名称
Delphi或c++ 的treeview控件,比如一个节点上面有个父节点,这个父节点上面还有一个父节点,如何获得这两个父节点的名字呢?请给出实现代码 先定义一个nodevarnode:TTreeNod ...