UVA 11440 - Help Tomisu
题意:给定n和m,求[2,n!]中,所有质因子个数都大于m的个数
思路:ϕ(m!)表示小于m!并与m!互质的个数,而与m!互质的个数,他的质因子肯定不包含1-m,因此就是满足条件的。然后对于这题而言,则是要求n!中,不与m!互质的个数,答案取模100000007
那么先看一个证明:
求kn中与n互质的个数,答案为kϕ(n)。
ϕ(n)表示1-n中与n互质的个数,那么由此考虑[n + 1, 2n], [2n + 1, 3n]...这每个区间中的每个数字都等于1-n中数字加上kn,对于原来就与n不互质的个数,加上n仍会有一个质因子重复,所以仍然不行,那么对于原来互质的数x,gcd(x, n) = 1,那么可知gcd(x + kn, n) = 1,仍然是互质的,所以每隔n的区间与n互质的个数是相同的,所以答案kϕ(n)
所以对于这道题目,答案就变成了n!/m!ϕ(m!),那么问题只剩下如何求ϕ(m!)。
已知ϕ(n)求法为n∗(1−1/p1)∗(1−1/p2)....(1−1/pn) (p为n的质因子),因此对于m!而言,分子为m!,分母为1 - m所有质数的(1−1/p)之乘积
到这里答案就可以求了,把m!消掉,得到n!/∏(1−1/pi)mod1000000007,先预处理那些表,每次去计算即可
代码:
#include <cstdio>
#include <cstring>
const long long N = 10000005;
const long long mod = 100000007;
long long ispri[N], fac[N], phi[N];
long long exgcd(long long a, long long b, long long &x, long long &y) {
if (!b) {x = 1; y = 0; return a;}
long long d = exgcd(b, a % b, y, x);
y -= a / b * x;
return d;
}
long long inv(long long a, long long n) {
long long x, y;
exgcd(a, n, x, y);
return (x + n) % n;
}
void get_table() {
fac[0] = fac[1] = 1; phi[0] = phi[1] = 1;
for (long long i = 2; i < N; i++) {
fac[i] = (fac[i - 1] * i) % mod;
if (ispri[i]) {
phi[i] = phi[i - 1];
continue;
}
phi[i] = phi[i - 1] * (i - 1) % mod * inv(i, mod) % mod;
for (long long j = i * i; j < N; j += i)
ispri[j] = 1;
}
}
int n, m;
int main() {
get_table();
while (~scanf("%d%d", &n, &m) && n) {
printf("%lld\n", ((fac[n] * phi[m] - 1) % mod + mod) % mod);
}
return 0;
}