sacgajcvs

Problem K

Nov 9th, 2019
266
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 4.18 KB | None | 0 0
  1. /*
  2. _____ _ _ _ _
  3. |_ _| |__ ___ / \ _ __ ___| |__ _ _| |
  4. | | | '_ \ / _ \ / _ \ | '_ \/ __| '_ \| | | | |
  5. | | | | | | __// ___ \| | | \__ \ | | | |_| | |
  6. |_| |_| |_|\___/_/ \_\_| |_|___/_| |_|\__,_|_|
  7.  
  8. */
  9. #include<bits/stdc++.h>
  10. #define ll long long
  11. #define pb push_back
  12. #define ppb pop_back
  13. #define endl '\n'
  14. #define mii map<ll int,ll int>
  15. #define msi map<string,ll int>
  16. #define mis map<ll int, string>
  17. #define rep(i,a,b) for(ll int i=a;i<b;i++)
  18. #define mpi map<pair<ll int,ll int>,ll int>
  19. #define pii pair<ll int,ll int>
  20. #define vi vector<ll int>
  21. #define vii vector<pair<ll int, ll int>>
  22. #define vs vector<string>
  23. #define all(a) (a).begin(),(a).end()
  24. #define F first
  25. #define S second
  26. #define sz(x) (ll int)x.size()
  27. #define hell 1000000007
  28. #define lbnd lower_bound
  29. #define ubnd upper_bound
  30. #define bs binary_search
  31. #define mp make_pair
  32. #define what_is(x) cerr << #x << " is " << x << endl;
  33. #define time cerr << "\nTime elapsed: " << 1000 * clock() / CLOCKS_PER_SEC << "ms\n";
  34. using namespace std;
  35. #define PI 3.141592653589793
  36. #define N 100005
  37. ll n,q;
  38. vector<vector<ll>> a[N],seg[4*N];
  39. bool lazy[4*N];
  40.  
  41. vector<vector<ll>> rev(vector<vector<ll>> p1)
  42. {
  43. vector<vector<ll>> p2(2,vector<ll>(2));
  44. rep(i,0,2)
  45. {
  46. rep(j,0,2)
  47. {
  48. p2[i][j]=p1[(!i)][(!j)];
  49. }
  50. }
  51. return p2;
  52. }
  53.  
  54. vector<vector<ll>> prod(vector<vector<ll>> p1 ,vector<vector<ll>> p2)
  55. {
  56. vector<vector<ll>> p3(2,vector<ll>(2,0));
  57. rep(i,0,2)
  58. {
  59. rep(j,0,2)
  60. {
  61. rep(k,0,2)
  62. {
  63. p3[i][j]+=(p1[i][k]*p2[k][j])%hell;
  64. p3[i][j]%=hell;
  65. }
  66. }
  67. }
  68. return p3;
  69. }
  70.  
  71. vi fun(ll x,ll y, vector<vector<ll>> p1)
  72. {
  73. vi ans(2);
  74. ans[0]=((p1[0][0]*x)%hell+(p1[1][0]*y)%hell)%hell;
  75. ans[1]=((p1[0][1]*x)%hell+(p1[1][1]*y)%hell)%hell;
  76. return ans;
  77. }
  78.  
  79. void build(ll cur,ll st,ll end)
  80. {
  81. if(st==end)
  82. {
  83. seg[cur]=a[st];
  84. return;
  85. }
  86. ll mid=(st+end)>>1;
  87. build(2*cur,st,mid);
  88. build(2*cur+1,mid+1,end);
  89. seg[cur]=prod(seg[2*cur],seg[2*cur+1]); /* 1-change */
  90. }
  91. void updaterange(ll cur,ll st,ll end,ll l,ll r)
  92. {
  93. if(lazy[cur]!=0) /* 2-change */
  94. {
  95. seg[cur]=rev(seg[cur]); /* 3-change */
  96. if(st!=end)
  97. {
  98. lazy[2*cur]!=lazy[2*cur]; /* 4-change */
  99. lazy[2*cur+1]!=lazy[2*cur+1]; /* 5-change */
  100. }
  101. lazy[cur]=0; /* 6-change */
  102. }
  103. if(st>end||st>r||end<l) return;
  104. if(l<=st&&end<=r)
  105. {
  106. seg[cur]=rev(seg[cur]); /* 7-change */
  107. if(st!=end)
  108. {
  109. lazy[2*cur]!=lazy[2*cur]; /* 4-change */
  110. lazy[2*cur+1]!=lazy[2*cur+1]; /* 5-change */
  111. }
  112. return;
  113. }
  114. ll mid=(st+end)>>1;
  115. updaterange(2*cur,st,mid,l,r);
  116. updaterange(2*cur+1,mid+1,end,l,r);
  117. seg[cur]=prod(seg[2*cur],seg[2*cur+1]); /* 10-change */
  118. }
  119. vector<vector<ll>> query(ll cur,ll st,ll end,ll l,ll r)
  120. {
  121. if(st>end||st>r||end<l) return {{1,0},{0,1}};
  122. if(lazy[cur]!=0) /* 2-change */
  123. {
  124. seg[cur]=rev(seg[cur]); /* 3-change */
  125. if(st!=end)
  126. {
  127. lazy[2*cur]!=lazy[2*cur]; /* 4-change */
  128. lazy[2*cur+1]!=lazy[2*cur+1]; /* 5-change */
  129. }
  130. lazy[cur]=0; /* 6-change */
  131. }
  132. if(st>=l&&end<=r) return seg[cur];
  133. ll mid=(st+end)>>1;
  134. vector<vector<ll>> ans1=query(2*cur,st,mid,l,r);
  135. vector<vector<ll>> ans2=query(2*cur+1,mid+1,end,l,r);
  136. return prod(ans1,ans2); /* 14-change */
  137. }
  138. void solve()
  139. {
  140. //ll n,q;
  141. cin>>n>>q;
  142. string s;
  143. cin>>s;
  144. vector<vector<ll>> A(2,vector<ll>(2,1)),B(2,vector<ll>(2,1));
  145. A[0][1]=0;
  146. B[1][0]=0;
  147. rep(i,0,n)
  148. {
  149. if(s[i]=='A')
  150. a[i+1]=A;
  151. else
  152. a[i+1]=B;
  153. }
  154. build(1,1,n);
  155. rep(i,0,q)
  156. {
  157. ll type;
  158. cin>>type;
  159. if(type==1)
  160. {
  161. ll l,r;
  162. cin>>l>>r;
  163. updaterange(1,1,n,l,r);
  164. }
  165. else
  166. {
  167. ll l,r,x,y;
  168. cin>>l>>r>>x>>y;
  169. vi ans=fun(x,y,query(1,1,n,l,r));
  170. cout<<ans[0]<<" "<<ans[1]<<endl;
  171. }
  172. }
  173. return;
  174. }
  175. int main()
  176. {
  177. ios_base::sync_with_stdio(false);
  178. cin.tie(0);
  179. cout.tie(0);
  180. int TESTS=1;
  181. // cin>>TESTS;
  182. while(TESTS--)
  183. {
  184. solve();
  185. }
  186. time
  187. return 0;
  188. }
Advertisement
Add Comment
Please, Sign In to add comment