Возведение большой матрицы, читаемой с файла, в степень N, и запись обратно в файл на C#

@МВлад · Регистрация: 07.01.2015

Author24 — интернет-сервис помощи студентам

Доброго времени. Очень нужно написать программу на C#, которая прочитает с файла огромную матрицу с вещественными числами (допустим, 1000 на 1000), возвести её в степень N, и результат записать в новый файл. Вот так... Думал, сам напишу, убил кучу времени, и теперь, от безысходности, пишу тут. Помогите пожалуйста. За ранее благодарен всем, кто уделит время.
И с Рождеством всех, ребята. =)
P.S.: это нужно делать в консольном приложении. Очень срочно.

Добавлено через 1 минуту
Пример матрицы тут можно взять (363х363):
http://rghost.ru/60192711

@rattrapper · 07.01.2015, 22:33

МВлад, алгоритм мудреный, но вроде работает

C#

static int[,] MatrixPower(int[,] matrix, int power)
{
    int l = matrix.GetLength(0);
    if (matrix.GetLength(1) != l)
        throw new ArgumentException("matrix");
    int[,] result = new int[l, l];
    for(int i = 0; i < l; i++)
        for(int j = 0; j < l; j++)
        {
            int sum = 0;
            var ind = new int[power - 1];
            while(ind[0] != l)
            {
                int el = matrix[i, ind[0]];
                for (int k = 1; k < power - 1; k++)
                    el *= matrix[ind[k - 1], ind[k]];
                el *= matrix[ind[power - 2], j];
                sum += el;
                int m = power - 2;
                while(++ind[m]==l && m > 0)
                {
                    ind[m--] = 0;
                }
            }
            result[i, j] = sum;
        }
    return result;
}

@chumich · 07.01.2015, 22:35

Посмотрите здесь, может быть подойдет:
http://easylab.net.ua/massivyi... tritsami-1
или здесь:
Возведение булевых матриц в степень

@rattrapper · 07.01.2015, 22:57

МВлад, вот еще вариант, который проще понять, рекурсивный

C#

static int[,] MatrixPower2(int[,] matrix, uint power)
{
    int l = matrix.GetLength(0);
    if (matrix.GetLength(1) != l)
        throw new ArgumentException();
    var result = new int[l, l];
    if(power == 0)
    {
        for (int i = 0; i < l; i++)
            result[i, i] = 1;
        return result;
    }
    else if (power == 1)
    {
        return matrix;
    }
    result = MatrixPower2(Multiply(matrix, matrix), power / 2);
    if (power % 2 != 0)
        result = Multiply(result, matrix);
    return result;
}
static int[,] Multiply(int[,] m1, int[,] m2)
{
    int l = m1.GetLength(1);
    if (m2.GetLength(0) != l)
        throw new ArgumentException();
    var result = new int[m1.GetLength(0), m2.GetLength(1)];
    for (int i = 0; i < m1.GetLength(0); i++)
        for (int j = 0; j < m2.GetLength(1); j++)
        {
            int sum = 0;
            for (int k = 0; k < l; k++)
                sum += m1[i, k] * m2[k, j];
            result[i, j] = sum;
        }
    return result;
}

@МВлад · 08.01.2015, 00:37 **[ТС]**

C#

...
Console.WriteLine("Введите размерность матрицы:");
int n = Convert.ToInt32(Console.ReadLine());
var a = new int[n, n];
int i, j;
for (i = 0; i < n; i++)
{
    for (j = 0; j < n; j++)
    {
        Console.Write("Введите элемент [" + i + "][" + j + "]: ");
        a[i, j] = Convert.ToInt32(Console.ReadLine());
     }
}
Console.WriteLine("\n");
int power = 2;
MatrixPower2(a, power);
...

Видимо, я что-то делаю не так, да? Неправильно массив заполнил? Или, скорее всего, неправильно передал? Ибо ответ неправильный он выдаёт. Не подскажете, в чём провинился?

@datey · 08.01.2015, 00:54

http://rghost.ru/60197862
Вот. Сохранение переделай, а том я с ним провозился, но сохраняет не так как нада (с тегами).

Добавлено через 1 минуту
Блюааа. Я сделал в форме... Не заметил про консоль...

@МВлад · 08.01.2015, 01:23 **[ТС]**

Ничего, я хоть на алгоритм посмотрю. Огромное Вам спасибо за такую работу большую!

Добавлено через 10 минут
Не, что-то там совсем не то...

Psilon · 08.01.2015, 03:01

C#

using System;
using System.Collections;
using System.Collections.Generic;
using System.Runtime.CompilerServices;
using System.Threading.Tasks;
 
namespace ConsoleApplication18
{
    public class Matrix : IEnumerable<double>
    {
        const string DimentionExceptionMessage = "Invalid matrix dimentions";
        public static Matrix GetE(int n)
        {
            var result = new Matrix(n, n);
            for (int i = 0; i < n; i++)
            {
                result[i, i] = 1.0;
            }
            return result;
        }
 
        private readonly double[,] _value;
 
        public Matrix(int rows, int cols)
            : this(new double[rows, cols])
        {
        }
 
        public Matrix(double[,] value)
        {
            _value = (double[,])value.Clone();
        }
 
        public static Matrix operator *(Matrix one, Matrix other)
        {
            return new Matrix(Multiply(one._value, other._value));
        }
 
        public static Matrix operator *(double scalar, Matrix matrix)
        {
            var result = (double[,])matrix._value.Clone();
            for (int i = 0; i < result.GetLength(0); i++)
            {
                for (int j = 0; j < result.GetLength(1); j++)
                {
                    result[i, j] *= scalar;
                }
            }
            return result;
        }
 
        public static Matrix operator *(Matrix matrix, double scalar)
        {
            return scalar * matrix;
        }
 
        public static Matrix operator +(Matrix one, Matrix another)
        {
            return GetMatrix(one, another, 1);
        }
 
        public static Matrix operator -(Matrix one, Matrix another)
        {
            return GetMatrix(one, another, -1);
        }
 
        private static Matrix GetMatrix(Matrix one, Matrix another, int sign)
        {
            var m = GetDim(one, another, 0);
            int n = GetDim(one, another, 1);
            var res = new double[m, n];
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    res[i, j] = one[i, j] + another[i, j] * sign;
                }
            }
            return new Matrix(res);
        }
 
        private static int GetDim(Matrix one, Matrix another, int dim)
        {
            int dimOne = one._value.GetLength(dim);
            int dimAnother = another._value.GetLength(dim);
 
            if (dimOne != dimAnother)
            {
 
                throw new ArgumentException(DimentionExceptionMessage);
            }
            return dimOne;
        }
 
        private int GetDimIfSquare()
        {
            int x = _value.GetLength(0);
            int y = _value.GetLength(1);
            if (x != y)
            {
                throw new ArgumentException(DimentionExceptionMessage);
            }
            return x;
        }
 
        public static Matrix operator /(Matrix matrix, double scalar)
        {
            return matrix * (1.0 / scalar);
        }
 
        public double this[int i, int j]
        {
            get { return _value[i, j]; }
            set { _value[i, j] = value; }
        }
 
        public IEnumerator<double> GetEnumerator()
        {
            // ReSharper disable once LoopCanBeConvertedToQuery
            foreach (double d in _value)
            {
                yield return d;
            }
        }
 
        IEnumerator IEnumerable.GetEnumerator()
        {
            return GetEnumerator();
        }
 
        public static implicit operator Matrix(double[,] doubles)
        {
            return new Matrix(doubles);
        }
 
        public static explicit operator double[,](Matrix matrix)
        {
            return matrix._value;
        }
 
        private static double[,] Multiply(double[,] a, double[,] b)
        {
            if (a.GetLength(1) != b.GetLength(0))
                throw new ArgumentException(DimentionExceptionMessage);
            var result = new double[a.GetLength(0), b.GetLength(1)];
            Parallel.For(0, a.GetLength(0), i =>
            {
                for (int j = 0; j < b.GetLength(1); j++)
                    result[i, j] = MultiplyR(a, b, i, j);
            });
            return result;
        }
 
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        private static double MultiplyR(double[,] a, double[,] b, int i, int j)
        {
            double result = 0;
            for (int k = 0; k < a.GetLength(1); k++)
                result += a[i, k] * b[k, j];
            return result;
        }
 
        public int GetLength(int dim)
        {
            return _value.GetLength(dim);
        }
 
        public Matrix Pow(int power)
        {
            int dim = GetDimIfSquare();
            var result = GetE(dim);
            var tmp = this;
            while (power != 0)
            {
                if (power % 2 != 0)
                {
                    result *= tmp;
                    power -= 1;
                }
                tmp *= tmp;
                power /= 2;
            }
            return result;
        }
    }
}

