Description

有一列元素，每一个元素有三个属性:标号、标识符、数值。这些元素按照标号从1~ n排列，标识符也是1~n的一个排列，初始时数值为0。当然我们可以把每个元素看成一个多维数字，那么这列元素就是一个数列。

现在请你维护这个数列，使其能支持以下两种操作：1.将标号为l~ r的所有元素的数值先乘上x，再加上y；2.将标识符为l~ r的所有元素的数值先乘上x，再加上y。当然你还得回答某些询问：1.标号为l~ r的所有元素的数值的和；2.标识符为l~r的所有元素的数值的和。

Input

第一行有两个正整数n、m，分别表示数列长度和操作与询问个数的总和。第二行有n个正整数，表示每个元素的标识符，保证这n个数是1~n的一个排列。接下来m行，每行的第一个数字为op。若op为0，则表示要进行第一个操作，接下去四个数字表示l,r,x,y;若op为1，则表示要进行第二个操作，接下去四个数字表示l,r,x,y;若op为2，则表示要回答第一个询问，接下去两个数字表示l,r;若op为3，则表示要回答第二个询问，接下去两个数字表示l,r。

Output

包含若干行，每行表示一个询问的答案。由于答案可能很大，只要请你输出答案对536870912取模后的值即可。

Sample Input

4 4 2 1 4 3 0 2 3 4 5 1 1 3 4 7 2 1 1 3 1 1

Sample Output

7 27

Solution

设第\(i\)个数的标识符为\(p_i\)，把每个数看做二维平面上的一个点\((i,p_i)\)，然后\(kd\_tree\)暴力维护就好了。

复杂度\(O(m\sqrt{n})\)。

关于卡常。。。

这题你只知道上面这些并不足以通过。。

卡了一个小时的常无果之后，，，

注意到模数是一个奇怪的数，并没有见过。

于是，，百度一下可以发现，模数\(=2^{29}\)。。。

然后就a掉了。。

#pragma GCC optimize(3)#include
    
     using namespace std;#define int unsigned intchar buf[1000000],*p1=buf,*p2=buf;#define getchar() ((p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin)),p1==p2)?EOF:*p1++)void read(int &x) {    x=0;int f=1;char ch=getchar();    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;    for(;isdigit(ch);ch=getchar()) x=x*10+ch-'0';x*=f;}char buf2[1000000],a[100];int p3=-1,P;inline void write(int x) {    do {a[++P]=x%10+48;} while(x/=10);    do {buf2[++p3]=a[P];} while(--P);    buf2[++p3]='\n';}const int maxn = 65535;int n,m,Dem,L,R,ans,D,mul,add;inline void chmin(int &x,int y) {if(y
     
      x) x=y;}struct node {    int add,mul,sum,val,l,r,cnt;    int mn[2],mx[2],p[2];};inline int cmp(node x,node y) {    return x.p[Dem]
      
       >1)int tot,rt;node t[maxn];void up(int x) {    register int l=t[x].l,r=t[x].r;    t[x].sum=(t[x].val+t[l].sum+t[r].sum);    t[x].cnt=t[l].cnt+t[r].cnt+1;    if(l) {        chmin(t[x].mn[0],t[l].mn[0]);        chmin(t[x].mn[1],t[l].mn[1]);        chmax(t[x].mx[0],t[l].mx[0]);        chmax(t[x].mx[1],t[l].mx[1]);    }    if(r) {        chmin(t[x].mn[0],t[r].mn[0]);        chmin(t[x].mn[1],t[r].mn[1]);        chmax(t[x].mx[0],t[r].mx[0]);        chmax(t[x].mx[1],t[r].mx[1]);    }}void update(int x) {t[x].sum=(t[x].val+t[t[x].l].sum+t[t[x].r].sum);}void push_tag(int x,int Add,int Mul) {    t[x].val=(t[x].val*Mul+Add);    t[x].add=(t[x].add*Mul+Add);    t[x].mul=t[x].mul*Mul;    t[x].sum=(t[x].sum*Mul+Add*t[x].cnt);}void pushdown(int x) {    if(t[x].l) push_tag(t[x].l,t[x].add,t[x].mul);    if(t[x].r) push_tag(t[x].r,t[x].add,t[x].mul);    t[x].add=0,t[x].mul=1;}int build(int De,int l,int r) {    Dem=De;nth_element(t+l+1,t+mid+1,t+r+1,cmp);    t[mid].mn[0]=t[mid].mx[0]=t[mid].p[0];    t[mid].mn[1]=t[mid].mx[1]=t[mid].p[1];    if(l!=mid) t[mid].l=build(De^1,l,mid-1);    if(r!=mid) t[mid].r=build(De^1,mid+1,r);    t[mid].cnt=1;up(mid);t[mid].mul=1;return mid;}void modify(int p) {    if(t[p].mn[D]>R||t[p].mx[D]
       
        R||t[p].mx[D]