Advertisement
Guest User

Untitled

a guest
Nov 15th, 2019
135
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.02 KB | None | 0 0
  1. #include <bits/stdc++.h>
  2.  
  3. using namespace std;
  4.  
  5. ifstream f("ecu.in");
  6. ofstream g("ecu.out");
  7.  
  8. const int NMAX = 35;
  9. int n,m;
  10. double x[NMAX], p[NMAX], y[NMAX], d[NMAX];
  11. double mat[NMAX][NMAX], c[NMAX][NMAX], ans[NMAX][NMAX];
  12.  
  13. void afisaremat(double mat[NMAX][NMAX]){
  14. int i,j;
  15. for(i = 1 ; i <= n + 1 ; i++){
  16. for(j = 1 ; j <= n + 1 ; j++)
  17. g << mat[i][j] << " ";
  18. g << "\n";
  19. }
  20. }
  21.  
  22. void inm(double a[NMAX][NMAX], double b[NMAX][NMAX]){
  23. int i,j,k;
  24.  
  25. for(i = 1 ; i <= n + 1 ; i++){
  26. for(j = 1 ; j <= n + 1 ; j++){
  27. c[i][j] = 0;
  28. for(k = 1 ; k <= n + 1 ; k++)
  29. c[i][j] += a[i][k] * b[k][j];
  30. }
  31. }
  32.  
  33. for(i = 1 ; i <= n + 1 ; i++)
  34. for(j = 1 ; j <= n + 1 ; j++)
  35. a[i][j] = c[i][j];
  36. }
  37.  
  38. void power(int p){
  39. int i,j;
  40.  
  41. for(i = 1 ; i <= n + 1 ; i++)
  42. for(j = 1 ; j <= n + 1 ; j++)
  43. ans[i][j] = mat[i][j];
  44.  
  45. while(p){
  46. if(p & 1){
  47. p--;
  48. inm(ans,mat);
  49. }
  50. inm(ans, ans);
  51. p /= 2;
  52. }
  53.  
  54. for(i = 1 ; i <= n + 1; i++)
  55. for(j = 1 ; j <= n + 1 ; j++)
  56. mat[i][j] = ans[i][j];
  57. }
  58.  
  59. void alh(){
  60. int i,j;
  61.  
  62. for(i = 1 ; i <= n ; i++){
  63. d[i] = 0;
  64.  
  65. for(j = 1 ; j <= n + 1 ; j++)
  66. d[i] += (x[j] * mat[j][i]);
  67. }
  68.  
  69. }
  70.  
  71. int main(){
  72. int i,j;
  73. f >> n >> m;
  74. for(i = 1 ; i <= n ; i++)
  75. f >> x[i];
  76. for(i = 1 ; i <= n ; i++)
  77. f >> p[i];
  78. for(i = 1 ; i <= n ; i++)
  79. f >> y[i];
  80.  
  81. for(i = 1 ; i < n ; i++){
  82. mat[i][i] = p[i];
  83. mat[i + 1][i] = 1 - p[i + 1];
  84. mat[n + 1][i] = y[i];
  85. //mat[i][n + 1] = 1;
  86. }
  87.  
  88. //mat[n][n + 1] = 1;
  89. mat[n + 1][n + 1] = 1;
  90. mat[1][n] = 1 - p[1];
  91. mat[n][n] = p[n];
  92. mat[n + 1][n] = y[n];
  93. x[n + 1] = 1;
  94.  
  95. power(2);
  96.  
  97. //afisaremat(mat);
  98.  
  99. alh();
  100.  
  101. for(i = 1 ; i <= n ; i++)
  102. g << d[i] << " ";
  103.  
  104. return 0;
  105. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement