Advertisement
myname0

Gauss

Nov 8th, 2016
95
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C# 2.50 KB | None | 0 0
  1. using System;
  2. using System.Collections.Generic;
  3. using System.Linq;
  4. using System.Text;
  5. using System.IO;
  6.  
  7.  
  8. namespace Gauss
  9. {
  10.     class Program
  11.     {
  12.         static List<List<double>> GetMatrix()
  13.         {
  14.             StreamReader fin = new StreamReader("input.txt");
  15.             int n = int.Parse(fin.ReadLine());
  16.             int m = int.Parse(fin.ReadLine());
  17.  
  18.             List<List<double>> matrix = new List<List<double>>();
  19.  
  20.             for (int i = 0; i < n; i++)
  21.             {
  22.                 List<double> tmp = new List<double>();
  23.                 string[] line = fin.ReadLine().Split(' ');
  24.                 for (int j = 0; j < m; j++)
  25.                     tmp.Add(double.Parse(line[j]));
  26.                 matrix.Add(tmp);
  27.             }
  28.  
  29.             fin.Close();
  30.  
  31.             return matrix;
  32.         }
  33.  
  34.         static void RightProcess(List<List<double>> matrix)
  35.         {
  36.             double diagonalElem = 0;
  37.             int count = matrix.Count + 1;
  38.             double firstElem = 0;
  39.  
  40.             for (int i = 0; i < matrix.Count; i++)
  41.             {
  42.                 diagonalElem = matrix[i][i];
  43.                 if (diagonalElem != 0)
  44.                 {
  45.                     for (int j = i; j < count; j++)
  46.                         matrix[i][j] /= diagonalElem;
  47.  
  48.                     for(int j = i + 1; j < matrix.Count; j++)
  49.                     {
  50.                         firstElem = matrix[j][i];
  51.                         for (int k = i; k < count; k++)
  52.                             matrix[j][k] -= matrix[i][k] * firstElem;
  53.                     }
  54.                 }
  55.                 else throw new Exception("Matrix have null diagonal element!");
  56.             }
  57.         }
  58.  
  59.         static void BackProcess(List<List<double>> matrix)
  60.         {
  61.             int current = matrix.Count - 1;
  62.            
  63.             for (int i = matrix.Count - 2; i >= 0; i--)
  64.             {
  65.                 for (int j = matrix.Count; j >= 0; j--)
  66.                 {
  67.                     for(int k = 0; k <= i; k++)
  68.  
  69.                         matrix[k][j] -= matrix[i + 1][j] * matrix[k][current];
  70.                 }
  71.                 current--;
  72.             }
  73.  
  74.             for (int i = matrix.Count - 1; i >= 0; i--)
  75.                 Console.Write("{0} ", matrix[i][matrix.Count]);
  76.             Console.WriteLine();
  77.  
  78.         }
  79.  
  80.         static void Main(string[] args)
  81.         {
  82.             List<List<double>> matrix = GetMatrix();
  83.             RightProcess(matrix);
  84.             BackProcess(matrix);
  85.         }
  86.     }
  87. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement