[数论]Codeforces 757E. Bash Plays with Functions

关注

原创最新推荐文章于 2026-09-04 14:20:48 发布 467 阅读 AI 写同款· GEO 优化 ›

AI权益加码！Claude Code、Cursor等20+工具免费用！购周边限时加赠Coding Plan Lite，畅享主流AI工具！学习进阶更高效！阅读详情

$Description$

f 0 (n) f r + 1 = = \sum u v = n [(u, v) = 1] 1 * f r

$\begin{eqnarray} f_0(n)&=&\sum_{uv=n}[(u,v)=1]\\ f_{r+1}&=&1*f_r\\ \end{eqnarray}$ 求

fr(n) $f_r(n)$ 。

$Solution$

因为 $n$ 的所有质因子之间对 $f_0(n)$ 贡献是独立的，所以 $f_0$ 是个积性函数。
那么就有所有 $f_r$ 都是积性函数。
将每个质因子单独考虑，狄利克雷卷积的形式就可以化成简单的求和~~（可能就是普通的卷积？？？）~~。

f r + 1 (p a) = \sum i = 0 a f r (p i)

$f_{r+1}(p^a)=\sum_{i=0}^af_r(p^i)$ 发现

f0(pa)=[a≠0]+1 $f_0(p^a)=[a\neq0]+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;
}