[CODEVS 3147] 矩阵乘法 2
来源:互联网 发布:百度文库离线文件json 编辑:程序博客网 时间:2024/05/21 14:59
描述
给出两个n*n的矩阵,m次询问它们的积中给定子矩阵的数值和。
http://codevs.cn/problem/3147/
分析
直接 n3 的矩阵乘法肯定超时, 要采用前缀和优化
row[s1] … row[t1]
col[s2] … col[t2](s1, s2) – (t1, t2)
row[x] * col[y] 表示用第 x 行的所有元素去乘第 y 行的所有元素.
==>
= row[s1] * col[s2] + row[s1] * col[s2+1] + … + row[s1] * col[t2] + row[s1+1] * col[s2] + … + row[s1+1] * col[t2] + … + row[t1] * col[t2]
// 分配律, 其实并没有看上去那么显然= row[s1] * (col[s2] + col[s2+1] + … + col[t2]) + row[s1+1] * (col[s2] + col[s2+1] + … + col[t2]) + … + row[t1] * (col[s2] + col[s2+1] + … + col[t2])
= (row[s1] + row[s1+1] + … + row[t1]) * (col[s2] + col[s2+1] + … + col[t2])
==> 用前缀和处理
注意: row 其实相当于一个 1 行 n 列的矩阵, 而 col 相当于一个 n 行 1 列的矩阵. 上式中 row[s1] + … + row[t1] 表示把 t1-s1+1 个这样的矩阵每个元素对应的加起来 (用前缀和优化) 得到一个新的 n 行 1 列的矩阵; col[s2] + … + col[t2] 采用同样方法得到一个新的 1 行 n 列的矩阵. 两个新矩阵再相乘就得到一个只有一个元素的矩阵了, 该元素就是最终答案.
代码
11809ms 49MB
#include<cstdio>#include<algorithm>using namespace std;const int maxn = 2000 + 10;typedef int Matrix[maxn][maxn];typedef long long LL;Matrix A, B;int main() { int n, m; scanf("%d %d", &n, &m); for(int x = 1; x <= n; x++) for(int y = 1; y <= n; y++) { scanf("%d", &A[x][y]); A[x][y] += A[x-1][y]; } // 前 x 行元素和 for(int x = 1; x <= n; x++) for(int y = 1; y <= n; y++) { scanf("%d", &B[x][y]); B[x][y] += B[x][y-1]; } // 前 y 列元素和 for(int i = 0; i < m; i++) { int x1, y1, x2, y2; scanf("%d %d %d %d", &x1, &y1, &x2, &y2); if(x1 > x2) swap(x1, x2); if(y1 > y2) swap(y1, y2); LL ans = 0; for(int i = 1; i <= n; i++) ans += (LL)(A[x2][i] - A[x1-1][i]) * (B[i][y2] - B[i][y1-1]); printf("%lld\n", ans); } return 0;}
主页
http://blog.csdn.net/qq_21110267
- [CODEVS 3147] 矩阵乘法 2
- codevs 3147 3147 矩阵乘法 2
- codevs#3147[题解]矩阵乘法2
- 【codevs 3147】矩阵乘法 2 2012年(乱搞)
- 【codevs 1287】 矩阵乘法
- Codevs 1287 矩阵乘法
- Codevs 3332 数列(矩阵乘法)
- 【矩阵乘法】【codevs 1250】Fibonacci数列
- codevs#1287[题解] Key:矩阵乘法
- Codevs 1070 普通递归关系(矩阵乘法)
- Codevs 1482 路线统计(矩阵乘法)
- Codevs 1574 广义斐波那契数列(矩阵乘法)
- 矩阵乘法-2
- 矩阵乘法(2)
- 矩阵乘法2(codevs3147)
- 矩阵乘法(2)【模板】
- codevs天梯高精度乘法
- codevs 高精度乘法 3117
- 抓取网络json数据并存入mongodb(2)
- Objective-C的对象等同性
- Axure Mobile原型发布工具说明
- No CPU/ABI system image available for this barget
- unity之台湾游戏教父王俊博:做永远的抬轿人
- [CODEVS 3147] 矩阵乘法 2
- 关于项目
- NYOJ--60 谁获得了最高奖学金【水题】
- FFmpeg示例程序合集-Git批量获取脚本
- unity游戏之年度最佳游戏奖项名单公布
- 便于理解的基本算法(常用的5种排序+打乱数组+两个范围内的一个随机数+1到100的不重复的随机数)
- 基于嵌入式Linux下GCC编译器__sync_sub_and_fetch_4错误的解决心得
- TypeScript 类(Classes)
- 土豪手游:[糖果传奇]累计收入破187亿元