智乃与无意义的题目(西工大校赛G题）

题目：（链接智乃与无意义的题目
）

思路:
本题可采用分治算法且题中多次询问，故想到用线段树来维护(原来线段树并不只是用来维护区间和的啊555

用线段树进行维护，对每一个叶子节点进行质因数分解，得到它以2,3,5,7做为因子的个数,向上合并的时候累加相同因子个数，最后查询的时候用唯一分解定理的因子个数公式计算得到。

在求数n的因子个数时采用 唯一分解定理 （那就先上背景知识（给我自己这个菜鸟看，qwq，大佬可跳过这一步

注：由于本题1<=ai<=10,10以内的质数只有2,3,4,7，故只需维护这4个数出现的次数即可

代码：

#include <cstdio>
#include <cmath>
#include <set>
#include <map>
#include <string>
#include <cstring>
#include <iostream>
#include <vector>
#include <queue>
#include <stack>
#include <algorithm>
#include <iterator>
#include<math.h>
#define debug() puts("what the fuck")
#define ll long long
#define INF 0x3f3f3f3f
#define maxn 2000000
using namespace std;
struct Node
{
  int p[4];
}s[maxn*4];
ll f[maxn];
int b[4]={2,3,5,7};
void countt(int m,int pos)
{
    for(int i=0;i<4;i++)
    {
        if(m%b[i]==0)
        {
            while(m%b[i]==0)
            {
                s[pos].p[i]++;
                m/=b[i];
            }
        }
    }
}
void pushup(ll v)
{
     for(int i=0;i<4;i++)
    {
        s[v].p[i]=s[v*2].p[i]+s[v*2+1].p[i];
    }
}
void build(ll l,ll r,ll v)
{
    if(l==r)
    {
        countt(f[l],v);
        return;
    }
    ll mid=(l+r)/2;
    build(l,mid,v*2);
    build(mid+1,r,v*2+1);
    pushup(v);
}
ll query(ll L,ll R,ll l,ll r,ll v,ll id)
{
    if(L<=l&&R>=r)
        return s[v].p[id];
    ll mid=(l+r)/2;
    ll ans=0;
    if(L<=mid)
        ans+=query(L,R,l,mid,v*2,id);
    if(R>mid)
        ans+=query(L,R,mid+1,r,v*2+1,id);
    return ans;
}
void update(ll L,ll R,ll l,ll r,ll v)
{
    if(l==r)
    {
        for(int i=0;i<4;i++)
            s[v].p[i]=0;
        countt(R,v);
        return;
    }
    ll mid=(l+r)/2;
    if(L<=mid)
        update(L,R,l,mid,v*2);
    else
        update(L,R,mid+1,r,v*2+1);
    pushup(v);
}
int main()
{
    ll n,m,a,b,c,w,ans,q;
    ll mod=998244353;
    scanf("%lld%lld",&n,&m);
    for(ll i=1; i<=n; i++)
        scanf("%lld",&f[i]);
    build(1,n,1);
    while(m--)
    {
        scanf("%lld%lld%lld",&c,&a,&b);
        if(c==2)
        {
            ll ans=1;
            for(int i=0;i<4;i++)
            {
                ans=((ans%mod)*((query(a,b,1,n,1,i)+1)%mod))%mod;
            }
            cout<<ans<<endl;
        }
        else if(c==1)
        {
            update(a,b,1,n,1);
        }
    }
    return 0;
}