AcWing 1294. 樱花数学推导+约数个数

最新推荐文章于 2022-08-01 17:36:04 发布

karshey

最新推荐文章于 2022-08-01 17:36:04 发布

阅读量182

点赞数

分类专栏：我的ACM之路文章标签：数论

本文链接：https://blog.csdn.net/karshey/article/details/121000250

版权

我的ACM之路专栏收录该内容

159 篇文章 15 订阅

订阅专栏

这篇博客探讨了一个编程问题，涉及计算小于给定数n的所有质数的分布。作者通过编写C++代码实现了求解过程，利用了质数筛法优化计算效率。在解决过程中遇到了整型溢出的问题，并进行了修正。最后，代码成功输出了质数分布并模运算的结果。

摘要由CSDN通过智能技术生成

题
 参考1
参考2

中途WA了一发，MOD错了，中间过程爆int了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> pii;
#define pb push_back
#define fi first
#define se second
#define mem(a,x) memset(a,x,sizeof(a));
#define db double 
#define fir(i,a,n) for(int i=a;i<=n;i++)
#define debug(x) cout<<#x<<" "<<x<<endl;
const int inf=0x3f3f3f3f;
const int MOD=1e9+7;
//======================
const int N=1e6+10;
int pri[N],cnt=0,v[N];
int n;
ll ans=1;
void prime()
{
	for(int i=2;i<=1e3;i++)
	{
		if(!v[i])
		{
			pri[cnt++]=i;
			for(int j=i*i;j<=1e6;j+=i)
			{
				v[j]=1;
			}
		}
	}
	for(int i=1e3+1;i<=1e6;i++)
	{
		if(!v[i]) pri[cnt++]=i;
	}
}

void solve()
{
	for(int i=0;i<cnt;i++)
	{
		if(pri[i]>n) break;
		
		ll t=pri[i];
		ll temp=0;
		while(t<=n)
		{
			temp+=n/t;
			//temp%=MOD;
			t*=pri[i];
		}
		ans=(ans*(2*temp+1))%MOD;
	}
	cout<<ans%MOD;
}
int main()
{
	prime();
	cin>>n;
	solve();
	return 0; 
}