C. Divisors of the Divisors of An Integer

题目

https://codeforces.com/gym/102040

分析

可以把 n! 分解质因数，成为 sigma( p[i] ^ a[i] ) 的形式，然后从空开始每次讨论加入一个质因数 p，次数为 a ，设上一次加了之后答案是 tmp，则在 n! 的所有约数中，所有最大质因子恰为 p 的数x的 d 函数和为 (2+3+...+(a+1)) * (1 + tmp)，左半部分意思是 x 中 p 的不同次数得到的贡献，右半部分是除了 p 以外其它因子得到的贡献（加一1代表无其他银子）。然后再加上 tmp 就得到了最大质因子不大于 p 的所有约数的 d 和。

1e6 以内有大概八万个质数。

怎么得出 n! 中质因子 p 的次数呢？考虑 1~n 中次数不小于 1 的数的个数，次数不小于 2 的数的个数。。。就会发现答案就是 n/p，n/p/p，n/p/p … 累加即可。

对每个输入，枚举质因子，然后依次算出 n! 中他们的次数并处理贡献就行了。

粗糙的时间复杂度大约是 Tlgnlgn 吧……

代码

#include <bits/stdc++.h>
using namespace std;
#define Ha 10000007
#define MAXN 1000005

long long P[800000];
long long mn[MAXN];
void get_p(long long *p, long long n)
{
	for (long long i=2; i<=n; ++i) {
		if (!mn[i]) p[mn[i]=++*p]=i;
		for (long long j=1; j<=mn[i] && i*p[j]<=n; ++j) mn[i*p[j]]=j;
	}
}

// n! 中 p 的次数 
long long get_cntP(long long n, long long p)
{
	long long ret=0;
	while (n>1)
		ret+= n/=p;
	return ret; 
}


long long fun(long long x)
{
	return (x+3)*x/2 %Ha;
}


int solve()
{
	long long n;
	scanf("%lld",&n);
	if (n==0) return 0;
	
	long long ans=0;
	for (int i=1; P[i]<=n; i++)
		ans= ( fun(get_cntP(n,P[i])) * (1+ans) %Ha +ans )%Ha;
		
	ans++;
	ans=(ans+Ha)%Ha;
	printf("%lld\n",ans);
	return 1;
	
}


int main()
{
	get_p(P,1000000);

	while (solve())
		continue;
	return 0;
}