zoukankan      html  css  js  c++  java
  • BZOJ.2440.[中山市选2011]完全平方数(莫比乌斯函数 二分)

    题目链接

    总感觉博客园的(Markdown)很。。(gouzhi),可以看这的


    题意即求第(k)个无平方因子数。

    无平方因子数(Square-Free Number),即分解之后所有质因数的次数都为1的数

    可以想到莫比乌斯函数,假设(n)是答案,那么有$$k=n-sum_{i=1}^n(1-|mu(i)|)$$
    (从这里能看出(x)的上界,后面的(sum)肯定是(<frac{n}{2})的,所以(nleq 2*k)
    二分一个(n),求([1,n])中有多少个无平方因子数。
    既然带着个平方就都开方。根据容斥,对于([1,sqrt{n}])中的质数,答案为$$[1,n]中0个质数平方倍数的个数-1个质数平方倍数的个数+2个质数平方倍数的个数-ldots$$
    即对于奇数个质数平方贡献为负,偶数个贡献为正,枚举这些质因子。若存在某个质因子的次数(>1),那么对答案没有影响(如(pi^2*pj^2)在计算(pi*pj)时统计了个数)。这也符合莫比乌斯函数的特点。那么答案可以写为:$$sum_{i=1}^{lfloorsqrt{n} floor}mu(i)*lfloorfrac{n}{i^2} floor$$

    也可以整除分块... 右端点是这样的:(r=sqrt{frac{n}{n/i^2}})

    (r)最大是(2e9),所以(l+r)可能爆int!
    就我被这个坑朝了吧。。

    //1308kb	2148ms
    #include <cmath>
    #include <cstdio>
    const int N=5e4;
    
    int cnt,P[10005],mu[N+3],pw[N+3];
    bool Not_P[N+3];
    
    void Init()
    {
    	mu[1]=1;
    	for(int i=2; i<N; ++i)
    	{
    		if(!Not_P[i]) P[++cnt]=i,mu[i]=-1;
    		for(int j=1; j<=cnt&&i*P[j]<N; ++j)
    		{
    			Not_P[i*P[j]]=1;
    			if(i%P[j]) mu[i*P[j]]=-mu[i];
    			else {mu[i*P[j]]=0; break;}
    		}
    	}
    	for(int i=1; i<N; ++i) pw[i]=i*i;
    }
    bool Check(long long n,int K)
    {
    	int res=0;//res=n
    	for(int i=1,lim=sqrt(n); i<=lim; ++i)
    		if(mu[i]/*这个?*/) res+=mu[i]*(n/pw[i]);
    	return res>=K;
    }
    
    int main()
    {
    	Init();
    	int T,K; scanf("%d",&T);
    	long long l,r,mid;//!
    	while(T--)
    	{
    		scanf("%d",&K), l=1, r=K<<1;
    		while(l<r)
    			if(Check(mid=(l+r)>>1,K)) r=mid;
    			else l=mid+1;
    		printf("%lld
    ",l);
    	}
    	return 0;
    }
    
  • 相关阅读:
    java+selenium自动化-IE浏览器搭建自动化环境
    python中的opencv
    随机森林参数说明
    剑指offer
    Python中常用的包--sklearn
    Anaconda安装,jupyter notebook 使用说明
    C++中的Public 、Private、Protected 区别
    C++类中的Static关键字二
    C++类中的Static关键字
    c语言二级指针内存模型
  • 原文地址:https://www.cnblogs.com/SovietPower/p/8696476.html
Copyright © 2011-2022 走看看