14
голосов
Принятый
Помогите понять код: алгоритм генерации матрицы со спиралью
# Создаётся функция, которая и будет всё делать
def zm(n):
# Переменной dx присваивается значение 1
# Переменной dy присваивается значение 0
# Обычно dx и dy - это некие приращения для ...
14
голосов
Какой алгоритм использовать для решения системы с приведенной матрицей коэффициентов?
Как правильно замечено MBo, это циклическая трехдиагональная система линейных уравнений.
Есть несколько вариантов решения:
Использование метода прогонки (также, называется алгоритмом Томаса в ...
12
голосов
Как повернуть матрицу (двумерный массив) на 90 градусов с помощью zip?
Нашел довольно элегантный способ сделать поворот матрицы в одну строчку без numpy и циклов. В рунете ничего толкового не смог найти, может кому то поможет.
Оригинал здесь:
Линк на оригинал
rotated = ...
Общий
12
голосов
Принятый
Какой алгоритм использовать для решения системы с приведенной матрицей коэффициентов?
Это циклическая трёхдиагональная система ЛУ.
Её можно решить, используя прогонку для чисто трёхдиагональной матрицы и формулу Шермана-Моррисона для модификации решения.
Код можно найти в Numerical ...
11
голосов
Принятый
Магический квадрат
Нет нужды сравнивать все строки со всеми, чтобы выяснить, что все суммы равны. Достаточно найти первую сумму и сравнивать остальные с ней.
first = sum(matrix[0])
for k in range(1, n):
if sum(...
10
голосов
Задача на сумму матриц
Сначала немного матричной математики...
Откуда совсем просто получается (I - единичная матрица)
и
так что
Ну, а дальше использовать быстрое возведение в степень, которое требует log(n) умножений. ...
9
голосов
Магический квадрат
Я бы, наверное, сделал так:
import numpy as np
def ismagic(a):
if np.array_equal(np.unique(a.sum(axis=1)),np.unique(a.sum(axis=0))):
return True
else:
return False
Проверяем:...
8
голосов
Принятый
Найти матрицу-множитель в произведении матриц при известном результате
Воспользуйтесь обратной (inverse) матрицей:
In [224]: b = np.linalg.inv(a) * c
In [225]: b
Out[225]:
matrix([[-3. ],
[ 2.5]])
Это будет работать для объектов типа numpy.matrix.
Если a и c - ...
8
голосов
Медленное умножение матриц 1024х1024
В комментарий не влезет, но это не ответ. Это именно комментарий.
Написал для проверки (см. ниже). Обалдел, ибо у меня, скомпилированное VC++ 2017, таки дало:
1023: 2106433
1024: 7664347
1025: ...
8
голосов
Принятый
Как эффективней вычислить определитель матрицы?
Как это можно оптимизировать?
воспользоваться модулем numpy:
In [32]: import numpy as np # pip install numpy
In [33]: a = np.random.rand(20, 20)
In [34]: res = np.linalg.det(a)
In [35]: res
Out[...
7
голосов
Деструктор удаляет дважды один и тот же адрес, из-за чего происходит сбой программы при перемножении матриц
В этом предложении
a = input1*input2;
используется созданный компилятором по умолчанию копирующий оператор присваивания, который результат перемножения двух матриц - временный объект этого же класса ...
6
голосов
Как правильно умножить матрицу на другую матрицу в нескольких потоках?
Вот накидал рабочий код:
public static class CalcThread extends Thread {
private int startRow, endRow;
private int[][] a, b, result;
private int n;
public CalcThread(int[][] a, int[][]...
6
голосов
Принятый
Как повернуть матрицу (двумерный массив) на 90 градусов с помощью zip?
Вполне себе рабочий вариант, если скорость не критична.
Но если есть требования по скорости обработки, то лучше все-таки воспользоваться NumPy.
Сравнение производительности для массива 100x100:
In [...
6
голосов
Принятый
Количество "особенных"елементов массива. С++
Чтобы не делать за вас домашку поясню на словах, без кода.
Просто подумайте логически: если элемент больше чем сумма других в колонке, то он точно наибольший из всех. Следовательно нужно искать ...
6
голосов
Принятый
Вывод матрицы на консоль в виде пирамиды
Задачка довольно тривиальная и немного математическая:
public static int[][] pyramideMatrix(int cells) {
int dim = cells * cells;
int[][] arr = new int[dim][dim];
for (int i = 0; i <...
6
голосов
Принятый
6
голосов
Создание цикла , складывающего два соседних числа из предыдущего элемента матрицы
pascal_triangle - бесконечный генератор строк треугольника Паскаля. Следующая строка получается суммированием соседних элементов предыдущей. Для этого к предыдущей строке прибавим нули слева и справа:
...
6
голосов
Принятый
Как найти максимальный и минимальный элемент из диапазона, не используя цикл?
Если речь идет о нахождении минимума/максимума в столбцах со второго по пятый:
a = np.array([[1,2,3,12,2,8], [15,25,32,75,57,35] , [18,22,32,53,57,46] , [13,23,35,53,55,83]])
min_, max_ = a[:, 2:5]....
6
голосов
Принятый
Как диагональ заполнить нулями, элементы выше главной диагонали заполнить -1 , ниже – 3?
Воспользуйтесь np.diag_indices(), np.triu_indices() и tril_indices():
N = a.shape[0]
a[np.diag_indices(N)] = 0
a[np.triu_indices(N, 1)] = -1
a[np.tril_indices(N, -1)] = -3
результат:
In [171]: a
...
5
голосов
Проблема вычисления обратной матрицы при типе float
Ну собственно что и следовало ожидать. Тип данных float занимает (обычно) 4 байта, а double — 8, есть ещё long double, он 10 байт. Соответственно чем больше места на тип, тем больше он хранит ...
5
голосов
Принятый
Как проще всего считать матрицу(f64) из файла на Rust?
Перевод моего кода с Си на Rust
use std::fs::File;
use std::io::{BufRead, BufReader};
fn main() {
// open the file
let mut f = BufReader::new(File::open("input.txt").unwrap());
// read ...
5
голосов
Магический квадрат
Решение с использованием модуля Numpy для проверки на настоящий магический квадрат (с проверкой сумм главной и побочной диагоналей):
import numpy as np
a = np.array([[2,7,6],[9,5,1],[4,3,8]])
In [...
5
голосов
Принятый
Матрица в C. Как это можно реализовать?
Мне все это видеться примерно так, если я конечно правильно понял вопрос
#include <stdio.h>
#include <math.h>
void somedo1(void *_x)
{
int *x = _x;
*x = *x * *x;
}
void somedo2(...
5
голосов
Принятый
Как заменить каждый элемент матрицы на среднее арифметическое?
Чтобы приведённый код хоть как-то заработал, стоит для начала вынести инициализацию суммы и количества элементов вне внутренних циклов
summator = arrayA[i][j];
num = 0;
for(int k = 0;k<formul;k++)...
5
голосов
Принятый
Как красиво прочитать двумерную матрицу?
мой вариант https://ideone.com/xwsY1A
import sys
from itertools import takewhile
def read_matrix():
return [[int(x) for x in l.split()] for l in takewhile(str.strip, sys.stdin)]
a = read_matrix(...
5
голосов
Принятый
Инициализация большого вектора в C++ из другой программы
Примерно так:
ostream& write(ostream& out, const vector<vector<double>>& d)
{
size_t sz = d.size();
out.write((char*)&sz,sizeof(sz));
for(const auto& v:d)
...
5
голосов
Задача Олимпиадная(Штурм)
Это решение находит только клетки с максимальным по карте количеством, а нужны локальные максимумы. Для их нахождения достаточно для каждой клетки проверить четырёх (если они есть) соседей.
if (a[i][j]...
5
голосов
Принятый
Задача Олимпиадная(Штурм)
Вы смотрите, сколько ячеек с максимальным (по всей матрице) значением, в то время как требуется только локальный максимум - больше, чем в соседних! Что-то типа
#include <vector>
#include <...
5
голосов
Ошибка транспонирования матрицы (замена столбцов на строки)
Идиоматический вариант транспонирования матрицы на "Vanilla Python":
arr_transposed = list(zip(*arr))
результат:
In [15]: arr_transposed
Out[15]:
[(2, 0, 0, 0, 6, 7, 0, 0, 1, 0),
(0, 1, 0, ...
Only top scored, non community-wiki answers of a minimum length are eligible
Связанные метки
матрицы × 702c++ × 205
python × 186
массивы × 185
c × 79
java × 62
алгоритм × 61
c# × 49
python-3.x × 49
математика × 49
numpy × 49
циклы × 26
list × 25
javascript × 23
сортировка × 22
динамические-массивы × 19
указатели × 17
графы × 14
память × 12
классы × 11
ооп × 10
функции × 10
файлы × 9
pascal × 9
оптимизация × 9