Добавлено через 44 секунды

C#

        static void Main(string[] args)
        {
            double[,] a =
            {
                {1, 2},
                {3, 4}
            };
            var matrix = new Matrix(a);
            var pow = matrix.Pow(5);
            for (int i = 0; i < pow.GetLength(0); i++)
            {
                for (int j = 0; j < pow.GetLength(1); j++)
                {
                    Console.Write(pow[i,j] + " ");
                }
                Console.WriteLine();
            }
        }

Добавлено через 10 минут
Возведение вашего файла в 256 степень например:

C#

using System;
using System.IO;
 
namespace ConsoleApplication18
{
    class Program
    {
        static void Main()
        {
            var a = GetMatrixFromFile("test1.A");
            
            var matrix = new Matrix(a);
            var pow = matrix.Pow(256);
            using (var sw = new StreamWriter("result.txt"))
            {
                for (int i = 0; i < pow.GetLength(0); i++)
                {
                    for (int j = 0; j < pow.GetLength(1); j++)
                    {
                        sw.Write(pow[i, j] + " ");
                    }
                    sw.WriteLine();
                }
            }
        }
 
        private static double[,] GetMatrixFromFile(string path)
        {
            string[] lines = File.ReadAllLines(path);
            double[,] a = new double[lines.Length, lines.Length];
            
            for (int i = 0; i < lines.Length; i++)
            {
                var split = lines[i].Split(new[] {' '}, StringSplitOptions.RemoveEmptyEntries);
                for (int j = 0; j < lines.Length; j++)
                {
                    a[i, j] = double.Parse(split[j]);
                }
            }
            return a;
        }
    }
}

Psilon · 08.01.2015, 04:14

После небольших оптимизаций лично на моем компе обсчитывает 256 степень этой матрицы за 4 с небольшим секунд:

C#

using System;
using System.Collections;
using System.Collections.Generic;
using System.Runtime.CompilerServices;
using System.Threading;
using System.Threading.Tasks;
 
namespace ConsoleApplication18
{
    public class Matrix : IEnumerable<double>
    {
        const string DimentionExceptionMessage = "Invalid matrix dimentions";
        public static Matrix GetE(int n)
        {
            var result = new Matrix(n, n);
            for (int i = 0; i < n; i++)
            {
                result[i, i] = 1.0;
            }
            return result;
        }
 
        private readonly double[,] _value;
 
        public Matrix(int rows, int cols)
            : this(new double[rows, cols])
        {
        }
 
        public Matrix(double[,] value)
        {
            _value = (double[,])value.Clone();
        }
 
        public static Matrix operator *(Matrix one, Matrix other)
        {
            return new Matrix(Multiply(one._value, other._value));
        }
 
        public static Matrix operator *(double scalar, Matrix matrix)
        {
            var result = (double[,])matrix._value.Clone();
            for (int i = 0; i < result.GetLength(0); i++)
            {
                for (int j = 0; j < result.GetLength(1); j++)
                {
                    result[i, j] *= scalar;
                }
            }
            return result;
        }
 
        public static Matrix operator *(Matrix matrix, double scalar)
        {
            return scalar * matrix;
        }
 
        public static Matrix operator +(Matrix one, Matrix another)
        {
            return GetMatrix(one, another, 1);
        }
 
        public static Matrix operator -(Matrix one, Matrix another)
        {
            return GetMatrix(one, another, -1);
        }
 
