Description
f0(n)fr+1==∑uv=n[(u,v)=1]1∗fr
求
fr(n)
。
Solution
因为
n
的所有质因子之间对
那么就有所有
fr
都是积性函数。
将每个质因子单独考虑,狄利克雷卷积的形式就可以化成简单的求和(可能就是普通的卷积???)。
fr+1(pa)=∑i=0afr(pi)
发现
f0(pa)=[a≠0]+1
。
直接DP就好啦。
#include <bits/stdc++.h>
using namespace std;
const int N = 1010101;
const int MOD = 1000000007;
typedef long long ll;
inline char get(void) {
static char buf[100000], *S = buf, *T = buf;
if (S == T) {
T = (S = buf) + fread(buf, 1, 100000, stdin);
if (S == T) return EOF;
}
return *S++;
}
template<typename T>
inline void read(T &x) {
static char c; x = 0;
for (c = get(); c < '0' || c > '9'; c = get());
for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0';
}
int n, r, q, x, sum, Pcnt, ans, lim;
vector<int> fact[N];
int f[N][30];
int prime[N], vis[N];
inline void Add(int &x, int a) {
x += a; while (x >= MOD) x -= MOD;
}
int main(void) {
freopen("1.in", "r", stdin);
read(q); f[0][0] = 1; lim = 1000000;
for (int i = 1; i <= 20; i++) f[0][i] = 2;
for (int i = 1; i <= lim; i++) {
sum = 0;
for (int j = 0; j <= 20; j++) {
Add(sum, f[i - 1][j]);
Add(f[i][j], sum);
}
}
for (int i = 2; i <= lim; i++) {
if (!vis[i]) prime[++Pcnt] = i;
for (int j = 1; j <= Pcnt && (x = prime[j] * i) <= lim; j++) {
vis[x] = 1; if (!(i % prime[j])) break;
}
}
for (int i = 1; i <= Pcnt; i++)
for (int j = prime[i]; j <= lim; j += prime[i])
fact[j].push_back(prime[i]);
while (q--) {
read(n); read(r); ans = 1;
for (int d: fact[r]) {
x = 0;
while (!(r % d)) {
r /= d; ++x;
}
ans = (ll)ans * f[n][x] % MOD;
}
printf("%d\n", ans);
}
return 0;
}