【BZOJ 2721】樱花

最新推荐文章于 2019-09-07 15:19:00 发布

forever_dreams

最新推荐文章于 2019-09-07 15:19:00 发布

阅读量302

点赞数 1

分类专栏： --------数学问题-------- # 线性筛文章标签：线性筛素数分解质因数

本文链接：https://blog.csdn.net/forever_dreams/article/details/83572641

版权

--------数学问题-------- 同时被 2 个专栏收录

141 篇文章 1 订阅

订阅专栏

线性筛

12 篇文章 0 订阅

订阅专栏

【题目】

传送门

Description
在这里插入图片描述
Input

Output

Sample Input

Sample Output

HINT
在这里插入图片描述

【分析】

令 $a = n!$ ，由于 $y > n!$ ，设 $y = a + b$

$\because \frac{1}{x}+\frac{1}{y}=\frac{1}{n!}$

$\therefore \frac{1}{x}+\frac{1}{a+b}=\frac{1}{a}$

$\therefore \frac{x+a+b}{(a+b)x}=\frac{1}{a}$

$\therefore ax+a^2+ab=ax+bx$

$\therefore x=\frac{a^2}{b}+a$

由于 $x$ 是正整数，所以 $\frac{a^2}{b}$ 必须是正整数，所以相当于是要找 $a^2$ 的因数个数

那么，题目就变成了：给定 $n$ ，找出 $n!)^2$ 的因数个数，直接线性筛+分解质因数就行了

【代码】

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1000005
#define Mod 1000000007
using namespace std;
int n,sum,num[N],prime[N],Min[N];
bool mark[N];
void linear_sieves()
{
	int i,j;
	memset(mark,true,sizeof(mark));
	mark[0]=mark[1]=false;
	for(i=2;i<=n;++i)
	{
		if(mark[i])  prime[++sum]=i,Min[i]=sum;
		for(j=1;j<=sum&&i*prime[j]<=n;++j)
		{
			Min[i*prime[j]]=j;
			mark[i*prime[j]]=false;
			if(i%prime[j]==0)  break;
		}
	}
}
void work(int x)
{
	while(x!=1)
	{
		num[Min[x]]++;
		x/=prime[Min[x]];
	}
}
int main()
{
	int i,ans=1;
	scanf("%d",&n);
	linear_sieves();
	for(i=1;i<=n;++i)  work(i);
	for(i=1;i<=sum;++i)  ans=1ll*ans*(num[i]*2+1)%Mod;
	printf("%d",ans);
	return 0;
}