bzoj2721[Violet 5]樱花 数论

Posted 2020-09-22 GXZlegend

题目描述

输入

输出

样例输入

2

样例输出

3

---

题解

数论

设1/x+1/y=1/m，那么xm+ym=xy，所以xy-xm-ym+m^2=m^2，所以(x-m)(y-m)=m^2.

所以解的数量就是m^2的约数个数。

所以只需要算出n!中每个素数的出现次数即可。

我们可以先快筛出1~n的素数，然后考虑每个素数出现的次数。

而p出现的次数为包含p^1的数的个数+包含p^2的数的个数+...+包含p^k的数的个数，我们可以迭代来求。

最后把它们乘2加1再乘到一起即可。

#include <cstdio>
#include <algorithm>
#define N 1000010
using namespace std;
int prime[N] , tot , cnt[N];
bool np[N];
int main()
{
	int n , i , j;
	long long ans = 1;
	scanf("%d" , &n);
	for(i = 2 ; i <= n ; i ++ )
	{
		if(!np[i]) prime[++tot] = i;
		for(j = 1 ; j <= tot && i * prime[j] <= n ; j ++ )
		{
			np[i * prime[j]] = 1;
			if(i % prime[j] == 0) break;
		}
	}
	for(i = 1 ; i <= tot ; i ++ )
	{
		for(j = n ; j ; j /= prime[i]) cnt[i] += j / prime[i];
		ans = ans * (2 * cnt[i] + 1) % 1000000007;
	}
	printf("%lld\n" , ans);
	return 0;
}