LOJ#6682. 梦中的数论（min_25筛）

题目描述：https://loj.ac/problem/6682

 

 

 

 

 

题解

比较简单的一道题

我们可以比较明显的看出，我们要求的式子是一个有实际意义的式子

如果把 i 分解质因数，把j+k换成k'（k'>j）

那么题目就是在求i的约数的无序对数，即

然后问题就转化为了求

和

直接上两个min_25筛

 

结果考试的时候把min_25写错了（各种错：预处理、递归s(a/pw,i+1)……），一直没调出来

思想不严谨，老是要犯低级错误（而且心态不好，最后比较慌。。。然后没调出来）

代码：

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
#define N 1000000
#define LL long long
const int mod=998244353;
const int inv2=499122177;
int tot,s[N+5];
bool vis[N+5];
int id1[N+5],id2[N+5];
LL prime[N+5],n,a[N+5];int m,blk;
inline int id(LL x){return x<=blk?id1[x]:id2[n/x];}
void min25()
{
	blk=sqrt(n+0.5);m=0;
	LL i,j;vis[1]=1;
	for(i=2;i<=blk;i++){
		if(!vis[i]) prime[++tot]=i;
		for(j=1;j<=tot;j++){
			LL tmp=1ll*i*prime[j];
			if(tmp>blk) break;
			vis[tmp]=1;
			if(i%prime[j]==0)break;
		}
	}
	for(i=1;i<=n;i=j+1){
		j=n/(n/i);
		a[++m]=n/i;
		if(a[m]<=blk)id1[a[m]]=m;
		else id2[n/a[m]]=m;
		s[m]=(a[m]-1)%mod;
	}
	for(j=1;j<=tot;j++)
		for(i=1;i<=m&&1ll*prime[j]*prime[j]<=a[i];i++)
			s[i]=((1ll*s[i]-(s[id(a[i]/prime[j])]-(j-1)))%mod+1ll*mod)%mod;
}
int solve0(LL a,int b)
{
	if(n<=1||a<prime[b])return 0;
	LL i,j,pw;
	int ret=2ll*(1ll*s[id(a)]-(b-1)+1ll*mod)%mod;
	for(i=b;i<=tot&&1ll*prime[i]*prime[i]<=a;i++)
		for(j=1,pw=prime[i];1ll*pw*prime[i]<=a;j++,pw*=prime[i])
			ret=(1ll*ret+1ll*solve0(a/pw,i+1)*(j+1)+1ll*(j+2))%mod;
	return ret;
}
int solve1(LL a,int b)
{
	if(n<=1||a<prime[b])return 0;
	LL i,j,pw;
	int ret=4ll*(1ll*s[id(a)]-(b-1)+1ll*mod)%mod;
	for(i=b;i<=tot&&1ll*prime[i]*prime[i]<=a;i++)
		for(j=1,pw=prime[i];1ll*pw*prime[i]<=a;j++,pw*=prime[i])
			ret=(1ll*ret+1ll*solve1(a/pw,i+1)*(j+1)%mod*(j+1)+1ll*(j+2)*(j+2))%mod;
	return ret;
}
int main()
{
	//freopen("5.in","r",stdin);
	//freopen("5.out","w",stdout);
	scanf("%lld",&n);
	min25();
	printf("%d",((1ll*solve1(n,1)-1ll*solve0(n,1))*inv2%mod+1ll*mod)%mod);
}

 

还有，思想不要僵化，不用把两个函数分别预处理，它们的 f 都是常函数，所以只需要筛一个质数个数就可以了。。。