A1339. JZPLCM(顾昱洲)|树状数组|hash表|逆元|分解质因数

各种细节实现

注意 树状数组 的下标不能为 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;
}