        private static Matrix GetMatrix(Matrix one, Matrix another, int sign)
        {
            var m = GetDim(one, another, 0);
            int n = GetDim(one, another, 1);
            var res = new double[m, n];
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    res[i, j] = one[i, j] + another[i, j] * sign;
                }
            }
            return new Matrix(res);
        }
 
        private static int GetDim(Matrix one, Matrix another, int dim)
        {
            int dimOne = one._value.GetLength(dim);
            int dimAnother = another._value.GetLength(dim);
 
            if (dimOne != dimAnother)
            {
 
                throw new ArgumentException(DimentionExceptionMessage);
            }
            return dimOne;
        }
 
        private int GetDimIfSquare()
        {
            int x = _value.GetLength(0);
            int y = _value.GetLength(1);
            if (x != y)
            {
                throw new ArgumentException(DimentionExceptionMessage);
            }
            return x;
        }
 
        public static Matrix operator /(Matrix matrix, double scalar)
        {
            return matrix * (1.0 / scalar);
        }
 
        public double this[int i, int j]
        {
            get { return _value[i, j]; }
            set { _value[i, j] = value; }
        }
 
        public IEnumerator<double> GetEnumerator()
        {
            // ReSharper disable once LoopCanBeConvertedToQuery
            foreach (double d in _value)
            {
                yield return d;
            }
        }
 
        IEnumerator IEnumerable.GetEnumerator()
        {
            return GetEnumerator();
        }
 
        public static implicit operator Matrix(double[,] doubles)
        {
            return new Matrix(doubles);
        }
 
        public static explicit operator double[,](Matrix matrix)
        {
            return matrix._value;
        }
 
        private static double[,] Multiply(double[,] a, double[,] b)
        {
            if (a.GetLength(1) != b.GetLength(0))
                throw new ArgumentException(DimentionExceptionMessage);
            return SafeMultiply(a, b);
        }
 
        private static double[,] SafeMultiply(double[,] a, double[,] b)
        {
            const int blocksize = 8;
            var imax = a.GetLength(0);
            int jmax = b.GetLength(1);
            int kmax = a.GetLength(1);
 
            var result = new double[imax, jmax];
            int completedTasks = 0, totalTasks = (int) Math.Ceiling(imax/(double) blocksize);
            var reset = new ManualResetEvent(false);
            for (int i0 = 0; i0 < imax; i0 += blocksize)
            {
                int ilocal = i0;
                ThreadPool.QueueUserWorkItem(state =>
                {
                    for (int j0 = 0; j0 < jmax; j0 += blocksize)
                        for (int k0 = 0; k0 < kmax; k0 += blocksize)
                            for (int i = ilocal; i < Math.Min(imax, ilocal + blocksize); i++)
                                for (int j = j0; j < Math.Min(jmax, j0 + blocksize); j++)
                                    for (int k = k0; k < Math.Min(kmax, k0 + blocksize); k++)
                                    {
                                        result[i, j] += a[i, k]*b[k, j];
                                    }
                    if (Interlocked.Increment(ref completedTasks) == totalTasks)
                    {
                        reset.Set();
                    }
                });
            }
            reset.WaitOne();
            return result;
        }
 
 
        public int GetLength(int dim)
        {
            return _value.GetLength(dim);
        }
 
        public Matrix Pow(int power)
        {
            int dim = GetDimIfSquare();
            var result = GetE(dim)._value;
            var tmp = _value;
            while (power != 0)
            {
                if (power % 2 != 0)
                {
                    result = SafeMultiply(tmp, result);
                    power -= 1;
                }
                tmp = SafeMultiply(tmp, tmp);
                power /= 2;
            }
            return result;
        }
    }
}

Psilon · 08.01.2015, 13:24

В общем после еще тучи всяких попыток на моей машине получилось выжать 3 с небольшим секунды для возведения в 256 степень:

C#

using System;
using System.Collections;
using System.Collections.Generic;
using System.Runtime.CompilerServices;
using System.Threading;
 
namespace ConsoleApplication18
{
    public class Matrix : IEnumerable<double>
    {
        const string DimentionExceptionMessage = "Invalid matrix dimensions";
        public static Matrix GetE(int n)
        {
            var result = new Matrix(n, n);
            for (int i = 0; i < n; i++)
            {
                result[i, i] = 1.0;
            }
            return result;
        }
 
        private readonly double[,] _value;
 
        public Matrix(int rows, int cols)
            : this(new double[rows, cols])
        {
        }
 
        public Matrix(double[,] value)
        {
            _value = (double[,])value.Clone();
        }
 
        public static Matrix operator *(Matrix one, Matrix other)
        {
            return new Matrix(Multiply(one._value, other._value));
        }
 
        public static Matrix operator *(double scalar, Matrix matrix)
        {
            var result = (double[,])matrix._value.Clone();
            for (int i = 0; i < result.GetLength(0); i++)
            {
                for (int j = 0; j < result.GetLength(1); j++)
                {
                    result[i, j] *= scalar;
                }
            }
            return result;
        }
 
        public static Matrix operator *(Matrix matrix, double scalar)
        {
            return scalar * matrix;
        }
 
        public static Matrix operator +(Matrix one, Matrix another)
        {
            return MatrixSum(one, another, 1);
        }
 
        public static Matrix operator -(Matrix one, Matrix another)
        {
            return MatrixSum(one, another, -1);
        }
 
        private static Matrix MatrixSum(Matrix one, Matrix another, int sign)
        {
            var m = GetDim(one, another, 0);
            int n = GetDim(one, another, 1);
            var res = new double[m, n];
            for (int i = 0; i < m; i++)
            {
                for (int j = 0; j < n; j++)
                {
                    res[i, j] = one[i, j] + another[i, j] * sign;
                }
            }
            return new Matrix(res);
        }
 
        private static int GetDim(Matrix one, Matrix another, int dim)
        {
            int dimOne = one._value.GetLength(dim);
            int dimAnother = another._value.GetLength(dim);
 
            if (dimOne != dimAnother)
                throw new ArgumentException(DimentionExceptionMessage);
            return dimOne;
        }
 
        private int GetDimIfSquare()
        {
            int x = _value.GetLength(0);
            int y = _value.GetLength(1);
            if (x != y)
                throw new ArgumentException(DimentionExceptionMessage);
            return x;
        }
 
        public static Matrix operator /(Matrix matrix, double scalar)
        {
            return matrix * (1.0 / scalar);
        }
 
        public double this[int i, int j]
        {
            get { return _value[i, j]; }
            set { _value[i, j] = value; }
        }
 
        public IEnumerator<double> GetEnumerator()
        {
            // ReSharper disable once LoopCanBeConvertedToQuery
            foreach (double d in _value)
            {
                yield return d;
            }
        }
 
        IEnumerator IEnumerable.GetEnumerator()
        {
            return GetEnumerator();
        }
 
        public static implicit operator Matrix(double[,] doubles)
        {
            return new Matrix(doubles);
        }
 
        public static explicit operator double[,](Matrix matrix)
        {
            return matrix._value;
        }
        public int GetLength(int dim)
        {
            return _value.GetLength(dim);
        }
 
        private static double[,] Multiply(double[,] a, double[,] b)
        {
            if (a.GetLength(1) != b.GetLength(0))
                throw new ArgumentException(DimentionExceptionMessage);
            return SafeMultiply(a, b, a.GetLength(0), b.GetLength(1), a.GetLength(1));
        }
 
 
        public Matrix Pow(int power)
        {
            int dim = GetDimIfSquare();
            var result = GetE(dim)._value;
            var tmp = _value;
            while (power != 0)
            {
                if (power % 2 != 0)
                {
                    result = SafeMultiply(tmp, result, dim, dim, dim);
                    power -= 1;
                }
                tmp = SafeMultiply(tmp, tmp, dim, dim, dim);
                power /= 2;
            }
            return result;
        }
 
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        static double[,] SafeMultiply(double[,] a, double[,] b, int imax, int jmax, int kmax)
        {
            double[,] c = new double[imax, jmax];
            int blockCount = Environment.ProcessorCount;
            int blockSize = (int) Math.Ceiling((double) jmax/blockCount);
            int completedBlocks = 0;
            var reset = new ManualResetEvent(false);
            for (int block = 0; block < blockCount; block++)
            {
                int from = block*blockSize;
                int to = Math.Min(from + blockSize, jmax);
                ThreadPool.QueueUserWorkItem(state =>
                {
                    var bcolj = new double[kmax];
                    for (int j = from; j < to; j++)
                    {
                        for (var k = 0; k < kmax; k++)
                            bcolj[k] = b[k, j];
 
                        for (var i = 0; i < imax; i++)
                        {
                            double s = 0;
                            for (var k = 0; k < kmax; k++)
                                s += a[i, k] * bcolj[k];
                            c[i, j] = s;
                        }
                    }
                    if (Interlocked.Increment(ref completedBlocks) == blockCount)
                    {
                        reset.Set();
                    }
                });
            }
            reset.WaitOne();
            return c;
        }
    }
}

@МВлад · 09.01.2015, 03:57 **[ТС]**

Воу!

Такого ответа я совсем не ожидал! Даже не знаю, как Вас отблагодарить за такую огромную и великолепную работу, но Вы меня очень здорово выручили. Все Ваши методы работают, и работают быстро и правильно. Думаю, решения лучше этого будет очень трудно найти!

Спасибо Вам огромное за Ваши старания!

Psilon · 09.01.2015, 04:28

МВлад, да не за что

Заодно потренировался в оптимизациях. Была бы 15 студия, впихнул бы AVX, возведение в 1024 степень происходило бы меньше чем за секунду скорее всего. Но, что поделать, пока что майкрософт её не выпустил.

Особенно приятно, когда работу ценят, а не "шо это за хрень, а где комментарии? че-та я преподу объяснить не смог, как работает"

Добавлено через 57 секунд
МВлад, и да, с днюхой, братан!

Не по теме:

и не надо на вы обращаться, мне и 22 еще нет хД

@МВлад · 09.01.2015, 08:19 **[ТС]**

Я очень ценю Ваш труд, правда, ибо работа сделана большая и долгая, причём, я не ожидал, что мне тут помогут (я первый раз за всю жизнь обратился за помощью на форум), да ещё и полную программу сделают. Не знал, как передать это словами, подумал, что Вы не поверите и не стал это писать... :-)

Спасибо большое за поздравление. Честно говоря, Вы первый, кто меня поздравил! :-)
Ещё раз спасибо за Вашу программу! Думаю, она ещё многим поможет. :-)

Psilon Master of Orion 6100 / 4956 / 905 Регистрация: 10.07.2011 Сообщений: 14,522 Записей в блоге: 5
	09.01.2015, 04:28	12
	МВлад, да не за что Заодно потренировался в оптимизациях. Была бы 15 студия, впихнул бы AVX, возведение в 1024 степень происходило бы меньше чем за секунду скорее всего. Но, что поделать, пока что майкрософт её не выпустил. Особенно приятно, когда работу ценят, а не "шо это за хрень, а где комментарии? че-та я преподу объяснить не смог, как работает" Добавлено через 57 секунд МВлад, и да, с днюхой, братан! Не по теме: и не надо на вы обращаться, мне и 22 еще нет хД 0

Блоги программистов
Как использовать GraphQL в C# с HotChocolate Programming 05.01.2025 GraphQL — это современный подход к разработке API, который позволяет клиентам запрашивать только те данные, которые им необходимы. Это делает взаимодействие с API более гибким и эффективным по. . .	Модель полного двоичного суматора с помощью логических операций (python) AlexSky-coder 04.01.2025 def binSum(x:list, y:list): s=^y] p=x and y for i in range(1,len(x)): s. append((x^y)^p) p=(x and y)or(p and (x or y)) return s x=list() y=list()	Это мы не проходили, это нам не задавали...(шестибитный асихронный счётчик с управляющим сигналом задержки). Hrethgir 04.01.2025 Асинхронный счётчик на сумматорах (шестиразрядный по числу диодов на плате, но наверное разрядов будет больше - восемь или шестнадцать, а диоды на старшие), так как триггеры прошли тестирование и. . .	Руководство по созданию бота для Телеграм на Python IT_Exp 04.01.2025 Боты для Телеграм представляют собой автоматизированные программы, которые выполняют различные задачи, взаимодействуя с пользователями через интерфейс мессенджера. В данной статье мы рассмотрим,. . .	Применение компонентов PrimeVue в Vue.js 3 на TypeScript BasicMan 04.01.2025 Введение в PrimeVue и настройка окружения PrimeVue представляет собой мощную библиотеку компонентов пользовательского интерфейса для Vue. js 3, которая предоставляет разработчикам богатый набор. . .	Как стать Senior developer cpp_developer 04.01.2025 В современной индустрии разработки программного обеспечения позиция Senior Developer представляет собой не просто следующую ступень карьерной лестницы, а качественно новый уровень профессионального. . .
Что известно о дате выхода Windows 12 и чего от нее ждать IT_Exp 04.01.2025 В мире технологий постоянно происходят изменения, и операционные системы не являются исключением. Windows 11, выпущенная в октябре 2021 года, принесла множество инноваций и улучшений, но. . .	Что новенького в .NET Core 9 Programming 04.01.2025 Обзор ключевых изменений в . NET Core 9 Платформа . NET Core продолжает активно развиваться, и версия 9 представляет собой значительный шаг вперед в эволюции этой технологии. Новый релиз. . .	Инструкция по установке python3.13.1 в Debian 12 AlexSky-coder 03.01.2025 sudo apt update sudo apt install build-essential zlib1g-dev libncurses5-dev libgdbm-dev libnss3-dev libssl-dev libreadline-dev libffi-dev wget. . .	Затестил триггеры. архив проекта прилагаю с GOA файлами в настройках архиватора проектов. Hrethgir 03.01.2025 В этот раз нет закольцованности, потому что от неё только глюки, как я понял, логика не вырезанная. Триггеры очень быстрые если верить измерениям с помощью анализатора от Gowin. Есть ещё регистры,. . .	Python в помощь DevOps IT_Exp 03.01.2025 Причины использования Python в работе DevOps Python стал неотъемлемой частью мира DevOps, и это не случайно. Этот язык программирования обладает множеством преимуществ, которые делают его. . .	Angular vs React vs Vue.js BasicMan 03.01.2025 О, друзья-разработчики и просто любопытные читатели! Сегодня мы отправимся в увлекательное путешествие по миру фронтенд-разработки, и первой остановкой станет Angular – этакий строгий немецкий. . .

Возведение большой матрицы, читаемой с файла, в степень N, и запись обратно в файл на C#

Решение

Решение

@МВлад 21 / 4 / 2 Регистрация: 07.01.2015 Сообщений: 45
		1
	Возведение большой матрицы, читаемой с файла, в степень N, и запись обратно в файл на C# 07.01.2015, 20:34. Показов 1340. Ответов 12 Метки нет (Все метки) Доброго времени. Очень нужно написать программу на C#, которая прочитает с файла огромную матрицу с вещественными числами (допустим, 1000 на 1000), возвести её в степень N, и результат записать в новый файл. Вот так... Думал, сам напишу, убил кучу времени, и теперь, от безысходности, пишу тут. Помогите пожалуйста. За ранее благодарен всем, кто уделит время. И с Рождеством всех, ребята. =) P.S.: это нужно делать в консольном приложении. Очень срочно. Добавлено через 1 минуту Пример матрицы тут можно взять (363х363): http://rghost.ru/60192711 0

@chumich 2081 / 1239 / 464 Регистрация: 20.12.2014 Сообщений: 3,237
	07.01.2015, 22:35	3
	Посмотрите здесь, может быть подойдет: http://easylab.net.ua/massivyi... tritsami-1 или здесь: Возведение булевых матриц в степень 0

@datey 25 / 25 / 10 Регистрация: 14.12.2014 Сообщений: 110
	08.01.2015, 00:54	6
	http://rghost.ru/60197862 Вот. Сохранение переделай, а том я с ним провозился, но сохраняет не так как нада (с тегами). Добавлено через 1 минуту Блюааа. Я сделал в форме... Не заметил про консоль... 1

@МВлад 21 / 4 / 2 Регистрация: 07.01.2015 Сообщений: 45
	08.01.2015, 01:23 [ТС]	7
	Ничего, я хоть на алгоритм посмотрю. Огромное Вам спасибо за такую работу большую! Добавлено через 10 минут Не, что-то там совсем не то... 0

@МВлад 21 / 4 / 2 Регистрация: 07.01.2015 Сообщений: 45
	09.01.2015, 03:57 [ТС]	11
	Воу! Такого ответа я совсем не ожидал! Даже не знаю, как Вас отблагодарить за такую огромную и великолепную работу, но Вы меня очень здорово выручили. Все Ваши методы работают, и работают быстро и правильно. Думаю, решения лучше этого будет очень трудно найти! Спасибо Вам огромное за Ваши старания! 0

@МВлад 21 / 4 / 2 Регистрация: 07.01.2015 Сообщений: 45
	09.01.2015, 08:19 [ТС]	13
	Я очень ценю Ваш труд, правда, ибо работа сделана большая и долгая, причём, я не ожидал, что мне тут помогут (я первый раз за всю жизнь обратился за помощью на форум), да ещё и полную программу сделают. Не знал, как передать это словами, подумал, что Вы не поверите и не стал это писать... :-) Спасибо большое за поздравление. Честно говоря, Вы первый, кто меня поздравил! :-) Ещё раз спасибо за Вашу программу! Думаю, она ещё многим поможет. :-) 0