[中山市选2011] 完全平方数

最新推荐文章于 2020-11-29 16:48:30 发布

转载最新推荐文章于 2020-11-29 16:48:30 发布 · 61 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/JYYHH/p/11320205.html

本文详细解析了莫比乌斯反演的基本原理及应用，通过具体题目讲解如何运用反演公式解决数学问题，提供了代码实现并分析了时间复杂度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题面

大概是比较基础的反演了，设 f(n) 为前n个数有多少个不是xxx，g(n)=n，则：

g(n) = ∑ f(n/i^2)，大概就是枚举 1~n中每个数的最大平方因子是多少了(注意极大也对，因为这里极大因子=最大因子)。

我们反演一下，可以得到 f(n) = ∑ μ(i) * g(n/i^2)，(这只是莫比乌斯反演的其中一种形式)，我可以来提示一下怎么推这个式子。

我们把 ∑ μ(i) * g(n/i^2) 中的所有g()用 g(n) = ∑ f(n/j^2) 展开，然后原式= ∑μ(i) * ∑ f(n/(i^2 * j^2))，我们把求和号交换，并让 i*j = d，可以得到原式= ∑f(n/d) * ∑μ(p) * [p|d]。

众所周知 e(x) = ∑ μ(y) * [y|x] = [x=1] (这个不管你用二项式定理还是数论函数卷积都可以证啦)，所以原式 = f(n) 得证。

然后这个题 f(n) 可以在 O(sqrt(n)) 的复杂度下求出之后就是个傻题了2333

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e6;

int zs[N/10+5],t,miu[N+5],T,n;
bool v[N+5];

inline void init(){
	miu[1]=1;
	for(int i=2;i<=N;i++){
		if(!v[i]) zs[++t]=i,miu[i]=-1;
		
		for(int j=1,u;j<=t&&(u=zs[j]*i)<=N;j++){
			v[u]=1;
			if(!(i%zs[j])) break;
			miu[u]=-miu[i];
		}
	}
}

inline ll Get(ll x){
	ll an=0,now=1;
	for(int i=1;now<=x;now+=(ll)(2*i+1),i++) an+=miu[i]*(ll)(x/now);
	return an; 
}

int main(){
	init();
	
	scanf("%d",&T);
	
	while(T--){
		scanf("%d",&n);
		
    	ll l=1,r=1644934081,mid,ans=0;
    	while(l<=r){
    		mid=l+r>>1;
    		if(Get(mid)>=n) ans=mid,r=mid-1;
    		else l=mid+1;
    	}	
	
    	printf("%lld\n",ans);
    }
	return 0;
}

转载于:https://www.cnblogs.com/JYYHH/p/11320205.html