14
Как правильно замечено MBo, это циклическая трехдиагональная система линейных уравнений.
Есть несколько вариантов решения:
Использование метода прогонки (также, называется алгоритмом Томаса в англоязычной литературе) с последующим применением метода Шермана-Моррисона.
Так как ваша матрица симметрична, то можно работать с Cyclic reduction.
Преобразование к ...
ответ дан 2 ноя '20 в 3:29
Anton Menshov
50411 золотой знак88 серебряных знаков2020 бронзовых знаков
13
# Создаётся функция, которая и будет всё делать
def zm(n):
# Переменной dx присваивается значение 1
# Переменной dy присваивается значение 0
# Обычно dx и dy - это некие приращения для переменных x и y
dx, dy = 1, 0
# Переменным x и y присваивается значение 0
x, y = 0, 0
# Создаётся список списков
# Это матрица n*n
# Пока ...
12
Это циклическая трёхдиагональная система ЛУ.
Её можно решить, используя прогонку для чисто трёхдиагональной матрицы и формулу Шермана-Моррисона для модификации решения.
Код можно найти в Numerical Recipes in C, раздел 2.7.2
11
Нет нужды сравнивать все строки со всеми, чтобы выяснить, что все суммы равны. Достаточно найти первую сумму и сравнивать остальные с ней.
first = sum(matrix[0])
for k in range(1, n):
if sum(matrix[0]) != first:
return False
Для столбцов - просто посчитать их суммы и так же сравнить с эталоном
for k in range(0, n):
if sum([row[k] for row ...
10
Сначала немного матричной математики...
Откуда совсем просто получается (I - единичная матрица)
и
так что
Ну, а дальше использовать быстрое возведение в степень, которое требует log(n) умножений. Каждое умножение - не будем умничать и звать на помощь Штрассена :) - выполняется за q3. Вычитание - за q2, обратная матрица - q3. Таким образом, в результате ...
ответ дан 25 апр '20 в 5:26
Harry
164k1212 золотых знаков9292 серебряных знака194194 бронзовых знака
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
Проверяем:
a = np.matrix([[1, 2], [4, 3]])
print(ismagic(a))
b = np.ones((5,5))
print(ismagic(b))
c = np.matrix([[2,7,6],[9,5,1],[4,3,8]])
print(ismagic(c))
d = np....
ответ дан 27 ноя '18 в 18:57
strawdog
11.5k22 золотых знака1313 серебряных знаков2727 бронзовых знаков
8
Нашел довольно элегантный способ сделать поворот матрицы в одну строчку без numpy и циклов. В рунете ничего толкового не смог найти, может кому то поможет.
Оригинал здесь:
Линк на оригинал
rotated = zip(*original[::-1]) # Python 2
rotated = tuple(zip(*original[::-1])) # Python 3
Как это работает.
original = [[1, 2],
[3, 4]]
Сначала работает ...
8
Воспользуйтесь обратной (inverse) матрицей:
In [224]: b = np.linalg.inv(a) * c
In [225]: b
Out[225]:
matrix([[-3. ],
[ 2.5]])
Это будет работать для объектов типа numpy.matrix.
Если a и c - объекты типа numpy.ndarray, то нужно использовать dot product (как в ответе @MarianD):
In [8]: np.linalg.inv(a).dot(c)
Out[8]:
matrix([[-3. ],
[ 2.5]])...
8
В комментарий не влезет, но это не ответ. Это именно комментарий.
Написал для проверки (см. ниже). Обалдел, ибо у меня, скомпилированное VC++ 2017, таки дало:
1023: 2106433
1024: 7664347
1025: 2106884
При отключенной оптимизации эффект выражен меньше:
1023: 9592402
1024: 11608342
1025: 9480406
Это для 64-разрядного приложения. 32-разрядное, впрочем, ...
ответ дан 28 дек '18 в 19:25
Harry
164k1212 золотых знаков9292 серебряных знака194194 бронзовых знака
7
В этом предложении
a = input1*input2;
используется созданный компилятором по умолчанию копирующий оператор присваивания, который результат перемножения двух матриц - временный объект этого же класса - присваивает объекту a.
Определенный компилятором копирующий оператор присваивания просто по-членно копирует данные одного объекта в другой. В результате в ...
ответ дан 25 окт '16 в 22:24
Vlad from Moscow
43.7k33 золотых знака3232 серебряных знака8484 бронзовых знака
7
Как это можно оптимизировать?
воспользоваться модулем 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[35]: 0.09252260373277807
время работы для матрицы 20x20:
In [36]: %timeit np.linalg.det(a)
27.6 µs ± 37 ns per loop (mean ± std. dev. of 7 runs, 10000 loops ...
6
Чтобы не делать за вас домашку поясню на словах, без кода.
Просто подумайте логически: если элемент больше чем сумма других в колонке, то он точно наибольший из всех. Следовательно нужно искать наибольший, попутно подсчитывая сумму других элементов. Как искать наибольший элемент, надеюсь знаете... Ну ладно, вот так это делается:
auto maximum = array[0];
...
6
Задачка довольно тривиальная и немного математическая:
public static int[][] pyramideMatrix(int cells) {
int dim = cells * cells;
int[][] arr = new int[dim][dim];
for (int i = 0; i < arr.length; i++) {
for (int j = 0; j < arr[i].length; j++) {
int ii;
int ij;
if (i < dim/2) {
...
ответ дан 9 мая '18 в 15:19
Peter Samokhin
6,01822 золотых знака1414 серебряных знаков4646 бронзовых знаков
6
for (int j = 0; j < 6; i++)
6
pascal_triangle - бесконечный генератор строк треугольника Паскаля. Следующая строка получается суммированием соседних элементов предыдущей. Для этого к предыдущей строке прибавим нули слева и справа:
1 3 3 1 -> 0 1 3 3 1
1 3 3 1 0
+ ---------
1 4 6 4 1
itertools.islice - инструмент для обрезания бесконечных ...
ответ дан 3 янв в 8:51
Stanislav Volodarskiy
4,05411 золотой знак22 серебряных знака1717 бронзовых знаков
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].min(), a[:, 2:5].max()
Чтобы найти минимум/максимум элементов матрицы в нулевой и третьей строках и учитывая столбцы только со 2-го по 5-й :
In [37]: min_, max_ ...
5
Перевод моего кода с Си на 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 the first line and extract the number from it
let mut num_line = String::new();
f.read_line(&mut num_line).unwrap();
let n: usize = num_line[1..]...
5
Ну собственно что и следовало ожидать. Тип данных float занимает (обычно) 4 байта, а double — 8, есть ещё long double, он 10 байт. Соответственно чем больше места на тип, тем больше он хранит информации и тем точнее результат вычисления.
Если вы вычисляете число обусловленности, то понимаете, что это по сути показатель, насколько ошибки вычислений скажутся ...
5
Вот накидал рабочий код:
public static class CalcThread extends Thread {
private int startRow, endRow;
private int[][] a, b, result;
private int n;
public CalcThread(int[][] a, int[][] b, int[][] result, int startRow, int endRow) {
this.a = a;
this.b = b;
this.result = result;
...
5
Решение с использованием модуля Numpy для проверки на настоящий магический квадрат (с проверкой сумм главной и побочной диагоналей):
import numpy as np
a = np.array([[2,7,6],[9,5,1],[4,3,8]])
In [90]: a
Out[90]:
array([[2, 7, 6],
[9, 5, 1],
[4, 3, 8]])
In [91]: a.sum(axis=0)
Out[91]: array([15, 15, 15])
In [92]: a.sum(axis=1)
Out[92]: ...
5
Мне все это видеться примерно так, если я конечно правильно понял вопрос
#include <stdio.h>
#include <math.h>
void somedo1(void *_x)
{
int *x = _x;
*x = *x * *x;
}
void somedo2(void *_x)
{
double *x = _x;
*x = log(*x);
}
void func(void *matrix, size_t nmemb, size_t size, void (*f)(void *))
{
char *m = matrix;
for (...
5
Чтобы приведённый код хоть как-то заработал, стоит для начала вынести инициализацию суммы и количества элементов вне внутренних циклов
summator = arrayA[i][j];
num = 0;
for(int k = 0;k<formul;k++){
...
Потом можно обдумать логику проверок. а лучше без проверок обходить только то, что нужно:
summator = arrayA[i][j];
num = 0;
int top = i>=n?i-n:0;
...
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()
print(a)
b = read_matrix()
print(b)
5
Примерно так:
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)
{
sz = v.size();
out.write((char*)&sz,sizeof(sz));
out.write((char*)v.data(),sz*sizeof(double));
}
return out;
}
...
5
Это решение находит только клетки с максимальным по карте количеством, а нужны локальные максимумы. Для их нахождения достаточно для каждой клетки проверить четырёх (если они есть) соседей.
if (a[i][j] > a[i-1][j] && a[i][j] > a[i][j-1] &&
a[i][j] > a[i+1][j] && a[i][j] > a[i][j+1]) ...
Для исключения особых случаев ...
5
Вы смотрите, сколько ячеек с максимальным (по всей матрице) значением, в то время как требуется только локальный максимум - больше, чем в соседних! Что-то типа
#include <vector>
#include <string>
#include <iostream>
#include <iomanip>
using namespace std;
int ok(const vector<vector<int>>& v, int i, int j)
{
if (i-...
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, 0, 0, 0, 8, 0, 1, 0),
(0, 1, 0, 1, 0, 1, 0, 1, 0, 0),
(0, 1, 0, 1, 0, 1, 0, 0, 0, 1),
(0, 0, 0, 1, 0, 1, 0, 1, 0, 1),
(0, 0, 0, 0, 0, 1, 0, 1, 0, 0),
(1, ...
5
Потому что вы делите на 0.
Запишите
assert(a[k][i] != 0);
перед
float s = a[k + 1][i] / a[k][i];
и убедитесь в этом сами...
Что и неудивительно - ведь матрица
10 23 17 44
15 35 26 69
25 57 42 108
30 69 51 133
вырожденная, ее определитель равен 0...
P.S. Что-то странный у вас вложенный цикл
for (int i = 0; i < n-1; i++) {
for (int i = 0; i < n -...
ответ дан 14 сен '20 в 4:02
Harry
164k1212 золотых знаков9292 серебряных знака194194 бронзовых знака
5
"Сдвинуть указатель" можно, хотя это будет несколько через... вобщем, не самым прямым путем.
В С++ проще всего воспользоваться ссылкой:
int m[5][5] =
{
{ 1,2,3,4,5},
{ 6,7,8,9,10},
{ 11,12,13,14,15},
{ 16,17,18,19,20},
{ 21,22,23,24,25}
};
typedef int reduced_matrix[4][5];
reduced_matrix& r = *(reduced_matrix*)(m+1);
for(...
ответ дан 23 сен '20 в 3:26
Harry
164k1212 золотых знаков9292 серебряных знака194194 бронзовых знака
Допускаются только превышающие минимальную длину ответы с наивысшим рейтингом, не являющиеся общими
Связанные метки
матрицы × 544c++ × 183
массивы × 147
python × 124
c × 59
алгоритм × 49
java × 48
математика × 40
python-3.x × 33
c# × 29
numpy × 28
циклы × 21
javascript × 20
сортировка × 18
динамические-массивы × 18
указатели × 16
память × 12
list × 11
классы × 9
графы × 9
линейная-алгебра × 9
pascal × 8
opengl × 8
графика × 8
вывод × 8