CodeChef Chef and Segments

hfj1975

　　题目
　　


　　求区间积%mod
　　因为每个数<=100,质因数分解之后,就只需要计算对应质数的次方.
　　树状数组&#43;快速幂
　　


　　#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==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[x]+=val;
x+=lowbit(x);
}
}
inline int getsum(int x,int i)
{
int ans=0;
while(x>0)
{
ans+=c[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(&quot;%d&quot;,&n))
{
memset(c,0,sizeof(c));
for(int i=1;i<=n;i++)
{
scanf(&quot;%d&quot;,&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(&quot;%d&quot;,&t);
while(t--)
{
scanf(&quot;%d%d%d&quot;,&a,&b,&m);
long long ans=1;
for(int i=0;i<prime.size();i++)
{
ans*=Pow(prime,getsum(b,i)-getsum(a-1,i),m);
ans%=m;
}
printf(&quot;%lld\n&quot;,ans%m);
}
}
}