选拔赛之 线段树

id=98

这道题自己陷入了一个坑，，

用了快速乘， 就会超时，，

快速乘只是用来防止中间结果溢出的，反而会比乘法慢

顺便复习下快速乘代码，，

这个地方还有一个坑就是数论里面的东西。。现在还不知道为什么。。

有个小坑就是ak≡ak mod (p−1)(mod p)a ^ k \equiv a ^ {k\ mod \ (p - 1)} (mod\ p)a​k​​≡a​k mod (p−1)​​(mod p),$p$是素数.

#include<cstdio>#include<cstring>#include<algorithm>#include<queue>#include<string>#include<cstring>#include<iomanip>#include<iostream>#include<stack>#include<cmath>#include<map>#include<vector>#define ll long long#define inf 0x3f3f3f3f#define INF 1000000000#define bug1 cout<<"bug1"<<endl;#define bug2 cout<<"bug2"<<endl;#define bug3 cout<<"bug3"<<endl;using namespace std;const int maxn=100005;const int mod=1e9+7;int a[maxn];ll sum[maxn*4],lazy[maxn*4];ll qpow(ll a,ll k){    ll ret=1;    while(k){        if(k&1){            ret=ret*a%mod;        }        a=a*a%mod;        k>>=1;    }    return ret;}void build(int rt,int l,int r){    if(l==r){sum[rt]=a[l];return;}    int mid=(l+r)>>1;    build(rt<<1,l,mid);    build(rt<<1|1,mid+1,r);    sum[rt]=sum[rt<<1]*sum[rt<<1|1]%mod;}void down(int rt){    sum[rt<<1]=qpow(sum[rt<<1],lazy[rt]);    sum[rt<<1|1]=qpow(sum[rt<<1|1],lazy[rt]);    lazy[rt<<1]=lazy[rt<<1]*lazy[rt]%(mod-1);    lazy[rt<<1|1]=lazy[rt<<1|1]*lazy[rt]%(mod-1);// 是mod   mod-1    lazy[rt]=1;}void update(int ql,int qr,int l,int r,int rt,int b){    if(ql<=l&&r<=qr){        sum[rt]=qpow(sum[rt],b);        lazy[rt]=lazy[rt]*b%(mod-1);        return;    }    if(lazy[rt]!=1)down(rt);    int mid=(l+r)>>1;    if(ql<=mid)update(ql,qr,l,mid,rt<<1,b);    if(mid<qr)update(ql,qr,mid+1,r,rt<<1|1,b);    sum[rt]=sum[rt<<1]*sum[rt<<1|1]%mod;}ll query(int ql,int qr,int l,int r,int rt){    if(ql<=l&&r<=qr){        return sum[rt];    }    if(lazy[rt]!=1)down(rt);    int mid=(l+r)>>1;    ll ret=1;    if(ql<=mid)ret=ret*query(ql,qr,l,mid,rt<<1)%mod;    if(mid<qr)ret=ret*query(ql,qr,mid+1,r,rt<<1|1)%mod;    return ret;}void init(){    for(int i=1;i<4*maxn;++i){        lazy[i]=1;        sum[i]=1;    }    memset(a,0,sizeof(a));}int main(){    int t;    scanf("%d",&t);    while(t--){        init();        int n,q;        scanf("%d%d",&n,&q);        for(int i=1;i<=n;++i){            scanf("%d",&a[i]);        }        build(1,1,n);        for(int i=1;i<=q;++i){            int op;            scanf("%d",&op);            if(op==1){                int l,r,b;                scanf("%d%d%d",&l,&r,&b);                update(l,r,1,n,1,b);            }            else{                int l,r;                scanf("%d%d",&l,&r);                cout<<query(l,r,1,n,1)<<'\n';            }        }    }}