题意:
有函数
f0,n=将n分解为两个互质的数的积的方案数
f
0
,
n
=
将
n
分
解
为
两
个
互
质
的
数
的
积
的
方
案
数
fr,n=∑fr−1,u+fr−1,v2(uv=n)
f
r
,
n
=
∑
f
r
−
1
,
u
+
f
r
−
1
,
v
2
(
u
v
=
n
)
m组询问,给出r,n,求
fr,n
f
r
,
n
对10^9+7取模的余数
m,r,n<=10^6
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<cmath>
#include<iostream>
#define N 1000010
#define M 22
#define mmod 1000000007
#define LL long long
using namespace std;
int m,f[N][M],p[N],pl,minp[N],ans,s[N][M];
bool b[N];
void get_p()
{
for(int i=2;i<N;i++)
{
if(b[i]==0) {p[++pl]=i;minp[i]=i;}
for(int j=1;j<=pl;j++)
{
if(i*p[j]>=N) break;
b[i*p[j]]=1;
minp[i*p[j]]=p[j];
if(i%p[j]==0) break;
}
}
}
int main()
{
get_p();
for(int i=1;i<M;i++) f[0][i]=2,s[0][i]=2*i+1;
f[0][0]=s[0][0]=1;
for(int i=1;i<N;i++)
{
for(int j=0;j<M;j++) f[i][j]=s[i-1][j];
s[i][0]=f[i][0];
for(int j=1;j<M;j++) s[i][j]=(s[i][j-1]+f[i][j])%mmod;
}
scanf("%d",&m);
while(m--)
{
int n,r;
scanf("%d%d",&r,&n);
ans=1;
while(n>1)
{
int t=minp[n],k=0;
while(minp[n]==t) {n/=t;k++;}
ans=(LL)ans*f[r][k]%mmod;
}
printf("%d\n",ans);
}
return 0;
}
题解:
设n有k种不同质因子。容易看出:
f0,n=2k
f
0
,
n
=
2
k
fr,n=∑d|nfr−1,d
f
r
,
n
=
∑
d
|
n
f
r
−
1
,
d
然后我就没什么想法了,看了一眼题解。。(捂脸)
注意
f0
f
0
是一个积性函数,而
fr=fr−1∗1
f
r
=
f
r
−
1
∗
1
也是积性函数,所以只需要求出
fr,pk
f
r
,
p
k
。
注意
fr,pk
f
r
,
p
k
的值只与r和k有关,暴力预处理,询问时分解质因数即可。
这题看起来不像数学题,但要用到数学题的分析方法。