Transformation

http://acm.hdu.edu.cn/showproblem.php?pid=4578

题意：

//长度为n的数组四个操作
//1 x y c [x,y]区间的数都加c
//2 x y c [x, y] 区间的数都乘以c
//3 x y c [x ,y] 区间的数都变为c
//4 x y p [x ,y] 求区间的数的p次方的和

题解：

这道题坑在有三种询问：set , add , mul。所以lazy标记要有三个，如果三个标记同时出现的处理方法——当更新set操作时，就把add标记和mul标记全部取消；当更新mul操作时，如果当前节点add标记存在，就把add标记改为：add * mul。这样的话就可以在PushDown()操作中先执行set，然后mul，最后add。

　　麻烦在有三种询问：和 , 平方和 , 立方和。对于set和mul操作来说，这三种询问都比较好弄。

　　对于add操作，和的话就比较好弄，按照正常方法就可以；

　　平方和这样来推：(a + c)2 = a2 + c2 + 2ac , 即sum2[rt] = sum2[rt] + (r - l + 1) * c * c + 2 * sum1[rt] * c；

　　立方和这样推：(a + c)3 = a3 + c3 + 3c(a2 + ac) , 即sum3[rt] = sum3[rt] + (r - l + 1) * c * c * c + 3 * c * (sum2[rt] + sum1[rt] * c)；

　　几个注意点：

add标记取消的时候是置0，mul标记取消的时候是置1；
在PushDown()中也也要注意取消标记，如set操作中取消add和mul，mul操作中更新add；
在add操作中要注意sum3 , sum2 , sum1的先后顺序，一定是先sum3 , 然后sum2 , 最后sum1；
int容易爆，还是用LL要保险一点；
最后就是运算较多，不要漏掉东西。

C++版本一

