C++前缀和算法:构造乘积矩阵

简介: C++前缀和算法:构造乘积矩阵

基础知识点

C++算法:前缀和基础

题目

给你一个下标从 0 开始、大小为 n * m 的二维整数矩阵 grid ,定义一个下标从 0 开始、大小为 n * m 的的二维矩阵 p。如果满足以下条件,则称 p 为 grid 的 乘积矩阵 :

对于每个元素 p[i][j] ,它的值等于除了 grid[i][j] 外所有元素的乘积。乘积对 12345 取余数。

返回 grid 的乘积矩阵。

示例 1:

输入:grid = [[1,2],[3,4]]

输出:[[24,12],[8,6]]

解释:p[0][0] = grid[0][1] * grid[1][0] * grid[1][1] = 2 * 3 * 4 = 24

p[0][1] = grid[0][0] * grid[1][0] * grid[1][1] = 1 * 3 * 4 = 12

p[1][0] = grid[0][0] * grid[0][1] * grid[1][1] = 1 * 2 * 4 = 8

p[1][1] = grid[0][0] * grid[0][1] * grid[1][0] = 1 * 2 * 3 = 6

所以答案是 [[24,12],[8,6]] 。

示例 2:

输入:grid = [[12345],[2],[1]]

输出:[[2],[0],[0]]

解释:p[0][0] = grid[0][1] * grid[0][2] = 2 * 1 = 2

p[0][1] = grid[0][0] * grid[0][2] = 12345 * 1 = 12345. 12345 % 12345 = 0 ,所以 p[0][1] = 0

p[0][2] = grid[0][0] * grid[0][1] = 12345 * 2 = 24690. 24690 % 12345 = 0 ,所以 p[0][2] = 0

所以答案是 [[2],[0],[0]] 。

感悟

原以为和MOD = 1000000007一样,直接使用封装好的此类,发现错误。赛场上时间紧急,来不及分析是两个不同的问题,还是我封装错误。只好使用笨办法。我记得1000000007是质数,才能转除为乘。考虑过12345是否是质数,当时觉判断烦恼,所以没判断。现在觉得很简单:以5结尾,就是5的倍数,不是质数。

分析

前缀和后缀和。第一轮的preRow记录[0,r) 所有元素的乘积,vLeft[r][c] 记录本行左边各元素的乘积,vRight[r][c]记录本行右边各元素的乘积。vRet[r][c]记录这三个的乘积。第二轮preRow记录[r+1,m_r)所有元素的乘积。第二轮vRet[r][c]乘以preRow就是结果。

测试用例

1 2 3
4 5 6
7 8 9

结果

当前数 前面行的乘积 前面行的乘积 左边乘积 右边乘积
1 1 4…9 1 6
2 1 4…9 1 3
3 1 4…9 2 1
4 6 7…9 1 30
5 6 7…9 4 6
6 6 7…9 20 1
7 1…6 1 1 72
8 1…6 1 7 9
9 1…6 1 56 1

解释

对{4,5,6}而言第一轮preRow是123=6,第二轮preRow是789。对4而言,left是1,right是30。对5而言,left是4,right是6。对6而言,left是20,right是1。

时间复杂度

O(n^2) 2轮,每轮2层循环,每层循环是O(n)。

代码

class Solution {
public:
vector<vector> constructProductMatrix(vector<vector>& grid) {
m_r = grid.size();
m_c = grid.front().size();
//vLeft记录当前行,左边的成绩
vector<vector> vLeft(m_r, vector(m_c)), vRight(m_r, vector(m_c)), vRet(m_r, vector(m_c));
int iPreRow = 1;
for (int r = 0; r < m_r; r++)
{
int pre = 1;
for (int c = 0; c < m_c; c++)
{
vLeft[r][c] = pre;
MulSelf(pre, grid[r][c]);
}
pre = 1;
for (int c = m_c-1 ; c >= 0 ; c-- )
{
vRight[r][c] = pre;
MulSelf(pre, grid[r][c]);
}
for (int c = 0; c < m_c; c++)
{
vRet[r][c] = 1;
MulSelf(vRet[r][c], iPreRow);
MulSelf(vRet[r][c], vLeft[r][c]);
MulSelf(vRet[r][c], vRight[r][c]);
}
MulSelf(iPreRow, pre);
}
iPreRow = 1;
for (int r = m_r-1; r >= 0 ; r-- )
{
int pre = 1;
for (int c = 0; c < m_c; c++)
{
MulSelf(vRet[r][c], iPreRow);
MulSelf(pre, grid[r][c]);
}
MulSelf(iPreRow, pre);
}
return vRet;
}
void MulSelf(int& self, int other)
{
const int MOD = 12345;
self = ((long long)self * other) % MOD;
}
int m_r, m_c;
};

一维化降低复杂度

分析

vLeft[r][c]记录 [0,r)行所有元素及r行[0,c)列元素的乘积,第二轮的pre记录(r,m_c)行所有元素及r行(c,m_c)列元素的乘积。

vLeft[1][1] = 1234 第二轮的pre = 9876

代码

class Solution {
public:
vector<vector> constructProductMatrix(vector<vector>& grid) {
m_r = grid.size();
m_c = grid.front().size();
vector < vector> vLeft(m_r, vector(m_c));
int pre = 1;
for (int r = 0; r < m_r; r++)
{
for (int c = 0; c < m_c; c++)
{
vLeft[r][c] = pre;
MulSelf(pre, grid[r][c]);
}
}
vector<vector> vRet(m_r, vector(m_c));
pre = 1;
for (int r = m_r-1 ; r >= 0 ;r–)
{
for (int c = m_c-1 ; c >= 0 ; c-- )
{
const int index = m_c * r + c;
vRet[r][c] = pre;
MulSelf(vRet[r][c], vLeft[r][c]);
MulSelf(pre, grid[r][c]);
}
}
return vRet;
}
void MulSelf(int& self, int other)
{
const int MOD = 12345;
self = ((long long)self * other) % MOD;
}
int m_r, m_c;
};

测试用例

template
void Assert(const vector& v1, const vector& v2)
{
if (v1.size() != v2.size())
{
assert(false);
return;
}
for (int i = 0; i < v1.size(); i++)
{
assert(v1[i] == v2[i]);
}
}
template
void Assert(const T& t1, const T& t2)
{
assert(t1 == t2);
}
int main()
{
vector<vector>grid = { {1,2,3},{4,5,6} };
vector<vector> ans = { {720,360,240},{180,144,120} };
auto res = Solution().constructProductMatrix(grid);
Assert(res, ans);
grid = { {1,2,},{3,4},{5,6 }};
ans = { {720,360},{240,180},{144,120} };
res = Solution().constructProductMatrix(grid);
Assert(res, ans);
grid = { { 1,2,3,4,5,6 } };
ans = { { 720,360,240,180,144,120} };
res = Solution().constructProductMatrix(grid);
Assert(res, ans);
grid = { { 1},{2},{3},{4},{5},{6} };
ans = { { 720},{360},{240},{180},{144},{120} };
res = Solution().constructProductMatrix(grid);
Assert(res, ans);
CConsole::Out(res);
}

其它

视频课程

要是你认为本篇难道较大,不好入手,推荐你先学习基础算法的课程,我已完成部分,余下部分持续更新中,就在CSDN学院。

https://edu.csdn.net/course/detail/38771

C#入职培训、C++入职培训等课程

https://edu.csdn.net/lecturer/6176

测试环境

操作系统:win7 开发环境: VS2019 C++17

相关下载

如果你想观其大略,建设下载《闻缺陷则喜算法册》doc版

https://download.csdn.net/download/he_zhidan/88348653


相关文章
|
5天前
|
存储 算法 安全
基于红黑树的局域网上网行为控制C++ 算法解析
在当今网络环境中,局域网上网行为控制对企业和学校至关重要。本文探讨了一种基于红黑树数据结构的高效算法,用于管理用户的上网行为,如IP地址、上网时长、访问网站类别和流量使用情况。通过红黑树的自平衡特性,确保了高效的查找、插入和删除操作。文中提供了C++代码示例,展示了如何实现该算法,并强调其在网络管理中的应用价值。
|
2天前
|
存储 算法 安全
基于哈希表的文件共享平台 C++ 算法实现与分析
在数字化时代,文件共享平台不可或缺。本文探讨哈希表在文件共享中的应用,包括原理、优势及C++实现。哈希表通过键值对快速访问文件元数据(如文件名、大小、位置等),查找时间复杂度为O(1),显著提升查找速度和用户体验。代码示例展示了文件上传和搜索功能,实际应用中需解决哈希冲突、动态扩容和线程安全等问题,以优化性能。
|
10天前
|
算法
|
10天前
|
算法 安全 C++
用 C++ 算法控制员工上网的软件,关键逻辑是啥?来深度解读下
在企业信息化管理中,控制员工上网的软件成为保障网络秩序与提升办公效率的关键工具。该软件基于C++语言,融合红黑树、令牌桶和滑动窗口等算法,实现网址精准过滤、流量均衡分配及异常连接监测。通过高效的数据结构与算法设计,确保企业网络资源优化配置与安全防护升级,同时尊重员工权益,助力企业数字化发展。
32 4
|
3月前
|
并行计算 算法 IDE
【灵码助力Cuda算法分析】分析共享内存的矩阵乘法优化
本文介绍了如何利用通义灵码在Visual Studio 2022中对基于CUDA的共享内存矩阵乘法优化代码进行深入分析。文章从整体程序结构入手,逐步深入到线程调度、矩阵分块、循环展开等关键细节,最后通过带入具体值的方式进一步解析复杂循环逻辑,展示了通义灵码在辅助理解和优化CUDA编程中的强大功能。
|
3月前
|
机器学习/深度学习 算法 搜索推荐
django调用矩阵分解推荐算法模型做推荐系统
django调用矩阵分解推荐算法模型做推荐系统
48 4
|
3月前
|
存储 算法 C++
高精度算法(加、减、乘、除,使用c++实现)
高精度算法(加、减、乘、除,使用c++实现)
779 0
高精度算法(加、减、乘、除,使用c++实现)
|
3月前
|
算法 数据处理 C++
c++ STL划分算法;partition()、partition_copy()、stable_partition()、partition_point()详解
这些算法是C++ STL中处理和组织数据的强大工具,能够高效地实现复杂的数据处理逻辑。理解它们的差异和应用场景,将有助于编写更加高效和清晰的C++代码。
54 0
|
3月前
|
存储 算法
动态规划算法学习一:DP的重要知识点、矩阵连乘算法
这篇文章是关于动态规划算法中矩阵连乘问题的详解,包括问题描述、最优子结构、重叠子问题、递归方法、备忘录方法和动态规划算法设计的步骤。
194 0
|
3月前
|
存储 算法 决策智能
【算法】博弈论(C/C++)
【算法】博弈论(C/C++)