各种细节实现
注意 树状数组 的下标不能为 0
注意 代码条理清晰
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
const int maxn=1e5+10;
const int mo=1000000007;
const int sqrn=31623;
#define LL long long
using namespace std;
int n,m;
//shai
int prime[sqrn],cnt_prime;
bool isprime[sqrn];
void shai()
{
for(int i=2;i<=sqrn;i++)
{
if(!isprime[i]) prime[++cnt_prime]=i;
for(int j=1;j<=cnt_prime&&prime[j]*i<=sqrn;j++)
{
isprime[prime[j]*i]=true;
if(i%prime[j]==0)break;
}
}
}
//num
int num1[maxn*20],num2[maxn*20],c[maxn*20],cnt_num;
int s[maxn],t[maxn];
void chai(int x)
{
for(int i=1;i<=cnt_prime;i++)
if(x%prime[i]==0)
{
int tmp=prime[i];
while(x%prime[i]==0) num1[++cnt_num]=tmp,num2[cnt_num]=prime[i],tmp=tmp*prime[i],x/=prime[i];
}
if(x!=1) num1[++cnt_num]=x,num2[cnt_num]=x;
}
//st
int st[maxn*20],st_loc[maxn*20],st_mo=maxn*20-1;
int st_q(int x)
{
int y=x%st_mo;
while(st[y]!=0&&st[y]!=x) y++,y%=st_mo;
return y;
}
//Q
int ans[maxn];
struct H {int l,r,id;}qq[maxn];
int tree[maxn*20];
bool cmp(H a,H b){return a.r<b.r;}
int lowbit(int x)
{
return -x&x;
}
void add(int x,int y)
{
x++;
while(x<=cnt_num)
{
tree[x]=(LL)tree[x]*y%mo;
x+=lowbit(x);
}
}
int q(int x)
{
x++;
int ans=1;
while(x>0)
{
ans=(LL)ans*tree[x]%mo;
x-=lowbit(x);
}
return ans;
}
//ni
int sum[maxn*20];
int f(int x,int y)
{
if(y==0) return 1;
int t=f(x,y/2);
t=(LL)t*t%mo;
if(y&1) t=(LL)t*x%mo;
return t;
}
//main
int main()
{
shai();
scanf("%d %d",&n,&m);
for(int i=1;i<=n;i++) {int x;scanf("%d",&x);s[i]=cnt_num+1;chai(x);t[i]=cnt_num;}
for(int i=1;i<=cnt_num;i++) {int x=st_q(num1[i]);c[i]=st_loc[x];st[x]=num1[i];st_loc[x]=i;}
for(int i=1;i<=m;i++) scanf("%d %d",&qq[i].l,&qq[i].r),qq[i].id=i;
sort(qq+1,qq+m+1,cmp);
int now=0;
for(int i=1;i<=cnt_num;i++) tree[i]=1;
sum[0]=1;for(int i=1;i<=cnt_num;i++) sum[i]=(LL)sum[i-1]*num2[i]%mo;
for(int i=1;i<=cnt_num;i++) sum[i]=f(sum[i],mo-2);
for(int i=1;i<=m;i++)
{
int r=qq[i].r,l=qq[i].l,id=qq[i].id;
while(now<r) for(int i=s[++now];i<=t[now];i++) add(c[i],num2[i]);
ans[id]=(LL)q(t[l-1])*sum[t[l-1]]%mo;
}
for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
return 0;
}