Advertisement
Guest User

Untitled

a guest
Mar 22nd, 2019
62
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.98 KB | None | 0 0
  1. #include <cstdio>
  2. using namespace std;
  3.  
  4. const int LEN = 2e5 + 9;
  5. const int INF = 2e9 + 9;
  6.  
  7. struct Matrix {
  8. // a b
  9. // c d
  10.  
  11. int a = INF, b = INF, c = INF, d = INF;
  12.  
  13. Matrix(int _a, int _b, int _c, int _d) {
  14. a = _a;
  15. b = _b;
  16. c = _c;
  17. d = _d;
  18. }
  19.  
  20. Matrix() {
  21. a = INF;
  22. b = INF;
  23. c = INF;
  24. d = INF;
  25. }
  26. };
  27.  
  28. bool operator == (Matrix a, int b) {
  29. return a.a == b && a.b == b && a.c == b && a.d == b;
  30. }
  31.  
  32. int r, n, m;
  33. Matrix tree[LEN * 4];
  34. int sz = 1;
  35.  
  36. int getMod(long long a, int b) {
  37. if (a < b) {
  38. return a;
  39. }
  40.  
  41. return a % b;
  42. }
  43.  
  44. Matrix operator * (Matrix A, Matrix B) {
  45. if (A == INF)
  46. return B;
  47. if (B == INF)
  48. return A;
  49.  
  50. int a = getMod(getMod((A.a * B.a), r) + getMod((A.b * B.c), r), r);
  51. int b = getMod(getMod(A.a * B.b, r) + getMod(A.b * B.d, r), r);
  52. int c = getMod(getMod(A.c * B.a, r) + getMod(A.d * B.c, r), r);
  53. int d = getMod(getMod(A.c * B.b, r) + getMod(A.d * B.d, r), r);
  54.  
  55. return Matrix(a, b, c, d);
  56. }
  57.  
  58. Matrix multiply(int v, int l, int r, int a, int b) {
  59. if (r < a || l > b)
  60. return Matrix();
  61.  
  62. if (a <= l && r <= b)
  63. return tree[v];
  64.  
  65. int m = (l + r) / 2;
  66. return multiply(v * 2, l, m, a, b) * multiply(v * 2 + 1, m + 1, r, a, b);
  67. }
  68.  
  69. int main() {
  70. freopen("crypto.in", "r", stdin);
  71. freopen("crypto.out", "w", stdout);
  72.  
  73. while (sz <= LEN) sz *= 2;
  74.  
  75. scanf("%d%d%d", &r, &n, &m);
  76.  
  77. int a, b, c, d;
  78. for (int i = 0; i < n; i++) {
  79. scanf("%d%d%d%d", &a, &b, &c, &d);
  80.  
  81. tree[sz + i] = Matrix(a, b, c, d);
  82. }
  83.  
  84. for (int i = sz - 1; i >= 1; i--) {
  85. tree[i] = tree[i * 2] * tree[i * 2 + 1];
  86. }
  87.  
  88. int l, r;
  89. for (int i = 0; i < m; i++) {
  90. scanf("%d%d", &l, &r);
  91.  
  92. Matrix ans = multiply(1, 1, sz, l, r);
  93.  
  94. printf("%d %d\n%d %d\n\n", ans.a, ans.b, ans.c, ans.d);
  95. }
  96.  
  97. return 0;
  98. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement