题目大意:
给定一个序列,每次询问一段区间的数的乘积的约数个数。
解题思路:
在太阳西斜的这个世界里,置身天上之森。等这场战争结束之后,不归之人与望眼欲穿的众人, 人人本着正义之名,长存不灭的过去、逐渐消逝的未来。我回来了,纵使日薄西山,即便看不到未来,此时此刻的光辉,盼君勿忘。————世界上最幸福的女孩
我永远喜欢珂朵莉。
---
\(10^9\)以内的数最多有10个不同的质因子。
考虑对其质因数分解。
由于值域范围过大,考虑使用Pollard-Rho算法。
这里普通的Pollard-Rho算法可能会TLE。如果你的代码能通过模板题,那基本上没问题(窝反正直接把以前写的板子拉过来然后调了调参)。
之后,你就会得到最多\(10n\)个不同的质因数。对其进行离散化,开桶记录。
然后上莫队,对于每次指针的偏移,把它所有的质因数加到桶里,同时维护约数个数即可。
这部分时间复杂度\(O(10n\sqrt n)\),加上上面的质因数分解的玄学期望复杂度,只能获得82分的好成绩。
---
我们考虑把每个数\(1000\)以内的质因子先取出来(\(1000\)以内共168个质数),然后,对其做前缀和,记录前缀的出现次数。
然后,由于\(1001^3>10^9\),所以每个数剩下最多不超过2个质因子。这部分用Pollard_Rho找即可。
然后莫队的时候,对于前面168个质数就可以不用维护,直接用前缀和。
而对于后面的大质因子,再离散化处理即可。由于每个数最多两个质因子,所以常数就小了很多。
而由于筛掉了很多小的质因子,Pollard_Rho的速度也会变快。然后就足以通过此题。
C++ Code:
#include<cctype>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cstring>
#define ctz __builtin_ctz
using namespace std;
#ifdef ONLINE_JUDGE
struct istream{
char buf[23333333],*s;
inline istream(){
buf[fread(s=buf,1,23333330,stdin)]='\n';
fclose(stdin);
}
inline istream&operator>>(int&d){
d=0;
for(;!isdigit(*s);++s);
while(isdigit(*s))
d=(d<<3)+(d<<1)+(*s++^'0');
return*this;
}
}cin;
struct ostream{
char buf[8000005],*s;
inline ostream(){s=buf;}
inline ostream&operator<<(int d){
if(!d){
*s++='0';
}else{
static int w;
for(w=1;w<=d;w*=10);
for(;w/=10;d%=w)*s++=d/w^'0';
}
return*this;
}
inline ostream&operator<<(const char&c){*s++=c;return*this;}
inline void flush(){
fwrite(buf,1,s-buf,stdout);
s=buf;
}
inline~ostream(){flush();}
}cout;
#else
#include<iostream>
#endif
int pri[170],cct=0,sum[100005][169],num[1005];
void sieve(){
for(int i=2;i<=1000;++i)num[i]=1;
for(int i=2;i<=1000;++i)
if(num[i]){
pri[num[i]=++cct]=i;
for(int j=i<<1;j<=1000;j+=i)num[j]=0;
}
}
using LoveLive=long long;
vector<int>tj;
const int pr[]={2,3,5,24251,61,19260817};
int gcd(int a,int b){
if(!a||!b)return a|b;
int t=ctz(a|b);
a>>=ctz(a);
do{
b>>=ctz(b);
if(a>b)swap(a,b);
b-=a;
}while(b);
return a<<t;
}
inline int power(int a,int b,const int&md){
int ans=1;
for(;b;b>>=1){
if(b&1)ans=(LoveLive)ans*a%md;
a=(LoveLive)a*a%md;
}
return ans;
}
int miller_rabin(int p){
if(p==2)return 1;
int b=p-1;
int t=0;
while(!(b&1))b>>=1,++t;
for(int i:pr){
int r=power(i%(p-2)+2,b,p);
if(r==1||r==p-1)continue;
int ok=1;
for(int j=1;j<=t&&ok;++j){
r=(LoveLive)r*r%p;
if(r==p-1)ok=0;
}
if(ok)return 0;
}
return 1;
}
void pollard_rho(int&n,int c){
int k=2,x=rand()%(n-1)+1,y=x,q=1,t=1;
for(;;k<<=1,y=x,q=1){
for(int i=1;i<=k;++i){
x=((LoveLive)x*x%n+c)%n;
q=(LoveLive)q*abs(x-y)%n;
if(!(i&63)){
t=gcd(q,n);
if(t>1)break;
}
}
if(t>1||(t=gcd(q,n))>1)break;
}
n=t;
}
void find(int n,int c){
if(n==1)return;
if(miller_rabin(n)){tj.push_back(n);return;}
int p=n;
while(p>=n)pollard_rho(p,c--);
n/=p;
tj.push_back(n),tj.push_back(p);
}
#define N 100005
const int md=19260817;
int n,m,a[N],inv[N],cnt[N],tot[N*2],now=1,ans[N];
int p[N][3];
vector<int>lr;
struct que{
static const int siz=317;
int l,r,id;
inline bool operator<(const que&rhs)const{
return((l/siz!=rhs.l/siz)?(l<rhs.l):r<rhs.r);
}
}q[N];
inline void add(int id){
for(register int i=1;i<=cnt[id];++i)
now=(LoveLive)now*inv[tot[p[id][i]]]%md*(tot[p[id][i]]+1)%md,++tot[p[id][i]];
}
inline void del(int id){
for(register int i=1;i<=cnt[id];++i)
now=(LoveLive)now*inv[tot[p[id][i]]]%md*(tot[p[id][i]]-1)%md,--tot[p[id][i]];
}
int main(){
srand(19260817);
sieve();
cin>>n>>m;
inv[1]=1;
for(int i=2;i<=n;++i)inv[i]=(md-md/i)*1LL*inv[md%i]%md;
for(int i=1;i<=n;++i){
cin>>a[i];
memcpy(sum[i],sum[i-1],sizeof*sum);
for(int j=1;j<=cct&&pri[j]*pri[j]<=a[i];++j)
while(!(a[i]%pri[j])){
++sum[i][j];
a[i]/=pri[j];
}
if(a[i]>1){
if(a[i]<=pri[cct]){
++sum[i][num[a[i]]];
continue;
}
tj.clear();
find(a[i],23333);
for(int it:tj)
p[i][++cnt[i]]=it,lr.push_back(it);
}
}
sort(lr.begin(),lr.end());
lr.erase(unique(lr.begin(),lr.end()),lr.end());
for(int i=1;i<=n;++i)
for(int j=1;j<=cnt[i];++j)p[i][j]=lower_bound(lr.begin(),lr.end(),p[i][j])-lr.begin();
for(int i=1;i<=m;++i)cin>>q[q[i].id=i].l>>q[i].r;
sort(q+1,q+m+1);
for(int i=0;i<n<<1;++i)tot[i]=1;
for(int i=1,l=1,r=0;i<=m;++i){
while(r<q[i].r)add(++r);
while(l>q[i].l)add(--l);
while(r>q[i].r)del(r--);
while(l<q[i].l)del(l++);
int&out=ans[q[i].id];out=now;
for(int j=1;j<=cct;++j)
out=(LoveLive)out*(sum[r][j]-sum[l-1][j]+1)%md;
}
for(int i=1;i<=m;++i)
cout<<ans[i]<<'\n';
return 0;
}
来源:oschina
链接:https://my.oschina.net/u/4347225/blog/3719066