CodeChef Chef and Segments

题目

求区间积%mod

因为每个数<=100,质因数分解之后,就只需要计算对应质数的次方.

树状数组+快速幂

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<vector>
using namespace std;
vector<int>prime;
bool vis[110];
inline void init()
{
    memset(vis,0,sizeof(vis));
    for(int i=2;i<=100;i++)
    {
        if(vis[i]==0)
        {
            prime.push_back(i);
            for(int j=i*2;j<=100;j+=i)
                vis[j]=1;
        }
    }
}
int c[25][100010],n,a,t,b,m;
inline int lowbit(int x)
{
    return x&(-x);
}
inline void add(int x,int i,int val)
{
    while(x<=n)
    {
        c[i][x]+=val;
        x+=lowbit(x);
    }
}
inline int getsum(int x,int i)
{
    int ans=0;
    while(x>0)
    {
        ans+=c[i][x];
        x-=lowbit(x);
    }
    return ans;
}
inline long long Pow(int a,int b,int m)
{
    long long d=1LL,t=a*1LL;
    while(b)
    {
        if(b&1) d=d*t%m;
        b>>=1;
        t=t*t%m;
    }
    return d;
}
int main()
{
    init();
    while(~scanf("%d",&n))
    {
        memset(c,0,sizeof(c));
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&a);
            for(int j=0;j<prime.size();j++)
            {
                int num=0;
                while(a%prime[j]==0)
                {
                    num++;
                    a/=prime[j];
                }
                if(num!=0) add(i,j,num);
                if(a==1) break;
            }
        }
        scanf("%d",&t);
        while(t--)
        {
            scanf("%d%d%d",&a,&b,&m);
            long long ans=1;
            for(int i=0;i<prime.size();i++)
            {
                ans*=Pow(prime[i],getsum(b,i)-getsum(a-1,i),m);
                ans%=m;
            }
            printf("%lld\n",ans%m);
        }
    }
}