众所周知lxl是个毒瘤，Ynoi道道都是神仙题，题面好评
原题传送门
一看这题没有修改操作就知道这是莫队题（我也只会莫队）
我博客里对莫队的简单介绍
一个数N可以分解成\(p_1^{c_1}p_2^{c_2}…p_m^{c_m}\)
它的约数个数就是\((c_1+1)(c_2+1)…(c_m+1)\)
我们考虑先把每一个数分解质因数
用试除法会使你tle到没救，所以我们要用pollard's Rho来解决问题
（用质因数分解是因为\(10^9<2*3*5*7*11*13*17*19*23*29\)，质因数数量较少）
分解质因数后离散化，否则mle飞天
进行莫队，正常计算贡献，你交一发，会发现只有82（73），因为这个算法需要扫描每个数的每个质因子
我们很容易想到，一百（五百）以下的质数很多，而且他们作为质因子的次数也很多，所耗时间巨大
所以我们考虑将一百（五百）以下的质因数建立前缀和
这样就珂以过此题了（话说我pollard's Rho写得好丑啊）
#include
#define N 100005
#define mod 19260817
#define ll long long
#define getchar nc
using namespace std;
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
register int x=0,f=1;register char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+ch-'0',ch=getchar();
return x*f;
}
inline void write(register int x)
{
if(!x)putchar('0');if(x<0)x=-x,putchar('-');
static int sta[20];register int tot=0;
while(x)sta[tot++]=x%10,x/=10;
while(tot)putchar(sta[--tot]+48);
}
inline int Max(register int x,register int y)
{
return x>y?x:y;
}
ll Test[10]={5,61,24251,19260817};
inline int mul(register int a,register int b,register int m)
{
int d=((long double)a/m*b+1e-8);
int r=a*b-d*m;
return r<0?r+m:r;
}
inline int Pow(register int a,register int b,register int m)
{
ll r=1;
for(;b;b>>=1,a=mul(a,a,m))
if(b&1)
r=mul(r,a,m);
return r;
}
inline bool Query(register int P)
{
if(P==1)
return 0;
if(P==2||P==3||P==5||P==7)
return 1;
if(!(P&1)||(P%3==0)||(P%5==0)||(P%7==0))
return 0;
int t=P-1;
int k=0;
while(!(t&1))
++k,t>>=1;
for(register int i=0;i<3;++i)
{
if(P==Test[i])
return 1;
int a=Pow(Test[i],t,P),nxt=a;
for(register int j=1;j<=k;++j)
{
nxt=mul(nxt,nxt,P);
if(nxt==1&&a!=1&&a!=P-1)
return 0;
a=nxt;
}
if(a!=1)
return 0;
}
return 1;
}
inline int gcd(register int a,register int b)
{
if(!a)
return b;
if(!b)
return a;
int t=__builtin_ctzll(a|b);
a>>=__builtin_ctzll(a);
do{
b>>=__builtin_ctzll(b);
if(a>b)
a^=b^=a^=b;
b-=a;
}while(b!=0);
return a<1)
break;
}
}
if(t>1||(t=gcd(q,n))>1)
break;
}
if (t==n)
{
t=1;
while(t==1)
t=gcd(abs((x=g(x, n))-y),n);
}
return t;
}
int f[N][105],np[N];
vector v;
inline void find(register int x,register int id)
{
if(x==1)
return;
if(Query(x))
{
f[id][++np[id]]=x;
v.push_back(x);
return;
}
int p=x;
while(p==x)
p=pollardrho(x);
find(p,id);
find(x/p,id);
}
struct query{
int l,r,id,bl;
}q[N];
int n,m,blocksize;
inline bool cmp(register query a,register query b)
{
return a.bl!=b.bl?a.lb.r);
}
int ans[N],qaq[N<<1],inv[N<<1];
int num;
inline void add(register int pos)
{
for(register int i=1;i<=np[pos];++i)
num=(ll)num*inv[qaq[f[pos][i]]]%mod*(qaq[f[pos][i]]+1)%mod,++qaq[f[pos][i]];
}
inline void del(register int pos)
{
for(register int i=1;i<=np[pos];++i)
num=(ll)num*inv[qaq[f[pos][i]]]%mod*(qaq[f[pos][i]]-1)%mod,--qaq[f[pos][i]];
}
int ispri[600],pri[600],cnt=0,sum[N][100];
inline void init()
{
for(register int i=2;i<=500;++i)
ispri[i]=1;
for(register int i=2;i<=500;++i)
if(ispri[i])
{
pri[ispri[i]=++cnt]=i;
for(register int j=i<<1;j<=500;j+=i)
ispri[j]=0;
}
}
int main()
{
srand(19260817);
init();
n=read(),m=read();
blocksize=sqrt(n);
inv[1]=1;
for(register int i=2;i<=N;++i)
inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
for(register int i=1;i<=n;++i)
{
int x=read();
memcpy(sum[i],sum[i-1],sizeof(sum[i-1]));
for(register int j=1;j<=cnt&&pri[j]*pri[j]<=x;++j)
while(!(x%pri[j]))
{
++sum[i][j];
x/=pri[j];
}
if(x>1)
{
if(x<=pri[cnt])
{
++sum[i][ispri[x]];
continue;
}
find(x,i);
}
}
sort(v.begin(),v.end());
v.erase(unique(v.begin(),v.end()),v.end());
for(register int i=1;i<=n;++i)
for(register int j=1;j<=np[i];++j)
f[i][j]=lower_bound(v.begin(),v.end(),f[i][j])-v.begin();
for(register int i=1;i<=m;++i)
q[i].l=read(),q[i].r=read(),q[i].id=i,q[i].bl=(q[i].l-1)/blocksize+1;
sort(q+1,q+1+m,cmp);
for(register int i=0;iq[i].l)
add(--l);
while(r>q[i].r)
del(r--);
while(l