/*
*@Author:   STZG
*@Language: C++
*/
#include <bits/stdc++.h>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<string>
#include<vector>
#include<bitset>
#include<queue>
#include<deque>
#include<stack>
#include<cmath>
#include<list>
#include<map>
#include<set>
//#define DEBUG
#define RI register int
using namespace std;
typedef long long ll;
//typedef __int128 lll;
const int N=100000+10;
const int M=100000+10;
const int MOD=10000+7;
const double PI = acos(-1.0);
const double EXP = 1E-8;
const int INF = 0x3f3f3f3f;
int t,n,m,k,q,p;
int ans,cnt,flag,temp,sum;
ll tree[N<<2];
ll tree2[N<<2];
ll tree3[N<<2];
ll lazy1[N<<2];
ll lazy2[N<<2];
ll lazy3[N<<2];
char str;
struct node{};
ll POW(ll a,int b){if (b==1) return a%MOD;if (b==2) return a*a%MOD;if (b==3) return a*a%MOD*a%MOD;return 1;
}
void pushup(int rt){tree[rt]=(tree[rt<<1]+tree[rt<<1|1])%MOD;tree2[rt]=(tree2[rt<<1]+tree2[rt<<1|1])%MOD;tree3[rt]=(tree3[rt<<1]+tree3[rt<<1|1])%MOD;
}
void pushdown(int l,int r,int rt){if(lazy3[rt]){int mid=(l+r)>>1;int llen=(mid-l+1);int rlen=(r-mid);tree[rt<<1]=(llen*lazy3[rt])%MOD;tree[rt<<1|1]=(rlen*lazy3[rt])%MOD;tree2[rt<<1]=(llen*POW(lazy3[rt],2))%MOD;tree2[rt<<1|1]=(rlen*POW(lazy3[rt],2))%MOD;tree3[rt<<1]=(llen*POW(lazy3[rt],3))%MOD;tree3[rt<<1|1]=(rlen*POW(lazy3[rt],3))%MOD;lazy3[rt<<1|1]=lazy3[rt<<1]=lazy3[rt];lazy3[rt]=0;lazy1[rt<<1|1]=lazy1[rt<<1]=1;lazy2[rt<<1|1]=lazy2[rt<<1]=0;}if(lazy1[rt]!=1){tree[rt<<1]=(tree[rt<<1]*lazy1[rt])%MOD;tree[rt<<1|1]=(tree[rt<<1|1]*lazy1[rt])%MOD;tree2[rt<<1]=(tree2[rt<<1]*POW(lazy1[rt],2))%MOD;tree2[rt<<1|1]=(tree2[rt<<1|1]*POW(lazy1[rt],2))%MOD;tree3[rt<<1]=(tree3[rt<<1]*POW(lazy1[rt],3))%MOD;tree3[rt<<1|1]=(tree3[rt<<1|1]*POW(lazy1[rt],3))%MOD;lazy2[rt<<1]=(lazy2[rt<<1]*lazy1[rt])%MOD;lazy2[rt<<1|1]=(lazy2[rt<<1|1]*lazy1[rt])%MOD;lazy1[rt<<1]=(lazy1[rt<<1]*lazy1[rt])%MOD;lazy1[rt<<1|1]=(lazy1[rt<<1|1]*lazy1[rt])%MOD;lazy1[rt]=1;}if(lazy2[rt]){int mid=(l+r)>>1;int llen=(mid-l+1);int rlen=(r-mid);tree3[rt<<1]=(tree3[rt<<1]+(llen*POW(lazy2[rt],3))%MOD+3*lazy2[rt]*tree2[rt<<1]%MOD+3*POW(lazy2[rt],2)*tree[rt<<1]%MOD)%MOD;tree3[rt<<1|1]=(tree3[rt<<1|1]+(rlen*POW(lazy2[rt],3))%MOD+3*lazy2[rt]*tree2[rt<<1|1]%MOD+3*POW(lazy2[rt],2)*tree[rt<<1|1]%MOD)%MOD;tree2[rt<<1]=(tree2[rt<<1]+(llen*POW(lazy2[rt],2))%MOD+2*tree[rt<<1]*lazy2[rt]%MOD)%MOD;tree2[rt<<1|1]=(tree2[rt<<1|1]+(rlen*POW(lazy2[rt],2))%MOD+2*tree[rt<<1|1]*lazy2[rt]%MOD)%MOD;tree[rt<<1]=(tree[rt<<1]+(llen*lazy2[rt])%MOD)%MOD;tree[rt<<1|1]=(tree[rt<<1|1]+(rlen*lazy2[rt])%MOD)%MOD;lazy2[rt<<1]=(lazy2[rt<<1]+lazy2[rt])%MOD;lazy2[rt<<1|1]=(lazy2[rt<<1|1]+lazy2[rt])%MOD;lazy2[rt]=0;}
}
void build(int l,int r,int rt){lazy1[rt]=1;lazy2[rt]=0;lazy3[rt]=0;if(l==r){//scanf("%lld",&tree[rt]);tree[rt]=0;tree2[rt]=0;tree3[rt]=0;return;}int mid=(l+r)>>1;build(l,mid,rt<<1);build(mid+1,r,rt<<1|1);pushup(rt);
}
void updata(int l,int r,int rt,int L,int R,ll C,int op){if(L<=l&&r<=R){int len=r-l+1;if(op==1){tree[rt]=(tree[rt]*C)%MOD;tree2[rt]=(tree2[rt]*POW(C,2))%MOD;tree3[rt]=(tree3[rt]*POW(C,3))%MOD;lazy2[rt]=(lazy2[rt]*C)%MOD;lazy1[rt]=(lazy1[rt]*C)%MOD;}else if(op==2){tree3[rt]=(tree3[rt]+(len*POW(C,3))%MOD+3*C*tree2[rt]%MOD+3*POW(C,2)*tree[rt]%MOD)%MOD;tree2[rt]=(tree2[rt]+(len*POW(C,2))%MOD+2*tree[rt]*C%MOD)%MOD;tree[rt]=(tree[rt]+len*C)%MOD;lazy2[rt]=(lazy2[rt]+C)%MOD;}else{tree[rt]=(len*C)%MOD;tree2[rt]=(len*POW(C,2))%MOD;tree3[rt]=(len*POW(C,3))%MOD;lazy1[rt]=1;lazy2[rt]=0;lazy3[rt]=C;}return;}int mid=(l+r)>>1;pushdown(l,r,rt);if(L<=mid)updata(l,mid,rt<<1,L,R,C,op);if(mid<R)updata(mid+1,r,rt<<1|1,L,R,C,op);pushup(rt);
}
ll query(int l,int r,int rt,int L,int R,int op){if(L<=l&&r<=R){if(op==1)return tree[rt];if(op==2)return tree2[rt];if(op==3)return tree3[rt];}int mid=(l+r)>>1;ll res=0;pushdown(l,r,rt);if(L<=mid)res=(res+query(l,mid,rt<<1,L,R,op))%MOD;if(mid<R)res=(res+query(mid+1,r,rt<<1|1,L,R,op))%MOD;return res;
}
void see(int l,int r,int rt){if(l==r){printf("%lld%c",tree3[rt]," \n"[r==n]);return;}int mid=(l+r)>>1;pushdown(l,r,rt);see(l,mid,rt<<1);see(mid+1,r,rt<<1|1);
}
int main()
{
#ifdef DEBUGfreopen("input.in", "r", stdin);//freopen("output.out", "w", stdout);
#endif//ios::sync_with_stdio(false);//cin.tie(0);//cout.tie(0);//scanf("%d",&t);//while(t--){while(~scanf("%d%d",&n,&m)&&n+m){build(1,n,1);int op;int l,r;ll c;for(int i=1;i<=m;i++){scanf("%d",&op);if(op==1){scanf("%d%d%lld",&l,&r,&c);updata(1,n,1,l,r,c,2);}else if(op==2){scanf("%d%d%lld",&l,&r,&c);updata(1,n,1,l,r,c,1);}else if(op==3){scanf("%d%d%lld",&l,&r,&c);updata(1,n,1,l,r,c,3);}else if(op==4){scanf("%d%d%d",&l,&r,&p);printf("%lld\n",query(1,n,1,l,r,p));}//see(1,n,1);}}//}#ifdef DEBUGprintf("Time cost : %lf s\n",(double)clock()/CLOCKS_PER_SEC);
#endif//cout << "Hello world!" << endl;return 0;
}

C++版本二

题解：

//用线段树维护里面的值都相等的区间 ,那么这个区间的所需答案为(r-l+1)*(tree[v].eq)^q
//对于懒惰操作 mul , eq , add的处理是
//对于每次eq操作可以直接操作tree[v].eq = eq ;tree[v].mul = 1 ; tree[v].add = 0;
//而对于mul ,和add的操作每次push_down(v)的时候保证下面一层的所有懒惰情况都是初始情况
//所以对于mul和add操作时就一直往下push_down知道找到一个可以

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn=100010;
#define left v<<1
#define right v<<1|1
#define mod 10007
struct node{int l,r,value;int eq,add,mul;
}tree[maxn<<2];
void build(int l,int r,int v){tree[v].l=l;tree[v].r=r;tree[v].add=0;tree[v].mul=1;tree[v].eq=-1;if(l==r){tree[v].eq=0;return;}int mid=(l+r)>>1;build(l,mid,left);build(mid+1,r,right);
}
void push_down(int v){if(tree[v].l==tree[v].r)return;if(tree[v].eq!=-1){//=tree[left].eq=tree[right].eq=tree[v].eq;tree[left].add=tree[right].add=0;tree[left].mul=tree[right].mul=1;tree[v].eq=-1;return;}if(tree[v].mul!=1){//*if(tree[left].eq!=-1)tree[left].eq=(tree[left].eq*tree[v].mul)%mod;else{push_down(left);tree[left].mul=(tree[left].mul*tree[v].mul)%mod;}if(tree[right].eq!=-1)tree[right].eq=(tree[right].eq*tree[v].mul)%mod;else{push_down(right);tree[right].mul=(tree[right].mul*tree[v].mul)%mod;}tree[v].mul=1;}if(tree[v].add){//+if(tree[left].eq!=-1)tree[left].eq=(tree[left].eq+tree[v].add)%mod;else{push_down(left);tree[left].add=(tree[left].add+tree[v].add)%mod;}if(tree[right].eq!=-1)tree[right].eq=(tree[right].eq+tree[v].add)%mod;else{push_down(right);tree[right].add=(tree[right].add+tree[v].add)%mod;}tree[v].add=0;}
}
void update(int l,int r,int v,int op,int c){if(l<=tree[v].l&&tree[v].r<=r){if(op==3){tree[v].add=0;tree[v].mul=1;tree[v].eq=c;return;}if(tree[v].eq!=-1){if(op==1)tree[v].eq=(tree[v].eq+c)%mod;else tree[v].eq=(tree[v].eq*c)%mod;}else{push_down(v);if(op==1)tree[v].add=(tree[v].add+c)%mod;else tree[v].mul=(tree[v].mul*c)%mod;}return;}push_down(v);int mid=(tree[v].l+tree[v].r)>>1;if(l<=mid)update(l,r,left,op,c);if(r>mid)update(l,r,right,op,c);
}
int query(int l,int r,int v,int q){if(tree[v].l>=l&&tree[v].r<=r&&tree[v].eq!=-1){int ans=1;for(int i=1;i<=q;i++)ans=(ans*tree[v].eq)%mod;return(ans*((tree[v].r-tree[v].l+1)%mod))%mod;}push_down(v);int mid=(tree[v].l+tree[v].r)>>1;if(l>mid)return query(l,r,right,q);else if(r<=mid)return query(l,r,left,q);elsereturn(query(l,mid,left,q)+query(mid+1,r,right,q))%mod;
}
int main()
{int n,m;while(scanf("%d%d",&n,&m)&&(n+m)){int op,x,y,c;build(1,n,1);while(m--){scanf("%d%d%d%d",&op,&x,&y,&c);if(op==4)printf("%d\n",(query(x,y,1,c)%mod));elseupdate(x,y,1,op,c);}}
return 0;
}

C++版本三

///先乘后加的原则，每次更新区间乘法，加法标记就*乘法标记#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long
#define lson rt<<1,l,m
#define rson rt<<1|1,m+1,r
#define mid int m=l+r>>1
#define debug(x) printf("----Line%s----\n",#x)
#define ls rt<<1
#define rs rt<<1|1
#define mod 10007
using namespace std;const int N = 1e5+5;ll add[N<<2],mtag[N<<2];
ll t[N<<2],tt[N<<2],ttt[N<<2];
ll isset[N<<2];ll Pow(ll a,int b){if (b==2) return a*a%mod;if (b==3) return a*a%mod*a%mod;
}void push_up(int rt)
{t[rt] = (t[ls] + t[rs])%mod;tt[rt] = (tt[ls] + tt[rs])%mod;ttt[rt] = (ttt[ls] + ttt[rs])%mod;
}void push_down(int rt,int len)
{if (isset[rt]){ll x = isset[rt];isset[ls] = isset[rs] = isset[rt];isset[rt] = 0;mtag[ls] = mtag[rs] = 1;add[ls] = add[rs] = 0;ttt[ls] = (len-len/2)*Pow(x,3)%mod;tt[ls] = (len-len/2)*Pow(x,2)%mod;t[ls] = (len-len/2)*x%mod;ttt[rs] = (len/2)*Pow(x,3)%mod;tt[rs] = (len/2)*Pow(x,2)%mod;t[rs] = (len/2)*x%mod;}if (mtag[rt]==1 && add[rt]==0) return ;ttt[ls] = (ttt[ls]*Pow(mtag[rt],3)%mod + (len-len/2)*Pow(add[rt],3)%mod + 3*tt[ls]*Pow(mtag[rt],2)%mod*add[rt]%mod + 3*t[ls]*mtag[rt]%mod*Pow(add[rt],2)%mod)%mod;ttt[rs] = (ttt[rs]*Pow(mtag[rt],3)%mod +  len/2     *Pow(add[rt],3)%mod + 3*tt[rs]*Pow(mtag[rt],2)%mod*add[rt]%mod + 3*t[rs]*mtag[rt]%mod*Pow(add[rt],2)%mod)%mod;tt[ls] = (tt[ls]*Pow(mtag[rt],2)%mod + (len-len/2)*Pow(add[rt],2)%mod + 2*t[ls]*add[rt]*mtag[rt])%mod;tt[rs] = (tt[rs]*Pow(mtag[rt],2)%mod + len/2*Pow(add[rt],2)%mod       + 2*t[rs]*add[rt]*mtag[rt])%mod;t[ls] = (t[ls]*mtag[rt] + (len-len/2)*add[rt])%mod;t[rs] = (t[rs]*mtag[rt] + (len/2)*add[rt])%mod;mtag[ls] = (mtag[ls]*mtag[rt])%mod;mtag[rs] = (mtag[rs]*mtag[rt])%mod;add[ls] = (add[ls]*mtag[rt]%mod+add[rt])%mod;///这个标记也要先乘后加add[rs] = (add[rs]*mtag[rt]%mod+add[rt])%mod;add[rt] = 0;mtag[rt] = 1;
}void update(int L,int R,int op,ll x,int rt,int l,int r)
{if (L<=l && r<=R){if (op==1){///+add[rt] = (add[rt]+x)%mod;ttt[rt] = (ttt[rt] + Pow(x,3)*(r-l+1)%mod + 3*tt[rt]*x%mod + 3*t[rt]*Pow(x,2)%mod)%mod;tt[rt] = (tt[rt] + Pow(x,2)*(r-l+1) + 2*t[rt]*x)%mod;t[rt] = (t[rt] + (r-l+1)*x)%mod;}if (op==2){///*mtag[rt] = (mtag[rt]*x)%mod;add[rt] = (add[rt]*x)%mod; ttt[rt] = (ttt[rt]*Pow(x,3))%mod;tt[rt] = (tt[rt]*Pow(x,2))%mod;t[rt] = (t[rt]*x)%mod;}if (op==3){///覆盖mtag[rt] = 1,add[rt] = 0;ttt[rt] = (r-l+1)*Pow(x,3)%mod;tt[rt] = (r-l+1)*Pow(x,2)%mod;t[rt] = (r-l+1)*x%mod;isset[rt] = x;}return ;}mid;push_down(rt,r-l+1);if (L<=m) update(L,R,op,x,lson);if (R>m)  update(L,R,op,x,rson);push_up(rt);
}ll query(int L,int R,ll x,int rt,int l,int r)
{if (L<=l && r<=R){if (x==1) return t[rt];if (x==2) return tt[rt];if (x==3) return ttt[rt];}push_down(rt,r-l+1);ll ans = 0;mid;if (L<=m) ans = (ans+query(L,R,x,lson))%mod;if (R>m)  ans = (ans+query(L,R,x,rson))%mod;return ans;
}void build(int rt,int l,int r)
{if (l==r){t[rt] = tt[rt] = ttt[rt] = add[rt] = isset[rt] = 0;mtag[rt] = 1;return ;}mid;build(lson);build(rson);push_up(rt);
}int main()
{int n,q;while (~scanf("%d %d",&n,&q),n||q){build(1,1,n);int op,l,r;ll x;while (q--){scanf("%d %d %d %lld",&op,&l,&r,&x);if (op!=4){update(l,r,op,x,1,1,n);}elseprintf("%lld\n",query(l,r,x,1,1,n)%mod);}}return 0;
}

C++版本四

原博客

#include <iostream>
#include <cstdio>
#include <vector>
#include <cmath>
#include <string>
#include <string.h>
#include <algorithm>
using namespace std;
#define LL __int64
#define eps 1e-8
#define INF INT_MAX
#define lson l , m , rt << 1
#define rson m + 1 , r , rt << 1 | 1
const int MOD = 10007;
const int maxn = 100000 + 5;
const int N = 12;
LL add[maxn << 2] , set[maxn << 2] , mul[maxn << 2];
LL sum1[maxn << 2] , sum2[maxn << 2] , sum3[maxn << 2];
void PushUp(int rt)
{sum1[rt] = (sum1[rt << 1] + sum1[rt << 1 | 1]) % MOD;sum2[rt] = (sum2[rt << 1] + sum2[rt << 1 | 1]) % MOD;sum3[rt] = (sum3[rt << 1] + sum3[rt << 1 | 1]) % MOD;
}
void build(int l , int r , int rt)
{add[rt] = set[rt] = 0;mul[rt] = 1;if(l == r) {sum1[rt] = sum2[rt] = sum3[rt] = 0;return;}int m = (l + r) >> 1;build(lson);build(rson);PushUp(rt);
}
void PushDown(int rt , int len)
{if(set[rt]) {set[rt << 1] = set[rt << 1 | 1] = set[rt];add[rt << 1] = add[rt << 1 | 1] = 0;    //注意这个也要下放mul[rt << 1] = mul[rt << 1 | 1] = 1;LL tmp = ((set[rt] * set[rt]) % MOD) * set[rt] % MOD;sum1[rt << 1] = ((len - (len >> 1)) % MOD) * (set[rt] % MOD) % MOD;sum1[rt << 1 | 1] = ((len >> 1) % MOD) * (set[rt] % MOD) % MOD;sum2[rt << 1] = ((len - (len >> 1)) % MOD) * ((set[rt] * set[rt]) % MOD) % MOD;sum2[rt << 1 | 1] = ((len >> 1) % MOD) * ((set[rt] * set[rt]) % MOD) % MOD;sum3[rt << 1] = ((len - (len >> 1)) % MOD) * tmp % MOD;sum3[rt << 1 | 1] = ((len >> 1) % MOD) * tmp % MOD;set[rt] = 0;}if(mul[rt] != 1) {    //这个就是mul[rt] != 1 ， 当时我这里没注意所以TLE了mul[rt << 1] = (mul[rt << 1] * mul[rt]) % MOD;mul[rt << 1 | 1] = (mul[rt << 1 | 1] * mul[rt]) % MOD;if(add[rt << 1])    //注意这个也要下放add[rt << 1] = (add[rt << 1] * mul[rt]) % MOD;if(add[rt << 1 | 1])add[rt << 1 | 1] = (add[rt << 1 | 1] * mul[rt]) % MOD;LL tmp = (((mul[rt] * mul[rt]) % MOD * mul[rt]) % MOD);sum1[rt << 1] = (sum1[rt << 1] * mul[rt]) % MOD;sum1[rt << 1 | 1] = (sum1[rt << 1 | 1] * mul[rt]) % MOD;sum2[rt << 1] = (sum2[rt << 1] % MOD) * ((mul[rt] * mul[rt]) % MOD) % MOD;sum2[rt << 1 | 1] = (sum2[rt << 1 | 1] % MOD) * ((mul[rt] * mul[rt]) % MOD) % MOD;sum3[rt << 1] = (sum3[rt << 1] % MOD) * tmp % MOD;sum3[rt << 1 | 1] = (sum3[rt << 1 | 1] % MOD) * tmp % MOD;mul[rt] = 1;}if(add[rt]) {add[rt << 1] += add[rt];    //add是+= ， mul是*=add[rt << 1 | 1] += add[rt];LL tmp = (add[rt] * add[rt] % MOD) * add[rt] % MOD;        //注意sum3 , sum2 , sum1的先后顺序sum3[rt << 1] = (sum3[rt << 1] + (tmp * (len - (len >> 1)) % MOD) + 3 * add[rt] * ((sum2[rt << 1] + sum1[rt << 1] * add[rt]) % MOD)) % MOD;sum3[rt << 1 | 1] = (sum3[rt << 1 | 1] + (tmp * (len >> 1) % MOD) + 3 * add[rt] * ((sum2[rt << 1 | 1] + sum1[rt << 1 | 1] * add[rt]) % MOD)) % MOD;sum2[rt << 1] = (sum2[rt << 1] + ((add[rt] * add[rt] % MOD) * (len - (len >> 1)) % MOD) + (2 * sum1[rt << 1] * add[rt] % MOD)) % MOD;sum2[rt << 1 | 1] = (sum2[rt << 1 | 1] + (((add[rt] * add[rt] % MOD) * (len >> 1)) % MOD) + (2 * sum1[rt << 1 | 1] * add[rt] % MOD)) % MOD;sum1[rt << 1] = (sum1[rt << 1] + (len - (len >> 1)) * add[rt]) % MOD;sum1[rt << 1 | 1] = (sum1[rt << 1 | 1] + (len >> 1) * add[rt]) % MOD;add[rt] = 0;}
}
void update(int L , int R , int c , int ch , int l , int r , int rt)
{if(L <= l && R >= r) {if(ch == 3) {set[rt] = c;add[rt] = 0;mul[rt] = 1;sum1[rt] = ((r - l + 1) * c) % MOD;sum2[rt] = ((r - l + 1) * ((c * c) % MOD)) % MOD;sum3[rt] = ((r - l + 1) * (((c * c) % MOD) * c % MOD)) % MOD;} else if(ch == 2) {mul[rt] = (mul[rt] * c) % MOD;if(add[rt])add[rt] = (add[rt] * c) % MOD;sum1[rt] = (sum1[rt] * c) % MOD;sum2[rt] = (sum2[rt] * (c * c % MOD)) % MOD;sum3[rt] = (sum3[rt] * ((c * c % MOD) * c % MOD)) % MOD;} else if(ch == 1) {add[rt] += c;LL tmp = (((c * c) % MOD * c) % MOD * (r - l + 1)) % MOD;    //(r - l + 1) * c^3sum3[rt] = (sum3[rt] + tmp + 3 * c * ((sum2[rt] + sum1[rt] * c) % MOD)) % MOD;sum2[rt] = (sum2[rt] + (c * c % MOD * (r - l + 1) % MOD) + 2 * sum1[rt] * c) % MOD;sum1[rt] = (sum1[rt] + (r - l + 1) * c) % MOD;}return;}PushDown(rt , r - l + 1);int m = (l + r) >> 1;if(L > m)update(L , R , c , ch , rson);else if(R <= m)update(L , R , c , ch , lson);else {update(L , R , c , ch , lson);update(L , R , c , ch , rson);}PushUp(rt);
}
LL query(int L , int R , int p , int l , int r , int rt)
{if(L <= l && R >= r) {if(p == 1)return sum1[rt] % MOD;else if(p == 2)return sum2[rt] % MOD;elsereturn sum3[rt] % MOD;}PushDown(rt , r - l + 1);int m = (l + r) >> 1;if(L > m)return query(L , R , p , rson);else if(R <= m)return query(L , R , p , lson);else return (query(L , R , p , lson) + query(L , R , p , rson)) % MOD;
}
int main()
{int n , m;int a , b , c , ch;while(~scanf("%d %d" , &n , &m)){if(n == 0 && m == 0)break;build(1 , n , 1);while(m--) {scanf("%d %d %d %d" , &ch , &a , &b , &c);if(ch != 4) {update(a , b , c , ch , 1 , n , 1);} else {printf("%I64d\n" , query(a , b , c , 1 , n , 1));}}}return 0;
}

C++版本五

题解：

我只能说这题很迷，迷之WA，这题我陆陆续续做了4天才AC，其中WA了3页，共计40多次吧，重写代码3次，用各种不同方法写。。最后无奈看着别人的博客照着一个大佬的写就AC了。。吐血，一开始我打算用3个v，3个tag来写，3个v表示区间的一二三次方值，tag分别表示加，乘等于的标记，先处理等，如果有等号标记，区间更新，清除加号和乘号标记，在处理乘和加之前看子区间有没有等号标记，有就先处理子区间的等号标记，然后处理乘，如果同时有加号标记把加的值乘上一个要乘的数，最后处理乘号标记，至于一次方二次方三次方公式是搞数学的老刘推的。。结果就是我样例过了，自己测的几组也过了还是WA，各种迷，怎么修改都是WA，后来看到一个大佬的没有推公式也没有储存三次方。。方法很神奇，而且速度比一般用储存3个值的那种快得多。。

说一下思路：

同样tag1表示加号标记，tag2表示乘，tag3表示等，无疑在向下更新的时候等号的优先级应该最高，先处理等号标记，并把加号和乘号标记置为0和1（初始值），然后重点来了，在处理乘号标记的时候，如果发现子区间有等号标记，直接修改等号标记，否则先将子区间pushdwon一下清空标记，再进行该子区间标记，加号同理，如果子区间有等号标记，就修改等号标记，否则将子区间pushdown清空标记，再将该子区间标记，所有操作完都要清空当前区间标记，然后修改的时候也差不多，先检查该区间有没有等号标记，有的话直接修改等号标记，否则pushdown清空该区间标记，再进行各种标记，查询的时候是整个算法的精髓，他是查询到要查询区间的子区间的等号标记不为初始标记，即为一段相同的数字的时候停止，然后返回该区间的值，进行各个子区间累加，得出答案。。。。一开始我还以为会超时，结果居然比一般的算法快得多，膜拜大佬

#include<algorithm>
#include<iostream>
#include<cstring>
#include<stdio.h>
#include<math.h>
#include<string>
#include<stdio.h>
#include<queue>
#include<stack>
#include<map>
#include<deque>
using namespace std;
#define left k*2
#define right k*2+1
const int N=10007;
struct node
{int l,r;int tag1,tag2,tag3;//分别表示加号，乘号，等号标记
}t[100005*4];
int n;
void Build(int l,int r,int k)
{t[k].l=l;t[k].r=r;t[k].tag1=0;t[k].tag2=1;t[k].tag3=-1;if(l==r){t[k].tag3=0;//最底层要赋值为0return;};int mid=(l+r)/2;Build(l,mid,left);Build(mid+1,r,right);
}
void pushdown(int k)
{if(t[k].l==t[k].r)//没有子区间了不用退了return;if(t[k].tag3!=-1)//处理等号{t[left].tag3=t[right].tag3=t[k].tag3;//更新子区间等号标记t[left].tag2=t[right].tag2=1;t[left].tag1=t[right].tag1=0;//清空子区间加乘标记t[k].tag3=-1;return;}if(t[k].tag2!=1)//处理乘号{if(t[left].tag3!=-1)//如果子区间有等号标记，直接修改等号标记t[left].tag3=(t[left].tag3*t[k].tag2)%N;else//否则清空该子区间标记，进行子区间标记{pushdown(left);t[left].tag2=(t[left].tag2*t[k].tag2)%N;}if(t[right].tag3!=-1)//同理处理右区间t[right].tag3=(t[right].tag3*t[k].tag2)%N;else{pushdown(right);t[right].tag2=(t[right].tag2*t[k].tag2)%N;}t[k].tag2=1;}if(t[k].tag1!=0)//处理加号标记，和上面同理处理{if(t[left].tag3!=-1)t[left].tag3=(t[left].tag3+t[k].tag1)%N;else{pushdown(left);t[left].tag1=(t[left].tag1+t[k].tag1)%N;}if(t[right].tag3!=-1)t[right].tag3=(t[right].tag3+t[k].tag1)%N;else{pushdown(right);t[right].tag1=(t[right].tag1+t[k].tag1)%N;}t[k].tag1=0;//记得还原}
}
void update(int l,int r,int v,int d,int k)
{if(t[k].l==l&&t[k].r==r){if(d==1){if(t[k].tag3!=-1)//如果有等号标记，就直接修改等号标记{t[k].tag3=(t[k].tag3+v%N)%N;}else{pushdown(k);//否则清空该区间，进行标记t[k].tag1=(t[k].tag1+v%N)%N;}}else if(d==2)//同理{if(t[k].tag3!=-1){t[k].tag3=(t[k].tag3*v%N)%N;}else{pushdown(k);t[k].tag2=(t[k].tag2*v%N)%N;}}else{t[k].tag3=v%N;t[k].tag1=0;t[k].tag2=1;}return;}pushdown(k);//向下更新int mid=(t[k].l+t[k].r)/2;if(r<=mid)update(l,r,v,d,left);else if(l>mid)update(l,r,v,d,right);else{update(l,mid,v,d,left);update(mid+1,r,v,d,right);}
}
int query(int l,int r,int p,int k)//查询
{if(t[k].l>=l&&t[k].r<=r&&t[k].tag3!=-1)//查到是查询区间的子区间且一段全为相同的数{int temp=1;for(int i=1;i<=p;i++)temp=(temp*t[k].tag3)%N;return ((t[k].r-t[k].l+1)%N*temp)%N;//注意要乘上长度}pushdown(k);int mid=(t[k].l+t[k].r)/2;if(r<=mid)return query(l,r,p,left)%N;else if(l>mid)return query(l,r,p,right)%N;else{return (query(l,mid,p,left)+query(mid+1,r,p,right))%N;}
}
int main()
{int i,j,k,m,d,x,y,c;while(scanf("%d%d",&n,&m)!=EOF){if(n==0&&m==0)break;Build(1,n,1);for(i=0;i<m;i++){scanf("%d%d%d%d",&d,&x,&y,&c);if(d>=1&&d<=3){update(x,y,c,d,1);}else{printf("%d\n",query(x,y,c,1)%N);}}}return 0;
}

Transformation相关推荐

2021年大数据Flink（十二）：流批一体API Transformation
目录 Transformation 官网API列表基本操作-略 map flatMap keyBy filter sum reduce 代码演示合并-拆分 union和connect split. ...
R语言使用caret包的preProcess函数进行数据预处理：对所有的数据列进行expoTrans指数变换（Exponential transformation、可以处理负数）
R语言使用caret包的preProcess函数进行数据预处理:对所有的数据列进行expoTrans指数变换(Exponential transformation.可以处理负数).设置method参数 ...
R语定义函数对宽分布(wide distribution )、有偏分布（skew distribution）的数据进行对数变换（符号对数变换函数、signed log transformation）
R语言自定义函数对宽分布(wide distribution ).有偏分布(skew distribution)的数据进行对数变换(符号对数变换函数.signed log transformation ...
R语言Box-Cox变换实战（Box-Cox Transformation）：将非正态分布数据转换为正态分布数据、计算最佳λ、变换后构建模型
R语言Box-Cox变换实战(Box-Cox Transformation):将非正态分布数据转换为正态分布数据.计算最佳λ.变换后构建模型目录
Spark的transformation和action算子简介
transformation算子 map(func) 返回一个新的分布式数据集,由每个原元素经过func函数处理后的新元素组成 filter(func) 返回一个新的数据集,由经过func函数处理后返 ...
T4 (Text Template Transformation Toolkit)实现简单实体代码生成
在很多场景下,我们都需要代码生成.你可以使用CodeSmith,不过它是商业软件.VisualStudio2008中自带也有代码生成功能.那就是T4 (Text Template Transforma ...
【图像处理】透视变换 Perspective Transformation
透视变换(Perspective Transformation)是将图片投影到一个新的视平面(Viewing Plane),也称作投影映射(Projective Mapping).通用的变换公式为: ...
《从0到1学习Flink》—— Flink Data transformation(转换)
前言在第一篇介绍 Flink 的文章 <<从0到1学习Flink>-- Apache Flink 介绍> 中就说过 Flink 程序的结构 Flink 应用程序结构就是如上图 ...
ASP.NET 4.0 新特性--Web.Config Transformation（原创）
. 系列文章 Web.Config Transformation详解,这部分内容比较简单,关键是用没有用过的问题,所以这里希望帮助大家实践一下. 一概述: 在VS2010中引入在Config 文件中 ...
使用Web.Config Transformation配置灵活的配置文件
发布Asp.net程序的时候,开发环境和发布环境的Web.Config往往不同,比如connectionstring等.如果常常有发布的需求,就需要常常修改web.config文件,这往往是一件非常麻 ...

Transformation

Transformation相关推荐

最新文章

热门文章