【51Nod 1769】クラークとMath2の
フェイス質問
問題の解決策
数論的関数のための\(F \)、\ (| N-F(D)=(F \タイムズ1)(N)\ \ sum_ {D}) 。
実際には、被験体が必要とされる\(G = F \回1 ^ K \) 。
考えてみましょう(1 ^ k個(N)\ \) ので、見つける方法\(1(n)は\)乗法関数であるので、\(^ 1 K(n)は \) にも乗法関数です。
我々は検討のために、\(N- \)各品質係数の\(P \)と回数それ\(R&LTの\)は、関数の対応する値を得ることができます。
各異なる等価\(I_ {J-1} \) と\(i_j \)で表さ品質係数挿入\(P \)回、また同じ場所に挿入することができます。
だから、順列と組み合わせの回答に基づいて、この問題がある\(K-1 +を選択してくださいRの\ \ r)は、その後、一緒に乗じ、各品質係数の寄与が答えです。
コード
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
const int Mod = 1e9 + 7;
inline int gi() {
register int data = 0, w = 1;
register char ch = 0;
while (!isdigit(ch) && ch != '-') ch = getchar();
if (ch == '-') w = -1, ch = getchar();
while (isdigit(ch)) data = (10ll * data + ch - '0') % Mod, ch = getchar();
return w == 1 ? data : (-data + Mod) % Mod;
}
int fpow(int x, int y) {
int res = 1;
while (y) {
if (y & 1) res = 1ll * res * x % Mod;
x = 1ll * x * x % Mod;
y >>= 1;
}
return res;
}
const int MAX_N = 5e5 + 5;
int fac[25], ifc[25], C[25];
int N, K, f[MAX_N], g[MAX_N], h[MAX_N], cur[MAX_N];
int main () {
#ifndef ONLINE_JUDGE
freopen("cpp.in", "r", stdin);
#endif
fac[0] = 1; for (int i = 1; i <= 20; i++) fac[i] = 1ll * i * fac[i - 1] % Mod;
ifc[20] = fpow(fac[20], Mod - 2);
for (int i = 19; ~i; i--) ifc[i] = 1ll * ifc[i + 1] * (i + 1) % Mod;
N = gi(), K = gi();
for (int i = 1; i <= N; i++) f[i] = gi(), g[i] = 1, cur[i] = i;
for (int i = 0; i <= 20; i++) {
int nw = (i + K - 1) % Mod; C[i] = ifc[i];
for (int j = 0; j < i; j++) C[i] = 1ll * C[i] * (nw - j + Mod) % Mod;
}
for (int i = 2; i <= N; i++) {
if (cur[i] == 1) continue;
for (int j = i; j <= N; j += i) {
int k = 0;
while (cur[j] % i == 0) ++k, cur[j] /= i;
g[j] = 1ll * g[j] * C[k] % Mod;
}
}
for (int i = 1; i <= N; i++)
for (int j = i; j <= N; j += i)
h[j] = (h[j] + 1ll * f[i] * g[j / i]) % Mod;
for (int i = 1; i <= N; i++) printf("%d ", h[i]);
putchar('\n');
return 0;
}