HDU_4920_Matrix multiplication_杭电多校1010_矩阵乘法的代码风格
2014-08-06 10:55
375 查看
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4920
http://vjudge.net/contest/view.action?cid=52719#problem/J
1.题意:
两个n*n(1<=n<=800)的矩阵A、B,求AxB。
2.题解
(1)取模运算耗时很长,要放到最后计算。
(2)循环的写法要写成:for i for j if(aij!=0) for k 更新cik+=aij*bjk这样的形式。
(3)和队友yy了一种算法,首先预测aij全都是1,所以把b的每一列之和sumj填入cij数组中,这样的复杂度就是n^2。然后扫描aij,如果aij为0,就说明cik中的预测值要比实际值大,就用cik-=bkj(j=1~n),如果aij的值为2,就说明cik中的预测值比实际值小,就用cik+=bkj(j=1~n),这样就可以避免乘法。
暴力过code:
yy过code:
http://vjudge.net/contest/view.action?cid=52719#problem/J
1.题意:
两个n*n(1<=n<=800)的矩阵A、B,求AxB。
2.题解
(1)取模运算耗时很长,要放到最后计算。
(2)循环的写法要写成:for i for j if(aij!=0) for k 更新cik+=aij*bjk这样的形式。
(3)和队友yy了一种算法,首先预测aij全都是1,所以把b的每一列之和sumj填入cij数组中,这样的复杂度就是n^2。然后扫描aij,如果aij为0,就说明cik中的预测值要比实际值大,就用cik-=bkj(j=1~n),如果aij的值为2,就说明cik中的预测值比实际值小,就用cik+=bkj(j=1~n),这样就可以避免乘法。
暴力过code:
#include <iostream> #include <cstdio> #include <cstring> using namespace std; const int MAXN=888; int a[MAXN][MAXN],b[MAXN][MAXN],c[MAXN][MAXN]; int main() { // freopen("data.in","r",stdin); // freopen("data.out","w",stdout); int n; while(scanf("%d",&n)!=EOF){ for(int i=0;i<n;i++) for(int j=0;j<n;j++){ scanf("%d",&a[i][j]); a[i][j]%=3; } for(int i=0;i<n;i++) for(int j=0;j<n;j++){ scanf("%d",&b[i][j]); b[i][j]%=3; } for(int i=0;i<n;i++) for(int j=0;j<n;j++) c[i][j]=0; for(int i=0;i<n;i++) for(int j=0;j<n;j++){ if(a[i][j]==0)continue; for(int k=0;k<n;k++) c[i][k]=c[i][k]+a[i][j]*b[j][k]; } for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ if(j!=0) printf(" "); printf("%d",c[i][j]%3); } puts(""); } } return 0; }
yy过code:
#include <iostream> #include <cstdio> #include <cstring> using namespace std; //typedef const int MAXN=888; typedef unsigned long long ULL; int a[MAXN][MAXN],b[MAXN][MAXN],c[MAXN][MAXN],sum[MAXN]; int main() { // freopen("data.in","r",stdin); // freopen("data.out","w",stdout); int n; while(scanf("%d",&n)!=EOF){ for(int i=0;i<n;i++) for(int j=0;j<n;j++){ scanf("%d",&a[i][j]); a[i][j]%=3; } memset(sum,0,sizeof sum); for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ scanf("%d",&b[i][j]); b[i][j]%=3; sum[j]+=b[i][j]; } } for(int i=0;i<n;i++) sum[i]%=3; for(int i=0;i<n;i++) for(int j=0;j<n;j++) c[i][j]=0; for(int i=0;i<n;i++) for(int j=0;j<n;j++){ c[i][j]=sum[j]; } for(int i=0;i<n;i++) for(int j=0;j<n;j++) if(a[i][j]==0){ for(int k=0;k<n;k++) c[i][k]-=b[j][k]; }else if(a[i][j]==2){ for(int k=0;k<n;k++) c[i][k]+=b[j][k]; } for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ if(j!=0) printf(" "); printf("%d",((c[i][j])%3+3)%3); } puts(""); } } return 0; }
相关文章推荐
- 2014多校第五场1010 || HDU 4920 Matrix multiplication(矩阵乘法优化)
- hdu 4920 Matrix multiplication(矩阵乘法)
- hdu 4920暑假多校5的1010 Matrix multiplication 矩阵
- 矩阵乘法 --- hdu 4920 : Matrix multiplication
- HDU 4920 坑爹的矩阵乘法
- HDU 4920 Matrix multiplication(稀疏矩阵的乘法)不优化很容易TLE
- HDU 4920(杭电多校训练#5 1010 题) Matrix multiplication(不知道该挂个什么帽子。。。)
- HDU 4920 Matrix multiplication(矩阵乘法小技巧)——2014 Multi-University Training Contest 5
- HDU 4920 Matrix multiplication(矩阵乘法)
- Hdu 4920 Matrix multiplication[矩阵乘法的时间复杂度优化]
- hdu 4920 Matrix multiplication(矩阵乘法)2014多培训学校5现场
- HDU 4920 Matrix multiplication(矩阵乘法,cin,cout超时呀)
- hdu 4920 Matrix multiplication (矩阵乘法)
- hdu--4920--原来一直写了速度慢的矩阵乘法
- hdu 4920 Matrix multiplication 多校第五场 稀疏矩阵乘法
- HDU 4920(Matrix multiplication-矩阵乘法优化)
- HDU 4920 Matrix multiplication(矩阵乘法)
- HDU 1082 矩阵乘法次数计算 写了半天搞定,很有成就感啊
- 【矩阵乘法实践】HDU 1575——Tr A
- hdu 1588 Gauss Fibonacci(矩阵乘法,二分)