树状数组区间修改求和

#include
using namespace std;
const int N=1e7+10;
const int mod=998244353;
bool p[N];
int prime[N],s[N],cp;
long long tree[4][N];
int lowbit(int x){return x&-x;}
void add(int x,int id,int p){while(id>n>>q;
    for(int i=1;i<=n;i++)
    {
        int x,sum=0;
        scanf("%d",&x);
        while(x!=1)
        {
            int tmp=s[x];
            while(x%tmp==0) sum+=tmp,x/=tmp;
        }
        if(sum!=0)
        {
            add(1,i,2);
            add(-1,i+1,2);
            add(i,i,3);
            add(-i-1,i+1,3);
        }
        add(sum,i,0);
        add(-sum,i+1,0);
        add(1ll*i*sum%mod,i,1);
        add(-1ll*(i+1)*sum%mod,i+1,1);
    }
    for(int i=1;i<=q;i++)
    {
        int op;
        scanf("%d",&op);
        if(op==1)
        {
            int l,r;
            scanf("%d%d",&l,&r);
            long long ans=1ll*(r+1)*query(0,r)%mod-query(1,r)-(1ll*l*query(0,l-1)%mod-query(1,l-1));
            if((1ll*(r+1)*query(2,r)%mod-query(3,r)-(1ll*l*query(2,l-1)%mod-query(3,l-1)))==0) ans=1;
            ans=(ans%mod+mod)%mod;
            printf("%lldn",ans);
        }
        else
        {
            int l,r,x,sum=0;
            scanf("%d%d%d",&l,&r,&x);
            while(x!=1)
            {
                int tmp=s[x];
                while(x%tmp==0) sum+=tmp,x/=tmp;
            }
            if(sum!=0)
            {
                add(1,l,2);
                add(-1,r+1,2);
                add(l,l,3);
                add(-r-1,r+1,3);
            }
            add(sum,l,0);
            add(-sum,r+1,0);
            add(1ll*l*sum%mod,l,1);
            add(-1ll*(r+1)*sum%mod,r+1,1);
        }
    }
}