“美登杯”上海市高校大学生程序设计邀请赛 (华东理工大学) E. 小花梨的数组 —— 线段树区间更新，乘除最小质因子运算

This way

题意：

给你n个数，每次给你一个操作，操作分为3种：
1 l r 这个区间所有的数乘上对应数的最小质因子（没有就乘1）
2 l r 这个区间所有的数除上对应数的最小质因子（没有就除1）
3 x这个数是多少

题解：

线段树无法区间更新不同的数，所以没次区间直接乘或者除是不现实的，那么我们就需要用一个东西去维护这个最小质因子。第一个想到的就是map，因为它自动排序而且查询方便。那么现在我们就可以用mi数组维护每一个位置的最小质因子数量。但是它肯定会小于0，怎么办？那肯定是变成0的时候改变最小质因子拉，这时候只需要erase(map.begin())，然后下一个map.begin()就是当前最小质因子。怎么找那些位置变成0？那么我们只需要维护区间最小值，用一个独立的方法在每一次操作结束之后看看哪个区间的最小值是0，已知追溯到l==r的时候，我们就可以更新mi数组了。
在3操作的时候注意我们之前1和2操作在线段树区间更新的时候并不能直接更新map里面的内容，因为那样复杂度可能过大，所以对于x位置的最小质因子数量我们是要用单点查询找出。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N=1e5+5;
const ll mod=1e9+7;
map<ll,int>mp[N];
int mi[N*4],flag[N*4];
void build(int l,int r,int root)
{
    if(l==r)
    {
        if(!mp[l].size())
            mi[root]=1e9;
        else
            mi[root]=mp[l].begin()->second;
        return ;
    }
    int mid=l+r>>1;
    build(l,mid,root<<1);
    build(mid+1,r,root<<1|1);
    mi[root]=min(mi[root<<1],mi[root<<1|1]);
}
void push_down(int root)
{
    if(!flag[root])
        return ;
    mi[root<<1]+=flag[root];
    mi[root<<1|1]+=flag[root];
    flag[root<<1]+=flag[root];
    flag[root<<1|1]+=flag[root];
    flag[root]=0;
}
void update(int l,int r,int root,int ql,int qr,int op)
{
    if(l>=ql&&r<=qr)
    {
        mi[root]+=op;
        flag[root]+=op;
        return ;
    }
    push_down(root);
    int mid=l+r>>1;
    if(mid>=ql)
        update(l,mid,root<<1,ql,qr,op);
    if(mid<qr)
        update(mid+1,r,root<<1|1,ql,qr,op);
    mi[root]=min(mi[root<<1],mi[root<<1|1]);
}
//查找最小质因子更新的地方：
void query(int l,int r,int root)
{
    if(l==r)
    {
        if(mi[root]==0)
        {
            mp[l].erase(mp[l].begin());
            if(!mp[l].size())
                mi[root]=1e9;
            else
                mi[root]=mp[l].begin()->second;
        }
        return ;
    }
    push_down(root);
    int mid=l+r>>1;
    if(mi[root<<1]==0)
        query(l,mid,root<<1);
    if(mi[root<<1|1]==0)
        query(mid+1,r,root<<1|1);
    mi[root]=min(mi[root<<1],mi[root<<1|1]);
}
ll qpow(ll a,ll b)
{
    ll ret=a,ans=1;
    while(b)
    {
        if(b&1)
            ans=ans*ret%mod;
        ret=ret*ret%mod;
        b>>=1;
    }
    return ans;
}
int findbegin(int l,int r,int root,int pos)
{
    if(l==r)
        return mi[root];
    push_down(root);
    int mid=l+r>>1;
    if(mid>=pos)
        return findbegin(l,mid,root<<1,pos);
    else
        return findbegin(mid+1,r,root<<1|1,pos);
}
int main()
{
    int n,m;
    scanf("%d%d",&n,&m);
    ll x;
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&x);
        for(ll j=2;j<=sqrt(x);j++)
        {
            while(x%j==0)
                x/=j,mp[i][j]++;
        }
        if(x>1)
            mp[i][x]++;
    }
    build(1,n,1);
    while(m--)
    {
        int op,l,r;
        scanf("%d",&op);
        if(op==1)
            scanf("%d%d",&l,&r),update(1,n,1,l,r,1);
        else if(op==2)
            scanf("%d%d",&l,&r),update(1,n,1,l,r,-1);
        else
        {
            scanf("%d",&l);
            ll ans=1;
            if(mp[l].size())
            {
                map<ll,int>::iterator it=mp[l].begin();
                ans=qpow(it->first,findbegin(1,n,1,l));
                it++;
                for(;it!=mp[l].end();it++)
                    ans=ans*qpow(it->first,it->second)%mod;
            }

            printf("%lld\n",ans);
        }
        query(1,n,1);
    }
    return 0;
}