Problem:
Matrix Operation

You are a student looking for a job.
Today you had an employment examination for an IT company.
They asked you to write an efficient program to perform several operations.
First, they showed you an N \times N square matrix and a list of operations.
All operations but one modify the matrix, and the last operation outputs the character in a specified cell.
Please remember that you need to output the final matrix after you finish all the operations.

Followings are the detail of the operations:

WR r c v

(Write operation) write a integer v into the cell (r,c) (1 \leq v \leq 1,000,000)

CP r1 c1 r2 c2

(Copy operation) copy a character in the cell (r1,c1) into the cell (r2,c2)

SR r1 r2

(Swap Row operation) swap the r1-th row and r2-th row

SC c1 c2

(Swap Column operation) swap the c1-th column and c2-th column

RL

(Rotate Left operation) rotate the whole matrix in counter-clockwise direction by 90 degrees

RR

(Rotate Right operation) rotate the whole matrix in clockwise direction by 90 degrees

RH

(Reflect Horizontal operation) reverse the order of the rows

RV

(Reflect Vertical operation) reverse the order of the columns

Input

First line of each testcase contains nine integers.
First two integers in the line, N and Q, indicate the size of matrix and the number of queries, respectively (1 \leq N,Q \leq 40,000).
Next three integers, A B, and C, are coefficients to calculate values in initial matrix (1 \leq A,B,C \leq 1,000,000),
and they are used as follows: A_{r,c} = (r * A + c * B) mod C where r and c are row and column indices, respectively (1\leq r,c\leq N).
Last four integers, D, E, F, and G, are coefficients to compute the final hash value mentioned in the next section (1 \leq D \leq E \leq N, 1 \leq F \leq G \leq N, E - D \leq 1,000, G - F \leq 1,000).
Each of next Q lines contains one operation in the format as described above.

Output

Output a hash value h computed from the final matrix B by using following pseudo source code.

h <- 314159265
for r = D...E
for c = F...G
h <- (31 * h + B_{r,c}) mod 1,000,000,007

where "<-" is a destructive assignment operator, "for i = S...T" indicates a loop for i from S to T (both inclusive), and "mod" is a remainder operation.

Sample Input 1

2 1 3 6 12 1 2 1 2
WR 1 1 1

Output for the Sample Input 1

676573821

Sample Input 2

2 1 3 6 12 1 2 1 2
RL

Output for the Sample Input 2

676636559

Sample Input 3

2 1 3 6 12 1 2 1 2
RH

Output for the Sample Input 3

676547189

Sample Input 4

39989 6 999983 999979 999961 1 1000 1 1000
SR 1 39989
SC 1 39989
RL
RH
RR
RV

Output for the Sample Input 4

458797120

Solution in C++ (judge verdict: Accepted)

#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
#include <cmath>
#include <algorithm>
#include <functional>
#include <map>
using namespace std;
#define rep(i,n) for(int i=0;i<(n);++i)
#define rep1(i,n) for(int i=1;i<=(n);++i)
#define repa(i,a,n) for(int i=a;i<a+n;++i)
#define repb(i,a,n) for(int i=a;i>a-n;--i)
#define all(c) (c).begin(),(c).end()
#define pb push_back
#define fs first
#define sc second
#define show(x) cout << #x << " = " << x << endl
typedef long long ll;
typedef pair<int,int> P;
map<P,ll> mp;
ll N,Q,A,B,C,D,E,F,G;
ll mod=1e9+7;
int row[40001],col[40001];
ll initial(int r,int c){
	return (A*r+B*c)%C;
}
ll value(int r,int c){
	if(mp.find(P(r,c))==mp.end()){
		return initial(r,c);
	}else{
		return mp[P(r,c)];
	}
}
int ro,inv;
void change(int &r,int &c){
	int nr,nc;
	if(inv==0){
		if(ro==0) nr=r,nc=c;
		if(ro==1) nr=c,nc=N+1-r;
		if(ro==2) nr=N+1-r,nc=N+1-c;
		if(ro==3) nr=N+1-c,nc=r;
	}else{
		if(ro==0) nr=N+1-r,nc=c;
		if(ro==1) nr=N+1-c,nc=N+1-r;
		if(ro==2) nr=r,nc=N+1-c;
		if(ro==3) nr=c,nc=r;
	}
	r=row[nr],c=col[nc];
}
void change2(int &r,int &c){
	int nr,nc;
	if(inv==0){
		if(ro==0) nr=r,nc=c;
		if(ro==1) nr=c,nc=N+1-r;
		if(ro==2) nr=N+1-r,nc=N+1-c;
		if(ro==3) nr=N+1-c,nc=r;
	}else{
		if(ro==0) nr=N+1-r,nc=c;
		if(ro==1) nr=N+1-c,nc=N+1-r;
		if(ro==2) nr=r,nc=N+1-c;
		if(ro==3) nr=c,nc=r;
	}
	r=nr,c=nc;
}
void showmat(){
	show(ro);
	show(inv);
	for(int r=D;r<=E;r++){
		for(int c=F;c<=G;c++){
			int tmpr=r,tmpc=c;
			change(tmpr,tmpc);
			cout<<value(tmpr,tmpc)<<" ";
		}
		cout<<endl;
	}
	cout<<endl;
}
int main(){
	cin>>N>>Q>>A>>B>>C>>D>>E>>F>>G;
	rep1(i,N) row[i]=i,col[i]=i;
	int X=E-D+1,Y=G-F+1;
	ro=0,inv=0;
	rep(i,Q){
		string type;
		cin>>type;
		if(type=="WR"){
			int r,c,v;
			cin>>r>>c>>v;
			change(r,c);
			mp[P(r,c)]=v;
		}
		if(type=="CP"){
			int r1,c1,r2,c2;
			cin>>r1>>c1>>r2>>c2;
			change(r1,c1);
			change(r2,c2);
			mp[P(r2,c2)]=value(r1,c1);
		}
		if(type=="SR"){
			int r1,r2,c1=0,c2=0;
			cin>>r1>>r2;
			change2(r1,c1);
			change2(r2,c2);
			// show(r1);
			// show(r2);
			// show(c1);
			// show(c2);
			if(c1==c2) swap(row[r1],row[r2]);
			if(r1==r2) swap(col[c1],col[c2]);
		}
		if(type=="SC"){
			int r1=0,r2=0,c1,c2;
			cin>>c1>>c2;
			change2(r1,c1);
			change2(r2,c2);
			if(c1==c2) swap(row[r1],row[r2]);
			if(r1==r2) swap(col[c1],col[c2]);
		}
		if(type=="RL"){
			ro=(ro+1)%4;
		}
		if(type=="RR"){
			ro=(ro+3)%4;
		}
		if(type=="RH"){
			if(ro%2==0){
				inv=1-inv;
			}else{
				inv=1-inv;
				ro=(ro+2)%4;
			}
		}
		if(type=="RV"){
			if(ro%2==0){
				inv=1-inv;
				ro=(ro+2)%4;
			}else{
				inv=1-inv;
			}
		}
//		showmat();
	}
	ll h=314159265;
	for(int r=D;r<=E;r++) for(int c=F;c<=G;c++){
		int tmpr=r,tmpc=c;
		change(tmpr,tmpc);
		h=h*31+value(tmpr,tmpc);
		h%=mod;
	}
	cout<<h<<endl;
}