fork download
  1. #include<bits/stdc++.h>
  2. using namespace std;
  3.  
  4. int main ()
  5. {
  6. int n; cin >> n;
  7. double a[101][101+1], x[100], ratio;
  8. for (int i=1; i<=n; i++){
  9. for (int j=1; j<=n+1; j++){
  10. cin >> a[i][j];
  11. }
  12. }
  13.  
  14. for (int i=1; i<n; i++){
  15. if (a[i][i]==0){
  16. cout << "Error!" << endl;
  17. break;
  18. }
  19. for (int j=i+1; j<=n; j++){
  20. ratio = a[j][i]/a[i][i];
  21. for (int k=1; k<=n+1; k++){
  22. a[j][k] = a[j][k] - a[i][k]*ratio;
  23. }
  24. }
  25. }
  26.  
  27. for (int i=1; i<=n; i++){
  28. for (int j=1; j<=n+1; j++){
  29. cout << a[i][j] << " ";
  30. }
  31. cout << endl;
  32. }
  33.  
  34. for (int i = n; i >= 2; i--) {
  35. for (int j = i - 1; j >= 1; j--) {
  36. double ratio = a[j][i] / a[i][i];
  37. for (int k = 1; k <= n + 1; k++) {
  38. a[j][k] = a[j][k] - ratio * a[i][k];
  39. }
  40. }
  41. }
  42.  
  43. for (int i=1; i<=n; i++){
  44. for (int j=1; j<=n+1; j++){
  45. cout << a[i][j] << " ";
  46. }
  47. cout << endl;
  48. }
  49.  
  50. for (int i=1; i<=n; i++){
  51. x[i]= a[i][n+1]/a[i][i];
  52. }
  53. for (int i=1; i<=n; i++){
  54. cout << x[i] <<endl;
  55. }
  56.  
  57.  
  58. }
Success #stdin #stdout 0.01s 5316KB
stdin
1 1 1 5
2 3 5 8
4 0 5 2
stdout
1 1 
1 1 
1